applications/solvers/combustion/reactingFoam/hsEqn.H

   1 {
   2     fvScalarMatrix hsEqn
   3     (
   4         fvm::ddt(rho, hs)
   5       + mvConvection->fvmDiv(phi, hs)
   6       - fvm::laplacian(turbulence->alphaEff(), hs)
   7 //      - fvm::laplacian(turbulence->muEff(), hs)  // unit lewis no.
   8      ==
   9         DpDt
  10       + chemistrySh
  11     );
  12
  13     hsEqn.relax();
  14     hsEqn.solve();
  15
  16     thermo.correct();
  17
  18     Info<< "T gas min/max   = " << min(T).value() << ", "
  19         << max(T).value() << endl;
  20 }