applications/solvers/combustion/fireFoam/UEqn.H

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