FreeFOAM The Cross-Platform CFD Toolkit
ignite.H
Go to the documentation of this file.
1  forAll (ign.sites(), i)
2  {
3  const ignitionSite& ignSite = ign.sites()[i];
4 
5  if (ignSite.igniting())
6  {
7  forAll(ignSite.cells(), icelli)
8  {
9  label ignCell = ignSite.cells()[icelli];
10  Info<< "Igniting cell " << ignCell;
11 
12  Info<< " state :"
13  << ' ' << b[ignCell]
14  << ' ' << Xi[ignCell]
15  << ' ' << Su[ignCell]
16  << ' ' << mgb[ignCell]
17  << endl;
18 
19  bEqn.diag()[ignSite.cells()[icelli]] +=
20  (
21  ignSite.strength()*ignSite.cellVolumes()[icelli]
22  *rhou[ignSite.cells()[icelli]]/ignSite.duration()
23  )/(b[ignSite.cells()[icelli]] + 0.001);
24  }
25  }
26  }
27 
28 // ************************ vim: set sw=4 sts=4 et: ************************ //