applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.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 nAlphaCorr
  15     (
  16         readLabel(piso.lookup("nAlphaCorr"))
  17     );
  18
  19     label nAlphaSubCycles
  20     (
  21         readLabel(piso.lookup("nAlphaSubCycles"))
  22     );
  23
  24     surfaceScalarField phic = mag(phi/mesh.magSf());
  25     phic = min(interface.cAlpha()*phic, max(phic));
  26
  27     volScalarField divU = fvc::div(phi);
  28
  29     dimensionedScalar totalDeltaT = runTime.deltaT();
  30
  31     if (nAlphaSubCycles > 1)
  32     {
  33         for
  34         (
  35             subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
  36             !(++alphaSubCycle).end();
  37         )
  38         {
  39 #           include "alphaEqn.H"
  40         }
  41     }
  42     else
  43     {
  44 #       include "alphaEqn.H"
  45     }
  46
  47     if (nOuterCorr == 1)
  48     {
  49         interface.correct();
  50     }
  51
  52     rho == alpha1*rho1 + (scalar(1) - alpha1)*rho2;
  53 }