search

AMReX-Fluids / IAMR

A parallel, adaptive mesh refinement (AMR) code that solves the variable-density incompressible Navier-Stokes equations.
https://amrex-fluids.github.io/IAMR/
81 stars 58 forks source link

Low Mach constraint options in IAMR #26

Closed tobyeh closed 4 years ago

tobyeh commented 6 years ago

Hello,

I am currently developing a model of atmospheric plume rise, for which I had previously implemented conservative projection method (Almgren et al 1998) in my own code, along with a variety of low mach constraints (incomp, bous, anelastic and psuedo-incomp). I'm keen now to take advantage of AMReX and develop IAMR for this purpose. My question is: are all these constraints currently implemented in IAMR or is there only an option for incompressible at present? I understand that MaestroeX does have the low mach psuedo incom constraint, but not the viscous terms - i intend to keep the diffusive parts available in IAMR to develop a suitable subgrid turbulence model, so my preference is to use and modify IAMR. Thanks in advance for your help

ajnonaka commented 6 years ago

IAMR allows for a non-zero divergence constraint you would need for low Mach number models. If you look at IAMR/Source/NavierStokesBase.cpp line 1794 there is a "have_divu" flag you can enable. I believe you specify what the RHS of the constraint equation div(u)=S using the calc_divu() function in NavierStokes.cpp.

I've never actually modified IAMR to try this but I know other people have. Does anybody know of any other modifications you have to make?

drummerdoc commented 6 years ago

There is a "hotspot" case that approximates a rising heated bubble, and the divu terms has the conduction bits that one needs. But I'm not 100% sure that all the models listed in the original question are supported and usable out-of-the-box.

asalmgren commented 6 years ago

There is also support there for anelastic but it hasn't been tested in quite a while so may be out of date.

Ann

On Wed, Oct 17, 2018 at 10:22 AM, Marc Day notifications@github.com wrote:

There is a "hotspot" case that approximates a rising heated bubble, and the divu terms has the conduction bits that one needs. But I'm not 100% sure that all the models listed in the original question are supported and usable out-of-the-box.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/AMReX-Codes/IAMR/issues/26#issuecomment-430715455, or mute the thread https://github.com/notifications/unsubscribe-auth/AKJPYvBhvlUGiMZ2qAjaP_X1zAqbvVbKks5ul2dwgaJpZM4XjaMq .