applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H

   1 // Solve the Momentum equation
   2
   3 tmp<fvVectorMatrix> UEqn
   4 (
   5     pZones.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 mrfZones.addCoriolis(rho, UEqn());
  27 pZones.addResistance(UEqn());
  28
  29 volScalarField rUA = 1.0/UEqn().A();
  30
  31 if (momentumPredictor)
  32 {
  33     if (oCorr == nOuterCorr-1)
  34     {
  35         solve(UEqn() == -fvc::grad(p), mesh.solutionDict().solver("UFinal"));
  36     }
  37     else
  38     {
  39         solve(UEqn() == -fvc::grad(p));
  40     }
  41 }
  42 else
  43 {
  44     U = rUA*(UEqn().H() - fvc::grad(p));
  45     U.correctBoundaryConditions();
  46 }