2 word scheme("div(phi,alpha)");
4 surfaceScalarField phir(phia - phib);
6 Info<< "Max Ur Courant Number = "
10 mesh.surfaceInterpolation::deltaCoeffs()*mag(phir)
16 for (int acorr=0; acorr<nAlphaCorr; acorr++)
18 fvScalarMatrix alphaEqn
21 + fvm::div(phi, alpha, scheme)
22 + fvm::div(-fvc::flux(-phir, beta, scheme), alpha, scheme)
28 fvScalarMatrix betaEqn
31 + fvm::div(phi, beta, scheme)
32 + fvm::div(-fvc::flux(phir, scalar(1) - beta, scheme), beta, scheme)
41 + sqr(scalar(1) - beta)
42 - sqr(scalar(1) - alpha)
46 beta = scalar(1) - alpha;
49 Info<< "Dispersed phase volume fraction = "
50 << alpha.weightedAverage(mesh.V()).value()
51 << " Min(alpha) = " << min(alpha).value()
52 << " Max(alpha) = " << max(alpha).value()
56 rho = alpha*rhoa + beta*rhob;