applications/solvers/multiphase/interFoam/gammaEqnSubCycle.H

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