36 namespace incompressible
55 (KsPlus - 2.25)/87.75 + Cs*KsPlus,
56 sin(0.4258*(
log(KsPlus) - 0.811))
61 return (1.0 + Cs*KsPlus);
68 const label patchI = patch().index();
74 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
76 const scalar Cmu25 =
pow(
Cmu_, 0.25);
83 label faceCellI = patch().faceCells()[faceI];
85 scalar uStar = Cmu25*
sqrt(k[faceCellI]);
86 scalar
yPlus = uStar*y[faceI]/nuw[faceI];
87 scalar KsPlus = uStar*
Ks_[faceI]/nuw[faceI];
98 scalar limitingNutw =
max(nutw[faceI], nuw[faceI]);
116 Info<<
"yPlus = " << yPlus
117 <<
", KsPlus = " << KsPlus
118 <<
", Edash = " << Edash
119 <<
", nutw = " << nutw[faceI]
151 Ks_(ptf.
Ks_, mapper),
164 Ks_(
"Ks", dict, p.
size()),
165 Cs_(
"Cs", dict, p.
size())
199 nutWallFunctionFvPatchScalarField::autoMap(m);
211 nutWallFunctionFvPatchScalarField::rmap(ptf, addr);
214 refCast<const nutRoughWallFunctionFvPatchScalarField>(ptf);
217 Cs_.rmap(nrwfpsf.
Cs_, addr);
227 writeEntry(
"value", os);