applications/solvers/compressible/rhoPimpleFoam/UEqn.H

   1 // Solve the Momentum equation
   2
   3 tmp<fvVectorMatrix> UEqn
   4 (
   5     fvm::ddt(rho, U)
   6   + fvm::div(phi, U)
   7   + turbulence->divDevRhoReff(U)
   8 );
   9
  10 UEqn().relax();
  11
  12 volScalarField rAU(1.0/UEqn().A());
  13
  14 if (pimple.momentumPredictor())
  15 {
  16     solve(UEqn() == -fvc::grad(p));
  17 }
  18 else
  19 {
  20     U = rAU*(UEqn().H() - fvc::grad(p));
  21     U.correctBoundaryConditions();
  22 }