35 namespace compressible
47 RNGkEpsilon::RNGkEpsilon
55 RASModel(typeName, rho, U, phi, thermophysicalModel),
188 mut_ = Cmu_*rho_*
sqr(k_)/(epsilon_ + epsilonSmall_);
189 mut_.correctBoundaryConditions();
192 alphat_.correctBoundaryConditions();
215 k_.boundaryField().types()
250 bool RNGkEpsilon::read()
252 if (RASModel::read())
254 Cmu_.readIfPresent(coeffDict());
255 C1_.readIfPresent(coeffDict());
256 C2_.readIfPresent(coeffDict());
257 C3_.readIfPresent(coeffDict());
258 Prt_.readIfPresent(coeffDict());
259 sigmaEps_.readIfPresent(coeffDict());
260 Prt_.readIfPresent(coeffDict());
261 eta0_.readIfPresent(coeffDict());
262 beta_.readIfPresent(coeffDict());
278 mut_ = rho_*Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_);
279 mut_.correctBoundaryConditions();
283 alphat_.correctBoundaryConditions();
307 ((eta*(-eta/eta0_ + scalar(1)))/(beta_*eta3 + scalar(1)));
310 epsilon_.boundaryField().updateCoeffs();
319 (C1_ - R)*G*epsilon_/k_
320 -
fvm::SuSp(((2.0/3.0)*C1_ + C3_)*rho_*divU, epsilon_)
321 -
fvm::Sp(C2_*rho_*epsilon_/k_, epsilon_)
326 epsEqn().boundaryManipulate(epsilon_.boundaryField());
329 bound(epsilon_, epsilon0_);
341 -
fvm::Sp(rho_*(epsilon_)/k_, k_)
350 mut_ = rho_*Cmu_*
sqr(k_)/epsilon_;
351 mut_.correctBoundaryConditions();
355 alphat_.correctBoundaryConditions();