applications/solvers/heatTransfer/buoyantPimpleFoam/UEqn.H

   1     // Solve the Momentum equation
   2
   3     fvVectorMatrix UEqn
   4     (
   5         fvm::ddt(rho, U)
   6       + fvm::div(phi, U)
   7       + turbulence->divDevRhoReff(U)
   8     );
   9
  10     UEqn.relax();
  11
  12     if (pimple.momentumPredictor())
  13     {
  14         solve
  15         (
  16             UEqn
  17          ==
  18             fvc::reconstruct
  19             (
  20                 (
  21                   - ghf*fvc::snGrad(rho)
  22                   - fvc::snGrad(p_rgh)
  23                 )*mesh.magSf()
  24             )
  25         );
  26     }