applications/solvers/combustion/dieselEngineFoam/hEqn.H

   1 {
   2     solve
   3     (
   4         fvm::ddt(rho, h)
   5       + mvConvection->fvmDiv(phi, h)
   6       - fvm::laplacian(turbulence->alphaEff(), h)
   7      ==
   8        DpDt
   9      + dieselSpray.heatTransferSource()
  10     );
  11
  12     thermo->correct();
  13
  14     forAll(dQ, i)
  15     {
  16         dQ[i] = 0.0;
  17     }
  18
  19     scalarField cp(dQ.size(), 0.0);
  20
  21     forAll(Y, i)
  22     {
  23         volScalarField RRi = chemistry.RR(i);
  24
  25         forAll(h, celli)
  26         {
  27             scalar Ti = T[celli];
  28             cp[celli] += Y[i][celli]*chemistry.specieThermo()[i].Cp(Ti);
  29             scalar hi = chemistry.specieThermo()[i].h(Ti);
  30             scalar RR = RRi[celli];
  31             dQ[celli] -= hi*RR;
  32         }
  33
  34     }
  35
  36     forAll(dQ, celli)
  37     {
  38         dQ[celli] /= cp[celli];
  39     }
  40 }