applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/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     fvc::makeAbsolute(phi, U);
  10     volScalarField divU(fvc::div(phi));
  11     fvc::makeRelative(phi, U);
  12
  13     if (nAlphaSubCycles > 1)
  14     {
  15         dimensionedScalar totalDeltaT = runTime.deltaT();
  16         surfaceScalarField rhoPhiSum(0.0*rhoPhi);
  17
  18         for
  19         (
  20             subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
  21             !(++alphaSubCycle).end();
  22         )
  23         {
  24             #include "alphaEqns.H"
  25             rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
  26         }
  27
  28         rhoPhi = rhoPhiSum;
  29     }
  30     else
  31     {
  32         #include "alphaEqns.H"
  33     }
  34
  35     if (pimple.corr() == 0)
  36     {
  37         interface.correct();
  38     }
  39 }