applications/solvers/heatTransfer/buoyantSimpleRadiationFoam/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)*fvc::interpolate(p))
   9       - p*fvc::div(phi/fvc::interpolate(rho))
  10       + radiation->Sh(thermo)
  11     );
  12
  13     hEqn.relax();
  14
  15     eqnResidual = hEqn.solve().initialResidual();
  16     maxResidual = max(eqnResidual, maxResidual);
  17
  18     thermo.correct();
  19
  20     radiation->correct();
  21 }