applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H

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