applications/solvers/compressible/rhoPimpleFoam/hEqn.H

   1 {
   2     fvScalarMatrix hEqn
   3     (
   4         fvm::ddt(rho, h)
   5       + fvm::div(phi, h)
   6       - fvm::laplacian(turbulence->alphaEff(), h)
   7      ==
   8         DpDt
   9     );
  10
  11     hEqn.relax();
  12     hEqn.solve();
  13
  14     thermo.correct();
  15 }