35 namespace incompressible
47 tmp<volScalarField> qZeta::fMu()
const
53 return exp((-scalar(2.5) + Rt/20.0)/
pow(scalar(1) + Rt/130.0, 3.0));
58 exp(-6.0/
sqr(scalar(1) + Rt/50.0))
59 *(scalar(1) + 3.0*
exp(-Rt/10.0));
64 tmp<volScalarField> qZeta::f2()
const
67 return scalar(1) - 0.3*
exp(-
sqr(Rt));
80 RASModel(typeName, U, phi, lamTransportModel),
165 k_.boundaryField().types()
179 epsilon_.boundaryField().types()
195 nut_ = Cmu_*fMu()*
sqr(k_)/(epsilon_ + epsilonSmall_);
196 nut_.correctBoundaryConditions();
257 if (RASModel::read())
297 (2.0*C1_ - 1)*G*zeta_/q_
327 epsilon_ = 2*q_*zeta_;
332 nut_ = Cmu_*fMu()*
sqr(k_)/epsilon_;