36 namespace incompressible
46 const label patchI = patch().index();
51 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
64 nuw[facei]*(yPlus[facei]*
kappa_/
log(
E_*yPlus[facei]) - 1.0);
78 const label patchI = patch().index();
82 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
89 scalar kappaRe = kappa_*magUp[facei]*y[facei]/nuw[facei];
91 scalar yp = yPlusLam_;
92 scalar ryPlusLam = 1.0/yp;
95 scalar yPlusLast = 0.0;
100 yp = (kappaRe + yp)/(1.0 +
log(E_*yp));
102 }
while(
mag(ryPlusLam*(yp - yPlusLast)) > 0.01 && ++iter < 10 );
104 yPlus[facei] =
max(0.0, yp);
175 const label patchI = patch().index();
190 writeLocalEntries(os);
191 writeEntry(
"value", os);