diff --git a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C index c440294bad..0932bb524e 100644 --- a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C +++ b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C @@ -155,13 +155,13 @@ SpalartAllmaras::SpalartAllmaras : LESModel(modelName, U, phi, transport), - alphaNut_ + sigmaNut_ ( dimensioned::lookupOrAddToDict ( - "alphaNut", + "sigmaNut", coeffDict_, - 1.5 + 0.66666 ) ), kappa_ @@ -169,7 +169,7 @@ SpalartAllmaras::SpalartAllmaras dimensioned::lookupOrAddToDict ( "kappa", - *this, + coeffDict_, 0.41 ) ), @@ -227,7 +227,7 @@ SpalartAllmaras::SpalartAllmaras 0.07 ) ), - Cw1_(Cb1_/sqr(kappa_) + alphaNut_*(1.0 + Cb2_)), + Cw1_(Cb1_/sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_), Cw2_ ( dimensioned::lookupOrAddToDict @@ -301,11 +301,11 @@ void SpalartAllmaras::correct(const tmp& gradU) + fvm::div(phi(), nuTilda_) - fvm::laplacian ( - alphaNut_*(nuTilda_ + nu()), + (nuTilda_ + nu())/sigmaNut_, nuTilda_, "laplacian(DnuTildaEff,nuTilda)" ) - - alphaNut_*Cb2_*magSqr(fvc::grad(nuTilda_)) + - Cb2_/sigmaNut_*magSqr(fvc::grad(nuTilda_)) == Cb1_*STilda*nuTilda_ - fvm::Sp(Cw1_*fw(STilda, dTilda)*nuTilda_/sqr(dTilda), nuTilda_) @@ -358,7 +358,7 @@ bool SpalartAllmaras::read() { if (LESModel::read()) { - alphaNut_.readIfPresent(coeffDict()); + sigmaNut_.readIfPresent(coeffDict()); kappa_.readIfPresent(*this); Cb1_.readIfPresent(coeffDict()); Cb2_.readIfPresent(coeffDict()); @@ -366,7 +366,7 @@ bool SpalartAllmaras::read() Cv2_.readIfPresent(coeffDict()); CDES_.readIfPresent(coeffDict()); ck_.readIfPresent(coeffDict()); - Cw1_ = Cb1_/sqr(kappa_) + alphaNut_*(1.0 + Cb2_); + Cw1_ = Cb1_/sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_; Cw2_.readIfPresent(coeffDict()); Cw3_.readIfPresent(coeffDict()); diff --git a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H index a1f951ac92..03058f5b68 100644 --- a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H +++ b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H @@ -71,7 +71,7 @@ protected: // Protected data - dimensionedScalar alphaNut_; + dimensionedScalar sigmaNut_; dimensionedScalar kappa_;