36 namespace compressible
44 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::calcYPlus
49 const label patchI = patch().index();
61 scalar kappaRe = kappa_*magUp[faceI]*y[faceI]/(muw[faceI]/rhow[faceI]);
63 scalar yp = yPlusLam_;
64 scalar ryPlusLam = 1.0/yp;
67 scalar yPlusLast = 0.0;
72 yp = (kappaRe + yp)/(1.0 +
log(E_*yp));
74 }
while (
mag(ryPlusLam*(yp - yPlusLast)) > 0.01 && ++iter < 10);
76 yPlus[faceI] =
max(0.0, yp);
84 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::calcMut()
const
86 const label patchI = patch().index();
101 if (yPlus[faceI] > yPlusLam_)
104 muw[faceI]*(yPlus[faceI]*kappa_/
log(E_*yPlus[faceI]) - 1.0);
114 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::
115 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField
125 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::
126 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField
138 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::
139 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField
150 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::
151 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField
160 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::
161 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField
174 mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::yPlus()
const
176 const label patchI = patch().index();
181 return calcYPlus(magUp);
185 void mutSpalartAllmarasStandardWallFunctionFvPatchScalarField::write
191 writeLocalEntries(os);
192 writeEntry(
"value", os);