applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/pEqn.H

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