2 word alphaScheme(
"div(phi,alpha)");
3 word alpharScheme(
"div(phirb,alpha)");
29 mesh.time().timeName(),
68 phiAlpha1 -= phiAlpha1BD;
105 phiAlpha2 -= phiAlpha2BD;
123 phiAlpha1 = phiAlpha1BD + lambda*phiAlpha1;
124 phiAlpha2 = phiAlpha2BD + lambda*phiAlpha2;
147 phiAlpha1*(
rho1 - rho3)
148 + (phiAlpha2 + alpha2Eqn.flux())*(
rho2 - rho3)
151 alpha3 = 1.0 - alpha1 - alpha2;
154 Info<<
"Air phase volume fraction = "
155 << alpha1.weightedAverage(mesh.V()).value()
156 <<
" Min(alpha1) = " <<
min(alpha1).value()
157 <<
" Max(alpha1) = " <<
max(alpha1).value()
160 Info<<
"Liquid phase volume fraction = "
161 << alpha2.weightedAverage(mesh.V()).value()
162 <<
" Min(alpha2) = " <<
min(alpha2).value()
163 <<
" Max(alpha2) = " <<
max(alpha2).value()