applications/solvers/combustion/reactingFoam/YEqn.H

   1 tmp<fv::convectionScheme<scalar> > mvConvection
   2 (
   3     fv::convectionScheme<scalar>::New
   4     (
   5         mesh,
   6         fields,
   7         phi,
   8         mesh.divScheme("div(phi,Yi_h)")
   9     )
  10 );
  11
  12 {
  13     label inertIndex = -1;
  14     volScalarField Yt = 0.0*Y[0];
  15
  16     for (label i=0; i<Y.size(); i++)
  17     {
  18         if (Y[i].name() != inertSpecie)
  19         {
  20             volScalarField& Yi = Y[i];
  21
  22             solve
  23             (
  24                 fvm::ddt(rho, Yi)
  25               + mvConvection->fvmDiv(phi, Yi)
  26               - fvm::laplacian(turbulence->muEff(), Yi)
  27              ==
  28                 kappa*chemistry.RR(i),
  29                 mesh.solver("Yi")
  30             );
  31
  32             Yi.max(0.0);
  33             Yt += Yi;
  34         }
  35         else
  36         {
  37             inertIndex = i;
  38         }
  39     }
  40
  41     Y[inertIndex] = scalar(1) - Yt;
  42     Y[inertIndex].max(0.0);
  43 }