AMReX-Fluids / incflo

Incompressible Navier-Stokes solver in complex geometries
https://amrex-fluids.github.io/incflo
BSD 3-Clause "New" or "Revised" License
46 stars 28 forks source link

Add InitialPressureProjection to enforce hydrostatic equilibrium #102

Closed cgilet closed 4 months ago

cgilet commented 4 months ago

in cases where that's not done through the background density and pressure.

Right now, this is a runtime option that is off by default. It does not check for any background state; it is up to the user to know when to use it and to ensure that the background density is turned off with incflo.ro_0 = 0

cgilet commented 4 months ago

@asalmgren we could also revert to having do_initial_pressure_projection as a runtime option instead of the last commit. We also still need to put in an abort for !constant_density && gravity transverse to outflow

asalmgren commented 4 months ago

I like making it a run-time option with default off for now