applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/hEqn.H

   1 {
   2     fvScalarMatrix hEqn
   3     (
   4         fvm::ddt(rho, h)
   5       + fvm::div(phi, h)
   6       - fvm::laplacian(turb.alphaEff(), h)
   7      ==
   8         DpDt
   9       + rad.Sh(thermo)
  10     );
  11
  12     hEqn.relax();
  13     hEqn.solve(mesh.solver(h.select(finalIter)));
  14
  15     thermo.correct();
  16
  17     rad.correct();
  18
  19     Info<< "Min/max T:" << min(thermo.T()).value() << ' '
  20         << max(thermo.T()).value() << endl;
  21 }