36 namespace incompressible
55 RASModel(typeName, U, phi, lamTransportModel),
188 nut_ = Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_);
189 nut_.correctBoundaryConditions();
191 if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
196 "(const volVectorField& U, const surfaceScalarField& phi,"
197 "transportModel& lamTransportModel)"
198 ) <<
"couplingFactor = " << couplingFactor_
199 <<
" is not in range 0 - 1" <<
nl
231 if (couplingFactor_.
value() > 0.0)
238 (1.0 - couplingFactor_)*nut_,
259 if (RASModel::read())
272 if (couplingFactor_.
value() < 0.0 || couplingFactor_.
value() > 1.0)
275 <<
"couplingFactor = " << couplingFactor_
276 <<
" is not in range 0 - 1"
314 -
fvm::Sp(C2_*epsilon_/k_, epsilon_)
333 if (isA<wallFvPatch>(curPatch))
337 label faceCelli = curPatch.
faceCells()[facei];
339 *=
min(G[faceCelli]/(0.5*
mag(
tr(P[faceCelli])) + SMALL), 1.0);
352 +
fvm::Sp(Clrr1_*epsilon_/k_, R_)
355 - (2.0/3.0*(1 - Clrr1_)*
I)*epsilon_
382 nut_ = Cmu_*
sqr(k_)/epsilon_;
392 if (isA<wallFvPatch>(curPatch))
410 = (faceAreas[facei]/magFaceAreas[facei])*snGradU[facei];
416 Rw[facei].xy() = tauw.
xy();
417 Rw[facei].xz() = tauw.
xz();
418 Rw[facei].yz() = tauw.
yz();