32 IOobject::groupName(
"rAUf",
phase1.name()),
36 + fvc::interpolate(U1Eqn.A())
43 IOobject::groupName(
"rAUf",
phase2.name()),
47 + fvc::interpolate(U2Eqn.A())
54tmp<surfaceScalarField>
Ff1;
55tmp<surfaceScalarField>
Ff2;
58 volScalarField
D(
fluid.D());
61 surfaceScalarField
Df1
63 fvc::interpolate(
D +
phase1.turbulence().pPrime())
67 surfaceScalarField
Df2
69 fvc::interpolate(
D +
phase2.turbulence().pPrime())
74 if (implicitPhasePressure)
98 volScalarField
rho(
"rho",
fluid.rho());
103 surfaceScalarField rhof1(fvc::interpolate(
rho1));
104 surfaceScalarField rhof2(fvc::interpolate(
rho2));
112 IOobject::groupName(
"alpharAUf",
phase1.name()),
118 IOobject::groupName(
"alpharAUf",
phase2.name()),
129 surfaceScalarField phig1
131 IOobject::groupName(
"phig",
phase1.name()),
140 surfaceScalarField phig2
142 IOobject::groupName(
"phig",
phase2.name()),
152 surfaceScalarField phiHbyA1
154 IOobject::groupName(
"phiHbyA",
phase1.name()),
163 + fvc::flux(U1Eqn.H())
170 surfaceScalarField phiHbyA2
172 IOobject::groupName(
"phiHbyA",
phase2.name()),
181 + fvc::flux(U2Eqn.H())
198 surfaceScalarField
rAUf
207 p_rgh.boundaryFieldRef(),
214 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
217 tmp<fvScalarMatrix> pEqnComp1;
218 tmp<fvScalarMatrix> pEqnComp2;
222 surfaceScalarField phid1
224 IOobject::groupName(
"phid",
phase1.name()),
227 surfaceScalarField phid2
229 IOobject::groupName(
"phid",
phase2.name()),
243 + fvm::div(phid1,
p_rgh) - fvm::Sp(fvc::div(phid1),
p_rgh)
246 pEqnComp1.ref().relax();
258 + fvm::div(phid2,
p_rgh) - fvm::Sp(fvc::div(phid2),
p_rgh)
261 pEqnComp2.ref().relax();
283 while (
pimple.correctNonOrthogonal())
285 fvScalarMatrix pEqnIncomp
293 pEqnComp1() + pEqnComp2() + pEqnIncomp,
297 if (
pimple.finalNonOrthogonalIter())
299 surfaceScalarField mSfGradp(
"mSfGradp", pEqnIncomp.flux()/
rAUf);
303 surfaceScalarField phi1s
310 surfaceScalarField phi2s
317 surfaceScalarField
phir
326 U1 = fvc::reconstruct(
MRF.absolute(
phi1));
327 U1.correctBoundaryConditions();
330 U2 = fvc::reconstruct(
MRF.absolute(
phi2));
331 U2.correctBoundaryConditions();
357 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()))
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())
surfaceScalarField snGradAlpha1(fvc::snGrad(alpha1) *mesh.magSf())
surfaceScalarField Df1(fvc::interpolate(rAU1 *(D+phase1.turbulence().pPrime())))
surfaceScalarField Df2(fvc::interpolate(rAU2 *(D+phase2.turbulence().pPrime())))
tmp< surfaceScalarField > Ff2
surfaceScalarField Vmf("Vmf", fluid.Vmf())
tmp< surfaceScalarField > Ff1
surfaceScalarField Ff(fluid.Ff())
const dimensionedScalar & D