37 namespace incompressible
46 if (!isA<wallFvPatch>(patch()))
48 FatalErrorIn(
"nutWallFunctionFvPatchScalarField::checkType()")
49 <<
"Invalid wall function specification" <<
nl
50 <<
" Patch type for patch " << patch().name()
51 <<
" must be wall" <<
nl
52 <<
" Current patch type is " << patch().type() <<
nl <<
endl
66 for (
int i=0; i<10; i++)
77 const label patchI = patch().index();
83 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
85 const scalar Cmu25 =
pow(
Cmu_, 0.25);
92 label faceCellI = patch().faceCells()[faceI];
94 scalar
yPlus = Cmu25*y[faceI]*
sqrt(k[faceCellI])/nuw[faceI];
98 nutw[faceI] = nuw[faceI]*(yPlus*
kappa_/
log(
E_*yPlus) - 1.0);
122 fixedValueFvPatchScalarField(p, iF),
126 yPlusLam_(calcYPlusLam(kappa_, E_))
140 fixedValueFvPatchScalarField(ptf, p, iF, mapper),
157 fixedValueFvPatchScalarField(p, iF, dict),
161 yPlusLam_(calcYPlusLam(kappa_, E_))
172 fixedValueFvPatchScalarField(wfpsf),
188 fixedValueFvPatchScalarField(wfpsf, iF),
210 const label patchI = patch().index();
218 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
228 writeEntry(
"value", os);