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

compressibleTwoPhaseEulerFoam: removed k contribution from the stress

Browse Source
This commit is contained in:
Henry
2012-03-15 15:10:14 +00:00
parent a70398ae2c
commit 12f2f0c09b
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H
View File

@ -18,7 +18,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
volTensorField Rc1 volTensorField Rc1
( (
"Rc1", "Rc1",
((2.0/3.0)*I)*(sqr(Ct)*k + nuEff1*tr(gradU1T)) - nuEff1*gradU1T (((2.0/3.0)*I)*nuEff1)*tr(gradU1T) - nuEff1*gradU1T
); );
if (kineticTheory.on()) if (kineticTheory.on())
@ -45,7 +45,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
- fvm::Sp(dragCoeff/rho1, U1) - fvm::Sp(dragCoeff/rho1, U1)
- alpha1*alpha2/rho1*(liftForce - Cvm*rho2*DDtU2) - alpha1*alpha2/rho1*(liftForce - Cvm*rho2*DDtU2)
); );
mrfZones.addCoriolis(alpha1, U1Eqn);
U1Eqn.relax(); U1Eqn.relax();
} }
@ -54,7 +54,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
volTensorField Rc2 volTensorField Rc2
( (
"Rc2", "Rc2",
((2.0/3.0)*I)*(k + nuEff2*tr(gradU2T)) - nuEff2*gradU2T (((2.0/3.0)*I)*nuEff2)*tr(gradU2T) - nuEff2*gradU2T
); );
U2Eqn = U2Eqn =
@ -76,7 +76,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
- fvm::Sp(dragCoeff/rho2, U2) - fvm::Sp(dragCoeff/rho2, U2)
+ alpha1*alpha2/rho2*(liftForce + Cvm*rho2*DDtU1) + alpha1*alpha2/rho2*(liftForce + Cvm*rho2*DDtU1)
); );
mrfZones.addCoriolis(alpha2, U2Eqn);
U2Eqn.relax(); U2Eqn.relax();
} }
} }