![]() |
|
00001 forAll (ign.sites(), i) 00002 { 00003 const ignitionSite& ignSite = ign.sites()[i]; 00004 00005 if (ignSite.igniting()) 00006 { 00007 forAll(ignSite.cells(), icelli) 00008 { 00009 label ignCell = ignSite.cells()[icelli]; 00010 Info<< "Igniting cell " << ignCell; 00011 00012 Info<< " state :" 00013 << ' ' << b[ignCell] 00014 << ' ' << Xi[ignCell] 00015 << ' ' << Su[ignCell] 00016 << ' ' << mgb[ignCell] 00017 << endl; 00018 00019 bEqn.diag()[ignSite.cells()[icelli]] += 00020 ( 00021 ignSite.strength()*ignSite.cellVolumes()[icelli] 00022 *rhou[ignSite.cells()[icelli]]/ignSite.duration() 00023 )/(b[ignSite.cells()[icelli]] + 0.001); 00024 } 00025 } 00026 }