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 (pimple.momentumPredictor())
  18     {
  19         solve
  20         (
  21             UEqn
  22           ==
  23             fvc::reconstruct
  24             (
  25                 (
  26                   - ghf*fvc::snGrad(rho)
  27                   - fvc::snGrad(p_rgh)
  28                 )*mesh.magSf()
  29             )
  30         );
  31     }