71 int main(
int argc,
char *argv[])
73 argList::validOptions.insert(
"ybl",
"scalar");
74 argList::validOptions.insert(
"Cbl",
"scalar");
75 argList::validOptions.insert(
"writenut",
"");
100 Info<<
"Calculating wall distance field" <<
endl;
119 <<
"Neither option 'ybl' or 'Cbl' have been provided to calculate"
120 " the boundary-layer thickness"
124 Info<<
"\nCreating boundary-layer for U of thickness "
125 << ybl.value() <<
" m" <<
nl <<
endl;
131 scalar yblv = ybl.value();
134 if (y[celli] <= yblv)
136 U[celli] *=
::pow(y[celli]/yblv, (1.0/7.0));
193 if (nuTildaHeader.headerOk())
195 Info<<
"Reading field nuTilda\n" <<
endl;
198 nuTilda.correctBoundaryConditions();
204 if (kHeader.headerOk() && epsilonHeader.headerOk())
209 Info<<
"Reading field epsilon\n" <<
endl;
213 k =
sqr(nut/(ck0*
min(y, ybl)));
214 k.correctBoundaryConditions();
218 epsilon.correctBoundaryConditions();
227 Info<<
nl <<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s"
228 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s"