applications/solvers/heatTransfer/buoyantPisoFoam/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 (momentumPredictor)
  13     {
  14         solve
  15         (
  16             UEqn
  17          ==
  18             fvc::reconstruct
  19             (
  20                 fvc::interpolate(rho)*(g & mesh.Sf())
  21               - fvc::snGrad(p)*mesh.magSf()
  22             )
  23         );
  24     }