applications/solvers/multiphase/settlingFoam/pmhEqn.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 pmhEqn
  23     (
  24         fvm::laplacian(rUAf, pmh) == fvc::ddt(rho) + fvc::div(phi)
  25     );
  26
  27     pmhEqn.setReference(pmhRefCell, pmhRefValue);
  28     pmhEqn.solve();
  29
  30     if (nonOrth == nNonOrthCorr)
  31     {
  32         phi -= pmhEqn.flux();
  33     }
  34 }
  35
  36 p == pmh + rho*gh;
  37
  38 if (pmh.needReference())
  39 {
  40     p += dimensionedScalar
  41     (
  42         "p",
  43         p.dimensions(),
  44         pRefValue - getRefCellValue(p, pmhRefCell)
  45     );
  46 }
  47
  48 #include "rhoEqn.H"
  49 #include "compressibleContinuityErrs.H"
  50
  51 U += rUA*fvc::reconstruct((phi - phiU)/rUAf);
  52 U.correctBoundaryConditions();