applications/solvers/compressible/sonicFoam/sonicDyMFoam/pEqn.H

   1 rho = thermo.rho();
   2
   3 volScalarField rAU(1.0/UEqn.A());
   4 U = UEqn.H()/UEqn.A();
   5
   6 surfaceScalarField phid
   7 (
   8     "phid",
   9     fvc::interpolate(psi)
  10    *(
  11         (fvc::interpolate(U) & mesh.Sf())
  12       - fvc::meshPhi(rho, U)
  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     phi = pEqn.flux();
  28 }
  29
  30 #include "compressibleContinuityErrs.H"
  31
  32 U -= rAU*fvc::grad(p);
  33 U.correctBoundaryConditions();