applications/solvers/incompressible/simpleFoam/pEqn.H

   1 {
   2     p.boundaryField().updateCoeffs();
   3
   4     volScalarField rAU(1.0/UEqn().A());
   5     U = rAU*UEqn().H();
   6     UEqn.clear();
   7
   8     phi = fvc::interpolate(U, "interpolate(HbyA)") & mesh.Sf();
   9     adjustPhi(phi, U, p);
  10
  11     // Non-orthogonal pressure corrector loop
  12     for (int nonOrth=0; nonOrth<=simple.nNonOrthCorr(); nonOrth++)
  13     {
  14         fvScalarMatrix pEqn
  15         (
  16             fvm::laplacian(rAU, p) == fvc::div(phi)
  17         );
  18
  19         pEqn.setReference(pRefCell, pRefValue);
  20
  21         pEqn.solve();
  22
  23         if (nonOrth == simple.nNonOrthCorr())
  24         {
  25             phi -= pEqn.flux();
  26         }
  27     }
  28
  29     #include "continuityErrs.H"
  30
  31     // Explicitly relax pressure for momentum corrector
  32     p.relax();
  33
  34     // Momentum corrector
  35     U -= rAU*fvc::grad(p);
  36     U.correctBoundaryConditions();
  37 }