applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H

   1 label nAlphaCorr(readLabel(pimple.dict().lookup("nAlphaCorr")));
   2
   3 label nAlphaSubCycles(readLabel(pimple.dict().lookup("nAlphaSubCycles")));
   4
   5 if (nAlphaSubCycles > 1)
   6 {
   7     dimensionedScalar totalDeltaT = runTime.deltaT();
   8     surfaceScalarField rhoPhiSum(0.0*rhoPhi);
   9
  10     for
  11     (
  12         subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
  13         !(++alphaSubCycle).end();
  14     )
  15     {
  16         #include "alphaEqn.H"
  17         rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
  18     }
  19
  20     rhoPhi = rhoPhiSum;
  21 }
  22 else
  23 {
  24     #include "alphaEqn.H"
  25 }
  26
  27 interface.correct();
  28
  29 rho == alpha1*rho1 + (scalar(1) - alpha1)*rho2;