applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H

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