applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H

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