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

reactingTwoPhaseEulerFoam: Completed mass-transfer terms in both forms of momentum equations

Browse Source
This commit is contained in:
Henry Weller
2015-07-03 23:18:25 +01:00
parent 9c9ac368fe
commit 7c6e1386f1
2 changed files with 7 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H
View File

@ -18,8 +18,7 @@ volScalarField Kd(fluid.Kd());
U1Eqn = U1Eqn =
( (
fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1) fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1)
- fvm::Sp(phase1.continuityError() - dmdt21, U1) - fvm::Sp(phase1.continuityError() + dmdt21, U1) + dmdt21*U2
- dmdt21*U2
+ MRF.DDt(alpha1*rho1, U1) + MRF.DDt(alpha1*rho1, U1)
+ phase1.turbulence().divDevRhoReff(U1) + phase1.turbulence().divDevRhoReff(U1)
== ==
@ -44,8 +43,7 @@ volScalarField Kd(fluid.Kd());
U2Eqn = U2Eqn =
( (
fvm::ddt(alpha2, rho2, U2) + fvm::div(alphaRhoPhi2, U2) fvm::ddt(alpha2, rho2, U2) + fvm::div(alphaRhoPhi2, U2)
- fvm::Sp(phase2.continuityError() + dmdt12, U2) - fvm::Sp(phase2.continuityError() - dmdt12, U2) - dmdt12*U1
+ dmdt12*U1
+ MRF.DDt(alpha2*rho2, U2) + MRF.DDt(alpha2*rho2, U2)
+ phase2.turbulence().divDevRhoReff(U2) + phase2.turbulence().divDevRhoReff(U2)
== ==

5
applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pUf/UEqns.H
View File

@ -21,10 +21,14 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
+ MRF.DDt(U2) + MRF.DDt(U2)
); );
const volScalarField dmdt12(posPart(fluid.dmdt()));
const volScalarField dmdt21(negPart(fluid.dmdt()));
{ {
U1Eqn = U1Eqn =
( (
fvm::div(alphaRhoPhi1, U1) - fvm::Sp(fvc::div(alphaRhoPhi1), U1) fvm::div(alphaRhoPhi1, U1) - fvm::Sp(fvc::div(alphaRhoPhi1), U1)
+ fvm::Sp(dmdt12, U1) + dmdt21*U2
+ MRF.DDt(alpha1*rho1, U1) + MRF.DDt(alpha1*rho1, U1)
+ phase1.turbulence().divDevRhoReff(U1) + phase1.turbulence().divDevRhoReff(U1)
+ Vm*(UgradU1 - (UgradU2 & U2)) + Vm*(UgradU1 - (UgradU2 & U2))
@ -39,6 +43,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
U2Eqn = U2Eqn =
( (
fvm::div(alphaRhoPhi2, U2) - fvm::Sp(fvc::div(alphaRhoPhi2), U2) fvm::div(alphaRhoPhi2, U2) - fvm::Sp(fvc::div(alphaRhoPhi2), U2)
- fvm::Sp(dmdt21, U2) - dmdt12*U1
+ MRF.DDt(alpha2*rho2, U2) + MRF.DDt(alpha2*rho2, U2)
+ phase2.turbulence().divDevRhoReff(U2) + phase2.turbulence().divDevRhoReff(U2)
+ Vm*(UgradU2 - (UgradU1 & U1)) + Vm*(UgradU2 - (UgradU1 & U1))