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

twoPhaseEulerFoam: Remove k from the Reynolds stress because there is only a single pressure in the system

Browse Source 
and to avoid the generation of flow from the gradient of k the pressure gradient must balance the k term in
each U equation
This commit is contained in:
Henry
2012-02-28 15:54:02 +00:00
parent 44b0daa71e
commit ef49aaf03d
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
View File

@ -18,7 +18,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
volTensorField Rca
(
"Rca",
((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
(((2.0/3.0)*I)*nuEffa)*tr(gradUaT) - nuEffa*gradUaT
);
if (kineticTheory.on())
@ -62,7 +62,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
volTensorField Rcb
(
"Rcb",
((2.0/3.0)*I)*(k + nuEffb*tr(gradUbT)) - nuEffb*gradUbT
(((2.0/3.0)*I)*nuEffb)*tr(gradUbT) - nuEffb*gradUbT
);
surfaceScalarField phiRb