414964 A Multi-Scale Modeling Study of Methyl Trans-3-Hexenoate Oxidation By HO2

Tuesday, November 10, 2015: 10:10 AM
355F (Salt Palace Convention Center)
Stefania Cagnina1, André Nicolle1, Theodorus de Bruin2, Yuri Georgievskii3 and Stephen J. Klippenstein3, (1)Engine and Vehicle Modeling Department, IFP Energies Nouvelles, Rueil-Malmaison Cedex, France, (2)Thermodynamics and Molecular Modeling Department, IFP Energies Nouvelles, Rueil-Malmaison Cedex, France, (3)Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, IL

Fatty acid methyl esters (FAME) represent a promising alternative to fossil middle distillates. The availability of predictive tools for the design of cutting-edge combustion technologies is strongly dependent upon a comprehensive understanding of ester oxidation kinetics. However, gas-phase oxidation of FAME at low temperatures (600-900 K) remains poorly understood. Due to the significant impact of chain length and unsaturation of esters on their low-temperature combustion and pollutant formation pathways [1], it is desirable to consider C5+ unsaturated esters as biodiesel surrogate compounds. Previous studies mainly focused on H-abstraction from the alkyl chain and the subsequent long-debated O2 addition [2], overlooking the impact of ester function on abstraction kinetics [3,4]. The reaction of esters with HO2 is expected to have a significant role in autoignition under low temperature conditions. This system has been widely studied either in the context of olefin oxidation [5-6] or as a product channel for alkylperoxy decomposition [7-8]. To our knowledge, no exhaustive kinetic modeling study on HO2 interaction with an unsaturated ester has been carried out.

The present study aims at unraveling the reaction mechanism involved in the interaction of a biodiesel surrogate, methyl trans-3-hexenoate, with hydroperoxy radicals using a bottom-up methodology, involving ab initio electronic structure calculations and a theoretical kinetics study. Using post-Hartree-Fock methods coupled with extended basis sets, an accurate potential energy surface was obtained for the major reaction channels identified. The corresponding kinetic rate constants were obtained by solving the master equation using an efficient eigenvalue based method [9]. The addition of HO2 to the carbonyl and oxygen of the ester and the impact of the ester function on the H abstraction kinetics were extensively studied. This study helps resolve the current controversy regarding the effect of stereoisomerism coupled with an ester moiety [10]. The obtained rate constants were finally implemented into several existing kinetic mechanisms, resulting in overall better agreement with experimental concentration profiles [11-12].

[1] Sarathy S. M.; Gaïl S.; Syed S. A.; Thomson M. J.; Dagaut P. Proc. Combust. Inst., 2007, 31, 1015-1022.

[2] Zádor J.; Taatjes C. A.; Fernandes R. X. Proc. Energ. Combust., 2011, 37, 371-421.   

[3] Metcalfe W. K.; Dooley S.; Curran H.J.; Simmie J.M.; El-Nahas A.M.; Navarro M.V. J. Phys. Chem. A, 2007, 111, 4001-4014.

[4] Glaude P. A.; Herbinet O.; Bax S.; Biet J.; Warth V.; Battin-Leclerc F. Combust. Flame, 2010, 157, 2035-2050. 

[5] Chen C-J.; Bozzelli J. W.; J. Phys. Chem. A, 2000, 104, 4997-5012.

[6] Villano S. M.; Carstensen H-H.; Dean A. M.; J. Phys. Chem. A, 2013, 117, 6458-6473.

[7] Miyoshi A. J. Phys. Chem. A, 2011, 115, 3301-3325.

[8] Villano S. M.; Huynh L. K.; Carstensen H-H.; Dean A. M.; J. Phys. Chem. A, 2012, 116, 5068- 5089.

[9] Georgievskii Y.; Miller J. A.; Burke M. P.; Klippenstein S. J. J. Phys. Chem. A, 2013, 117, 12146- 12154.

[10] Fridlyand A.; Goldsborough S. S.; Brezinsky K.; Merchant S.S.; Green, W. H. Proc. Combust. Inst., 2015, 35, 333-340.  

[11] Zhang K.; Togbé C.; Dayma G.; Dagaut P. Combust. Flame, 2014, 161, 818-825.

[12] Wagnon S. W.; Karwat D. M. A.; Wooldridge M. S.; Westbrook C. K. Energy Fuels, 2014, 28, 7227-7234.

Extended Abstract: File Not Uploaded