applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H

   1     // Solve the Momentum equation
   2
   3     tmp<fvVectorMatrix> UEqn
   4     (
   5         fvm::div(phi, U)
   6       - fvm::Sp(fvc::div(phi), U)
   7       + turbulence->divDevRhoReff(U)
   8     );
   9
  10     UEqn().relax();
  11
  12     eqnResidual = solve
  13     (
  14         UEqn() == -fvc::grad(pd) - fvc::grad(rho)*gh
  15     ).initialResidual();
  16
  17     maxResidual = max(eqnResidual, maxResidual);
  18