applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/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
  10 if (p.needReference())
  11 {
  12     fvc::makeRelative(phi, U);
  13     adjustPhi(phi, U, p);
  14     fvc::makeAbsolute(phi, U);
  15 }
  16
  17 for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  18 {
  19     fvScalarMatrix pEqn
  20     (
  21         fvm::laplacian(rAU, p) == fvc::div(phi)
  22     );
  23
  24     pEqn.setReference(pRefCell, pRefValue);
  25
  26     pEqn.solve
  27     (
  28         mesh.solver(p.select(pimple.finalInnerIter(corr, nonOrth)))
  29     );
  30
  31     if (nonOrth == pimple.nNonOrthCorr())
  32     {
  33         phi -= pEqn.flux();
  34     }
  35 }
  36
  37 #include "continuityErrs.H"
  38
  39 // Explicitly relax pressure for momentum corrector
  40 p.relax();
  41
  42 // Make the fluxes relative to the mesh motion
  43 fvc::makeRelative(phi, U);
  44
  45 U -= rAU*fvc::grad(p);
  46 U.correctBoundaryConditions();