applications/solvers/lagrangian/reactingParcelFilmFoam/pEqn.H

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