2010 CEFRC Conference Chemical Kinetics of Combustion Processes Hai Wang B. Yang, J. Camacho, S. Lieb, S. Memarzadeh, S.-K. Gao and S. Koumlis University of Southern California
Benzene + O( 3 P) Products Overall rate coefficient extensively studied, but the products and branching ratios not well known. Theoretical challenges in dealing with spin-state crossing. Product distributions determined over 300 to 1000 K and 1 to 10 Torr (with Taatjes/Sandia). Cyclopentadiene was directly identified for the first time, in addition to phenol and phenoxy. ab initio calculations (with Krylov/ USC) and master equation/rrkm modeling for data extrapolation. Taatjes, C. A., Osborn, D. L., Selby, T. M., Meloni, G., Trevitt, A. J., Epifanivskii, E. Krylov, A. I., Sirjean, B., Dames, E., Wang, H. Products of the benzene + O(3P) reaction, Journal of Physical Chemistry A, 114, 3355-3370 (2010) Figure 1. Branching ratios observed as a function of temperature at 4 Torr for (a) m/z=94/93 and m/z=66/65 by magnetic sector mass spectrometer (filled symbols) and time-offlight mass spectrometer (open symbols), and (b) branching ratios determined by time-of-flight mass spectrometer. Lines are drawn to guide the eye. Error bars are 2σ-standard deviations.
Benzene + O( 3 P) Products 3 A, 1 A, 3 A, 1 A 20 within 8 kcal/mol 1 2 O H + O ( 3 P) 3 O +H+CO Relative Energy ( kcal / mol ) 0 20 40 60 80 6 O C 4 O ISC + H 5 O 7 8 O O + HCO 10 O 9 OH 11 + CO 100
Benzene + O( 3 P) Products Table 6. Branching ratios experimentally determined and computed for reaction R1. Branching ratio Reaction channel P (Torr) T (K) Experimental Computational C 6 H 5 OH (9) 4 700 0.58 ± 0.08 a 0.56 (0.65) d 4 800 0.47 ± 0.08 a 0.43 (0.50) d 4 900 0.33 ± 0.08 a 0.29 (0.37) d 10 800 0.59 ± 0.08 a 0.51 (0.61) d 10 900 0.41 ± 0.09 a 0.40 (0.49) d C 6 H 5 O (3) + H 4 700 0.18 ± 0.08 a 0.18 4 800 0.24 ± 0.10 a 0.21 4 900 0.33 ± 0.13 a 0.23 10 800 0.19 ± 0.09 a 0.21 10 900 0.28 ± 0.12 a 0.23 C 5 H 6 (11) + CO 4 700 0.22 ± 0.05 a 0.17 4 800 0.27 ± 0.06 a 0.29 4 900 0.33 ± 0.08 a 0.40 10 800 0.21 ± 0.05 a 0.18 10 900 0.27 ± 0.07 a 0.28
Kinetic Modeling Studies of Iso-butane and Iso-butene Combustion (work-in-progress) Chemical kinetic submodel of iso-butanol combustion. In collaboration with N. Hansen/Sandia through Bin Yang (roving postdoc). A total of 61 datasets considered, including two burner stabilized flames (30 Torr) probed by Synchrotron Photoionization Molecular Beam Mass Spectrometry (Hansen/Yang). Preliminary results show that USC Mech II behaves rather well against the data, but improvements are needed for the initial pyrolysis and oxidation of isobutane and isobutene.
Kinetic Modeling Studies of Iso-butane and Iso-butene Combustion (work-in-progress) Chemical kinetic submodel of iso-butanol combustion. In collaboration with N. Hansen/Sandia through Bin Yang. A total of 61 datasets considered, including two burner stabilized flames (30 Torr) probed by Synchrotron Photoionization Molecular Beam Mass Spectrometry. Preliminary results show that USC Mech II behaves rather well against the data, but improvements are needed for the initial pyrolysis and oxidation of isobutane and isobutene. Oehlschlaeger et al. (2004) Koshi and coworkers (2007)
Sooting Behaviors of n- and iso-butanol Flames (Work-in-Progress) Fuel-rich chemistry and sooting behaviors of butanol isomers are unknown. Soot formation characterized in burner-stabilized stagnation flames. Burner-stabilized stagnation flame approach carrier N 2 H p cooling assembly stagnation plate/sample probe (T s ) x u v T b to SMPS 1cm r 1.2 1.2 1.0 1.0 Temperature, T (K) 1500 1000 500 0.0 0.2 0.4 0.6 0.8 Distance, H (cm) 1.0 1.2 0.8 0.7 0.6 H p = 0.55 cm dn/dlogd p 10 10 10 8 3 5 10 20 50 Diameter, D p (nm) 0.8 0.7 0.6 H p = 0.55 cm Abid et al. (2009)
Sooting Behaviors of n- and iso-butanol Flames (Work-in-Progress) Fuel-rich chemistry and sooting behaviors of butanol isomers are unknown. Soot formation characterized in burner-stabilized stagnation flames. No significant difference between n-butanol and isobutanol in particle inception and soot mass growth. Soot formed show a strong and persistent nucleation behavior well into the post flame. In comparison to comparable ethylene flames, soot inception in butanol flames takes substantially longer reaction time, but once nucleation starts, surface reaction is as fast as the ethylene flames. dn/dlogd p (cm -3 ) 10 11 10 10 10 9 10 8 10 7 H p = 0.90 cm 10 100 H p = 1.10 cm H p = 1.20 cm 10 100 10 100 Diameter, D p (nm) H p = 1.40 cm 10 100 H p = 1.50 cm 10 100
PAH Precursor Chemistry (1) PAH formation is sensitive to a multitude of elementary reactions and local flame conditions. Accurate prediction of PAH formation may require a precision in main flame chemistry currently unavailable. PAH formation can be highly sensitive to fuel structures. Possibly a large number of pathways to PAHs have yet been considered. Spectral sensitivity of pyrene concentration 90 Torr burner stabilized C 2 H 2 /O 2 /Ar flame (Bockhorn), H = 0.55 cm H+O 2 =O+OH HO 2+H=OH+OH HO 2 +OH=O 2 +H 2 O HCO+O 2 =CO+HO 2 CH+H 2 =CH 2 +H CH 2 +O 2 =CO 2 +H+H CH 2 *+H 2 =CH 3 +H C 2 H 2 +O=CH 2 +CO C 2H 2 +OH=C 2 H+H 2 O C 2 H 2 +H(+M)=C 2 H 3 (+M) HCCO+H=CH 2 *+CO HCCO+CH 3 =C 2 H4+CO C 2H 3 +O 2 =C 2 H 3 O+O C 2 H 2 +CH 2 *=C 3 H 3 +H C 2 H 2 +CH 2 =C 3 H 3 +H C 3 H 3 +OH=C 3 H 2 +H 2 O C 3 H 3 +OH=C 2 H 3 +HCO C 3 H 3 +C 3 H 3 =>A 1 c-c 6H 4+H=A 1- C 4 H 2 +H=n-C4H 3 A 1 +H=A 1 -+H 2 A 1 +OH=A 1-+H 2O A 1 -+H(+M)=A 1(+M) n-a 1 C 2 H 2 +C 2 H 2 =A 2 +H A 2 +OH=A 2-1+H 2 O A 2-1+H(+M)=A 2(+M) A 2 C 2 HA*+H(+M)=A 2 C 2 HA A 2 C 2 HA+OH=A 2C 2 HA*+H 2 A 2 C 2 HA*+C 2 H 2=A 3-4 P 2 +H=P 2 -+H 2 P 2 -+C 2 H 2 =A 3 +H A 3 +H=A 3-4+H 2 A 3-4+H(+M)=A 3(+M) A 3-4+C 2H 2=A 4 +H Main flame chemistry First aromatic ring Aromatic growth chemistry -0.5 0.0 0.5 1.0 Logarithmic Sensitivity Coefficient Wang & Frenklach, C&F (1997)
PAH Precursor Chemistry (2) Thermokinetic Origin of PAH Formation/Growth Beyond HACA G (kcal/mol-c) r 0.5 0.0 0.5 1.0 +H( H 2 ) +H (+M) +C 2 H 2 ( H) +H( H 2 ) +H (+M) +C 2 H 2 ( H) +C 2 H 2 +H (+M) +H( H 2 ) +H (+M) +C 2 H 2 +C 2 H 2 ( H) +H (+M) C 2 H 2 5 5 4 4 3 2 2 1 0 H 1 0 1 0 0 1 0 0 1 H 2 0 1 1 2 2 2 3 3 3 While HACA captures the thermokinetic requirements for PAH formation, its reversibility opens it to competitions from other pathways
Soot Nucleation (1) C B A C/H 1.2-2 ~2 ~10 Violi/D Anna Frenklach/Wang Miller Homann
Soot Nucleation (2) Binding Energy (ckal/mol) 100 10 circumcoronene ovalene coronene chrysene benzo[ghi]perylene pyrene anthracene & phenanthrene naphthalene 10 100 Number of C atoms 800 Herdman and Miller (2008) Binding energy of coronene = 25 kcal/mol Relative Concentration 10 0 200 10-1 10-2 10-3 10-4 600 400 Boiling/Sublimation Temperature (K) 10-5 1 2 3 4 5 6 7 Number of Aromatic Rings
Soot Nucleation (3) Is 25.4 kcal/mol enough to bind a pair of coronene together? 6 1 H E0+ 1 hν 4 i 1 i kt = + B exp hν i ( kt B ) 1 32 3 0 2 S 2B h p σ 1 ln 2 4 Ru m π ( ek ) σ BT 2 6 hν i kt B + = 1 ln 1 i hν i kt B e 1 Assumptions: v i = 200 cm -1 B (cm -1 ) = 1510 x MW -2.12 σ 2 = 1 2 coronene (coronene) 2 hν i kt B ( e ) K p (atm -1 ) Ovalene E 0 = 35 kcal/mol Circumcoronene E 0 = 63 kcal/mol Only circumcoronenes bind > 1600 K. G too positive to allow binding above 1000 K Entropy tears the dimer apart. 10 0 10-2 10-4 10-6 10-8 coronene ovalene circumcoronene (b) 0 500 1000 1500 2000 2500 T (K)
Soot Nucleation (1) C Possible, but not general C/H 1.2-2 Violi/D Anna B ~2 Frenklach/ Miller A ~10 Homann
Soot Nucleation (4) Polyacenes are singlet diradicals (though arguable). Ground-state polyacenes are close-shell singlets, but the adiabatic S0-T1 excitation energy is only 13 kcal/mol for heptacene - Hajgató et al. (2009). Applications in organic light emitting diodes and organic semiconductors and capacitors.
Soot Nucleation (5) Zigzag edges of graphene have localized π-electronic states Kobayashi 1993; Klein 1994 Zigzag edges have an open-shell singlet ground state e.g., Fujita et al. 1996; Nakada 1996 Finite-sized graphenes have radical or even multiradical characteristics. e.g., Nakano et al. 2008, Nagai 2010 Side chain can induce π-radical characteristics Nakano et al. 2007 Nonlinear optics applications. j = 1 2 3 zigzag edge i = 1 2 3 4 i j y i j y 1 1 0.000 1 3 0.037 2 1 0.050 2 3 0.217 3 1 0.149 3 3 0.510 4 1 0.281 4 3 0.806 Nagai et al. 2010 UBHandHLYP/6-31G(D) calculations 0 y 1 y = 0: close shell singlet y = 1: open shell singlet (diradical) Armchair edge
Soot Nucleation - Summary If PAHs with π-radicals do play a role in soot nucleation, we need to Understand the nature and structures of these PAH species, Determine their binding energies with relevant species, including aromatics, Probe them in flames (however small their concentrations may be), Account for the mechanism of their formation.
Work Summary Products of benzene + O( 3 P) reaction Kinetic modeling studies of iso-butane and iso-butene combustion Sooting behaviors of isobutanol and n-butanol flames Aromatic π-diradicals in soot nucleation and mass growth