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     if (oCorr == nOuterCorr-1)
  12     {
  13         hEqn.relax();
  14         hEqn.solve(mesh.solver("hFinal"));
  15     }
  16     else
  17     {
  18         hEqn.relax();
  19         hEqn.solve();
  20     }
  21
  22     thermo.correct();
  23 }