35 namespace incompressible
47 tmp<volScalarField> kOmegaSST::F1(
const volScalarField& CDkOmega)
const
61 (scalar(1)/betaStar_)*
sqrt(k_)/(omega_*y_),
62 scalar(500)*
nu()/(
sqr(y_)*omega_)
64 (4*alphaOmega2_)*k_/(CDkOmegaPlus*
sqr(y_))
72 tmp<volScalarField> kOmegaSST::F2()
const
78 (scalar(2)/betaStar_)*
sqrt(k_)/(omega_*y_),
79 scalar(500)*
nu()/(
sqr(y_)*omega_)
97 RASModel(typeName, U, phi, lamTransportModel),
238 bound(omega_, omega0_);
241 nut_.correctBoundaryConditions();
302 if (RASModel::read())
343 omega_.boundaryField().updateCoeffs();
359 -
fvm::Sp(beta(F1)*omega_, omega_)
362 (F1 - scalar(1))*CDkOmega/omega_,
369 omegaEqn().boundaryManipulate(omega_.boundaryField());
382 min(G, c1_*betaStar_*k_*omega_)
383 -
fvm::Sp(betaStar_*omega_, k_)
392 nut_ = a1_*k_/
max(a1_*omega_, F2()*
sqrt(2*S2));