Revert back to just having do_initial_pressure_proj runtime option

cgilet · Jun 20, 2024 · 55f5576 · 55f5576
1 parent 1ddaf78
commit 55f5576
Show file tree

Hide file tree

Showing 4 changed files with 5 additions and 4 deletions.
diff --git a/src/incflo.H b/src/incflo.H
@@ -360,7 +360,8 @@ private:
     amrex::Real m_dt_change_max = amrex::Real(1.1);
 
     // Initial projection / iterations
-    bool m_do_initial_proj    = true;
+    bool m_do_initial_proj          = true;
+    bool m_do_initial_pressure_proj = false;
     int  m_initial_iterations = 3;
 
     // Use Boussinesq approximation for buoyancy?

diff --git a/src/incflo.cpp b/src/incflo.cpp
@@ -76,9 +76,7 @@ void incflo::InitData ()
         if (m_do_initial_proj) {
             InitialProjection();
 
-            if (!m_constant_density
-                && std::abs(m_gravity.sum()) > std::numeric_limits<Real>::epsilon()
-                && !(m_ro_0 > std::numeric_limits<Real>::epsilon()) ) {
+            if (m_do_initial_pressure_proj) {
                 InitialPressureProjection();
             }
         }

diff --git a/src/setup/init.cpp b/src/setup/init.cpp
@@ -38,6 +38,7 @@ void incflo::ReadParameters ()
         pp.query("steady_state_tol", m_steady_state_tol);
         pp.query("initial_iterations", m_initial_iterations);
         pp.query("do_initial_proj", m_do_initial_proj);
+        pp.query("do_initial_pressure_proj", m_do_initial_pressure_proj);
 
         pp.query("fixed_dt", m_fixed_dt);
         pp.query("cfl", m_cfl);

diff --git a/test_no_eb_2d/benchmark.rayleigh_taylor b/test_no_eb_2d/benchmark.rayleigh_taylor
@@ -38,6 +38,7 @@ incflo.advection_type   = "Godunov"
 
 incflo.constant_density = false
 incflo.ro_0             = 0             # Turn off background state
+incflo.do_initial_pressure_proj = 1     # to enforce hydrostatic equilibrium
 
 incflo.advect_tracer    = true