applications/solvers/lagrangian/coalChemistryFoam/YEqn.H

   1 tmp<fv::convectionScheme<scalar> > mvConvection
   2 (
   3     fv::convectionScheme<scalar>::New
   4     (
   5         mesh,
   6         fields,
   7         phi,
   8         mesh.schemesDict().divScheme("div(phi,Yi_h)")
   9     )
  10 );
  11
  12
  13 {
  14     label inertIndex = -1;
  15     volScalarField Yt("Yt", 0.0*Y[0]);
  16
  17     forAll(Y, i)
  18     {
  19         if (Y[i].name() != inertSpecie)
  20         {
  21             volScalarField& Yi = Y[i];
  22             solve
  23             (
  24                 fvm::ddt(rho, Yi)
  25               + mvConvection->fvmDiv(phi, Yi)
  26               - fvm::laplacian(turbulence->muEff(), Yi)
  27               ==
  28                 coalParcels.Srho(i)
  29               + kappa*chemistry.RR(i)().dimensionedInternalField()
  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 }