applications/solvers/lagrangian/reactingParcelFoam/UEqn.H

   1     fvVectorMatrix UEqn
   2     (
   3         fvm::ddt(rho, U)
   4       + fvm::div(phi, U)
   5       + turbulence->divDevRhoReff(U)
   6      ==
   7         rho.dimensionedInternalField()*g
   8       + parcels.SU(U)
   9     );
  10
  11     UEqn.relax();
  12
  13     if (pimple.momentumPredictor())
  14     {
  15         solve(UEqn == -fvc::grad(p));
  16     }