applications/solvers/lagrangian/reactingParcelFilmFoam/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       + parcels.Sh(hs)
  10       + surfaceFilm.Sh()
  11       + radiation->Shs(thermo)
  12       + chemistrySh
  13     );
  14
  15     hsEqn.relax();
  16
  17     hsEqn.solve();
  18
  19     thermo.correct();
  20
  21     radiation->correct();
  22
  23     Info<< "min/max(T) = " << min(T).value() << ", " << max(T).value() << endl;
  24 }