Go to the documentation of this file.00001 {
00002 word alphaScheme("div(phi,alpha)");
00003 word alpharScheme("div(phirb,alpha)");
00004
00005 surfaceScalarField phic = mag(phi/mesh.magSf());
00006 phic = min(interface.cAlpha()*phic, max(phic));
00007 surfaceScalarField phir = phic*interface.nHatf();
00008
00009 for (int aCorr=0; aCorr<nAlphaCorr; aCorr++)
00010 {
00011 surfaceScalarField phiAlpha =
00012 fvc::flux
00013 (
00014 phi,
00015 alpha1,
00016 alphaScheme
00017 )
00018 + fvc::flux
00019 (
00020 -fvc::flux(-phir, scalar(1) - alpha1, alpharScheme),
00021 alpha1,
00022 alpharScheme
00023 );
00024
00025 MULES::explicitSolve(alpha1, phi, phiAlpha, 1, 0);
00026
00027 rhoPhi = phiAlpha*(rho1 - rho2) + phi*rho2;
00028 }
00029
00030 Info<< "Liquid phase volume fraction = "
00031 << alpha1.weightedAverage(mesh.V()).value()
00032 << " Min(alpha1) = " << min(alpha1).value()
00033 << " Max(alpha1) = " << max(alpha1).value()
00034 << endl;
00035 }
