diff --git a/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.C b/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.C index 1b0d5feb31fa25f55017bdd4cadf47e2d2caecdd..8db844d9a8fde8fc70e34df0cb3de5cb135a8393 100644 --- a/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.C +++ b/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.C @@ -132,8 +132,7 @@ kOmegaSSTSato<BasicTurbulenceModel>::gasTurbulence() const template<class BasicTurbulenceModel> void kOmegaSSTSato<BasicTurbulenceModel>::correctNut ( - const volScalarField& S2, - const volScalarField& F2 + const volScalarField& S2 ) { const PhaseCompressibleTurbulenceModel<transportModel>& gasTurbulence = @@ -145,7 +144,12 @@ void kOmegaSSTSato<BasicTurbulenceModel>::correctNut ); this->nut_ = - this->a1_*this->k_/max(this->a1_*this->omega_, this->b1_*F2*sqrt(S2)) + this->a1_*this->k_ + /max + ( + this->a1_*this->omega_, + this->b1_*this->F23()*sqrt(S2) + ) + sqr(1 - exp(-yPlus/16.0)) *Cmub_*gasTurbulence.transport().d()*gasTurbulence.alpha() *(mag(this->U_ - gasTurbulence.U())); diff --git a/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.H b/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.H index 37fd2a8b4c705f55b7dde7b8f329e3b4db0f87e5..ac42f2a77770f128c28ac9af4c85984332b940eb 100644 --- a/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.H +++ b/src/TurbulenceModels/phaseCompressible/RAS/kOmegaSSTSato/kOmegaSSTSato.H @@ -156,11 +156,7 @@ protected: // Protected Member Functions - virtual void correctNut - ( - const volScalarField& S2, - const volScalarField& F2 - ); + virtual void correctNut(const volScalarField& S2); public: