applications/solvers/incompressible/pimpleFoam/pEqn.H

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