applications/solvers/lagrangian/LTSReactingParcelFoam/chemistry.H

   1 if (chemistry.chemistry())
   2 {
   3     Info<< "Solving chemistry" << endl;
   4
   5     // update reaction rates
   6     chemistry.calculate();
   7
   8     // turbulent time scale
   9     if (turbulentReaction)
  10     {
  11         typedef DimensionedField<scalar, volMesh> dsfType;
  12
  13         const dimensionedScalar e0("e0", sqr(dimLength)/pow3(dimTime), SMALL);
  14
  15         const dsfType tk =
  16             Cmix*sqrt(turbulence->muEff()/rho/(turbulence->epsilon() + e0));
  17
  18         const dsfType tc = chemistry.tc()().dimensionedInternalField();
  19
  20         kappa = tc/(tc + tk);
  21     }
  22     else
  23     {
  24         kappa = 1.0;
  25     }
  26
  27     chemistrySh = kappa*chemistry.Sh()();
  28 }