applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H

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