Baptiste Sirjean

Towards cleaner combustion engines through groundbreaking detailed chemical kinetic models

In the context of limiting the environmental impact of transportation, this critical review discusses new directions which are being followed in the development of more predictive and more accurate detailed chemical kinetic models for the combustion of fuels. In the first part, the performance of current models, especially in terms of the prediction of pollutant formation, is evaluated. In the next parts, recent methods and ways to improve these models are described. An emphasis is given on the development of detailed models based on elementary reactions, on the production of the related thermochemical and kinetic parameters, and on the experimental techniques available to produce the data necessary to evaluate model predictions under well defined conditions (212 references).

Combustion chemistry and flame structure of furan group biofuels using molecular-beam mass spectrometry and gas chromatography – Part I: Furan

Fuels of the furan family, i.e. furan itself, 2-methylfuran (MF), and 2,5-dimethylfuran (DMF) are being proposed as alternatives to hydrocarbon fuels and are potentially accessible from cellulosic biomass. While some experiments and modeling results are becoming available for each of these fuels, a comprehensive experimental and modeling analysis of the three fuels under the same conditions, simulated using the same chemical reaction model, has - to the best of our knowledge - not been attempted before. The present series of three papers, detailing the results obtained in flat flames for each of the three fuels separately, reports experimental data and explores their combustion chemistry using kinetic modeling. The first part of this series focuses on the chemistry of low-pressure furan flames. Two laminar premixed low-pressure (20 and 40 mbar) flat argon-diluted (50%) flames of furan were studied at two equivalence ratios (φ=1.0 and 1.7) using an analytical combination of high-resolution electron-ionization molecular-beam mass spectrometry (EI-MBMS) in Bielefeld and gas chromatography (GC) in Nancy. The time-of-flight MBMS with its high mass resolution enables the detection of both stable and reactive species, while the gas chromatograph permits the separation of isomers. Mole fractions of reactants, products, and stable and radical intermediates were measured as a function of the distance to the burner. A single kinetic model was used to predict the flame structure of the three fuels: furan (in this paper), 2-methylfuran (in Part II), and 2,5-dimethylfuran (in Part III). A refined sub-mechanism for furan combustion, based on the work of Tian et al. [Combustion and Flame 158 (2011) 756-773] was developed which was then compared to the present experimental results. Overall, the agreement is encouraging. The main reaction pathways involved in furan combustion were delineated computing the rates of formation and consumption of all species. It is seen that the predominant furan consumption pathway is initiated by H-addition on the carbon atom neighboring the O-atom with acetylene as one of the dominant products.

Shock Tube and Chemical Kinetic Modeling Study of the Oxidation of 2,5-Dimethylfuran

A detailed kinetic model describing the oxidation of 2,5-dimethylfuran (DMF), a potential second-generation biofuel, is proposed. The kinetic model is based upon quantum chemical calculations for the initial DMF consumption reactions and important reactions of intermediates. The model is validated by comparison to new DMF shock tube ignition delay time measurements (over the temperature range 1300-1831 K and at nominal pressures of 1 and 4 bar) and the DMF pyrolysis speciation measurements of Lifshitz et al. [ J. Phys. Chem. A 1998 , 102 ( 52 ), 10655 - 10670 ]. Globally, modeling predictions are in good agreement with the considered experimental targets. In particular, ignition delay times are predicted well by the new model, with model-experiment deviations of at most a factor of 2, and DMF pyrolysis conversion is predicted well, to within experimental scatter of the Lifshitz et al. data. Additionally, comparisons of measured and model predicted pyrolysis speciation provides validation of theoretically calculated channels for the oxidation of DMF. Sensitivity and reaction flux analyses highlight important reactions as well as the primary reaction pathways responsible for the decomposition of DMF and formation and destruction of key intermediate and product species.

Theoretical Kinetic Study of Thermal Unimolecular Decomposition of Cyclic Alkyl Radicals

Whereas many studies have been reported on the reactions of aliphatic hydrocarbons, the chemistry of cyclic hydrocarbons has not been explored extensively. In the present work, a theoretical study of the gas-phase unimolecular decomposition of cyclic alkyl radicals was performed by means of quantum chemical calculations at the CBS-QB3 level of theory. Energy barriers and high-pressure-limit rate constants were calculated systematically. Thermochemical data were obtained from isodesmic reactions, and the contribution of hindered rotors was taken into account. Classical transition state theory was used to calculate rate constants. The effect of tunneling was taken into account in the case of CH bond breaking. Three-parameter Arrhenius expressions were derived in the temperature range of 500-2000 K at atmospheric pressure, and the CC and CH bond breaking reactions were studied for cyclic alkyl radicals with a ring size ranging from three to seven carbon atoms, with and without a lateral alkyl chain. For the ring-opening reactions, the results clearly show an increase of the activation energy as the pi bond is being formed in the ring (endo ring opening) in contrast to the cases in which the pi bond is formed on the side chain (exo ring opening). These results are supported by analyses of the electronic charge density that were performed with Atoms in Molecules (AIM) theory. For all cycloalkyl radicals considered, CH bond breaking exhibits larger activation energies than CC bond breaking, except for cyclopentyl for which the ring-opening and H-loss reactions are competitive over the range of temperatures studied. The theoretical results compare rather well with the experimental data available in the literature. Evans-Polanyi correlations for CC and CH beta-scissions in alkyl and cycloalkyl free radicals were derived. The results highlight two different types of behavior depending on the strain energy in the reactant.

Products of the benzene + O(3P) reaction.

The gas-phase reaction of benzene with O((3)P) is of considerable interest for modeling of aromatic oxidation, and also because there exist fundamental questions concerning the prominence of intersystem crossing in the reaction. While its overall rate constant has been studied extensively, there are still significant uncertainties in the product distribution. The reaction proceeds mainly through the addition of the O atom to benzene, forming an initial triplet diradical adduct, which can either dissociate to form the phenoxy radical and H atom or undergo intersystem crossing onto a singlet surface, followed by a multiplicity of internal isomerizations, leading to several possible reaction products. In this work, we examined the product branching ratios of the reaction between benzene and O((3)P) over the temperature range 300-1000 K and pressure range 1-10 Torr. The reactions were initiated by pulsed-laser photolysis of NO(2) in the presence of benzene and helium buffer in a slow-flow reactor, and reaction products were identified by using the multiplexed chemical kinetics photoionization mass spectrometer operating at the Advanced Light Source (ALS) of Lawrence Berkeley National Laboratory. Phenol and phenoxy radical were detected and quantified. Cyclopentadiene and cyclopentadienyl radical were directly identified for the first time. Finally, ab initio calculations and master equation/RRKM modeling were used to reproduce the experimental branching ratios, yielding pressure-dependent rate expressions for the reaction channels, including phenoxy + H, phenol, cyclopentadiene + CO, which are proposed for kinetic modeling of benzene oxidation.

Improvement of the modeling of the low-temperature oxidation of n-butane: study of the primary reactions.

This paper revisits the primary reactions involved in the oxidation of n-butane from low to intermediate temperatures (550-800 K) including the negative temperature coefficient (NTC) zone. A model that was automatically generated is used as a starting point and a large number of thermochemical and kinetic data are then re-estimated. The kinetic data of the isomerization of alkylperoxy radicals giving (•)QOOH radicals and the subsequent decomposition to give cyclic ethers has been calculated at the CBS-QB3 level of theory. The newly obtained model allows a satisfactory prediction of experimental data recently obtained in a jet-stirred reactor and in rapid compression machines. A considerable improvement of the prediction of the selectivity of cyclic ethers is especially obtained compared to previous models. Linear and global sensitivity analyses have been performed to better understand which reactions are of influence in the NTC zone.

Tunneling in hydrogen-transfer isomerization of n-alkyl radicals.

The role of quantum tunneling in hydrogen shift in linear heptyl radicals is explored using multidimensional, small-curvature tunneling method for the transmission coefficients and a potential energy surface computed at the CBS-QB3 level of theory. Several one-dimensional approximations (Wigner, Skodje and Truhlar, and Eckart methods) were compared to the multidimensional results. The Eckart method was found to be sufficiently accurate in comparison to the small-curvature tunneling results for a wide range of temperature, but this agreement is in fact fortuitous and caused by error cancellations. High-pressure limit rate constants were calculated using the transition state theory with treatment of hindered rotations and Eckart transmission coefficients for all hydrogen-transfer isomerizations in n-pentyl to n-octyl radicals. Rate constants are found in good agreement with experimental kinetic data available for n-pentyl and n-hexyl radicals. In the case of n-heptyl and n-octyl, our calculated rates agree well with limited experimentally derived data. Several conclusions made in the experimental studies of Tsang et al. (Tsang, W.; McGivern, W. S.; Manion, J. A. Proc. Combust. Inst. 2009, 32, 131-138) are confirmed theoretically: older low-temperature experimental data, characterized by small pre-exponential factors and activation energies, can be reconciled with high-temperature data by taking into account tunneling; at low temperatures, transmission coefficients are substantially larger for H-atom transfers through a five-membered ring transition state than those with six-membered rings; channels with transition ring structures involving greater than 8 atoms can be neglected because of entropic effects that inhibit such transitions. The set of computational kinetic rates were used to derive a general rate rule that explicitly accounts for tunneling. The rate rule is shown to reproduce closely the theoretical rate constants.

Theoretical study of the thermal decomposition of the 5-methyl-2-furanylmethyl radical.

The thermal decomposition of the 5-methyl-2-furanylmethyl radical (R(1)), the most important primary radical formed during the combustion and thermal decomposition of 2,5-dimethylfuran (a promising next-generation biofuel), was studied using CBS-QB3 calculations and master equation (ME)/RRKM modeling. Because very little information is available in the literature, the detailed potential energy surface (PES) was investigated thoroughly. Only the main pathways, having a kinetic influence on the decomposition of R(1), were retained in the final ME/RRKM model. Among all the channels studied, the ring-opening of the 5-methyl-2-furanylmethyl radical, followed by ring enlargement to form cyclohexadienone molecules is predicted to be the easiest decomposition channel of R(1). The C(6) cyclic species formed can undergo unimolecular reactions to yield phenol and to a lesser extent cyclopentadiene and CO. Our calculations predict that these species are important products formed during the pyrolysis of 2,5-dimethylfuran (DMF). Other channels involved in the decomposition of R(1) lead directly to the formation of linear and cyclic unsaturated C(5) species and constitute an additional source of cyclopentadiene and CO. High-pressure limit rate constants were computed as well as thermochemical properties for important species. ME/RRKM analysis was performed to probe the influence of pressure on the rate coefficients and pressure dependent rate coefficients were proposed for pressures and temperatures ranging, respectively, from 10(-2) bar to 10 bar and 1000 to 2000 K.

Study of the Low Temperature Oxidation of Propane

The low-temperature oxidation of propane was investigated using a jet-stirred reactor at atmospheric pressure and two methods of analysis: gas chromatography and synchrotron vacuum ultraviolet photoionization mass spectrometry (SVUV-PIMS) with direct sampling through a molecular jet. The second method allowed the identification of products, such as molecules with hydroperoxy functions, which are not stable enough to be detected by gas chromatography. Mole fractions of the reactants and reaction products were measured as a function of the temperature (530-730 K), with a particular attention to reaction products involved in the low temperature oxidation, such as cyclic ethers, aldehydes, alcohols, ketones, and hydroperoxides. A new model has been obtained from an automatically generated one, which was used as a starting point, with a large number of re-estimated thermochemical and kinetic data. The kinetic data of the most sensitive reactions, i.e., isomerizations of alkylperoxy radicals and the subsequent decompositions, have been calculated at the CBS-QB3 level of theory. The model allows a satisfactory prediction of the experimental data. A flow rate analysis has allowed highlighting the important reaction channels.

Theoretical kinetic study of the reactions of cycloalkylperoxy radicals.

Reactions of alkyl radicals with oxygen are key reactions in the low-temperature oxidation of hydrocarbons, but they have not been extensively studied yet in the case of cycloalkanes. Isomerizations of cycloalkylperoxy radicals and formation of cyclic ethers are especially important. In the present work, a theoretical study of the gas-phase reactions of cyclopentylperoxy and cyclohexylperoxy radicals has been carried out by means of quantum chemical calculations at the CBS-QB3 level. Computations on cyclopentylperoxy decomposition pathways are reported here for the first time. Thermochemical data have been obtained by means of isodesmic reactions, and the contribution of hindered rotors has been explicitly taken into account. Transition state theory has been used to calculate rate constants for all the elementary reactions. Three-parameter Arrhenius expressions have been derived in the temperature range 300-1000 K. Tunneling effects have been accounted for in the case of H-atom transfers. Our results compare well with experimental data and previous calculations available in the literature. In particular, the predicted rate constants for processes involving cyclohexylperoxy radicals, which have been introduced in a reaction mechanism scheme proposed before, exhibit excellent agreement with experiments at low and intermediate temperatures.