33 namespace incompressible
45 tmp<volScalarField> SpalartAllmaras::chi()
const
51 tmp<volScalarField> SpalartAllmaras::fv1(
const volScalarField& chi)
const
54 return chi3/(chi3 +
pow3(Cv1_));
58 tmp<volScalarField> SpalartAllmaras::fv2
64 return 1.0/
pow3(scalar(1) + chi/Cv2_);
68 tmp<volScalarField> SpalartAllmaras::fv3
79 *(3*(scalar(1) + chiByCv2) +
sqr(chiByCv2))
80 /
pow3(scalar(1) + chiByCv2);
84 tmp<volScalarField> SpalartAllmaras::fw(
const volScalarField& Stilda)
const
95 r.boundaryField() == 0.0;
112 RASModel(typeName, U, phi, lamTransportModel),
151 Cw1_(Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_),
324 if (RASModel::read())
330 Cw1_ = Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_;
364 + fv2(chi, fv1)*nuTilda_/
sqr(kappa_*d_);
375 -
fvm::Sp(Cw1_*fw(Stilda)*nuTilda_/
sqr(d_), nuTilda_)
378 nuTildaEqn().relax();