alphatJayatillekeWallFunction: More robust stabilisation

Resolves bug report https://bugs.openfoam.org/view.php?id=3888
OpenFOAM · Sep 20, 2022 · a23937d · a23937d
1 parent 34d8167
commit a23937d
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/...Functions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.C b/...Functions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.C
@@ -222,7 +222,7 @@ void alphatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs()
 
             const scalar C = Pr*0.5*rhow[facei]*uTau*sqr(magUp[facei]);
 
-            alphaEff = (A - C)/(B + sign(B)*vSmall);
+            alphaEff = (A - C)/(B + sign(B)*rootVSmall);
         }
         else
         {
@@ -238,7 +238,7 @@ void alphatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs()
                 0.5*rhow[facei]*uTau
                *(Prt_*sqr(magUp[facei]) + (Pr - Prt_)*sqr(magUc));
 
-            alphaEff = (A - C)/(B + sign(B)*vSmall);
+            alphaEff = (A - C)/(B + sign(B)*rootVSmall);
         }
 
         // Bounds on turbulent thermal diffusivity