32 IOobject::groupName(
"rAUf",
phase1.name()),
43 IOobject::groupName(
"rAUf",
phase2.name()),
54 tmp<surfaceScalarField>
Ff1;
55 tmp<surfaceScalarField>
Ff2;
74 if (implicitPhasePressure)
112 IOobject::groupName(
"alpharAUf",
phase1.name()),
118 IOobject::groupName(
"alpharAUf",
phase2.name()),
131 IOobject::groupName(
"phig",
phase1.name()),
142 IOobject::groupName(
"phig",
phase2.name()),
154 IOobject::groupName(
"phiHbyA",
phase1.name()),
172 IOobject::groupName(
"phiHbyA",
phase2.name()),
207 p_rgh.boundaryFieldRef(),
214 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
217 tmp<fvScalarMatrix> pEqnComp1;
218 tmp<fvScalarMatrix> pEqnComp2;
224 IOobject::groupName(
"phid",
phase1.name()),
229 IOobject::groupName(
"phid",
phase2.name()),
246 pEqnComp1.ref().relax();
261 pEqnComp2.ref().relax();
283 while (
pimple.correctNonOrthogonal())
293 pEqnComp1() + pEqnComp2() + pEqnIncomp,
297 if (
pimple.finalNonOrthogonalIter())
327 U1.correctBoundaryConditions();
331 U2.correctBoundaryConditions();
357 p_rgh.correctBoundaryConditions();
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix by subtracting the matrix multiplied by the current fi...
surfaceScalarField Vmf("Vmf", fluid.Vmf())
fvMatrix< scalar > fvScalarMatrix
const surfaceScalarField alphaf1("alphaf1", fvc::interpolate(alpha1))
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar & pMin
surfaceScalarField Df1(fvc::interpolate(rAU1 *(D+phase1.turbulence().pPrime())))
surfaceScalarField alphaRhof10("alphaRhof10", fvc::interpolate(max(alpha1.oldTime(), phase1.residualAlpha()) *rho1.oldTime()))
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
const volScalarField & alpha2
surfaceScalarField alphaRhof20("alphaRhof20", fvc::interpolate(max(alpha2.oldTime(), phase2.residualAlpha()) *rho2.oldTime()))
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
surfaceScalarField & phi1
const surfaceScalarField & alphaPhi2
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const volScalarField & psi2
surfaceScalarField phir(fvc::flux(UdmModel.Udm()))
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
const surfaceScalarField & rAUf2
surfaceScalarField & phi2
const uniformDimensionedVectorField & g
const surfaceScalarField Kdf("Kdf", fluid.Kdf())
tmp< surfaceScalarField > Ff1
const surfaceScalarField alpharAUf2(fvc::interpolate(max(alpha2, phase2.residualAlpha()) *rAU2))
const volScalarField & psi1
mixture thermo2().correctRho(psi2 *(p_rgh - p_rgh_0))
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
surfaceScalarField Df2(fvc::interpolate(rAU2 *(D+phase2.turbulence().pPrime())))
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
surfaceScalarField ghSnGradRho(ghf *fvc::snGrad(rho) *mesh.magSf())
surfaceScalarField snGradAlpha1(fvc::snGrad(alpha1) *mesh.magSf())
setSnGrad< fixedFluxPressureFvPatchScalarField >(p_rgh.boundaryFieldRef(),(phiHbyA.boundaryField() - MRF.relative(phib))/(mesh.magSf().boundaryField() *rAUf.boundaryField()))
tmp< surfaceScalarField > Ff2
const volScalarField & gh
volScalarField p_rgh_0(p_rgh)
const surfaceScalarField alphaf2("alphaf2", scalar(1) - alphaf1)
const surfaceScalarField alpharAUf1(fvc::interpolate(max(alpha1, phase1.residualAlpha()) *rAU1))
const dimensionedScalar & D
const surfaceScalarField & rAUf1
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const surfaceScalarField & alphaPhi1
surfaceScalarField Ff(fluid.Ff())
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
mixture thermo1().correctRho(psi1 *(p_rgh - p_rgh_0))
const volScalarField & alpha1