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

MRFtwoPhaseEulerFoam: improved handling of the flux BCs for MRF

Browse Source
This commit is contained in:
Henry
2012-03-16 12:26:15 +00:00
parent 438dd02098
commit 82be4453f4
1 changed files with 5 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
applications/solvers/multiphase/twoPhaseEulerFoam/MRFtwoPhaseEulerFoam/pEqn.H
View File

@ -69,6 +69,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);
surfaceScalarField Dp
@ -92,18 +97,8 @@
{
surfaceScalarField SfGradp(pEqn.flux()/Dp);
phi1.boundaryField() ==
(fvc::interpolate(U1) & mesh.Sf())().boundaryField();
mrfZones.relativeFlux(phi1);
phi1 = phiHbyA1 - rAU1f*SfGradp/rho1;
mrfZones.relativeFlux(phi1.oldTime());
phi2.boundaryField() ==
(fvc::interpolate(U2) & mesh.Sf())().boundaryField();
mrfZones.relativeFlux(phi2);
phi2 = phiHbyA2 - rAU2f*SfGradp/rho2;
mrfZones.relativeFlux(phi2.oldTime());
phi = alpha1f*phi1 + alpha2f*phi2;
p.relax();