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/dynamicLagrangian: Added support for fvOptions

Browse Source
This commit is contained in:
Henry Weller
2015-12-03 14:52:48 +00:00
parent 2485269283
commit 1ae82d6e0c
1 changed files with 11 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C
View File

@ -24,6 +24,7 @@ License
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "dynamicLagrangian.H" #include "dynamicLagrangian.H"
#include "fvOptions.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -42,6 +43,9 @@ void dynamicLagrangian<BasicTurbulenceModel>::correctNut
{ {
this->nut_ = (flm_/fmm_)*sqr(this->delta())*mag(dev(symm(gradU))); this->nut_ = (flm_/fmm_)*sqr(this->delta())*mag(dev(symm(gradU)));
this->nut_.correctBoundaryConditions(); this->nut_.correctBoundaryConditions();
fv::options::New(this->mesh_).correct(this->nut_);
BasicTurbulenceModel::correctNut();
} }
@ -157,6 +161,7 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
// Local references // Local references
const surfaceScalarField& phi = this->phi_; const surfaceScalarField& phi = this->phi_;
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
fv::options& fvOptions(fv::options::New(this->mesh_));
LESeddyViscosity<BasicTurbulenceModel>::correct(); LESeddyViscosity<BasicTurbulenceModel>::correct();
@ -190,11 +195,13 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
== ==
invT*LM invT*LM
- fvm::Sp(invT, flm_) - fvm::Sp(invT, flm_)
+ fvOptions(flm_)
); );
flmEqn.relax(); flmEqn.relax();
fvOptions.constrain(flmEqn);
flmEqn.solve(); flmEqn.solve();
fvOptions.correct(flm_);
bound(flm_, flm0_); bound(flm_, flm0_);
volScalarField MM(M && M); volScalarField MM(M && M);
@ -206,11 +213,13 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
== ==
invT*MM invT*MM
- fvm::Sp(invT, fmm_) - fvm::Sp(invT, fmm_)
+ fvOptions(fmm_)
); );
fmmEqn.relax(); fmmEqn.relax();
fvOptions.constrain(fmmEqn);
fmmEqn.solve(); fmmEqn.solve();
fvOptions.correct(fmm_);
bound(fmm_, fmm0_); bound(fmm_, fmm0_);
correctNut(gradU); correctNut(gradU);