applications/solvers/multiphase/settlingFoam/alphaEqn.H

   1 {
   2     surfaceScalarField phiAlpha
   3     (
   4         IOobject
   5         (
   6             "phiAlpha",
   7             runTime.timeName(),
   8             mesh
   9         ),
  10         phi + rhoc*(mesh.Sf() & fvc::interpolate(Vdj))
  11     );
  12
  13     solve
  14     (
  15           fvm::ddt(rho, Alpha)
  16         + fvm::div(phiAlpha, Alpha)
  17         - fvm::laplacian(mut, Alpha)
  18     );
  19
  20     Info<< "Solid phase fraction = "
  21         << Alpha.weightedAverage(mesh.V()).value()
  22         << "  Min(Alpha) = " << min(Alpha).value()
  23         << "  Max(Alpha) = " << max(Alpha).value()
  24         << endl;
  25
  26     Alpha.min(1.0);
  27     Alpha.max(0.0);
  28
  29     rho == rhoc/(scalar(1) + (rhoc/rhod - 1.0)*Alpha);
  30     alpha == rho*Alpha/rhod;
  31 }