applications/solvers/combustion/fireFoam/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)")
   9     )
  10 );
  11
  12 {
  13     fvScalarMatrix ftEqn
  14     (
  15         fvm::ddt(rho, ft)
  16       + mvConvection->fvmDiv(phi, ft)
  17       - fvm::laplacian(turbulence->alphaEff(), ft)
  18     );
  19
  20     ftEqn.relax();
  21     ftEqn.solve();
  22 }
  23
  24 Info<< "max(ft) = " << max(ft).value() << endl;
  25 Info<< "min(ft) = " << min(ft).value() << endl;