34 namespace incompressible
45 void RASModel::printCoeffs()
78 turbulence_(lookup(
"turbulence")),
79 printCoeffs_(lookupOrDefault<Switch>(
"printCoeffs",
false)),
80 coeffDict_(subOrEmptyDict(type +
"Coeffs")),
83 epsilon0_(
"epsilon0", k0_.dimensions()/
dimTime, SMALL),
84 epsilonSmall_(
"epsilonSmall", epsilon0_.dimensions(), SMALL),
86 omegaSmall_(
"omegaSmall", omega0_.dimensions(), SMALL),
91 epsilon0_.readIfPresent(*
this);
92 epsilonSmall_.readIfPresent(*
this);
93 omega0_.readIfPresent(*
this);
94 omegaSmall_.readIfPresent(*
this);
129 dict.
lookup(
"RASModel") >> modelName;
132 Info<<
"Selecting RAS turbulence model " << modelName <<
endl;
134 dictionaryConstructorTable::iterator cstrIter =
135 dictionaryConstructorTablePtr_->find(modelName);
137 if (cstrIter == dictionaryConstructorTablePtr_->end())
141 "RASModel::New(const volVectorField&, "
142 "const surfaceScalarField&, transportModel&)"
143 ) <<
"Unknown RASModel type " << modelName
145 <<
"Valid RASModel types are :" << endl
146 << dictionaryConstructorTablePtr_->sortedToc()
157 scalar RASModel::yPlusLam(
const scalar
kappa,
const scalar
E)
const
161 for (
int i=0; i<10; i++)
172 const fvPatch& curPatch = mesh_.boundary()[patchNo];
177 if (isA<wallFvPatch>(curPatch))
182 /
nu().boundaryField()[patchNo];
188 "tmp<scalarField> RASModel::yPlus(const label patchNo) const"
189 ) <<
"Patch " << patchNo <<
" is not a wall. Returning null field"
203 if (turbulence_ && mesh_.changing())
210 bool RASModel::read()
214 lookup(
"turbulence") >> turbulence_;
218 coeffDict_ <<= *dictPtr;
221 k0_.readIfPresent(*
this);
222 epsilon0_.readIfPresent(*
this);
223 epsilonSmall_.readIfPresent(*
this);
224 omega0_.readIfPresent(*
this);
225 omegaSmall_.readIfPresent(*
this);