research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

reactingEulerFoam, twoPhaseEulerFoam: Added fvOption sources to the face-based momentum equations

Browse source
This commit is contained in:
Henry Weller 2016-07-19 16:56:08 +01:00
parent d5da93c477
commit eab6e3199d
2 changed files with 5 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/pUf/UEqns.H
View file

@ -29,6 +29,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
+ 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))
- fvOptions(alpha1, rho1, U1)
); );
U1Eqn.relax(); U1Eqn.relax();
fvOptions.constrain(U1Eqn); fvOptions.constrain(U1Eqn);
@ -44,6 +45,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
+ 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))
- fvOptions(alpha2, rho2, U2)
); );
U2Eqn.relax(); U2Eqn.relax();
fvOptions.constrain(U2Eqn); fvOptions.constrain(U2Eqn);

3
applications/solvers/multiphase/twoPhaseEulerFoam/pUf/UEqns.H
View file

@ -3,6 +3,7 @@ Info<< "Constructing face momentum equations" << endl;
MRF.correctBoundaryVelocity(U1); MRF.correctBoundaryVelocity(U1);
MRF.correctBoundaryVelocity(U2); MRF.correctBoundaryVelocity(U2);
MRF.correctBoundaryVelocity(U); MRF.correctBoundaryVelocity(U);
Info<< "Constructing face momentum equations" << endl;
fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime); fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime);
fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime); fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
@ -29,6 +30,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
+ 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))
- fvOptions(alpha1, rho1, U1)
); );
U1Eqn.relax(); U1Eqn.relax();
fvOptions.constrain(U1Eqn); fvOptions.constrain(U1Eqn);
@ -43,6 +45,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
+ 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))
- fvOptions(alpha2, rho2, U2)
); );
U2Eqn.relax(); U2Eqn.relax();
fvOptions.constrain(U2Eqn); fvOptions.constrain(U2Eqn);