zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: kOmegaSSTLM improved robustness

Browse Source
This commit is contained in:
Andrew Heather
2018-01-18 12:09:06 +00:00
parent 8716ca25c0
commit f37b3c479a
1 changed files with 9 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/TurbulenceModels/turbulenceModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C
View File

@ -81,7 +81,12 @@ tmp<volScalarField::Internal> kOmegaSSTLM<BasicTurbulenceModel>::Fthetat
const volScalarField::Internal& omega = this->omega_();
const volScalarField::Internal& y = this->y_();
const volScalarField::Internal delta(375*Omega*nu*ReThetat_()*y/sqr(Us));
dimensionedScalar deltaMin("deltaMin", dimLength, SMALL);
volScalarField::Internal delta
(
max(375*Omega*nu*ReThetat_()*y/sqr(Us), deltaMin)
);
const volScalarField::Internal ReOmega(sqr(y)*omega/nu);
const volScalarField::Internal Fwake(exp(-sqr(ReOmega/1e5)));
@ -617,11 +622,11 @@ void kOmegaSSTLM<BasicTurbulenceModel>::correct()
return;
}
// Correct ReThetat and gammaInt
correctReThetatGammaInt();
// Correct k and omega
kOmegaSST<BasicTurbulenceModel>::correct();
// Correct ReThetat and gammaInt
correctReThetatGammaInt();
}