Merge pull request #174 from AMReX-FHD/fix_giantflux_3d

Fix inputs_giantfluct_3d to work by ensuring bc YK's sum to 1
AMReX-FHD · Dec 5, 2024 · d61978a · d61978a
2 parents 9cd3bd9 + c3773dd
commit d61978a
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 2 deletions.
diff --git a/exec/compressible/inputs_giantfluct_3d b/exec/compressible/inputs_giantfluct_3d
@@ -98,7 +98,7 @@
 
   # Xk and Yk at the wall for Dirichlet (concentrations) - set one to zero
   # Ordering: (species 1, x-dir), (species 2, x-dir), ... (species 1, y-dir), ...
-  bc_Yk_y_lo = 0.2 0.09316672 0.70683296     # lo BC
-  bc_Yk_y_hi = 0.3 0.40683328 0.29316704     # hi BC
+  bc_Yk_y_lo = 0.2 0.09316672 0.70683328     # lo BC
+  bc_Yk_y_hi = 0.3 0.40683328 0.29316672     # hi BC
 
 
diff --git a/src_compressible_stag/boundaryStag.cpp b/src_compressible_stag/boundaryStag.cpp
@@ -32,6 +32,9 @@ void SetupCWallStag() {
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
            GetMolfrac(bc_Yk_x_lo,bc_Xk_x_lo);
        }
+       else {
+           Abort("SetupCWallStag: lo-x; mass or mole fractions do not sum to 1");
+       }
     }
 
     if (bc_mass_lo[0] >= 3) {
@@ -69,6 +72,8 @@ void SetupCWallStag() {
           GetMassfrac(bc_Xk_x_hi,bc_Yk_x_hi);
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
           GetMolfrac(bc_Yk_x_hi,bc_Xk_x_hi);
+       } else {
+           Abort("SetupCWallStag: hi-x; mass or mole fractions do not sum to 1");
        }
     }
 
@@ -108,6 +113,8 @@ void SetupCWallStag() {
           GetMassfrac(bc_Xk_y_lo,bc_Yk_y_lo);
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
           GetMolfrac(bc_Yk_y_lo,bc_Xk_y_lo);
+       } else {
+           Abort("SetupCWallStag: lo-y; mass or mole fractions do not sum to 1");
        }
     }
 
@@ -146,6 +153,8 @@ void SetupCWallStag() {
           GetMassfrac(bc_Xk_y_hi,bc_Yk_y_hi);
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
           GetMolfrac(bc_Yk_y_hi,bc_Xk_y_hi);
+       } else {
+           Abort("SetupCWallStag: hi-y; mass or mole fractions do not sum to 1");
        }
     }
 
@@ -185,6 +194,8 @@ void SetupCWallStag() {
           GetMassfrac(bc_Xk_z_lo,bc_Yk_z_lo);
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
           GetMolfrac(bc_Yk_z_lo,bc_Xk_z_lo);
+       } else {
+           Abort("SetupCWallStag: lo-z; mass or mole fractions do not sum to 1");
        }
     }
 
@@ -223,6 +234,8 @@ void SetupCWallStag() {
           GetMassfrac(bc_Xk_z_hi,bc_Yk_z_hi);
        } else if (amrex::Math::abs(sumy-1) < 1.e-10) {
           GetMolfrac(bc_Yk_z_hi,bc_Xk_z_hi);
+       } else {
+           Abort("SetupCWallStag: hi-z; mass or mole fractions do not sum to 1");
        }
     }