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

kOmegaSST: Added support for fvOptions

Browse Source
This commit is contained in:
Henry Weller
2015-12-02 20:33:43 +00:00
parent 1a71623b46
commit c005f6d68c
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/kOmegaSST/kOmegaSST.C
View File

@ -24,6 +24,7 @@ License
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "kOmegaSST.H" #include "kOmegaSST.H"
#include "fvOptions.H"
#include "bound.H" #include "bound.H"
#include "wallDist.H" #include "wallDist.H"
@ -112,6 +113,7 @@ void kOmegaSST<BasicTurbulenceModel>::correctNut(const volScalarField& S2)
{ {
this->nut_ = a1_*k_/max(a1_*omega_, b1_*F23()*sqrt(S2)); this->nut_ = a1_*k_/max(a1_*omega_, b1_*F23()*sqrt(S2));
this->nut_.correctBoundaryConditions(); this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut(); BasicTurbulenceModel::correctNut();
} }
@ -399,6 +401,7 @@ void kOmegaSST<BasicTurbulenceModel>::correct()
const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_; const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_;
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
volScalarField& nut = this->nut_; volScalarField& nut = this->nut_;
fv::options& fvOptions(fv::options::New(this->mesh_));
eddyViscosity<RASModel<BasicTurbulenceModel> >::correct(); eddyViscosity<RASModel<BasicTurbulenceModel> >::correct();
@ -446,13 +449,14 @@ void kOmegaSST<BasicTurbulenceModel>::correct()
) )
+ Qsas(S2, gamma, beta) + Qsas(S2, gamma, beta)
+ omegaSource() + omegaSource()
+ fvOptions(alpha, rho, omega_)
); );
omegaEqn().relax(); omegaEqn().relax();
fvOptions.constrain(omegaEqn());
omegaEqn().boundaryManipulate(omega_.boundaryField()); omegaEqn().boundaryManipulate(omega_.boundaryField());
solve(omegaEqn); solve(omegaEqn);
fvOptions.correct(omega_);
bound(omega_, this->omegaMin_); bound(omega_, this->omegaMin_);
} }
@ -467,10 +471,13 @@ void kOmegaSST<BasicTurbulenceModel>::correct()
- fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_) - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_)
- fvm::Sp(alpha*rho*betaStar_*omega_, k_) - fvm::Sp(alpha*rho*betaStar_*omega_, k_)
+ kSource() + kSource()
+ fvOptions(alpha, rho, k_)
); );
kEqn().relax(); kEqn().relax();
fvOptions.constrain(kEqn());
solve(kEqn); solve(kEqn);
fvOptions.correct(k_);
bound(k_, this->kMin_); bound(k_, this->kMin_);
correctNut(S2); correctNut(S2);