35 namespace compressible
61 (scalar(1)/betaStar_)*
sqrt(k_)/(omega_*y_),
62 scalar(500)*(
mu()/rho_)/(
sqr(y_)*omega_)
64 (4*alphaOmega2_)*k_/(CDkOmegaPlus*
sqr(y_))
78 (scalar(2)/betaStar_)*
sqrt(k_)/(omega_*y_),
79 scalar(500)*(
mu()/rho_)/(
sqr(y_)*omega_)
98 RASModel(typeName, rho, U, phi, thermophysicalModel),
260 bound(omega_, omega0_);
271 mut_.correctBoundaryConditions();
274 alphat_.correctBoundaryConditions();
297 k_.boundaryField().types()
332 bool kOmegaSST::read()
334 if (RASModel::read())
336 alphaK1_.readIfPresent(coeffDict());
337 alphaK2_.readIfPresent(coeffDict());
338 alphaOmega1_.readIfPresent(coeffDict());
339 alphaOmega2_.readIfPresent(coeffDict());
340 Prt_.readIfPresent(coeffDict());
341 gamma1_.readIfPresent(coeffDict());
342 gamma2_.readIfPresent(coeffDict());
343 beta1_.readIfPresent(coeffDict());
344 beta2_.readIfPresent(coeffDict());
345 betaStar_.readIfPresent(coeffDict());
346 a1_.readIfPresent(coeffDict());
347 c1_.readIfPresent(coeffDict());
366 mut_.correctBoundaryConditions();
370 alphat_.correctBoundaryConditions();
379 if (mesh_.changing())
396 omega_.boundaryField().updateCoeffs();
412 -
fvm::SuSp((2.0/3.0)*rhoGammaF1*divU, omega_)
416 rho_*(F1 - scalar(1))*CDkOmega/omega_,
423 omegaEqn().boundaryManipulate(omega_.boundaryField());
426 bound(omega_, omega0_);
435 min(G, (c1_*betaStar_)*rho_*k_*omega_)
437 -
fvm::Sp(rho_*betaStar_*omega_, k_)
446 mut_ = a1_*rho_*k_/
max(a1_*omega_,
F2()*
sqrt(2.0*S2));
447 mut_.correctBoundaryConditions();
451 alphat_.correctBoundaryConditions();