applications/solvers/combustion/engineFoam/pEqn.H

   1 rho = thermo.rho();
   2
   3 volScalarField rAU(1.0/UEqn.A());
   4 U = rAU*UEqn.H();
   5
   6 if (pimple.transonic())
   7 {
   8     surfaceScalarField phid
   9     (
  10         "phid",
  11         fvc::interpolate(psi)
  12        *((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U))
  13     );
  14
  15     for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  16     {
  17         fvScalarMatrix pEqn
  18         (
  19             fvm::ddt(psi, p)
  20           + fvm::div(phid, p)
  21           - fvm::laplacian(rho*rAU, p)
  22         );
  23
  24         pEqn.solve
  25         (
  26             mesh.solver(p.select(pimple.finalInnerIter(corr, nonOrth)))
  27         );
  28
  29         if (nonOrth == pimple.nNonOrthCorr())
  30         {
  31             phi == pEqn.flux();
  32         }
  33     }
  34 }
  35 else
  36 {
  37     phi = fvc::interpolate(rho)
  38          *((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U));
  39
  40     for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  41     {
  42         fvScalarMatrix pEqn
  43         (
  44             fvm::ddt(psi, p)
  45           + fvc::div(phi)
  46           - fvm::laplacian(rho*rAU, p)
  47         );
  48
  49         pEqn.solve
  50         (
  51             mesh.solver(p.select(pimple.finalInnerIter(corr, nonOrth)))
  52         );
  53
  54         if (nonOrth == pimple.nNonOrthCorr())
  55         {
  56             phi += pEqn.flux();
  57         }
  58     }
  59 }
  60
  61 #include "rhoEqn.H"
  62 #include "compressibleContinuityErrs.H"
  63
  64 U -= rAU*fvc::grad(p);
  65 U.correctBoundaryConditions();
  66
  67 DpDt = fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);