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