applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/UEqn.H

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