applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H

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