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

dynamicLagragian implemetation from dev. Adding phase support.

Browse Source
This commit is contained in:
sergio
2015-12-11 10:26:31 -08:00
parent 8179509780
commit 06a1ff719c
1 changed files with 11 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -155,6 +155,9 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
}
// Local references
const alphaField& alpha = this->alpha_;
const rhoField& rho = this->rho_;
const surfaceScalarField& alphaRhoPhi = this->alphaRhoPhi_;
const surfaceScalarField& phi = this->phi_;
const volVectorField& U = this->U_;
@ -185,11 +188,11 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
fvScalarMatrix flmEqn
(
fvm::ddt(flm_)
+ fvm::div(phi, flm_)
fvm::ddt(alpha, rho, flm_)
+ fvm::div(alphaRhoPhi, flm_)
==
invT*LM
- fvm::Sp(invT, flm_)
rho*invT*LM
- fvm::Sp(rho*invT, flm_)
);
flmEqn.relax();
@ -201,11 +204,11 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
fvScalarMatrix fmmEqn
(
fvm::ddt(fmm_)
+ fvm::div(phi, fmm_)
fvm::ddt(alpha, rho, fmm_)
+ fvm::div(alphaRhoPhi, fmm_)
==
invT*MM
- fvm::Sp(invT, fmm_)
rho*invT*MM
- fvm::Sp(rho*invT, fmm_)
);
fmmEqn.relax();