applications/solvers/multiphase/interPhaseChangeFoam/gammaEqnSubCycle.H

   1 surfaceScalarField rhoPhi
   2 (
   3     IOobject
   4     (
   5         "rhoPhi",
   6         runTime.timeName(),
   7         mesh
   8     ),
   9     mesh,
  10     dimensionedScalar("0", dimensionSet(1, 0, -1, 0, 0), 0)
  11 );
  12
  13 {
  14     label nGammaCorr
  15     (
  16         readLabel(piso.lookup("nGammaCorr"))
  17     );
  18
  19     label nGammaSubCycles
  20     (
  21         readLabel(piso.lookup("nGammaSubCycles"))
  22     );
  23
  24     surfaceScalarField phic = mag(phi/mesh.magSf());
  25     phic = min(interface.cGamma()*phic, max(phic));
  26
  27     volScalarField divU = fvc::div(phi);
  28
  29     dimensionedScalar totalDeltaT = runTime.deltaT();
  30
  31     if (nGammaSubCycles > 1)
  32     {
  33         for
  34         (
  35             subCycle<volScalarField> gammaSubCycle(gamma, nGammaSubCycles);
  36             !(++gammaSubCycle).end();
  37         )
  38         {
  39 #           include "gammaEqn.H"
  40         }
  41     }
  42     else
  43     {
  44 #       include "gammaEqn.H"
  45     }
  46
  47     if (nOuterCorr == 1)
  48     {
  49         interface.correct();
  50     }
  51
  52     rho == gamma*rho1 + (scalar(1) - gamma)*rho2;
  53 }