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 if (oCorr == nOuterCorr - 1)
  11 {
  12     if (mesh.solutionDict().relax("UFinal"))
  13     {
  14         UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal"));
  15     }
  16     else
  17     {
  18         UEqn().relax(1);
  19     }
  20 }
  21 else
  22 {
  23     UEqn().relax();
  24 }
  25
  26 volScalarField rUA = 1.0/UEqn().A();
  27
  28 if (momentumPredictor)
  29 {
  30     if (oCorr == nOuterCorr-1)
  31     {
  32         solve(UEqn() == -fvc::grad(p), mesh.solutionDict().solver("UFinal"));
  33     }
  34     else
  35     {
  36         solve(UEqn() == -fvc::grad(p));
  37     }
  38 }
  39 else
  40 {
  41     U = rUA*(UEqn().H() - fvc::grad(p));
  42     U.correctBoundaryConditions();
  43 }