applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/UEqn.H

   1     // Solve the momentum equation
   2
   3     tmp<fvVectorMatrix> UEqn
   4     (
   5         fvm::div(phi, U)
   6       + turbulence->divDevReff(U)
   7     );
   8
   9     UEqn().relax();
  10
  11     eqnResidual = solve
  12     (
  13         UEqn()
  14       ==
  15         fvc::reconstruct
  16         (
  17             (
  18                 fvc::interpolate(rhok)*(g & mesh.Sf())
  19               - fvc::snGrad(p)*mesh.magSf()
  20             )
  21         )
  22     ).initialResidual();
  23
  24     maxResidual = max(eqnResidual, maxResidual);