applications/solvers/compressible/rhoSimpleFoam/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     eqnResidual = solve
  12     (
  13         UEqn() == -fvc::grad(p)
  14     ).initialResidual();
  15
  16     maxResidual = max(eqnResidual, maxResidual);