applications/solvers/combustion/XiFoam/ftEqn.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,ft_b_h_hu)")
   9     )
  10 );
  11
  12 if (composition.contains("ft"))
  13 {
  14     volScalarField& ft = composition.Y("ft");
  15
  16     solve
  17     (
  18         fvm::ddt(rho, ft)
  19       + mvConvection->fvmDiv(phi, ft)
  20       - fvm::laplacian(turbulence->alphaEff(), ft)
  21     );
  22 }