35 namespace compressible
47 tmp<volScalarField> realizableKE::rCmu
54 tmp<volSymmTensorField> tS =
dev(
symm(gradU));
58 (2*
sqrt(2.0))*((S&S)&&S)
71 return 1.0/(A0_ + As*Us*k_/(epsilon_ + epsilonSmall_));
75 tmp<volScalarField> realizableKE::rCmu
82 return rCmu(gradU, S2, magS);
88 realizableKE::realizableKE
96 RASModel(typeName, rho, U, phi, thermophysicalModel),
203 bound(epsilon_, epsilon0_);
205 mut_ = rCmu(
fvc::grad(U_))*rho_*
sqr(k_)/(epsilon_ + epsilonSmall_);
206 mut_.correctBoundaryConditions();
209 alphat_.correctBoundaryConditions();
232 k_.boundaryField().types()
267 bool realizableKE::read()
269 if (RASModel::read())
271 Cmu_.readIfPresent(coeffDict());
272 A0_.readIfPresent(coeffDict());
273 C2_.readIfPresent(coeffDict());
274 sigmak_.readIfPresent(coeffDict());
275 sigmaEps_.readIfPresent(coeffDict());
276 Prt_.readIfPresent(coeffDict());
293 mut_.correctBoundaryConditions();
297 alphat_.correctBoundaryConditions();
321 epsilon_.boundaryField().updateCoeffs();
330 C1*rho_*magS*epsilon_
333 C2_*rho_*epsilon_/(k_ +
sqrt((
mu()/rho_)*epsilon_)),
340 epsEqn().boundaryManipulate(epsilon_.boundaryField());
343 bound(epsilon_, epsilon0_);
355 -
fvm::Sp(rho_*(epsilon_)/k_, k_)
363 mut_ = rCmu(gradU, S2, magS)*rho_*
sqr(k_)/epsilon_;
364 mut_.correctBoundaryConditions();
368 alphat_.correctBoundaryConditions();