applications/solvers/lagrangian/coalChemistryFoam/hEqn.H

   1 {
   2     fvScalarMatrix hEqn
   3     (
   4         fvm::ddt(rho, h)
   5       + mvConvection->fvmDiv(phi, h)
   6       - fvm::laplacian(turbulence->alphaEff(), h)
   7      ==
   8         DpDt
   9      +  coalParcels.Sh()
  10      +  limestoneParcels.Sh()
  11      +  enthalpySource.Su()
  12      +  radiation->Sh(thermo)
  13     );
  14
  15     hEqn.relax();
  16
  17     hEqn.solve();
  18
  19     thermo.correct();
  20
  21     radiation->correct();
  22 }