applications/solvers/compressible/sonicFoam/pEqn.H

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