applications/solvers/multiphase/twoLiquidMixingFoam/gammaEqn.H

   1 {
   2     fvScalarMatrix gammaEqn
   3     (
   4         fvm::ddt(gamma)
   5       + fvm::div(phi, gamma)
   6       - fvm::laplacian(Dab, gamma)
   7     );
   8
   9     gammaEqn.solve();
  10
  11     rhoPhi = gammaEqn.flux()*(rho1 - rho2) + phi*rho2;
  12     rho = gamma*rho1 + (scalar(1) - gamma)*rho2;
  13
  14     Info<< "Phase 1 volume fraction = "
  15         << gamma.weightedAverage(mesh.V()).value()
  16         << "  Min(gamma) = " << min(gamma).value()
  17         << "  Max(gamma) = " << max(gamma).value()
  18         << endl;
  19 }