exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Updated cross-coupling terms in coupledKEpsilon turbulence model. Robert Keser

Browse source
This commit is contained in:
Hrvoje Jasak 2017-09-21 13:45:06 +01:00
commit 64a0b7bf57
1 changed files with 7 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C
View file

@ -253,9 +253,9 @@ void coupledKEpsilon::correct()
+ fvm::div(phi_, epsilon_) + fvm::div(phi_, epsilon_)
+ fvm::SuSp(-fvc::div(phi_), epsilon_) + fvm::SuSp(-fvc::div(phi_), epsilon_)
- fvm::laplacian(DepsilonEff(), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_)
+ fvm::Sp(2*C2_*epsilon_/k_, epsilon_)
== ==
2*C1_*Cmu_*magSqr(symm(fvc::grad(U_)))*k_ C1_*G*epsilon_/k_
- fvm::Sp(2*C2_*epsilon_/k_, epsilon_)
+ C2_*sqr(epsilon_)/k_ + C2_*sqr(epsilon_)/k_
); );
@ -268,7 +268,6 @@ void coupledKEpsilon::correct()
volScalarField coupling volScalarField coupling
( (
"coupling", "coupling",
-2*C1_*Cmu_*magSqr(symm(fvc::grad(U_)))
-C2_*sqr(epsilon_/k_) -C2_*sqr(epsilon_/k_)
); );
scalarField& couplingIn = coupling.internalField(); scalarField& couplingIn = coupling.internalField();
@ -295,8 +294,10 @@ void coupledKEpsilon::correct()
+ fvm::div(phi_, k_) + fvm::div(phi_, k_)
+ fvm::SuSp(-fvc::div(phi_), k_) + fvm::SuSp(-fvc::div(phi_), k_)
- fvm::laplacian(DkEff(), k_) - fvm::laplacian(DkEff(), k_)
+ fvm::Sp(Cmu_*k_/(nut_ + nutSmall), k_) ==
- G G
+ Cmu_*sqr(k_)/(nut_+nutSmall)
- fvm::Sp(2*Cmu_*k_/(nut_+nutSmall), k_)
); );
kEqn.relax(); kEqn.relax();