applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H

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