applications/solvers/multiphase/settlingFoam/UEqn.H

   1     // Solve the Momentum equation
   2
   3     fvVectorMatrix UEqn
   4     (
   5         fvm::ddt(rho, U)
   6       + fvm::div(phi, U)
   7       + fvc::div
   8         (
   9             (Alpha/(scalar(1.001) - Alpha))*(sqr(rhoc)/rho)*Vdj*Vdj,
  10             "div(phiVdj,Vdj)"
  11         )
  12       - fvm::laplacian(mu, U, "laplacian(muEff,U)")
  13     );
  14
  15     UEqn.relax();
  16
  17     if (momentumPredictor)
  18     {
  19         solve
  20         (
  21             UEqn
  22           ==
  23             fvc::reconstruct
  24             (
  25                 fvc::interpolate(rho)*(g & mesh.Sf())
  26               - fvc::snGrad(p)*mesh.magSf()
  27             )
  28         );
  29     }