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/phaseCompressible/RAS/mixtureKEpsilon: Added support for fvOptions

Browse Source
This commit is contained in:
Henry Weller
2015-12-03 15:43:27 +00:00
parent fde2866a8b
commit 934aa987d9
1 changed files with 12 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/TurbulenceModels/phaseCompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C
View File

@ -24,9 +24,9 @@ License
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "mixtureKEpsilon.H" #include "mixtureKEpsilon.H"
#include "fvOptions.H"
#include "bound.H" #include "bound.H"
#include "twoPhaseSystem.H" #include "twoPhaseSystem.H"
#include "dragModel.H"
#include "virtualMassModel.H" #include "virtualMassModel.H"
#include "fixedValueFvPatchFields.H" #include "fixedValueFvPatchFields.H"
#include "inletOutletFvPatchFields.H" #include "inletOutletFvPatchFields.H"
@ -335,6 +335,9 @@ void mixtureKEpsilon<BasicTurbulenceModel>::correctNut()
{ {
this->nut_ = Cmu_*sqr(k_)/epsilon_; this->nut_ = Cmu_*sqr(k_)/epsilon_;
this->nut_.correctBoundaryConditions(); this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut();
} }
@ -575,6 +578,8 @@ void mixtureKEpsilon<BasicTurbulenceModel>::correct()
volScalarField& km = km_(); volScalarField& km = km_();
volScalarField& epsilonm = epsilonm_(); volScalarField& epsilonm = epsilonm_();
fv::options& fvOptions(fv::options::New(this->mesh_));
eddyViscosity<RASModel<BasicTurbulenceModel> >::correct(); eddyViscosity<RASModel<BasicTurbulenceModel> >::correct();
// Update the effective mixture density // Update the effective mixture density
@ -657,13 +662,14 @@ void mixtureKEpsilon<BasicTurbulenceModel>::correct()
- fvm::SuSp(((2.0/3.0)*C1_)*divUm, epsilonm) - fvm::SuSp(((2.0/3.0)*C1_)*divUm, epsilonm)
- fvm::Sp(C2_*epsilonm/km, epsilonm) - fvm::Sp(C2_*epsilonm/km, epsilonm)
+ epsilonSource() + epsilonSource()
+ fvOptions(epsilonm)
); );
epsEqn().relax(); epsEqn().relax();
fvOptions.constrain(epsEqn());
epsEqn().boundaryManipulate(epsilonm.boundaryField()); epsEqn().boundaryManipulate(epsilonm.boundaryField());
solve(epsEqn); solve(epsEqn);
fvOptions.correct(epsilonm);
bound(epsilonm, this->epsilonMin_); bound(epsilonm, this->epsilonMin_);
@ -679,10 +685,13 @@ void mixtureKEpsilon<BasicTurbulenceModel>::correct()
- fvm::SuSp((2.0/3.0)*divUm, km) - fvm::SuSp((2.0/3.0)*divUm, km)
- fvm::Sp(epsilonm/km, km) - fvm::Sp(epsilonm/km, km)
+ kSource() + kSource()
+ fvOptions(km)
); );
kmEqn().relax(); kmEqn().relax();
fvOptions.constrain(kmEqn());
solve(kmEqn); solve(kmEqn);
fvOptions.correct(km);
bound(km, this->kMin_); bound(km, this->kMin_);
km.correctBoundaryConditions(); km.correctBoundaryConditions();