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

turbulenceModels/LES/SpalartAllmarasDES: Added support for fvOptions

Browse Source
This commit is contained in:
Henry Weller
2015-12-03 14:53:12 +00:00
parent 1ae82d6e0c
commit 154893db0f
1 changed files with 7 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/TurbulenceModels/turbulenceModels/LES/SpalartAllmarasDES/SpalartAllmarasDES.C
View File

@ -24,6 +24,7 @@ License
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "SpalartAllmarasDES.H" #include "SpalartAllmarasDES.H"
#include "fvOptions.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -170,6 +171,7 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correctNut
{ {
this->nut_ = nuTilda_*fv1; this->nut_ = nuTilda_*fv1;
this->nut_.correctBoundaryConditions(); this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut(); BasicTurbulenceModel::correctNut();
} }
@ -412,6 +414,7 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correct()
const rhoField& rho = this->rho_; const rhoField& rho = this->rho_;
const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_; const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_;
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
fv::options& fvOptions(fv::options::New(this->mesh_));
LESeddyViscosity<BasicTurbulenceModel>::correct(); LESeddyViscosity<BasicTurbulenceModel>::correct();
@ -436,11 +439,14 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correct()
Cw1_*alpha*rho*fw(Stilda, dTilda)*nuTilda_/sqr(dTilda), Cw1_*alpha*rho*fw(Stilda, dTilda)*nuTilda_/sqr(dTilda),
nuTilda_ nuTilda_
) )
+ fvOptions(alpha, rho, nuTilda_)
); );
nuTildaEqn().relax(); nuTildaEqn().relax();
fvOptions.constrain(nuTildaEqn());
solve(nuTildaEqn); solve(nuTildaEqn);
bound(nuTilda_, dimensionedScalar("zero", nuTilda_.dimensions(), 0.0)); fvOptions.correct(nuTilda_);
bound(nuTilda_, dimensionedScalar("0", nuTilda_.dimensions(), 0.0));
nuTilda_.correctBoundaryConditions(); nuTilda_.correctBoundaryConditions();
correctNut(fv1); correctNut(fv1);