applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H

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