applications/solvers/coupled/conjugateHeatFoam/UEqn.H

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