multiphaseEulerFoam: Updated MRF support

This commit is contained in:
Henry
2012-03-19 11:54:52 +00:00
parent aa7066dad9
commit 7c7fc6e02b

View File

@ -132,6 +132,17 @@
phasei++; phasei++;
} }
phasei = 0;
forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
{
phaseModel& phase = iter();
mrfZones.relativeFlux(phase.phi().oldTime());
mrfZones.relativeFlux(phase.phi());
phasei++;
}
surfaceScalarField Dp surfaceScalarField Dp
( (
IOobject IOobject
@ -185,7 +196,6 @@
phase.phi() = phase.phi() =
phiHbyAs[phasei] phiHbyAs[phasei]
+ rAlphaAUfs[phasei]*mSfGradp/phase.rho(); + rAlphaAUfs[phasei]*mSfGradp/phase.rho();
mrfZones.relativeFlux(phase.phi().oldTime());
phi += alphafs[phasei]*phase.phi(); phi += alphafs[phasei]*phase.phi();
phasei++; phasei++;