initial commits for bulk viscosity in compressible FHD

AMReX-FHD · Sep 20, 2023 · b2f3450 · b2f3450
1 parent c39702f
commit b2f3450
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 0 deletions.
diff --git a/src_compressible/compressible_functions.cpp b/src_compressible/compressible_functions.cpp
@@ -9,6 +9,7 @@ AMREX_GPU_MANAGED int compressible::do_2D;
 AMREX_GPU_MANAGED int compressible::all_correl;
 AMREX_GPU_MANAGED int compressible::nspec_surfcov = 0;
 AMREX_GPU_MANAGED bool compressible::do_reservoir = false;
+AMREX_GPU_MANAGED amrex::Real compressible::zeta_ratio = -1.0;
 
 void InitializeCompressibleNamespace()
 {
@@ -70,6 +71,11 @@ void InitializeCompressibleNamespace()
         do_reservoir = true;
     }
 
+    // bulk viscosity ratio
+    pp.query("zeta_ratio",zeta_ratio);
+    if ((amrex::Math::abs(visc_type) == 3) and (zeta_ratio < 0.0)) amrex::Abort("need non-negative zeta_ratio (ratio of bulk to shear viscosity) for visc_type = 3 (use bulk viscosity)");
+    if ((amrex::Math::abs(visc_type) == 3) and (zeta_ratio >= 0.0)) amrex::Print() << "bulk viscosity model selected; bulk viscosity ratio is: " << zeta_ratio << "\n";
+
     return;
 }
 

diff --git a/src_compressible/compressible_namespace.H b/src_compressible/compressible_namespace.H
@@ -8,6 +8,7 @@ namespace compressible {
     extern AMREX_GPU_MANAGED int all_correl;
     extern AMREX_GPU_MANAGED int nspec_surfcov;
     extern AMREX_GPU_MANAGED bool do_reservoir;
+    extern AMREX_GPU_MANAGED amrex::Real zeta_ratio;
 
 }
 
diff --git a/src_compressible/transCoeffs.cpp b/src_compressible/transCoeffs.cpp
@@ -75,6 +75,12 @@ void calculateTransportCoeffs(const MultiFab& prim_in,
                 }
             }
 
+            // set bulk viscosity
+            if (amrex::Math::abs(visc_type) == 3) {
+                zeta(i,j,k) = zeta_ratio * eta(i,j,k);
+            }
+
+
         });
     }
 }