applications/solvers/lagrangian/coalChemistryFoam/hsEqn.H

   1 {
   2     fvScalarMatrix hsEqn
   3     (
   4         fvm::ddt(rho, hs)
   5       + mvConvection->fvmDiv(phi, hs)
   6       - fvm::laplacian(turbulence->alphaEff(), hs)
   7      ==
   8         DpDt
   9      +  coalParcels.Sh()
  10      +  limestoneParcels.Sh()
  11      +  enthalpySource.Su()
  12      +  radiation->Shs(thermo)
  13      +  chemistrySh
  14     );
  15
  16     hsEqn.relax();
  17
  18     hsEqn.solve();
  19
  20     thermo.correct();
  21
  22     radiation->correct();
  23
  24     Info<< "T gas min/max   = " << min(T).value() << ", "
  25         << max(T).value() << endl;
  26 }