research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

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 b84ae0ecac
commit a1f9be94ed
1 changed files with 7 additions and 1 deletions
Show stats 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 "fvOptions.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -170,6 +171,7 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correctNut
{
this->nut_ = nuTilda_*fv1;
this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut();
}
@ -412,6 +414,7 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correct()
const rhoField& rho = this->rho_;
const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_;
const volVectorField& U = this->U_;
fv::options& fvOptions(fv::options::New(this->mesh_));
LESeddyViscosity<BasicTurbulenceModel>::correct();
@ -436,11 +439,14 @@ void SpalartAllmarasDES<BasicTurbulenceModel>::correct()
Cw1_*alpha*rho*fw(Stilda, dTilda)*nuTilda_/sqr(dTilda),
nuTilda_
)
+ fvOptions(alpha, rho, nuTilda_)
);
nuTildaEqn().relax();
fvOptions.constrain(nuTildaEqn());
solve(nuTildaEqn);
bound(nuTilda_, dimensionedScalar("zero", nuTilda_.dimensions(), 0.0));
fvOptions.correct(nuTilda_);
bound(nuTilda_, dimensionedScalar("0", nuTilda_.dimensions(), 0.0));
nuTilda_.correctBoundaryConditions();
correctNut(fv1);