applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H

   1 {
   2     rho = thermo.rho();
   3
   4     // Thermodynamic density needs to be updated by psi*d(p) after the
   5     // pressure solution - done in 2 parts. Part 1:
   6     thermo.rho() -= psi*p_rgh;
   7
   8     volScalarField rAU(1.0/UEqn.A());
   9     surfaceScalarField rhorAUf("(rho*(1|A(U)))", fvc::interpolate(rho*rAU));
  10
  11     U = rAU*UEqn.H();
  12
  13     phi = fvc::interpolate(rho)*
  14     (
  15         (fvc::interpolate(U) & mesh.Sf())
  16       + fvc::ddtPhiCorr(rAU, rho, U, phi)
  17     );
  18
  19     surfaceScalarField buoyancyPhi(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf());
  20     phi += buoyancyPhi;
  21
  22     fvScalarMatrix p_rghDDtEqn
  23     (
  24         fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh))
  25       + fvc::div(phi)
  26     );
  27
  28     for (int nonOrth=0; nonOrth<=pimple.nNonOrthCorr(); nonOrth++)
  29     {
  30         fvScalarMatrix p_rghEqn
  31         (
  32             p_rghDDtEqn
  33           - fvm::laplacian(rhorAUf, p_rgh)
  34         );
  35
  36         p_rghEqn.solve
  37         (
  38             mesh.solver(p_rgh.select(pimple.finalInnerIter(corr, nonOrth)))
  39         );
  40
  41         if (nonOrth == pimple.nNonOrthCorr())
  42         {
  43             // Calculate the conservative fluxes
  44             phi += p_rghEqn.flux();
  45
  46             // Explicitly relax pressure for momentum corrector
  47             p_rgh.relax();
  48
  49             // Correct the momentum source with the pressure gradient flux
  50             // calculated from the relaxed pressure
  51             U += rAU*fvc::reconstruct((buoyancyPhi + p_rghEqn.flux())/rhorAUf);
  52             U.correctBoundaryConditions();
  53         }
  54     }
  55
  56     p = p_rgh + rho*gh;
  57
  58     // Second part of thermodynamic density update
  59     thermo.rho() += psi*p_rgh;
  60
  61     DpDt = fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);
  62
  63     #include "rhoEqn.H"
  64     #include "compressibleContinuityErrs.H"
  65 }