34 namespace compressible
45 void RASModel::printCoeffs()
79 turbulence_(lookup(
"turbulence")),
80 printCoeffs_(lookupOrDefault<Switch>(
"printCoeffs",
false)),
81 coeffDict_(subOrEmptyDict(type +
"Coeffs")),
84 epsilon0_(
"epsilon0", k0_.dimensions()/
dimTime, SMALL),
85 epsilonSmall_(
"epsilonSmall", epsilon0_.dimensions(), SMALL),
87 omegaSmall_(
"omegaSmall", omega0_.dimensions(), SMALL),
92 epsilon0_.readIfPresent(*
this);
93 epsilonSmall_.readIfPresent(*
this);
94 omega0_.readIfPresent(*
this);
95 omegaSmall_.readIfPresent(*
this);
131 dict.
lookup(
"RASModel") >> modelName;
134 Info<<
"Selecting RAS turbulence model " << modelName <<
endl;
136 dictionaryConstructorTable::iterator cstrIter =
137 dictionaryConstructorTablePtr_->find(modelName);
139 if (cstrIter == dictionaryConstructorTablePtr_->end())
143 "RASModel::New(const volScalarField&, "
144 "const volVectorField&, const surfaceScalarField&, "
146 ) <<
"Unknown RASModel type " << modelName
148 <<
"Valid RASModel types are :" << endl
149 << dictionaryConstructorTablePtr_->sortedToc()
155 cstrIter()(
rho,
U,
phi, thermophysicalModel)
162 scalar RASModel::yPlusLam(
const scalar
kappa,
const scalar
E)
const
166 for (
int i=0; i<10; i++)
177 const fvPatch& curPatch = mesh_.boundary()[patchNo];
182 if (isA<wallFvPatch>(curPatch))
188 mu().boundaryField()[patchNo].patchInternalField()
189 /rho_.boundaryField()[patchNo]
196 "tmp<scalarField> RASModel::yPlus(const label patchNo) const"
197 ) <<
"Patch " << patchNo <<
" is not a wall. Returning null field"
209 if (mesh_.changing())
216 bool RASModel::read()
220 lookup(
"turbulence") >> turbulence_;
224 coeffDict_ <<= *dictPtr;
227 k0_.readIfPresent(*
this);
228 epsilon0_.readIfPresent(*
this);
229 epsilonSmall_.readIfPresent(*
this);
230 omega0_.readIfPresent(*
this);
231 omegaSmall_.readIfPresent(*
this);