36 namespace incompressible
55 RASModel(typeName, U, phi, lamTransportModel),
226 nut_ = Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_);
227 nut_.correctBoundaryConditions();
229 if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
233 "LaunderGibsonRSTM::LaunderGibsonRSTM"
234 "(const volVectorField& U, const surfaceScalarField& phi,"
235 "transportModel& lamTransportModel)"
236 ) <<
"couplingFactor = " << couplingFactor_
237 <<
" is not in range 0 - 1" <<
nl
269 if (couplingFactor_.
value() > 0.0)
274 +
fvc::laplacian((1.0-couplingFactor_)*nut_, U,
"laplacian(nuEff,U)")
292 if (RASModel::read())
308 if (couplingFactor_.
value() < 0.0 || couplingFactor_.
value() > 1.0)
311 <<
"couplingFactor = " << couplingFactor_
312 <<
" is not in range 0 - 1"
355 -
fvm::Sp(C2_*epsilon_/k_, epsilon_)
374 if (isA<wallFvPatch>(curPatch))
378 label faceCelli = curPatch.
faceCells()[facei];
380 min(G[faceCelli]/(0.5*
mag(
tr(P[faceCelli])) + SMALL), 1.0);
394 +
fvm::Sp(Clg1_*epsilon_/k_, R_)
397 + (2.0/3.0*(Clg1_ - 1)*
I)*epsilon_
403 I*((yr_.
n() & reflect) & yr_.
n())
404 - 1.5*(yr_.
n()*(reflect & yr_.
n())
405 + (yr_.
n() & reflect)*yr_.
n())
406 )*
pow(Cmu_, 0.75)*
pow(k_, 1.5)/(kappa_*yr_*epsilon_)
432 nut_ = Cmu_*
sqr(k_)/epsilon_;
442 if (isA<wallFvPatch>(curPatch))
460 = (faceAreas[facei]/magFaceAreas[facei])*snGradU[facei];
466 Rw[facei].xy() = tauw.
xy();
467 Rw[facei].xz() = tauw.
xz();
468 Rw[facei].yz() = tauw.
yz();