8 #include "alphaControls.H" 10 surfaceScalarField phic(mag(phi/mesh.magSf()));
11 phic = min(interface.cAlpha()*phic, max(phic));
13 volScalarField divU(fvc::div(fvc::absolute(phi, U)));
15 if (nAlphaSubCycles > 1)
17 dimensionedScalar totalDeltaT = runTime.deltaT();
18 surfaceScalarField rhoPhiSum
27 dimensionedScalar(
"0", rhoPhi.dimensions(), 0)
32 subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
33 !(++alphaSubCycle).end();
37 rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
builds and solves the phase volume fraction equation.