applications/solvers/combustion/dieselEngineFoam/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
  14     label inertIndex = -1;
  15     volScalarField Yt = 0.0*Y[0];
  16
  17     for(label i=0; i<Y.size(); i++)
  18     {
  19         if (Y[i].name() != inertSpecie)
  20         {
  21             volScalarField& Yi = Y[i];
  22
  23             solve
  24             (
  25                 fvm::ddt(rho, Yi)
  26               + mvConvection->fvmDiv(phi, Yi)
  27               - fvm::laplacian(turbulence->muEff(), Yi)
  28               ==
  29                 dieselSpray.evaporationSource(i)
  30               + kappa*chemistry.RR(i),
  31                 mesh.solver("Yi")
  32             );
  33
  34             Yi.max(0.0);
  35             Yt += Yi;
  36         }
  37         else
  38         {
  39             inertIndex = i;
  40         }
  41     }
  42
  43     Y[inertIndex] = scalar(1) - Yt;
  44     Y[inertIndex].max(0.0);
  45
  46 }