songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

compressibleTwoPhaseEulerFoam: Corrected boundary relative fluxes for MRF

Browse Source
This commit is contained in:
Henry 2012-03-15 20:01:49 +00:00
parent 0a50fdabe0
commit 7f2ff7c8ff
1 changed files with 5 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
View File

@ -42,6 +42,11 @@
);
mrfZones.relativeFlux(phiHbyA2);
mrfZones.relativeFlux(phi1.oldTime());
mrfZones.relativeFlux(phi1);
mrfZones.relativeFlux(phi2.oldTime());
mrfZones.relativeFlux(phi2);
surfaceScalarField phiHbyA("phiHbyA", alpha1f*phiHbyA1 + alpha2f*phiHbyA2);
HbyA1 += (1.0/rho1)*rAU1*dragCoeff*U2;
@ -101,16 +106,8 @@
{
surfaceScalarField mSfGradp = pEqnIncomp.flux()/Dp;
phi1.boundaryField() ==
(fvc::interpolate(U1) & mesh.Sf())().boundaryField();
phi1 = phiHbyA1 + rAlphaAU1f*mSfGradp/fvc::interpolate(rho1);
mrfZones.relativeFlux(phi1.oldTime());
phi2.boundaryField() ==
(fvc::interpolate(U2) & mesh.Sf())().boundaryField();
phi2 = phiHbyA2 + rAlphaAU2f*mSfGradp/fvc::interpolate(rho2);
mrfZones.relativeFlux(phi2.oldTime());
phi = alpha1f*phi1 + alpha2f*phi2;
dgdt =