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