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

   1 {
   2     fvScalarMatrix hEqn
   3     (
   4         fvm::div(phi, h)
   5       - fvm::Sp(fvc::div(phi), h)
   6       - fvm::laplacian(turb.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       + rad.Sh(thermo)
  11     );
  12
  13     hEqn.relax();
  14
  15     hEqn.solve();
  16
  17     thermo.correct();
  18
  19     rad.correct();
  20
  21     Info<< "Min/max T:" << min(thermo.T()).value() << ' '
  22         << max(thermo.T()).value() << endl;
  23 }