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

SSG: Added support for fvOptions

Browse Source
This commit is contained in:
Henry Weller
2015-12-03 13:12:27 +00:00
parent a6f4bf60b8
commit 05fa38f5e3
1 changed files with 9 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/TurbulenceModels/turbulenceModels/RAS/SSG/SSG.C
View File

@ -24,6 +24,7 @@ License
\*---------------------------------------------------------------------------*/
#include "SSG.H"
#include "fvOptions.H"
#include "wallFvPatch.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -40,6 +41,7 @@ void SSG<BasicTurbulenceModel>::correctNut()
{
this->nut_ = this->Cmu_*sqr(k_)/epsilon_;
this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut();
}
@ -290,6 +292,7 @@ void SSG<BasicTurbulenceModel>::correct()
const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_;
const volVectorField& U = this->U_;
volSymmTensorField& R = this->R_;
fv::options& fvOptions(fv::options::New(this->mesh_));
ReynoldsStress<RASModel<BasicTurbulenceModel> >::correct();
@ -311,13 +314,14 @@ void SSG<BasicTurbulenceModel>::correct()
==
Ceps1_*alpha*rho*G*epsilon_/k_
- fvm::Sp(Ceps2_*alpha*rho*epsilon_/k_, epsilon_)
+ fvOptions(alpha, rho, epsilon_)
);
epsEqn().relax();
fvOptions.constrain(epsEqn());
epsEqn().boundaryManipulate(epsilon_.boundaryField());
solve(epsEqn);
fvOptions.correct(epsilon_);
bound(epsilon_, this->epsilonMin_);
@ -364,10 +368,13 @@ void SSG<BasicTurbulenceModel>::correct()
+ C4_*dev(twoSymm(b&S))
+ C5_*twoSymm(b&Omega)
)
+ fvOptions(alpha, rho, R)
);
REqn().relax();
fvOptions.constrain(REqn());
solve(REqn);
fvOptions.correct(R);
this->boundNormalStress(R);