36 namespace compressible
56 RASModel(typeName, rho, U, phi, thermophysicalModel),
219 if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
224 "( const volScalarField&, const volVectorField&"
225 ", const surfaceScalarField&, incompressibleTransportModel&)"
226 ) <<
"couplingFactor = " << couplingFactor_
227 <<
" is not in range 0 - 1" <<
nl
231 mut_ = Cmu_*rho_*
sqr(k_)/(epsilon_ + epsilonSmall_);
232 mut_.correctBoundaryConditions();
235 alphat_.correctBoundaryConditions();
265 if (couplingFactor_.value() > 0.0)
290 if (RASModel::read())
292 Cmu_.readIfPresent(coeffDict());
293 Clrr1_.readIfPresent(coeffDict());
294 Clrr2_.readIfPresent(coeffDict());
295 C1_.readIfPresent(coeffDict());
296 C2_.readIfPresent(coeffDict());
297 Cs_.readIfPresent(coeffDict());
298 Ceps_.readIfPresent(coeffDict());
299 sigmaR_.readIfPresent(coeffDict());
300 sigmaEps_.readIfPresent(coeffDict());
301 Prt_.readIfPresent(coeffDict());
302 couplingFactor_.readIfPresent(coeffDict());
304 if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
307 <<
"couplingFactor = " << couplingFactor_
308 <<
" is not in range 0 - 1" <<
nl
326 mut_ = rho_*Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_);
327 mut_.correctBoundaryConditions();
331 alphat_.correctBoundaryConditions();
342 epsilon_.boundaryField().updateCoeffs();
352 C1_*rho_*G*epsilon_/k_
353 -
fvm::Sp(C2_*rho_*epsilon_/k_, epsilon_)
358 epsEqn().boundaryManipulate(epsilon_.boundaryField());
361 bound(epsilon_, epsilon0_);
372 if (isA<wallFvPatch>(curPatch))
376 label faceCelli = curPatch.
faceCells()[facei];
378 *=
min(G[faceCelli]/(0.5*
mag(
tr(P[faceCelli])) + SMALL), 100.0);
390 +
fvm::Sp(Clrr1_*rho_*epsilon_/k_, R_)
393 - (2.0/3.0*(1 - Clrr1_)*
I)*rho_*epsilon_
408 k0_.value(), -GREAT, -GREAT,
420 mut_ = rho_*Cmu_*
sqr(k_)/epsilon_;
421 mut_.correctBoundaryConditions();
425 alphat_.correctBoundaryConditions();
434 if (isA<wallFvPatch>(curPatch))
444 = mesh_.Sf().boundaryField()[
patchi];
447 = mesh_.magSf().boundaryField()[
patchi];
453 = (faceAreas[facei]/magFaceAreas[facei])*snGradU[facei];
456 tensor tauw = -(mutw[facei]/rhow[facei])*2*
dev(
symm(gradUw));
459 Rw[facei].xy() = tauw.
xy();
460 Rw[facei].xz() = tauw.
xz();
461 Rw[facei].yz() = tauw.
yz();