25const surfaceScalarField
Kdf(
"Kdf",
fluid.Kdf());
30PtrList<surfaceScalarField>
rAUfs;
33 new surfaceScalarField
35 IOobject::groupName(
"rAUf",
phase1.name()),
39 + fvc::interpolate(U1Eqn.A())
46 new surfaceScalarField
48 IOobject::groupName(
"rAUf",
phase2.name()),
52 + fvc::interpolate(U2Eqn.A())
70 volScalarField
rho(
"rho",
fluid.rho());
75 surfaceScalarField rhof1(fvc::interpolate(
rho1));
76 surfaceScalarField rhof2(fvc::interpolate(
rho2));
84 IOobject::groupName(
"alpharAUf",
phase1.name()),
90 IOobject::groupName(
"alpharAUf",
phase2.name()),
101 const surfaceScalarField phigF1
111 const surfaceScalarField phigF2
123 surfaceScalarField phiHbyA1
125 IOobject::groupName(
"phiHbyA",
phase1.name()),
133 + fvc::flux(U1Eqn.H())
137 surfaceScalarField phiHbyA2
139 IOobject::groupName(
"phiHbyA",
phase2.name()),
147 + fvc::flux(U2Eqn.H())
152 PtrList<surfaceScalarField> phiKdPhifs(
fluid.phiKdPhifs(
rAUfs));
153 const surfaceScalarField& phiKdPhif1 = phiKdPhifs[0];
154 const surfaceScalarField& phiKdPhif2 = phiKdPhifs[1];
167 const surfaceScalarField
rAUf
177 p_rgh.boundaryFieldRef(),
184 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
188 tmp<fvScalarMatrix> pEqnComp1, pEqnComp2;
190 if (
phase1.compressible())
194 surfaceScalarField phid1
196 IOobject::groupName(
"phid",
phase1.name()),
210 + fvm::div(phid1,
p_rgh) - fvm::Sp(fvc::div(phid1),
p_rgh)
214 pEqnComp1.ref().relax();
227 if (
phase2.compressible())
231 surfaceScalarField phid2
233 IOobject::groupName(
"phid",
phase2.name()),
247 + fvm::div(phid2,
p_rgh) - fvm::Sp(fvc::div(phid2),
p_rgh)
251 pEqnComp2.ref().relax();
269 if (pEqnComp1.valid())
271 pEqnComp1.ref() -= (optEqn1 &
rho1)/
rho1;
281 if (pEqnComp2.valid())
283 pEqnComp2.ref() -= (optEqn2 &
rho2)/
rho2;
295 PtrList<volScalarField> dmdts(
fluid.dmdts());
298 if (pEqnComp1.valid())
300 pEqnComp1.ref() -= dmdts[0]/
rho1;
304 pEqnComp1 = fvm::Su(- dmdts[0]/
rho1,
p_rgh);
309 if (pEqnComp2.valid())
311 pEqnComp2.ref() -= dmdts[1]/
rho2;
315 pEqnComp2 = fvm::Su(- dmdts[1]/
rho2,
p_rgh);
323 while (
pimple.correctNonOrthogonal())
325 fvScalarMatrix pEqnIncomp
332 fvScalarMatrix pEqn(pEqnIncomp);
334 if (pEqnComp1.valid())
339 if (pEqnComp2.valid())
351 if (
pimple.finalNonOrthogonalIter())
353 surfaceScalarField mSfGradp(
"mSfGradp", pEqnIncomp.flux()/
rAUf);
357 const surfaceScalarField phi1s
362 const surfaceScalarField phi2s
367 surfaceScalarField
phir
376 U1 = fvc::reconstruct(
MRF.absolute(
phi1));
377 U1.correctBoundaryConditions();
380 U2 = fvc::reconstruct(
MRF.absolute(
phi2));
381 U2.correctBoundaryConditions();
385 if (pEqnComp1.valid())
389 if (pEqnComp2.valid())
409 p_rgh.correctBoundaryConditions();
const uniformDimensionedVectorField & g
const surfaceScalarField & ghf
const volScalarField & gh
const volScalarField & alpha1
const volScalarField & psi2
surfaceScalarField & phi2
const volScalarField & alpha2
const surfaceScalarField & alphaPhi1
surfaceScalarField & phi1
const volScalarField & psi1
const dimensionedScalar & pMin
const surfaceScalarField & alphaPhi2
surfaceScalarField phir(fvc::flux(UdmModel.Udm()))
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
volScalarField p_rgh_0(p_rgh)
surfaceScalarField ghSnGradRho(ghf *fvc::snGrad(rho) *mesh.magSf())
setSnGrad< fixedFluxPressureFvPatchScalarField >(p_rgh.boundaryFieldRef(),(phiHbyA.boundaryField() - MRF.relative(phib))/(mesh.magSf().boundaryField() *rAUf.boundaryField()))
PtrList< surfaceScalarField > phiFfs(fluid.phiFfs(rAUfs))
const surfaceScalarField alphaf2("alphaf2", scalar(1) - alphaf1)
const surfaceScalarField alpharAUf1(fvc::interpolate(max(alpha1, phase1.residualAlpha()) *rAU1))
const surfaceScalarField alpharAUf2(fvc::interpolate(max(alpha2, phase2.residualAlpha()) *rAU2))
const surfaceScalarField alphaf1("alphaf1", fvc::interpolate(alpha1))
const surfaceScalarField & rAUf1
const surfaceScalarField & rAUf2
surfaceScalarField alphaRhof20("alphaRhof20", fvc::interpolate(max(alpha2.oldTime(), phase2.residualAlpha()) *rho2.oldTime()))
surfaceScalarField alphaRhof10("alphaRhof10", fvc::interpolate(max(alpha1.oldTime(), phase1.residualAlpha()) *rho1.oldTime()))
const surfaceScalarField Kdf("Kdf", fluid.Kdf())
const surfaceScalarField & phiFf2
PtrList< surfaceScalarField > rAUfs
const surfaceScalarField & phiFf1
PtrList< surfaceScalarField > AFfs(fluid.AFfs())