applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/hEqn.H

   1 {
   2     fvScalarMatrix hEqn
   3     (
   4         fvm::div(phi, h)
   5       - fvm::Sp(fvc::div(phi), h)
   6       - fvm::laplacian(turbulence->alphaEff(), h)
   7      ==
   8         fvc::div(phi/fvc::interpolate(rho), rho/psi, "div(U,p)")
   9       - (rho/psi)*fvc::div(phi/fvc::interpolate(rho))
  10     );
  11
  12     pZones.addEnthalpySource(thermo, rho, hEqn);
  13
  14     hEqn.relax();
  15
  16     hEqn.solve();
  17
  18     thermo.correct();
  19 }