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 } 00019 00020 // ************************ vim: set sw=4 sts=4 et: ************************ //