applications/solvers/lagrangian/coalChemistryFoam/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         DimensionedField<scalar, volMesh> tk =
  14             Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());
  15         DimensionedField<scalar, volMesh> tc =
  16             chemistry.tc()().dimensionedInternalField();
  17
  18         // Chalmers PaSR model
  19         kappa = (runTime.deltaT() + tc)/(runTime.deltaT() + tc + tk);
  20     }
  21     else
  22     {
  23         kappa = 1.0;
  24     }
  25
  26     chemistrySh = kappa*chemistry.Sh()();
  27 }