applications/solvers/multiphase/compressibleInterFoam/alphaEqnsSubCycle.H

   1 {
   2     label nAlphaCorr(readLabel(pimple.dict().lookup("nAlphaCorr")));
   3
   4     label nAlphaSubCycles(readLabel(pimple.dict().lookup("nAlphaSubCycles")));
   5
   6     surfaceScalarField phic(mag(phi/mesh.magSf()));
   7     phic = min(interface.cAlpha()*phic, max(phic));
   8
   9     volScalarField divU(fvc::div(phi));
  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 "alphaEqns.H"
  23             rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
  24         }
  25
  26         rhoPhi = rhoPhiSum;
  27     }
  28     else
  29     {
  30         #include "alphaEqns.H"
  31     }
  32
  33     if (pimple.corr() == 0)
  34     {
  35         interface.correct();
  36     }
  37 }