2 word alphaScheme(
"div(phi,alpha)");
3 word alpharScheme(
"div(phirb,alpha)");
9 volScalarField::DimensionedInternalField
Sp
21 volScalarField::DimensionedInternalField
Su
31 divU*
min(alpha1, scalar(1))
36 if (dgdt[celli] > 0.0 && alpha1[celli] > 0.0)
38 Sp[celli] -= dgdt[celli]*alpha1[celli];
39 Su[celli] += dgdt[celli]*alpha1[celli];
41 else if (dgdt[celli] < 0.0 && alpha1[celli] < 1.0)
43 Sp[celli] += dgdt[celli]*(1.0 - alpha1[celli]);
66 rhoPhi = phiAlpha1*(rho1f - rho2f) +
phi*rho2f;
68 alpha2 = scalar(1) - alpha1;
71 Info<<
"Liquid phase volume fraction = "
72 << alpha1.weightedAverage(
mesh.V()).value()
73 <<
" Min(alpha1) = " <<
min(alpha1).value()
74 <<
" Min(alpha2) = " <<
min(alpha2).value()