High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method.

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔHf(0 K) values with 2σ uncertainties of 0.2-0.5 kcal mol-1. This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔHf(0 K), in which the reference species may be up to five heavy atoms. The high-level ANL0 and ANL1 reference ΔHf(0 K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol-1) and the propagation of those uncertainties when used in the calculation of ΔHf(0 K) for the 194 target species. We examine the sensitivity of the predicted ΔHf(0 K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

Hierarchical Study of the Reactions of Hydrogen Atoms with Alkenes: A Theoretical Study of the Reactions of Hydrogen Atoms with C2-C4 Alkenes.

The present study complements our previous studies of the reactions of hydrogen atoms with C5 alkene species including 1- and 2-pentene and the branched isomers (2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene), by studying the reactions of hydrogen atoms with C2-C4 alkenes (ethylene, propene, 1- and 2-butene, and isobutene). The aim of the current work is to develop a hierarchical set of rate constants for H atom addition reactions to C2-C5 alkenes, both linear and branched, which can be used in the development of chemical kinetic models. High-pressure limiting and pressure-dependent rate constants are calculated using the Rice-Ramsperger-Kassel-Marcus (RRKM) theory and a one-dimensional master equation (ME). Rate constant recommendations for H atom addition and abstraction reactions in addition to alkyl radical decomposition reactions are also proposed and provide a useful tool for use in mechanisms of larger alkenes for which calculations do not exist. Additionally, validation of our theoretical results with single-pulse shock-tube pyrolysis experiments is carried out. An improvement in species mole fraction predictions for alkene pyrolysis is observed, showing the relevance of the present study.

Theoretical Study of the Reaction of Hydrogen Atoms with Three Pentene Isomers: 2-Methyl-1-butene, 2-Methyl-2-butene, and 3-Methyl-1-butene.

This paper presents a comprehensive potential energy surface (PES) for hydrogen atom addition to and abstraction from 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene and the subsequent ß-scission and H atom transfer reactions. Thermochemical parameters for species on the C5H11 potential energy surface (PES) were calculated as a function of temperature (298-2000 K), using a series of isodesmic reactions to determine the formation enthalpies. High-pressure limiting and pressure-dependent rate constants were calculated using Rice-Ramsperger-Kassel-Marcus theory with a one-dimensional master equation. A number of studies have highlighted the fact that C5 intermediate species play a role in polyaromatic hydrocarbon formation and that a fuel's chemical structure can be key in understanding the intermediate species formed during fuel decomposition. Rate constant recommendations for both H atom addition to, and H-atom abstraction by H atoms from, linear and branched alkenes have subsequently been proposed by incorporating our earlier work on 1- and 2-pentene, and these can be used in mechanisms of larger alkenes for which calculations do not exist. The current set of rate constants for the reactions of H atoms with both linear and branched C5 alkenes, including their chemically activated pathways, are the first available in the literature of any reasonable fidelity for combustion modeling and are important for gasoline mechanisms. Validation of our theoretical results with pyrolysis experiments of 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene at 2 bar in a single pulse shock tube (SPST) were carried out, with satisfactory agreement observed.

Hindered rotor benchmarks for the transition states of free radical additions to unsaturated hydrocarbons.

The hindered internal rotors of 32 transition states (TSs) formed through four free radicals, namely methyl, vinyl, ethyl, methoxy (CH3, C2H3, C2H5, CH3) additions to acetylene, ethylene, allene, propyne, and propene (C2H2/C2H4/C3H4-a/C3H4-p/C3H6) are studied. To validate the uncertainties of rate constants that stem from the use of different electronic structure methods to treat hindered rotors, the rotations of the newly formed C-C and/or C-O rotors in the transition states are calculated using commonly used DFT methods (B3LYP, M06-2X, ωB97X-D and B2PLYP-D3 with two Pople basis sets (6-31+G(d,p), 6-311++G(d,p)) and cc-pVTZ). The hindrance potential energies V(χ) calculated using the M06-2X/6-311++G(d,p) method are benchmarked at the CCSD(T), CCSD(T)-F12, DLPNO-CCSD(T) levels of theory with cc-pVTZ-F12 and cc-pVXZ (X = T, Q) basis sets and are extrapolated to the complete basis set (CBS) limit. The DLPNO-CCSD(T)/CBS method is proven to reproduce the CCSD(T)/CBS energies within 0.5 kJ mol-1 and this method is selected as the benchmark for all of the rotors in this study. Rotational constants B(χ) are computed for each method based on the optimized geometries for the hindrance potential via the I(2,3) approximation. Thereafter, the V(χ) and B(χ) values are used to compute hindered internal rotation partition functions, QHR, as a function of temperature. The uncertainties in the V(χ), B(χ) and QHR calculations stem from the use of different DFT methods for the internal rotor treatment are discussed for these newly formed rotors. For rotors formed by + C2 alkenes/alkynes, the V(χ) and QHR values calculated using DFT methods are compared with the DLPNO-CCSD(T)/CBS results and analysed according to reaction types. Based on comparisons of the DFT methods with the benchmarking method, reliable DFT methods are recommended for the treatment of internal rotors for different reaction types considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark hindrance potentials which can be used to estimate uncertainties in theoretically derived rate constants arising from the choice of different electronic structure methods.

An ab Initio/Transition State Theory Study of the Reactions of C5H9 Species of Relevance to 1,3-Pentadiene, Part II: Pressure Dependent Rate Constants and Implications for Combustion Modeling.

The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C5H8/1,4-C5H8/cyC5H8 + H, C2H4 + C3H5-a, C3H6 + C2H3) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. C5H9 radicals can be produced via H atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive C5H9 potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work ( J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C5H8 molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified C5H9 PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified C5H9 PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for H atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement.

Ab Initio/Transition-State Theory Study of the Reactions of C5H9 Species of Relevance to 1,3-Pentadiene, Part I: Potential Energy Surfaces, Thermochemistry, and High-Pressure Limiting Rate Constants.

In this study, the reactions of C5H9 radicals are theoretically investigated, with a particular emphasis on hydrogen atom addition reactions to 1,3-pentadiene (C5H8) to form C5H9 radicals, although the subsequent isomerization and decomposition reactions of the C5H9 radicals are also of direct relevance to the radicals formed from the pyrolysis and oxidation of species including pentene and cyclopentane. Moreover, H-atom abstraction reactions by hydrogen atoms from 1,3-pentadiene are also investigated. The geometries and frequencies of 63 potential energy surface (PES) minima and 88 transition states are optimized at the ωB97XD/aug-cc-pVTZ level of theory. Spin-unrestricted open-shell single-point energies for all the species are calculated at the CCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). A one-dimensional hindered rotor treatment is employed for torsional modes, with the M06-2X/6-311++G(D,P) method used to compute the potential energy as a function of the dihedral angle. The high-pressure limiting rate constants and the thermochemical properties for C5 species are calculated using the Master Equation System Solver (MESS) with conventional transition-state theory and comparisons made with existing available literature data. A hydrogen atom can add to the terminal carbon atom of 1,3-pentadiene to form the 2,4-C5H9 radical and/or the internal carbon atoms to form 2,5-C5H9, 1,4-C5H9, and 1,3-C5H9 radicals. Among the four entrance channels for H atom addition reactions, the formation of 2,4-C5H9 and 1,3-C5H9 radicals is more exothermic in comparison to the other C5H9 isomers (2,5-C5H9, 1,4-C5H9) because of the resonantly stabilized allylic structure. Consequently, the formation of the former is generally dominant in terms of barrier heights. H atom addition reactions to 1,3-pentadiene are compared to available C3-C5 alkenes and dienes, with external addition calculated to be kinetically favored over internal addition. However, the correlation between heats of formation and energy barriers for H atom addition to 1,2-dienes is different from that for 1,3- and 1,4-dienes. Hydrogen atom addition and abstraction rate constants are also compared for 1,3-pentadiene, with addition found to be dominant. The subsequent unimolecular reactions on the C5H9 PES are found to be highly complex with reactions taking place on a multiple-well multiple-channel PES. For clarity, the reaction mechanism and kinetics of each C5H9 radical are discussed individually in terms of the computed enthalpy of the reaction and activation, the transition-state structure/reaction class, and also in terms of the combustion species for which the reactions are of potential importance. The reactions on the C5H9 PES are divided into three reaction classes (H-shift isomerization, cycloaddition, and ß-scission reactions), and the reactivity-structure-based estimation rules for energy barriers are derived for these three reaction classes and compared to literature results for alkyl radicals.

Theoretical, Experimental, and Modeling Study of the Reaction of Hydrogen Atoms with 1- and 2-Pentene.

Alkyl radicals are prominent in combustion chemistry as they are formed by hydrocarbon decomposition or from a radical attack on hydrocarbons. Accurate determinations of the thermochemistry and kinetics of their unimolecular isomerization and decomposition reactions and related addition reactions of alkenes are therefore important in simulating the combustion chemistry of virtually all hydrocarbon fuels. In this work, a comprehensive potential energy surface (PES) for H-atom addition to and abstraction from 1- and 2-pentene, and the subsequent C-C and C-H ß-scission reactions, and H-atom transfer reactions has been considered. Thermochemical values for the species on the C5H11 PES were calculated as a function of temperature (298-2000 K), with enthalpies of formation determined using a network of isodesmic reactions. High-pressure limiting and pressure-dependent rate constants were calculated using the Rice-Ramsperger-Kassel-Marcus theory coupled with a one-dimensional master equation. As a validation of our theoretical results, hydrogen atomic resonance absorption spectrometry experiments were performed on the H-atom addition and abstraction reactions of 1- and 2-pentene. By incorporating our calculations into a detailed chemical kinetic model (AramcoMech 3.0), excellent agreement with these experiments is observed. The theoretical results are further validated via a comprehensive series of simulations of literature data. Our a priori model is found to reproduce important absolute species concentrations and product ratios reported therein.

Chemical Kinetics of Hydrogen Atom Abstraction from Allylic Sites by 3O2; Implications for Combustion Modeling and Simulation.

Hydrogen atom abstraction from allylic C-H bonds by molecular oxygen plays a very important role in determining the reactivity of fuel molecules having allylic hydrogen atoms. Rate constants for hydrogen atom abstraction by molecular oxygen from molecules with allylic sites have been calculated. A series of molecules with primary, secondary, tertiary, and super secondary allylic hydrogen atoms of alkene, furan, and alkylbenzene families are taken into consideration. Those molecules include propene, 2-butene, isobutene, 2-methylfuran, and toluene containing the primary allylic hydrogen atom; 1-butene, 1-pentene, 2-ethylfuran, ethylbenzene, and n-propylbenzene containing the secondary allylic hydrogen atom; 3-methyl-1-butene, 2-isopropylfuran, and isopropylbenzene containing tertiary allylic hydrogen atom; and 1-4-pentadiene containing super allylic secondary hydrogen atoms. The M06-2X/6-311++G(d,p) level of theory was used to optimize the geometries of all of the reactants, transition states, products and also the hinder rotation treatments for lower frequency modes. The G4 level of theory was used to calculate the electronic single point energies for those species to determine the 0 K barriers to reaction. Conventional transition state theory with Eckart tunnelling corrections was used to calculate the rate constants. The comparison between our calculated rate constants with the available experimental results from the literature shows good agreement for the reactions of propene and isobutene with molecular oxygen. The rate constant for toluene with O2 is about an order magnitude slower than that experimentally derived from a comprehensive model proposed by Oehlschlaeger and coauthors. The results clearly indicate the need for a more detailed investigation of the combustion kinetics of toluene oxidation and its key pyrolysis and oxidation intermediates. Despite this, our computed barriers and rate constants retain an important internal consistency. Rate constants calculated in this work have also been used in predicting the reactivity of the target fuels of 1-butene, 2-butene, isobutene, 2-methylfuran, 2,5-dimethylfuran, and toluene, and the results show that the ignition delay times for those fuels have been increased by a factor of 1.5-3. This work provides a first systematic study of one of the key initiation reaction for compounds containing allylic hydrogen atoms.

Modeling Nitrogen Species as Pollutants: Thermochemical Influences.

To simulate emissions of nitrogen-containing compounds in practical combustion environments, it is necessary to have accurate values for their thermochemical parameters, as well as accurate kinetic values to describe the rates of their formation and decomposition. Significant disparity is observed in the literature for the former, and we therefore present herein high-accuracy ab initio gas-phase thermochemistry for 60 nitrogenous compounds, many of which are important in the formation and consumption chemistry of NOx species. Several quantum-chemical composite methods (CBS-APNO, G3, and G4) were utilized to derive enthalpies of formation via the atomization method. Entropies and heat capacities were calculated from traditional statistical thermodynamics, with oscillators treated as anharmonic based on ro-vibrational property analyses carried out at the B3LYP/cc-pVTZ level of theory. The use of quantum chemical methods, along with the treatments of anharmonicities and hindered rotors, ensures accurate enthalpy of formation, entropy, and heat capacity values across the temperature range 298.15-3000 K. The implications of these results for atmospheric and combustion modeling are discussed.

Benchmarking Compound Methods (CBS-QB3, CBS-APNO, G3, G4, W1BD) against the Active Thermochemical Tables: Formation Enthalpies of Radicals.

The 298.15 K formation enthalpies of 38 radicals with molecular formula CxHyOz have been computed via the atomization procedure using the five title methods. The computed formation enthalpies are then benchmarked against the values recommended in the Active Thermochemical Tables (ATcT). The accuracy of the methods have been interpreted in terms of descriptive statistics, including the mean-signed error, mean-unsigned error, maximum average deviation, 2σ uncertainties, and 2×root-mean-square-deviations (2RMSD). The results highlight the following rank order of accuracy for the methods studied G4 > G3 > W1BD > CBS-APNO > CBS-QB3. The findings of this work are also considered in light of a recent companion study, which took an identical approach to quantifying the accuracies of these methods for 48 closed-shell singlet CxHyOz compounds. A similar order of accuracies and precisions were observed therein: G3 > G4 > W1BD > CBS-APNO > CBS-QB3. Both studies highlight systematic biases/deviations from the ATcT for the methods investigated, which are discussed in some detail, with methods having clear tendencies to over- or underpredict the recommended formation enthalpies for radical and/or closed-shell CxHyOz compounds. We show that one can improve the accuracy of their computation, and simultaneously reduce the uncertainty, by taking unweighted average formation enthalpies from various combinations of methods used. The reader should note that the statistical analyses preceding these conclusions also highlight that these error cancellation effects are unique for closed-shell and radical species. By extension, these error-cancellation effects can be expected to be different for various homologous series and chemical functionalities and their closed- and open-shell subgroups. Hence, further benchmarking studies are advised for other homologous series, such that the scientists and engineers (e.g., combustion/atmospheric/astrochemical) who frequently use these methods can assign reasonable uncertainties to their computations, while simultaneously optimizing their computational costs. For CxHyOz compounds, a combination of the CBS-APNO/G3/G4 methods is shown to be quite powerful when the atomization method is employed and is capable of reproducing the ATcT to within "near-chemical-accuracy", with 2RMSD (≈95% confidence interval) values of 0.0 ± 4.34 kJ mol(-1) computed for CxHyOz radicals, 0.0 ± 4.22 kJ mol(-1) for closed-shell CxHyOz compounds, with a total uncertainty of 0.0 ± 4.27 kJ mol(-1) subsequently computed considering all 85 CxHyOz compounds. Given the performance of these methods for determination of formation enthalpies when the atomization procedure is employed, we expect isodesmic reactions involving these methods to be capable of achieving chemical accuracy, as illustrated for the case of the tert-butyl radical. We also highlight that there is still disagreement between experiment and theory for this radical, despite its significance in gas-phase chemistry. Kineticists, thermodynamicists, and chemical kinetic modellers alike are warned that the popular CBS-QB3 method is found to have particularly poor performance, with a computed 2RMSD of 0.0 ± 12.51 kJ mol(-1), indicating that one should not apply this method in isolation for formation enthalpy determination unless other error-cancellation strategies are employed.

Revisiting the Kinetics and Thermodynamics of the Low-Temperature Oxidation Pathways of Alkanes: A Case Study of the Three Pentane Isomers.

This paper describes our developing understanding of low-temperature oxidation kinetics. We have investigated the ignition of the three pentane isomers in a rapid compression machine over a wide range of temperatures and pressures, including conditions of negative temperature coefficient behavior. The pentane isomers are small alkanes, yet have structures that are complex enough to allow for the application of their kinetic and thermochemical rules to larger molecules. Updates to the thermochemistry of the species important in the low-temperature oxidation of hydrocarbons have been made based on a thorough literature review. An evaluation of recent quantum-chemically derived rate coefficients from the literature pertinent to important low-temperature oxidation reaction classes has been performed, and new rate rules are recommended for these classes. Several reaction classes have also been included to determine their importance with regard to simulation results, and we have found that they should be included when developing future chemical kinetic mechanisms. A comparison of the model simulations with pressure-time histories from experiments in a rapid compression machine shows very good agreement for both ignition delay time and pressure rise for both the first- and second-stage ignition events. We show that revisions to both the thermochemistry and the kinetics are required in order to replicate experiments well. A broader validation of the models with ignition delay times from shock tubes and a rapid compression machine is presented in an accompanying paper. The results of this study enhance our understanding of the combustion of straight- and branched-chained alkanes.

Benchmarking Compound Methods (CBS-QB3, CBS-APNO, G3, G4, W1BD) against the Active Thermochemical Tables: A Litmus Test for Cost-Effective Molecular Formation Enthalpies.

The theoretical atomization energies of some 45 CxHyOz molecules present in the Active Thermochemical Tables compilation and of particular interest to the combustion chemistry community have been computed using five composite model chemistries as titled. The species contain between 1-8 "heavy" atoms, and a few are conformationally diverse with up to nine conformers. The enthalpies of formation at 0 and 298.15 K are then derived via the atomization method and compared against the recommended values. In general, there is very good agreement between our averaged computed values and those in the ATcT; those for 1,3-cyclopentadiene exceptionally differ considerably, and we show from isodesmic reactions that the true value for 1,3-cyclopentadiene is closer to 134 kJ mol(-1) than the reported 101 kJ mol(-1). If one is restricted to using a single method, statistical measures indicate that the best methods are in the rank order G3 ≈ G4 > W1BD > CBS-APNO > CBS-QB3. The CBS-x methods do on average predict ΔfH(⊖)(298.15 K) within ≈5 kJ mol(-1) but are prone to occasional lapses. There are statistical advantages to be gained from using a number of methods in tandem, and all possible combinations have been tested. We find that the average formation enthalpy coming from using CBS-APNO/G4, CBS-APNO/G3, and G3/G4 show lower mean signed and mean unsigned errors, and lower standard and root-mean-squared deviations, than any of these methods in isolation. Combining these methods also leads to the added benefit of providing an uncertainty rooted in the chemical species under investigation. In general, CBS-APNO and W1BD tend to underestimate the formation enthalpies of target species, whereas CBS-QB3, G3, and G4 have a tendency to overestimate the same. Thus, combining CBS-APNO with a G3/G4 combination leads to an improvement in all statistical measures of accuracy and precision, predicting the ATcT values to within 0.14 ± 4.21 kJ mol(-1), thus rivalling "chemical accuracy" (±4.184 kJ mol(-1)) without the excessive cost associated with higher-level methods such as W1BD.

The pyrolysis of 2-methylfuran: a quantum chemical, statistical rate theory and kinetic modelling study.

Due to the rapidly growing interest in the use of biomass derived furanic compounds as potential platform chemicals and fossil fuel replacements, there is a simultaneous need to understand the pyrolysis and combustion properties of such molecules. To this end, the potential energy surfaces for the pyrolysis relevant reactions of the biofuel candidate 2-methylfuran have been characterized using quantum chemical methods (CBS-QB3, CBS-APNO and G3). Canonical transition state theory is employed to determine the high-pressure limiting kinetics, k(T), of elementary reactions. Rice-Ramsperger-Kassel-Marcus theory with an energy grained master equation is used to compute pressure-dependent rate constants, k(T,p), and product branching fractions for the multiple-well, multiple-channel reaction pathways which typify the pyrolysis reactions of the title species. The unimolecular decomposition of 2-methylfuran is shown to proceed via hydrogen atom transfer reactions through singlet carbene intermediates which readily undergo ring opening to form collisionally stabilised acyclic C5H6O isomers before further decomposition to C1-C4 species. Rate constants for abstraction by the hydrogen atom and methyl radical are reported, with abstraction from the alkyl side chain calculated to dominate. The fate of the primary abstraction product, 2-furanylmethyl radical, is shown to be thermal decomposition to the n-butadienyl radical and carbon monoxide through a series of ring opening and hydrogen atom transfer reactions. The dominant bimolecular products of hydrogen atom addition reactions are found to be furan and methyl radical, 1-butene-1-yl radical and carbon monoxide and vinyl ketene and methyl radical. A kinetic mechanism is assembled with computer simulations in good agreement with shock tube speciation profiles taken from the literature. The kinetic mechanism developed herein can be used in future chemical kinetic modelling studies on the pyrolysis and oxidation of 2-methylfuran, or the larger molecular structures for which it is a known pyrolysis/combustion intermediate (e.g. cellulose, coals, 2,5-dimethylfuran).

A comprehensive experimental and detailed chemical kinetic modelling study of 2,5-dimethylfuran pyrolysis and oxidation.

The pyrolytic and oxidative behaviour of the biofuel 2,5-dimethylfuran (25DMF) has been studied in a range of experimental facilities in order to investigate the relatively unexplored combustion chemistry of the title species and to provide combustor relevant experimental data. The pyrolysis of 25DMF has been re-investigated in a shock tube using the single-pulse method for mixtures of 3% 25DMF in argon, at temperatures from 1200-1350 K, pressures from 2-2.5 atm and residence times of approximately 2 ms. Ignition delay times for mixtures of 0.75% 25DMF in argon have been measured at atmospheric pressure, temperatures of 1350-1800 K at equivalence ratios (ϕ) of 0.5, 1.0 and 2.0 along with auto-ignition measurements for stoichiometric fuel in air mixtures of 25DMF at 20 and 80 bar, from 820-1210 K. This is supplemented with an oxidative speciation study of 25DMF in a jet-stirred reactor (JSR) from 770-1220 K, at 10.0 atm, residence times of 0.7 s and at ϕ = 0.5, 1.0 and 2.0. Laminar burning velocities for 25DMF-air mixtures have been measured using the heat-flux method at unburnt gas temperatures of 298 and 358 K, at atmospheric pressure from ϕ = 0.6-1.6. These laminar burning velocity measurements highlight inconsistencies in the current literature data and provide a validation target for kinetic mechanisms. A detailed chemical kinetic mechanism containing 2768 reactions and 545 species has been simultaneously developed to describe the combustion of 25DMF under the experimental conditions described above. Numerical modelling results based on the mechanism can accurately reproduce the majority of experimental data. At high temperatures, a hydrogen atom transfer reaction is found to be the dominant unimolecular decomposition pathway of 25DMF. The reactions of hydrogen atom with the fuel are also found to be important in predicting pyrolysis and ignition delay time experiments. Numerous proposals are made on the mechanism and kinetics of the previously unexplored intermediate temperature combustion pathways of 25DMF. Hydroxyl radical addition to the furan ring is highlighted as an important fuel consuming reaction, leading to the formation of methyl vinyl ketone and acetyl radical. The chemically activated recombination of HȮ2 or CH3Ȯ2 with the 5-methyl-2-furanylmethyl radical, forming a 5-methyl-2-furylmethanoxy radical and ȮH or CH3Ȯ radical is also found to exhibit significant control over ignition delay times, as well as being important reactions in the prediction of species profiles in a JSR. Kinetics for the abstraction of a hydrogen atom from the alkyl side-chain of the fuel by molecular oxygen and HȮ2 radical are found to be sensitive in the estimation of ignition delay times for fuel-air mixtures from temperatures of 820-1200 K. At intermediate temperatures, the resonantly stabilised 5-methyl-2-furanylmethyl radical is found to predominantly undergo bimolecular reactions, and as a result sub-mechanisms for 5-methyl-2-formylfuran and 5-methyl-2-ethylfuran, and their derivatives, have also been developed with consumption pathways proposed. This study is the first to attempt to simulate the combustion of these species in any detail, although future refinements are likely necessary. The current study illustrates both quantitatively and qualitatively the complex chemical behavior of what is a high potential biofuel. Whilst the current work is the most comprehensive study on the oxidation of 25DMF in the literature to date, the mechanism cannot accurately reproduce laminar burning velocity measurements over a suitable range of unburnt gas temperatures, pressures and equivalence ratios, although discrepancies in the experimental literature data are highlighted. Resolving this issue should remain a focus of future work.

A high temperature and atmospheric pressure experimental and detailed chemical kinetic modelling study of 2-methyl furan oxidation.

An experimental ignition delay time study for the promising biofuel 2-methyl furan (2MF) was performed at equivalence ratios of 0.5, 1.0 and 2.0 for mixtures of 1% fuel in argon in the temperature range 1200-1800 K at atmospheric pressure. Laminar burning velocities were determined using the heat-flux method for mixtures of 2MF in air at equivalence ratios of 0.55-1.65, initial temperatures of 298-398 K and atmospheric pressure. A detailed chemical kinetic mechanism consisting of 2059 reactions and 391 species has been constructed to describe the oxidation of 2MF and is used to simulate experiment. Accurate reproduction of the experimental data has been obtained over all conditions with the developed mechanism. Rate of production and sensitivity analyses have been carried out to identify important consumption pathways of the fuel and key kinetic parameters under these conditions. The reactions of hydrogen atom with the fuel are highlighted as important under all experimental conditions studied, with abstraction by the hydrogen atom promoting reactivity and hydrogen atom addition to the furan ring inhibiting reactivity. This work, to the authors knowledge, is the first to combine theoretical and experimental work to describe the oxidation of any of the alkylated furans. The mechanism developed herein to describe 2MF combustion should also function as a sub-mechanism to describe the oxidation of 2,5-dimethyl furan whilst also providing key insights into the oxidation of this similar biofuel candidate.