huEqn.H

Go to the documentation of this file.
00001 if (ign.ignited())
00002 {
00003     solve
00004     (
00005         betav*fvm::ddt(rho, hu)
00006       + mvConvection->fvmDiv(phi, hu)
00007       - fvm::laplacian(Db, hu)
00008 
00009     // These terms cannot be used in partially-premixed combustion due to
00010     // the resultant inconsistency between ft and hu transport.
00011     // A possible solution would be to solve for ftu as well as ft.
00012     //- fvm::div(muEff*fvc::grad(b)/(b + 0.001), hu)
00013     //+ fvm::Sp(fvc::div(muEff*fvc::grad(b)/(b + 0.001)), hu)
00014 
00015      ==
00016         betav*DpDt*rho/thermo.rhou()
00017     );
00018 }