84 void calcIncompressibleYPlus
93 wallFunctionPatchField;
104 const volScalarField::GeometricBoundaryField nutPatches =
107 bool foundNutPatch =
false;
110 if (isA<wallFunctionPatchField>(nutPatches[
patchi]))
112 foundNutPatch =
true;
114 const wallFunctionPatchField& nutPw =
115 dynamic_cast<const wallFunctionPatchField&
>
121 Info<<
"Patch " << patchi
122 <<
" named " << nutPw.patch().name()
123 <<
" y+ : min: " <<
min(Yp) <<
" max: " <<
max(Yp)
130 Info<<
" no " << wallFunctionPatchField::typeName <<
" patches"
136 void calcCompressibleYPlus
144 typedef compressible::RASModels::mutWallFunctionFvPatchScalarField
145 wallFunctionPatchField;
156 if (!rhoHeader.headerOk())
162 Info <<
"Reading field rho\n" <<
endl;
169 basicPsiThermo::New(mesh)
175 compressible::RASModel::New
184 const volScalarField::GeometricBoundaryField mutPatches =
187 bool foundMutPatch =
false;
188 forAll(mutPatches, patchi)
190 if (isA<wallFunctionPatchField>(mutPatches[patchi]))
192 foundMutPatch =
true;
194 const wallFunctionPatchField& mutPw =
195 dynamic_cast<const wallFunctionPatchField&
>
201 Info<<
"Patch " << patchi
202 <<
" named " << mutPw.patch().name()
203 <<
" y+ : min: " <<
min(Yp) <<
" max: " <<
max(Yp)
210 Info<<
" no " << wallFunctionPatchField::typeName <<
" patches"
216 int main(
int argc,
char *argv[])
218 timeSelector::addOptions();
222 argList::validOptions.insert(
"compressible",
"");
233 runTime.
setTime(timeDirs[timeI], timeI);
235 fvMesh::readUpdateState state = mesh.
readUpdate();
238 if (timeI == 0 || state != fvMesh::UNCHANGED)
240 Info<<
"Calculating wall distance\n" <<
endl;
242 Info<<
"Writing wall distance to field "
270 if (UHeader.headerOk())
272 Info <<
"Reading field U\n" <<
endl;
277 calcCompressibleYPlus(mesh, runTime, U, yPlus);
281 calcIncompressibleYPlus(mesh, runTime, U, yPlus);