applications/solvers/combustion/XiFoam/huEqn.H

   1 if (ign.ignited())
   2 {
   3     solve
   4     (
   5         fvm::ddt(rho, hu)
   6       + mvConvection->fvmDiv(phi, hu)
   7       - fvm::laplacian(turbulence->alphaEff(), hu)
   8
   9     // These terms cannot be used in partially-premixed combustion due to
  10     // the resultant inconsistency between ft and hu transport.
  11     // A possible solution would be to solve for ftu as well as ft.
  12     //- fvm::div(muEff*fvc::grad(b)/(b + 0.001), hu)
  13     //+ fvm::Sp(fvc::div(muEff*fvc::grad(b)/(b + 0.001)), hu)
  14
  15      ==
  16         DpDt*rho/thermo.rhou()
  17     );
  18 }