9 word alphaScheme(
"div(phi,alpha)");
10 word alpharScheme(
"div(phirb,alpha)");
12 surfaceScalarField phir(phic*interface.nHatf());
14 for (
int gCorr=0; gCorr<nAlphaCorr; gCorr++)
16 volScalarField::Internal Sp
25 dimensionedScalar(
"Sp", dgdt.dimensions(), 0.0)
28 volScalarField::Internal Su
38 divU*min(alpha1, scalar(1))
43 if (dgdt[celli] > 0.0 && alpha1[celli] > 0.0)
45 Sp[celli] -= dgdt[celli]*alpha1[celli];
46 Su[celli] += dgdt[celli]*alpha1[celli];
48 else if (dgdt[celli] < 0.0 && alpha1[celli] < 1.0)
50 Sp[celli] += dgdt[celli]*(1.0 - alpha1[celli]);
55 surfaceScalarField phiAlpha1
65 -fvc::flux(-phir, alpha2, alpharScheme),
83 surfaceScalarField rho1f(fvc::interpolate(rho_gas));
84 surfaceScalarField rho2f(fvc::interpolate(rho_foam));
85 rhoPhi = phiAlpha1*(rho1f - rho2f) + phi*rho2f;
86 alpha2 = scalar(1) - alpha1;
89 Info<<
"Liquid phase volume fraction = " 90 << alpha1.weightedAverage(mesh.V()).value()
91 <<
" Min(alpha1) = " << min(alpha1).value()
92 <<
" Min(alpha2) = " << min(alpha2).value()