applications/solvers/lagrangian/coalChemistryFoam/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       + coalParcels.SU()
   9       + limestoneParcels.SU()
  10     );
  11
  12     UEqn.relax();
  13
  14     if (momentumPredictor)
  15     {
  16         solve(UEqn == -fvc::grad(p));
  17     }