4 PtrList<volScalarField>
rAUs;
9 IOobject::groupName(
"rAU",
phase1.name()),
21 IOobject::groupName(
"rAU",
phase2.name()),
95 IOobject::groupName(
"HbyA",
phase1.name()),
108 IOobject::groupName(
"HbyA",
phase2.name()),
120 PtrList<surfaceScalarField> ddtCorrByAs(
fluid.ddtCorrByAs(
rAUs));
125 IOobject::groupName(
"phiHbyA",
phase1.name()),
126 fvc::flux(HbyA1) - phigF1 - ddtCorrByAs[0]
131 IOobject::groupName(
"phiHbyA",
phase2.name()),
132 fvc::flux(HbyA2) - phigF2 - ddtCorrByAs[1]
138 PtrList<surfaceScalarField> phiKdPhis(
fluid.phiKdPhis(
rAUs));
163 p_rgh.boundaryFieldRef(),
170 )/(
mesh.magSf().boundaryField()*
rAUf.boundaryField())
174 tmp<fvScalarMatrix> pEqnComp1, pEqnComp2;
175 if (
phase1.compressible())
181 IOobject::groupName(
"phid",
phase1.name()),
199 pEqnComp1.ref().relax();
211 if (
phase2.compressible())
217 IOobject::groupName(
"phid",
phase2.name()),
235 pEqnComp2.ref().relax();
253 if (pEqnComp1.valid())
255 pEqnComp1.ref() -= (optEqn1 &
rho1)/
rho1;
265 if (pEqnComp2.valid())
267 pEqnComp2.ref() -= (optEqn2 &
rho2)/
rho2;
278 PtrList<volScalarField> dmdts(
fluid.dmdts());
281 if (pEqnComp1.valid())
283 pEqnComp1.ref() -= dmdts[0]/
rho1;
292 if (pEqnComp2.valid())
294 pEqnComp2.ref() -= dmdts[1]/
rho2;
307 while (
pimple.correctNonOrthogonal())
319 if (pEqnComp1.valid())
324 if (pEqnComp2.valid())
333 if (
pimple.finalNonOrthogonalIter())
340 pEqnIncomp.flux()/
rAUf 366 if (pEqnComp1.valid())
370 if (pEqnComp2.valid())
378 mSfGradp = pEqnIncomp.flux()/
rAUf;
405 U1.correctBoundaryConditions();
409 U2.correctBoundaryConditions();
428 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...
fvMatrix< scalar > fvScalarMatrix
const surfaceScalarField & phiF2
const surfaceScalarField alphaf1("alphaf1", fvc::interpolate(alpha1))
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar & pMin
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.
const volScalarField Kd(fluid.Kd())
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
const volScalarField & alpha2
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
const volScalarField & rAU1
surfaceScalarField & phi1
GeometricField< vector, fvPatchField, volMesh > volVectorField
const surfaceScalarField & phiF1
const surfaceScalarField & alphaPhi2
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const volScalarField & psi2
surfaceScalarField phir(fvc::flux(UdmModel.Udm()))
void clear()
Clear the list, i.e. set size to zero.
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
PtrList< surfaceScalarField > phiFs(fluid.phiFs(rAUs))
tmp< volScalarField > byDt(const volScalarField &vf)
surfaceScalarField & phi2
const uniformDimensionedVectorField & g
const surfaceScalarField alpharAUf2(fvc::interpolate(max(alpha2, phase2.residualAlpha()) *rAU2))
const volScalarField & psi1
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
const surfaceScalarField phiKd1(fvc::interpolate(rAUKd1))
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
PtrList< volScalarField > rAUs
surfaceScalarField ghSnGradRho(ghf *fvc::snGrad(rho) *mesh.magSf())
setSnGrad< fixedFluxPressureFvPatchScalarField >(p_rgh.boundaryFieldRef(),(phiHbyA.boundaryField() - MRF.relative(phib))/(mesh.magSf().boundaryField() *rAUf.boundaryField()))
const volScalarField & gh
volScalarField p_rgh_0(p_rgh)
const surfaceScalarField alphaf2("alphaf2", scalar(1) - alphaf1)
const volScalarField rAUKd1(rAU1 *Kd)
const surfaceScalarField alpharAUf1(fvc::interpolate(max(alpha1, phase1.residualAlpha()) *rAU1))
const volScalarField rAUKd2(rAU2 *Kd)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const surfaceScalarField & alphaPhi1
const volScalarField & rAU2
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
const surfaceScalarField phiKd2(fvc::interpolate(rAUKd2))
const volScalarField & alpha1