A wide range kinetic modelling study of laminar flame speeds of reference fuels and their mixtures A. Frassoldati, R. Grana, A. Cuoci, T. Faravelli, E. Ranzi Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico di Milano Italy Abstract Aim of this work is to further validate a general detailed mechanism of pyrolysis and oxidation for the high temperature combustion of a large variety of different fuels. The attention is first placed on the validation of the chemical mechanism in respect of the laminar flame speed of H 2 and H 2 /CO mixtures as well as the burning velocity of small hydrocarbon fuels such as methane and C2-C4 species. Alkanes, alkenes and alkynes are considered, together with butadiene. Reference fuels for gasoline and jet fuel surrogates are also analysed. More specifically, three hydrocarbon classes are considered: alkanes (n-heptane: nc 7 H 16 and iso-octane: ic 8 H 18, and higher n-alkanes: nc 1 H 22 and nc 12 H 26 ), cyclo-alkanes (cyclo-hexane) and aromatics (benzene: C 6 H 6, toluene: C 7 H 8 ). Finally, burning velocity of C1-C4 alcohols are also compared with experimental data. 1. Introduction Since several years now, new burners and improved experimental methods produced reliable measurements for the flame speeds [1-2]. Flame speed is one of the key properties which are very useful to validate and optimize chemical reaction mechanisms at high temperatures. Burning velocity of different fuels are compared with experimental measurements, moving from small hydrocarbons up to reference species of gasoline and diesel fuels. Alkanes, alkenes and alkynes are considered, together with butadiene. Flame speeds of cyclo-alkanes and aromatics as well as burning velocity of alcohol fuels are also included in this kinetic study. This work takes advantage of a similar and recent kinetic study of Blanquart et al. [3], used to validate the Stanford high temperature kinetic mechanism. According to the conclusion of Kumar and Sung [4], comparison of experimental and predicted results indicates some deficiencies of the existing mechanisms. Sensitivity analysis shows that the most sensitive reactions pertinent to laminar flame speed are very similar for different hydrocarbons of the same family. 2. Kinetic model The semi-detailed oxidation mechanism of hydrocarbon fuels up to C16 already discussed elsewhere [5] is here adopted. The overall kinetic scheme is based on hierarchical modularity and is made up of more than 1 reactions among 35 species. Thermochemical data for most species was obtained from the CHEMKIN thermodynamic database [6]. For those species for which thermodynamic data is not available in the literature, the group additive method was used to estimate these properties. The complete mechanism, with thermodynamic and transport properties, is available in CHEMKIN format [http://www.chem.polimi.it/creckmodeling]. 1
Processes and Technologies for a Sustainable Energy 3. Comparisons of laminar flame speeds and model predictions The hierarchical structure of the kinetic model of hydrocarbon pyrolysis and combustion suggests analysing laminar flame speeds moving from H 2 and syngas before the investigation of laminar flame speeds of heavier fuel species. Recently, Frassoldati et al. [7] already revised and discussed burning rates of H 2 and syngas mixtures. This work will then sequentially analyse the following systems: 3.1 Flame speed of small hydrocarbon fuels (Akanes, alkanes, alkynes and dienes) 3.2 Flame speed of Heavy hydrocarbon species: Linear and branched alkanes (n-heptane, n-decane, n-dodecane and iso-octane); Cyclo-alkanes (cyclo-hexane); Aromatics. 3.3 Flame speed of alcohols 3.1. Flame speed of small hydrocarbon fuels 3.1.1 Alkanes (CH 4, C 2 H 6, C 3 H 8, n-c 4 H 1, iso-c 4 H 1 ) Babushok and Tsang [8] observed strong the similarities amongst alkane/air flames and they also reported the relative sensitivity coefficients for the burning velocity of these flames. For small alkane flames, the kinetic model shows large positive coefficients for the following reactions: H + O 2 OH + O HCO + [M] CO + H + [M] CO+OH CO 2 +H CH 3 +OH CH 2 S+ H 2 O, while a large negative coefficient is shown by the reaction: H+OH+[M] H 2 O + [M]. The reaction CH 3 +H+[M] CH 4 +[M] indicates a negative coefficient of both k and k, mainly for methane flame. The reaction CH 3 +CH 3 C 2 H 5 +H promotes the reactivity of methane, while has a larger negative impact on ethane flame. Model predictions at stoichiometric conditions, slightly underestimate methane speed, while overestimate ethane burning velocity (see Fig. 1a [1,9-13]). Fig. 1b shows a good agreement between model predictions and experimental measurements of the burning velocities of n- and iso-butane [14, 15]. Model predictions confirm the highest reactivity of ethane, the lowest of methane and also the lower burning rate of branched species when compared to the normal ones. 5 45 4 35 3 25 2 15 1 5.4.6.8 1. 1.2 1.4 1.6 1.8 Vagelopoulos 1998 Halter et al. 25 methane Egolfopoulos 199 Bosschaart 24 ethane Vagelopoulos et al. 1994 Jomaas et al. 25 propane 45 4 35 3 25 2 15 1 5.6.8 1. 1.2 1.4 1.6 1.8 (normal) Hirasawa et al. 22 (normal) (iso) model nc 4 H 1 model ic 4 H 1 Fig. 1 Burning velocity of small alkanes. Fig. 1a Flame velocity of CH 4, C 2 H 6 and C 3 H 8 Fig. 1b Flame velocity of nc 4 H 1 and isoc 4 H 1 2
Ischia, June, 27-3 - 21 3.1.2 Alkenes (C 2 H4, C 3 H 6, 1-C 4 H 8, iso-c 4 H 8 ) 8 7 6 5 4 3 2 1 Hirasawa et al. 22 Jomaas et al. 25 ethylene Jomaas et al. 25 Propylene.4.6.8 1. 1.2 1.4 1.6 1.8 5 4 3 2 1.6.8 1. 1.2 1.4 1.6 1.8 n butene Davis and Law 1998 isobutene Davis and Law 1998 model nc 4 H 8 model ic 4 H 8 Fig. 2 Burning velocity of alkenes. Fig. 2a Flame velocity of C 2 H 4, C 3 H 6 Fig. 2b Flame velocity of nc 4 H 8 and isoc 4 H 8 Contrary to the Stanford model [3], acetylene isomerization into singlet vinylidene C 2 H 2 H 2 CC: and successive reactions with molecular oxygen are not significant in our kinetic model, neither in acetylene nor in ethylene flames. Fig. 2a shows the flame velocity of C 2 H 4 and C 3 H 6. The resonance stability of allyl radicals from propylene explains the lower burning rates of propylene in respect of ethylene. Sensitivity analysis reveals the importance of the reaction of vinyl radical with oxygen CH 2 CH+O 2 CH 2 CHO + O in ethylene flame. As already mentioned, the H and OH abstractions on propylene to form allyl radicals reduce the propylene flame velocity. Fig. 2b compares the flame velocity of nc 4 H 8 and isoc 4 H 8. Again the model properly reproduces the lower reactivity of branched species. While model predictions are very good in lean conditions, it seems necessary to improve the kinetics in order to increase the reactivity of n-butene in rich conditions, by favouring dehydrogenation reactions of methyl-allyl radicals in respect of the dealkylation ones. 3.1.3 Alkynes and Dienes (C 2 H 2, C 3 H 4, 1-C 4 H 6 ) 16 14 12 1 8 6 4 2 Egelfopoulos et al. 199 Jomaas et al. 25 acetylene Hirasawa et al. 22 Jomaas et al. 25 ethylene Egolfopoulos et al. 199 Bosschaart et al. 24 ethane.4.6.8 1. 1.2 1.4 1.6 1.8 2. 2.2 7 6 5 4 3 2 1.6.8 1. 1.2 1.4 1.6 1.8 (nc4h1) Hirasawa et al. 22 (nc4h1) model nc4h1 ndavis (n butene Davis and Law 1998) model nc4h8 Laskin et al. 2 (linear) Laskin et al. 2 (NonLinear) model C4H6 Fig. 3 Burning velocity of alkynes and butadiene. Fig. 3a Flame velocity of C 2 H 2, C 2 H 4, C 2 H 6 Fig. 3b Flame velocity of nc 4 H 1, nc 4 H 8 and 1,3C 4 H 6 Fig. 3a shows the flame velocity of C 2 H 2, C 2 H 4 and C 2 H 6. The model properly captures the different reactivity of these C2 species. Flame velocity of acetylene is more than the double of the ethylene one and three times the velocity of ethane flame. Similarly, Fig. 3b shows that the flame velocity of 1,3C 4 H 6 is higher than the one of nc 4 H 8. nbutane velocity is the lowest. The interest in 1,3-butadiene (C 4 H 6 ) flames is due to its role in the formation of butadienyl radical and successive aromatic rings. Model predictions well confirm that the reactivity of 3
Processes and Technologies for a Sustainable Energy unsaturated species is larger than the one of saturated species with the same carbon number. 3.2. Flame speed of Heavy hydrocarbon species 3.2.1 Primary reference fuels (n-heptane and iso-octane) Babushok and Tsang [8] already studied the burning velocity of n-heptane and they observed that sensitivity of burning velocity to reactions of C5-C7 species is relatively small, while the sensitive reactions are similar to those determining burning velocity of C1-C4 hydrocarbons. Fig. 4 shows that model predictions well agree with the experiments of [14,17], both for n- heptane and iso-octane, further confirming the lower velocities of branched species. 3.2.2 Cyclo-alkanes (cyclo-hexane) Fig. 5 shows the flame velocity of cyclo-hexane. This comparison of cyclo-alkanes flame velocity clearly confirms not only the very good predictions of the model but also the similar reactivity of cyclo and normal- alkanes. 5 5 4 3 2 1 Kumar et al. 27 model C 7 H 16 Kumar et al. 27 model C 8 H 18.6.8 1. 1.2 1.4 1.6 1.8 4 3 2 1.6.8 1. 1.2 1.4 1.6 model Fig. 4 Flame speed of n-heptane and iso-octane Fig. 5 Flame speed of cyclo-hexane 3.2.3 Heavy normal alkanes (n-decane and n-dodecane) Fig. 6 shows the comparisons of burning velocities of n-decane and n-dodecane at different initial temperatures [4]. Note that these burning rates of heavy alkanes are very close to the ones of n-heptane, at the same temperature conditions. 9 8 7 6 5 4 3 2 1 Decane Kumar and Sung 27 (36 K) Kumar and Sung 27 (4 K) Kumar and Sung 27 (47 K) model 36 K model 4 K model 47 K.6.7.8.9 1. 1.1 1.2 1.3 1.4 1.5 Laminar Flame Speed [cm/s] 9 75 6 45 3 15 Dodecane Kumar and Sung 27 (4 K) model 4 K Kumar and Sung 27 (47 K) model 47 K.6.7.8.9 1. 1.1 1.2 1.3 1.4 1.5 Fig. 6 Burning velocity of large normal-alkanes in air at different initial temperatures 3.2.4 Aromatics (benzene and toluene) Fig. 7 shows the preliminary comparisons with the benzene and toluene flame speeds at 1 and 3 atm and different temperatures [14,18]. Model predictions capture the lower reactivity of toluene, due to benzyl radicals, but the kinetic model requires a further refinement in both the 4
Ischia, June, 27-3 - 21 systems. 8 7 6 5 4 3 2 1.6.8 1. 1.2 1.4 1.6 1.8 Davis and Law 1998 (1 atm 298 K) model 1 atm 298 K Johnston and Farrell 25 (3 atm 45 K) model 3 atm 45 K 6 5 4 3 2 1.6.8 1. 1.2 1.4 1.6 Davis and Law 1998 (1 atm 298 K) model (1 atm 298 K) Johnston and Farrell 25 (3 atm 45 K) model (3 atm 45 K) Fig. 7 Burning velocity of aromatics at different initial temperatures and pressures 3.3. Flame speed of oxygenated hydrocarbon species The kinetic mechanism of alcohol fuels, from ethanol to butanol isomers, was recently discussed [19-2-21]. Fig 8 show the very good agreement for all the C1-C4 alcohols. Note that methanol is faster than the remaining alcohols. This is mainly due to the lower propensity of methanol to form methyl radicals. Once again, while ethanol, n-propanol and n-butanol show very similar reactivity, the branched iso-propanol has a lower flame speed. 4. Conclusions This kinetic study not only contribute to the validation of the kinetic study, but also confirms Davis and Law conclusions [14]. Laminar flame speeds of various fuel hydrocarbons, including alkanes, alkenes, alkynes, aromatics, and alcohols, compared over an extensive range of equivalence ratios, temperatures and pressures, show that the laminar flame speeds of normal alkanes are close in the whole field of operating conditions, except for the lower flame speed of CH 4 and the higher ones of C 2 H 6. For fuels with the same carbon number, unsaturated species have higher flame speeds: alkanes < alkenes < alkynes. Branched species show lower burning rates than the corresponding linear components. The flame speed of cyclo-hexane is similar to the one of n-heptane [14]. Fall-out of this kinetic analysis is further validation of the kinetic scheme at high temperature conditions. Future activities will further analyze real fuels mixtures and the effect of high pressure conditions. 7 6 5 4 3 2 1.6.8 1. 1.2 1.4 1.6 Veloo et al. 21 methanol model methanol Veloo et al. 21 ethanol model Ethanol Laminar Burning Velocity (cm/s) 6 5 4 3 2 1 n-c 4 H 9 OH i-c 3 H 7 OH n-c 3 H 7 OH.6.8 1. 1.2 1.4 1.6 1.8 Fig. 8 Burning velocity of alcohol fuels. Fig. 8a Flame velocity of methanol and ethanol. Fig. 8b Flame velocity of n- propanol, iso-propanol and n-butanol (experimental data [2]) 5
Processes and Technologies for a Sustainable Energy 5. References 1. Bosschaart, K.J., De Goey, L.P.H.,: Combustion and Flame, 136: 261 (24). 2. Veloo, P.S., Wang, Y.L., Egolfopoulos, F.N., Westbrook, C.K.: Combustion and Flame, Article in Press (21). 3. Blanquart, G., Pepiot-Desjardins, P., Pitsch, H.: Combustion and Flame, 156: 588 (29). 4. Kumar, K., Sung, C.J.: Combustion and Flame, 151: 29 (27). 5. Ranzi, E., Dente, M., Goldaniga, A., Bozzano, G., Faravelli, T.: Progress in Energy and Combustion Science 27: 99 (21) 6. Kee, R. J., Rupley, F., Miller, J. A.: Tech. Rep., Sandia National Laboratories, Livermore, California, 1989. 7. Frassoldati, A., Faravelli, T., Ranzi, E.: International Journal of Hydrogen Energy, 31: 231 (26). 8. Babushok, V.I., Tsang, W.: Journal of Propulsion and Power, 2: 43 (24). 9. Vagelopoulos, Christine M., Egolfopoulos, Fokion N.: Symposium (International) on Combustion, 1 : 513 (1998). 1. Halter, F., Chauveau, C., Djebaïli-Chaumeix, N., Gökalp, I.: Proceedings of the Combustion Institute, 3: 21 (25). 11. Egolfopoulos, F.N., Law, C.K.: Symposium (International) on Combustion, 23: 333 (1991). 12. Vagelopoulos, C.M., Egolfopoulos, F.N. Symposium (International) on Combustion: 25, 1317 (1994). 13. Jomaas, G., Zheng, X.L., Zhu, D.L., Law, C.K. Proceedings of the Combustion Institute, 3: 193 (25). 14. Davis, S.G., Law, C.K.: Combustion Science and Technology, 14: 427 (1998). 15. Hirasawa, T., Sung, C.J., Joshi, A., Yang, Z., Wang, H., Law, C.K.: International Symposium on Combustion Abstracts of Accepted Papers, p. 13 (22). 16. Laskin, A., Wang, H., Law C.K.: International Journal of Chemical Kinetics, 32: 589 (2). 17. Kumar, K., Freeh, J.E., Sung, C.J., Huang, Y.: Journal of Propulsion and Power, 23: 428 (27). 18. Johnston R.J., Farrell, J.T.: Proceedings of the Combustion Institute, 3: 217 (25). 19. Frassoldati, A., Cuoci, A., Faravelli, T., Niemann, U., Ranzi, E., Seiser, R., Seshadri, K.: Combustion and Flame, 157: 2 (21). 2. Grana, R., Frassoldati, A., Faravelli, T, Niemann, U., Ranzi, E., Seiser, R., Cattolica, R., Seshadri, K.: Combustion and Flame, in press (21). 21. Frassoldati, A., Cuoci, A., Faravelli, T., Ranzi, E.: Combustion Science and Technology, in press (21). 6