87 bool writeResults = !args.optionFound(
"noWrite");
97 if (phiHeader.headerOk())
99 autoPtr<surfaceScalarField> PePtr;
116 IOobject RASPropertiesHeader
125 IOobject LESPropertiesHeader
136 if (
phi.dimensions() == dimensionSet(0, 3, -1, 0, 0))
138 if (RASPropertiesHeader.headerOk())
144 autoPtr<incompressible::RASModel>
RASModel
146 incompressible::RASModel::New
168 * mesh.surfaceInterpolation::deltaCoeffs()
174 else if (LESPropertiesHeader.headerOk())
176 IOdictionary LESProperties(LESPropertiesHeader);
180 autoPtr<incompressible::LESModel>
sgsModel
199 * mesh.surfaceInterpolation::deltaCoeffs()
211 "transportProperties",
232 mesh.surfaceInterpolation::deltaCoeffs()
233 * (
mag(
phi)/mesh.magSf())*(runTime.deltaT()/
nu)
238 else if (
phi.dimensions() == dimensionSet(1, 0, -1, 0, 0))
240 if (RASPropertiesHeader.headerOk())
244 autoPtr<basicPsiThermo>
thermo(basicPsiThermo::New(mesh));
257 autoPtr<compressible::RASModel>
RASModel
259 compressible::RASModel::New
282 * mesh.surfaceInterpolation::deltaCoeffs()
288 else if (LESPropertiesHeader.headerOk())
290 IOdictionary LESProperties(LESPropertiesHeader);
292 autoPtr<basicPsiThermo>
thermo(basicPsiThermo::New(mesh));
305 autoPtr<compressible::LESModel>
sgsModel
324 * mesh.surfaceInterpolation::deltaCoeffs()
336 "transportProperties",
357 mesh.surfaceInterpolation::deltaCoeffs()
358 * (
mag(
phi)/(mesh.magSf()))*(runTime.deltaT()/
mu)
366 <<
"Incorrect dimensions of phi: " <<
phi.dimensions()
390 Info <<
"Pe max : " <<
max(PePtr()).value() <<
endl;