applications/solvers/combustion/PDRFoam/UEqn.H

   1     fvVectorMatrix UEqn
   2     (
   3         betav*fvm::ddt(rho, U)
   4       + fvm::div(phi, U)
   5       + turbulence->divDevRhoReff(U)
   6      ==
   7         betav*rho*g
   8     );
   9
  10     volSymmTensorField invA = inv(I*UEqn.A() + drag->Dcu());
  11
  12     if (momentumPredictor)
  13     {
  14         U = invA & (UEqn.H() - betav*fvc::grad(p));
  15         U.correctBoundaryConditions();
  16     }