Replace nuFraExp to nuFra in UEqns.H

SedFoam · Oct 20, 2023 · d3c6f40 · d3c6f40
1 parent 3940686
commit d3c6f40
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 15 deletions.
diff --git a/solver/pU/UEqns.H b/solver/pU/UEqns.H
@@ -49,19 +49,14 @@ if (SUSaniso)
     susb = SUS*(K/rhob)*turbulenceb->nut()
                    *((SUS_I*iso-SUS_A*aniso) & fvc::grad(alpha));
 }
-// limiter for the explicit stress terms to avoid numerical instability and
-// too strong limitation on the time step
-volScalarField nuFraExp(nuFra);
-nuFraExp.min(nuMaxExp);
-nuFraExp.correctBoundaryConditions();
 
 // Compute explicit stress term for phase a
 volTensorField Rca
 (
     "Rca",
-    ((2.0/3.0)*I)*((nuEffa+nuFraExp/(alpha+alphaSmall))*tr(gradUaT)
+    ((2.0/3.0)*I)*((nuEffa+nuFra/(alpha+alphaSmall))*tr(gradUaT)
     + spherSigmaSGSa)
-    - (nuEffa+nuFraExp/(alpha+alphaSmall))*gradUaT
+    - (nuEffa+nuFra/(alpha+alphaSmall))*gradUaT
 );
 
 // Compute concentration gradient stress contribution on cell faces for phase a

diff --git a/solver/pUf/UEqns.H b/solver/pUf/UEqns.H
@@ -50,19 +50,13 @@ if (SUSaniso)
                    *((SUS_I*iso-SUS_A*aniso) & fvc::grad(alpha));
 }
 
-// limiter for the explicit stress terms to avoid numerical instability and
-// too strong limitation on the time step
-volScalarField nuFraExp(nuFra);
-nuFraExp.min(nuMaxExp);
-nuFraExp.correctBoundaryConditions();
-
 // Compute explicit stress term for phase a
 volTensorField Rca
 (
     "Rca",
-    ((2.0/3.0)*I)*((nuEffa+nuFraExp/(alpha+alphaSmall))*tr(gradUaT)
+    ((2.0/3.0)*I)*((nuEffa+nuFra/(alpha+alphaSmall))*tr(gradUaT)
     + spherSigmaSGSa)
-    - (nuEffa+nuFraExp/(alpha+alphaSmall))*gradUaT
+    - (nuEffa+nuFra/(alpha+alphaSmall))*gradUaT
 );
 
 // Compute concentration gradient stress contribution on cell faces for phase a