applications/solvers/multiphase/interMixingFoam/alphaEqnsSubCycle.H

   1 const dictionary& pimpleDict = pimple.dict();
   2 label nAlphaCorr(readLabel(pimpleDict.lookup("nAlphaCorr")));
   3 label nAlphaSubCycles(readLabel(pimpleDict.lookup("nAlphaSubCycles")));
   4
   5 if (nAlphaSubCycles > 1)
   6 {
   7     surfaceScalarField rhoPhiSum(0.0*rhoPhi);
   8     dimensionedScalar totalDeltaT = runTime.deltaT();
   9
  10     for
  11     (
  12         subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
  13         !(++alphaSubCycle).end();
  14     )
  15     {
  16         #include "alphaEqns.H"
  17         rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
  18     }
  19
  20     rhoPhi = rhoPhiSum;
  21 }
  22 else
  23 {
  24     #include "alphaEqns.H"
  25 }
  26
  27 interface.correct();
  28
  29 {
  30     volScalarField rhoNew(alpha1*rho1 + alpha2*rho2 + alpha3*rho3);
  31
  32     //solve(fvm::ddt(rho) + fvc::div(rhoPhi));
  33     //Info<< "density error = "
  34     //    << max((mag(rho - rhoNew)/mag(rhoNew))().internalField()) << endl;
  35
  36     rho == rhoNew;
  37 }