applications/solvers/multiphase/interFoam/gammaEqn.H

   1 {
   2     word gammaScheme("div(phi,gamma)");
   3     word gammarScheme("div(phirb,gamma)");
   4
   5     surfaceScalarField phic = mag(phi/mesh.magSf());
   6     phic = min(interface.cGamma()*phic, max(phic));
   7     surfaceScalarField phir = phic*interface.nHatf();
   8
   9     for (int gCorr=0; gCorr<nGammaCorr; gCorr++)
  10     {
  11         surfaceScalarField phiGamma =
  12             fvc::flux
  13             (
  14                 phi,
  15                 gamma,
  16                 gammaScheme
  17             )
  18           + fvc::flux
  19             (
  20                 -fvc::flux(-phir, scalar(1) - gamma, gammarScheme),
  21                 gamma,
  22                 gammarScheme
  23             );
  24
  25         MULES::explicitSolve(gamma, phi, phiGamma, 1, 0);
  26
  27         rhoPhi = phiGamma*(rho1 - rho2) + phi*rho2;
  28     }
  29
  30     Info<< "Liquid phase volume fraction = "
  31         << gamma.weightedAverage(mesh.V()).value()
  32         << "  Min(gamma) = " << min(gamma).value()
  33         << "  Max(gamma) = " << max(gamma).value()
  34         << endl;
  35 }