applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H

   1 {
   2     volScalarField rAU("rAU", 1.0/UEqn().A());
   3     surfaceScalarField rAUf("(1|A(U))", fvc::interpolate(rAU));
   4
   5     U = rAU*UEqn().H();
   6     UEqn.clear();
   7
   8     phi = fvc::interpolate(U) & mesh.Sf();
   9     adjustPhi(phi, U, p_rgh);
  10
  11     surfaceScalarField buoyancyPhi(rAUf*ghf*fvc::snGrad(rhok)*mesh.magSf());
  12     phi -= buoyancyPhi;
  13
  14     for (int nonOrth=0; nonOrth<=simple.nNonOrthCorr(); nonOrth++)
  15     {
  16         fvScalarMatrix p_rghEqn
  17         (
  18             fvm::laplacian(rAUf, p_rgh) == fvc::div(phi)
  19         );
  20
  21         p_rghEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));
  22
  23         p_rghEqn.solve();
  24
  25         if (nonOrth == simple.nNonOrthCorr())
  26         {
  27             // Calculate the conservative fluxes
  28             phi -= p_rghEqn.flux();
  29
  30             // Explicitly relax pressure for momentum corrector
  31             p_rgh.relax();
  32
  33             // Correct the momentum source with the pressure gradient flux
  34             // calculated from the relaxed pressure
  35             U -= rAU*fvc::reconstruct((buoyancyPhi + p_rghEqn.flux())/rAUf);
  36             U.correctBoundaryConditions();
  37         }
  38     }
  39
  40     #include "continuityErrs.H"
  41
  42     p = p_rgh + rhok*gh;
  43
  44     if (p_rgh.needReference())
  45     {
  46         p += dimensionedScalar
  47         (
  48             "p",
  49             p.dimensions(),
  50             pRefValue - getRefCellValue(p, pRefCell)
  51         );
  52         p_rgh = p - rhok*gh;
  53     }
  54 }