applications/solvers/combustion/rhoReactingFoam/chemistry.H

   1 {
   2     Info << "Solving chemistry" << endl;
   3
   4     chemistry.solve
   5     (
   6         runTime.value() - runTime.deltaT().value(),
   7         runTime.deltaT().value()
   8     );
   9
  10     // turbulent time scale
  11     if (turbulentReaction)
  12     {
  13         volScalarField tk =
  14                 Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());
  15         volScalarField tc = chemistry.tc();
  16
  17         // Chalmers PaSR model
  18         kappa = (runTime.deltaT() + tc)/(runTime.deltaT() + tc + tk);
  19     }
  20     else
  21     {
  22         kappa = 1.0;
  23     }
  24 }