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

   1     // Solve the Momentum equation
   2     tmp<fvVectorMatrix> UEqn
   3     (
   4         fvm::ddt(rho, U)
   5       + fvm::div(phi, U)
   6       + turb.divDevRhoReff(U)
   7     );
   8
   9     UEqn().relax();
  10
  11     if (momentumPredictor)
  12     {
  13          solve
  14          (
  15             UEqn()
  16          ==
  17             fvc::reconstruct
  18             (
  19                 (
  20                   - ghf*fvc::snGrad(rho)
  21                   - fvc::snGrad(p_rgh)
  22                 )*mesh.magSf()
  23             ),
  24             mesh.solver(U.select(finalIter))
  25         );
  26     }