applications/solvers/multiphase/settlingFoam/pEqn.H

   1 volScalarField rUA = 1.0/UEqn.A();
   2
   3 surfaceScalarField rUAf
   4 (
   5     "(rho*(1|A(U)))",
   6     fvc::interpolate(rho)*fvc::interpolate(rUA)
   7 );
   8
   9 U = rUA*UEqn.H();
  10 phi =
  11     fvc::interpolate(rho)
  12    *(
  13        (fvc::interpolate(U) & mesh.Sf())
  14      + fvc::ddtPhiCorr(rUA, rho, U, phi)
  15     );
  16
  17 surfaceScalarField phiU("phiU", phi);
  18 phi -= ghf*fvc::snGrad(rho)*rUAf*mesh.magSf();
  19
  20 for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
  21 {
  22     fvScalarMatrix pEqn
  23     (
  24         fvm::laplacian(rUAf, p) == fvc::ddt(rho) + fvc::div(phi)
  25     );
  26
  27     pEqn.setReference(pRefCell, pRefValue);
  28     pEqn.solve();
  29
  30     if (nonOrth == nNonOrthCorr)
  31     {
  32         phi -= pEqn.flux();
  33     }
  34 }
  35
  36 #include "rhoEqn.H"
  37 #include "compressibleContinuityErrs.H"
  38
  39 U += rUA*fvc::reconstruct((phi - phiU)/rUAf);
  40 U.correctBoundaryConditions();