applications/solvers/multiphase/settlingFoam/pEqn.H

   1 volScalarField rAU(1.0/UEqn.A());
   2
   3 surfaceScalarField rAUf
   4 (
   5     "(rho*(1|A(U)))",
   6     fvc::interpolate(rho)*fvc::interpolate(rAU)
   7 );
   8
   9 U = rAU*UEqn.H();
  10 phi =
  11     fvc::interpolate(rho)
  12    *(
  13        (fvc::interpolate(U) & mesh.Sf())
  14      + fvc::ddtPhiCorr(rAU, rho, U, phi)
  15     );
  16
  17 surfaceScalarField phiU("phiU", phi);
  18 phi -= ghf*fvc::snGrad(rho)*rAUf*mesh.magSf();
  19
  20 for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  21 {
  22     fvScalarMatrix p_rghEqn
  23     (
  24         fvm::laplacian(rAUf, p_rgh) == fvc::ddt(rho) + fvc::div(phi)
  25     );
  26
  27     p_rghEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));
  28     p_rghEqn.solve();
  29
  30     if (nonOrth == pimple.nNonOrthCorr())
  31     {
  32         phi -= p_rghEqn.flux();
  33     }
  34 }
  35
  36 p == p_rgh + rho*gh;
  37
  38 if (p_rgh.needReference())
  39 {
  40     p += dimensionedScalar
  41     (
  42         "p",
  43         p.dimensions(),
  44         pRefValue - getRefCellValue(p, pRefCell)
  45     );
  46     p_rgh = p - rho*gh;
  47 }
  48
  49 #include "rhoEqn.H"
  50 #include "compressibleContinuityErrs.H"
  51
  52 U += rAU*fvc::reconstruct((phi - phiU)/rAUf);
  53 U.correctBoundaryConditions();