Detailed and simplified kinetic models of n-dodecane oxidation: The role of fuel cracking in aliphatic hydrocarbon combustion

A detailed kinetic model is proposed for the combustion of normal alkanes up to n-dodecane above 850 K. The model was validated against experimental data, including fuel pyrolysis in plug flow and jet-stirred reactors, laminar flame speeds, and ignition delay times behind reflected shock waves, with n-dodecane being the emphasis. Analysis of the computational results reveal that for a wide range of combustion conditions, the kinetics of fuel cracking to form smaller molecular fragments is fast and may be decoupled from the oxidation kinetics of the fragments. Subsequently, a simplified model containing a minimal set of 4 species and 20 reaction steps was developed to predict the fuel pyrolysis rate and product distribution. Combined with the base C₁-C₄ model, the simplified model predicts fuel pyrolysis rate and product distribution, laminar flame speeds, and ignition delays as close as the detailed reaction model.     

Combustion chemistry probed by synchrotron VUV photoionization mass spectrometry

Combustion is directly related to energy conversion and the environment. Gas-phase chemical reactions such as thermal decomposition, oxidation and recombination play a critical role in combustion processes. Here we review six applications of synchrotron vacuum-ultraviolet (VUV) photoionization mass spectrometry (PIMS) in fundamental studies of combustion chemistry. These applications range from the use of flow reactors to probe elementary reaction kinetics, studies of pyrolysis in plug-flow reactors and oxidation in jet-stirred reactors, studies of spatial evolution of species concentrations in premixed and non-premixed flames, product distributions in pyrolysis of biomass, and analysis of polycyclic aromatic hydrocarbon (PAH) formation. These experiments provide valuable data for the development and validation of detailed chemical kinetic models. Furthermore, some additional potential applications are proposed.     

Effects of benzene,cyclo-hexane and n-hexane addition to methane on soot yields in high-pressure laminar diffusion flames

Effects of doping high pressure methane diffusion flames with benzene, cyclo-hexane and n-hexane were investigated to assess the sooting propensity of three hydrocarbons with six carbons at elevated pressures. Amount of liquid hydrocarbons added to methane constituted 7.5% of the total carbon content of the fuel stream. The pressure range investigated extended up to 10 bar and the experiments were carried out in a high pressure combustion chamber capable of establishing stable laminar diffusion flames with various fuels at elevated pressures and was used in similar experiments previously. Temperatures and soot volume fractions were measured using the spectral soot emission technique capturing spectrally-resolved line-of-sight intensities which were subsequently inverted using an Abel type algorithm to obtain radial distributions assuming that the flames are axisymmetric. The total mass carbon flow of the fuel stream was kept constant at 0.524 mg/s in neat methane, benzene-doped methane, cyclo-hexane-doped methane, and n-hexane-doped methane flames to have tractable measurements at all pressures. Measured maximum soot volume fractions and evaluated maximum soot yields showed that benzene-doped methane flame had the higher values than cyclo-hexane doped methane flames which in turn had higher values than n-hexane doped methane flames at all pressures. Sooting propensity dependence of the three hydrocarbons on pressure can be ranked as, in descending order, n-hexane, cyclo-hexane, and benzene; however, the difference between pressure dependencies of n-hexane and cyclo-hexane was within the measurement error margins. Ratio of soot yields of benzene to n-hexane doped flames changed from about 2 at 2 bar to 1.2 at 10 bar; the ratio of benzene to cyclo-hexane doped flames showed similar trends.     

Mechanisms of the oxidation and combustion of normal paraffin hydrocarbons: Transition from C1–C10 to C11–C16

Recently, detailed kinetic mechanisms of the oxidation and combustion of higher hydrocarbons, composed of hundreds of components and thousands of elementary reactions, have been proposed. Despite the undoubtful advantages of such detailed mechanisms, their application to simulations of turbulent combustion and gas dynamic phenomena is difficult because of their complexity. At the same time, to some extent limited, they cannot be considered exhaustive. This work applies previously proposed algorithm for constructing an optimal mechanism of the high- and low-temperature oxidation and combustion of normal paraffin hydrocarbons, which takes into account the main processes determining the reaction rate and the formation of key intermediates and final products. The mechanism has the status of a nonempirical detailed mechanism, since all the constituent elementary reactions have a kinetic substantiation. The mechanism has two specific features: (1) it does not include reactions of so-called double oxygen addition (first to the peroxide radical, and then to its isomeric form), i.e., the first addition turns out to be sufficient; (2) it does not include isomeric compounds and their derivatives as intermediates, since this oxidation pathway is slower than the oxidation of molecules and radicals with normal structure. Application of the algorithm makes it possible to compile a compact mechanism, which is important for modeling chemical processes involving paraffin hydrocarbons C _n with large n. Previously, based on this algorithm, compact mechanisms of the oxidation and combustion of propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane have been constructed. In this work, we constructed a nonempirical detailed mechanism of the oxidation and combustion of hydrocarbons from n-undecane to n-hexadecane. The most important feature of the new mechanism is its staged nature, which manifests itself through the emergence of cool and blue flames during low-temperature autoignition. The calculation results are compared with experimental data.     

二甲醚/空气对向流扩散燃烧反应历程的数值研究

本文采用详细化学反应动力学模型对二甲醚对向流扩散燃烧火焰进行了数值模拟,通过分析二甲醚燃烧过程中基元反应速度、关键中间产物和自由基,得到了二甲醚在对向流扩散燃烧中氧化的主要反应途径.结果表明,反应主要发生在高温区域(大于 800 K),CH2O、H2 和 CH4是重要的中间产物;OH、H 和 CH3是重要的自由基,在生成 CO、CO2和 H2O 的过程中起到关键性作用.     

Skeletal mechanism generation and analysis for n-heptane with CSP

We use a procedure based on the decomposition into fast and slow dynamical components offered by the Computational Singular Perturbation (CSP) method to generate automatically skeletal kinetic mechanisms for the simplification of the kinetics of n-heptane oxidation. The detailed mechanism of the n-heptane oxidation here considered has been proposed by Curran et al. and involves 561 species and 2538 reactions. After carrying out a critical assessment of important aspects of this procedure, we show that the comprehensive skeletal kinetic mechanisms so generated are able to reproduce the main features of n-heptane ignition at various initial pressures and temperatures and equivalence ratios. A by-product of the algorithm that generates the skeletal mechanisms is the identification of the network of important species and reactions at a given state of the kinetic system. The analysis of this network is carried out by resorting to a visual representation of the pathways at selected time instants of the ignition process. Visual inspection of the pathways enables the identification and comparison of the relevant kinetic processes as obtained at different ignition regimes. The graphs are generated by interfacing the model reduction procedure with the open-source package graphviz.     

Effect of pressure on soot formation in the pyrolysis of <Emphasis Type="Italic">n</Emphasis>-hexane and the oxidation of fuel-rich mixtures of <Emphasis Type="Italic">n</Emphasis>-heptane behind reflected shock waves

The results of detailed kinetic simulations of the formation of soot particles in the pyrolysis of n-hexane–argon mixtures and in the oxidation of fuel-rich (φ = 5) n-heptane–oxygen–argon mixtures behind reflected shock waves at pressures of 20–100 bar and a constant concentration of carbon atoms or a constant fraction of argon in the initial mixture within the framework of a modified reaction mechanism are reported. The choice of n-hexane and n-heptane for examining the effect of pressure on the process of soot formation was motivated by the availability for these hydrocarbons of experimental measurements in reflected shock waves at high pressures (up to ~100 bar). The temperature dependences of the yield of soot particles formed in the pyrolysis of n-hexane are found to be very weakly dependent on pressure and slightly shifting to lower temperatures with increasing pressure. In general, pressure produces a very weak effect on the soot formation in the pyrolysis of n-hexane. The effect of pressure and concentration of carbon atoms in the initial mixture on the process of soot formation during the oxidation of fuel-rich n-heptane mixtures behind reflected shock waves is studied. The results of our kinetic simulations show that, for both the pyrolysis of n-hexane and the oxidation of fuel-rich n-heptane–oxygen mixtures, the influence of pressure on the process of soot formation is negligible. By contrast, the concentration of carbon atoms in the initial reaction mixture produces a much more pronounced effect.     

Modeling high pressure ethane oxidation and pyrolysis

A previously developed detailed chemical kinetic model for ethane oxidation and pyrolysis has been extended by considering the impact of much higher pressures on the chemically- and thermally-activated reactions in the mechanism. The resulting mechanism is now compared to data at much higher pressures. These data include both ethane oxidation and pyrolysis shock tube experiments up to pressures of 986 atm as well as autoignition and ethane dissociation data near 10 atm. The temperature and pressure dependencies of the rate coefficients in the model are represented by Chebyshev polynomials. This model, with no adjustments, describes the varied experimental data very well. A first order sensitivity analysis identified the most important reactions in each of the kinetic regimes. These results, coupled with the earlier validation studies at lower pressures, suggest this model is broadly applicable to analysis of ethane pyrolysis and oxidation over very wide temperature and pressure ranges.     

A wide-range experimental and kinetic modeling study of the pyrolysis and oxidation of 2-butyne

To reduce particulate emissions leading to a cleaner environment, it is important to understand how polycyclic-aromatic hydrocarbons (PAHs) and their precursors are formed during combustion. 2-butyne can decompose to propargyl and allyl radicals. These radicals can produce benzene and other PAHs, leading to the formation of soot. In the present study, pyrolysis, oxidation, and laminar flame speed experiments were performed for 2-butyne. The pyrolysis experiments were conducted in a single-pulse shock tube at 2 bar in the temperature range 1000 – 1500 K. Ignition delay times for 2-butyne/‘air’ mixtures were measured in the pressures range 1 – 50 bar, over the temperature range 660 – 1630 K, at equivalence ratios of 0.5, 1.0, and 2.0 using rapid compression machines and shock tubes. Moreover, laminar flame speed (LFS) experiments were performed at ambient temperature, at p = 1 – 3 atm, over an equivalence ratio range of 0.6 – 1.8. A new, detailed chemical kinetic model for 2-butyne has been developed and widely validated against the data measured in this study and those available in the literature. The significant reactions for 2-butyne pyrolysis, ignition, and oxidation are identified and discussed using flux and sensitivity analyses.     

Development of a reduced biodiesel surrogate model for compression ignition engine modeling

Methylbutanoate (MB), a C₄ methyl ester, represents the simplest surrogate that captures the chemical effects of the ester moiety in biodiesel and biodiesel surrogates. An updated chemical kinetic model has been developed to characterize the ignition and flame characteristics of MB. The mechanistic elements within this model that relate to the MB and smaller ester/oxygenate sub-mechanisms are drawn from the prototypical Fisher et al. model and from more recent theory and modeling efforts. The MB model development which is based on an iterative procedure involving global sensitivity analyses to identify elementary reactions that govern ignition and subsequent high level ab initio based theoretical updates to these reaction rates are presented. The MB model makes reasonable predictions of ignition delays and laminar flame speeds.The C₅–C₇ submechanisms from the LLNL n-heptane (NH) model were merged with the present MB model to obtain a detailed chemical kinetics model for a surrogate blend representing biodiesel. The detailed MB-NH model (661 species) was reduced using graph based techniques. The robust reduction techniques employed result in a reduced model (145 species) that is in good agreement with the detailed model over a wide range of conditions. 3-D compression ignition (CI) engine simulations utilizing this reduced chemistry model for MB-NH blends as a surrogate for biodiesel show good agreement with the experimental data suggesting the utility of this model for predictions of combustion and emission characteristics of biodiesel in realistic CI engine simulations.     

二甲醚HCCI燃烧高温反应动力学分析

应用单区燃烧模型对二甲醚均质压燃燃烧的化学反应动力学过程进行了数值模拟研究。通过分析在内燃机压燃燃烧边界条件下二甲醚高温氧化反应过程中的关键基元反应速度、关键中间产物以及自由基的浓度随曲轴转角的变化,得到了二甲醚燃烧氧化的高温反应途经。结果表明,二甲醚均质压燃燃烧具有明显的两阶段放热特性,即低温反应放热和高温反应放热。高温反应阶段又可分为蓝焰反应阶段和热焰反应阶段,其中蓝焰反应阶段是甲醛氧化成CO的过程,热焰反应主要是CO氧化成CO2的过程。二甲醚氧化产物之一甲酸(HOCHO)在蓝焰反应阶段分解生成CO2。     

A detailed chemical kinetic model of high-temperature ethylene glycol gasification

In recent experimental investigations, ethylene glycol is used as a model substance for biomass-based pyrolysis oil in an entrained flow gasifier. In order to gain a deeper insight into process sequences and to conduct parametric analysis, this study describes the development and validation of a detailed chemical kinetic model of high-temperature ethylene glycol gasification. A detailed reaction mechanism based on elementary reactions has been developed considering 80 species and 1243 reactions for application in CFD software. In addition to mechanism validation based on ignition delay times, laminar flame speeds and concentration profiles, simulation results are compared to experimental data of ethylene glycol gasification under complex turbulent reactive flow conditions.     

Isomer-specific influences on ignition and intermediates of two C5 ketones in an RCM

Ketones have been considered as potential biofuels and main components of blend stock for internal engines. To better understand the chemical kinetics of ketones, ignition delay times of 2-pentanone (propyl methyl ketone, PMK) and 3-pentanone (diethyl ketone, DEK) were measured at temperatures of 895–1128 K under 10 and 20 bar, at equivalence ratios (?) of 0.5 and 1.0 in a rapid compression machine (RCM). To explore the impact of carbonyl functionality and resonance stabilized structures of fuel radicals on their combustion kinetics, high-temperature pyrolysis at 1130 K and relatively low-temperature oxidation at 950 K studies were performed in an RCM, and the time-resolved species concentration profiles under these two conditions were quantified using a fast sampling system and gas chromatography (GC). A new kinetic model containing low-temperature reactions was built aiming at predicting the pyrolysis and oxidation behaviors of both ketones. The consumption pathways of the resonance stabilization fuel radicals through oxygen addition and following reactions are promoted since the decomposition rates of these radicals are about 4 orders magnitudes lower than regular fuel radicals. The occurrences of the so-called “addition-dissociation reactions”, i.e., ketones reacting with a hydrogen yielding aldehyde or reacting with a methyl radical yielding shorter-chain-length ketones, are verified in pyrolysis experiments. Based on experiments and model analysis, the carbonyl functionality in both ketones is preserved during the process of β-scissions of fuel radicals and α-scissions of fuel-related acyl radicals, resulting in the direct formation of CO and ketene. However, the position of carbonyl functionality has a significant impact on the species pools.     

Mechanism of the oxidation and combustion of normal paraffin hydrocarbons: Transition from C<Subscript>1</Subscript>–C<Subscript>6</Subscript> to C<Subscript>7</Subscript>H<Subscript>16</Subscript>

A previously proposed algorithm for constructing an optimal mechanism of the high- and low-temperature oxidation and combustion of normal paraffin hydrocarbons was used, which includes the major processes that determine the rate of reaction and the formation of the main intermediate and final products. The mechanism has the status of a nonempirical detailed mechanism, since all the constituent elementary reactions have a kinetic substantiation. The mechanism has two specific features: it included no reactions of so-called double addition of oxygen and no isomeric compounds and derivatives thereof as intermediate species. Realization of this algorithm leads to fairly compact models, a circumstance important for studies of chemical processes involving paraffin hydrocarbons C _n with large n. Previously, based on this algorithm, compact mechanisms of oxidation and combustion of propane, n-butane, n-pentane, and n-hexane were constructed. In this paper, we develop a nonempirical detailed mechanism of oxidation and combustion of n-heptane. The most important feature of the new mechanism is its ability to predict the staging of the process in the form of cool and blue flames at low autoignition temperatures. A comparison of the simulation results with the available experimental data is conducted.     

Insights into the oxidation of propylene oxide through the analysis of experiments and kinetic modeling

Cyclic ethers are important intermediates in the oxidation of hydrocarbons and biofuels. Studying the oxidation and pyrolysis of cyclic ethers will help in improving our understanding of this functional group and provide consistency to the base mechanism where they play an important role. In this aspect, propylene oxide has been investigated in this study by obtaining ignition delay time measurements in the rapid compression machine and shock tube. The experiments were performed in a range of pressures varying from 10 to 40 bar at different equivalence ratios (0.5–2.0) and dilution percentages. Additionally, speciation measurements in the shock tube at pyrolysis conditions have been performed at a pressure of 40 bar to explore the isomerization pathways. A detailed kinetic mechanism was developed to describe both the oxidation and pyrolysis chemistry of propylene oxide. The mechanism is not only able to predict the data obtained from this study but also reproduces the data from the literature in a consistent trend. For a better understanding of the oxidation and pyrolysis chemistry of propylene oxide, the kinetic analyses were performed using the developed mechanism to comprehend the important reaction pathways and sensitive reactions. At the investigated regime, the consumption of propylene oxide through its isomerization channels is the critical pathway that controls the reactivity of the fuel.     

A functional-group-based approach to modeling real-fuel combustion chemistry – III: Application to biodiesels

Real biodiesel fuels are mixtures comprising many high molecular weight components, making it a challenge to predict their combustion chemistry with detailed kinetic models. Our group previously proposed a functional-group approach (FGMech) to model the combustion chemistry of real gasoline and jet fuels; this approach has now been extended to model real biodiesel combustion and mixtures with petroleum fuels. As in our previous work, a decoupling philosophy is adopted for construction of the model. A lumped reaction mechanism describes the (oxidative) pyrolysis of fuels, while a detailed base chemistry model represents the oxidation of key pyrolysis intermediates. However, due to the presence of the ester group, several oxygenated species are identified as additional primary products and incorporated into the lumped reaction steps. In addition to the lumped reactions initiated by unimolecular decomposition and H-atom abstraction reactions, a lumped H-atom addition-elimination reaction is also incorporated as a new reaction class to account for the presence of double bonds. Stoichiometric parameters are obtained based on a multiple linear regression (MLR) model, which establishes relationships between the fuel's functional group distributions and the stoichiometric parameters of the lumped reactions. Global rate constants are developed from consistent rate rules obtained from pure fuels. New pyrolysis experimental data for methyl pentanoate/methyl nonanoate and methyl heptanoate/n-heptane mixtures (50%/50% in mol) are obtained in a jet-stirred reactor at atmospheric pressure. In general, kinetic models developed using the FGMech approach can reasonably reproduce all the validation targets obtained in this work, as well as those in the literature, confirming that functional-group-modeling is a promising approach to simulate combustion behavior of diesel/biodiesel surrogate fuels and real biodiesels.     

The decisive role of pericyclic reactions in the thermal decomposition of organophosphorus compounds

The understanding of the thermal decomposition chemistry of chemical warfare nerve agents is largely limited by the scarcity of kinetic data. Because of the high toxicity of these molecules, experimental determination of their chemical properties is very difficult. In the present work, a comprehensive detailed kinetic model for the decomposition of sarin and some simulants, i.e. di-isopropyl methyl phosphonate (DIMP), diethyl methylphosphonate (DEMP), and triethyl phosphate (TEP) were developed, containing possible molecular and radical pathways. The importance of unimolecular pericyclic decomposition led to evaluate precisely the rate constants of these reactions with high level theoretical calculations. The QCISD(T)/cc-PV∞QZ//B2PLYPD3/6–311+G(2d,d,p) level of theory was selected after a benchmark. The contribution of hindered rotors was included with the 1D-HR-U approach. Tunneling was taken into account for H-atom transfer. Transition state theory was used to calculate high-pressure limit rate constants and pressure dependent rate constants were calculated using Master Equation modeling. The model was validated against experimental pyrolysis and oxidation experimental data available in literature. Flux analyses showed that whatever the conditions are, the first step of decomposition of the studied phosphorus compounds are pericyclic eliminations leading to successive decompositions, whereas bond-breaking or H-atom abstraction remain negligible, even at high temperature.     

Experimental and semi-detailed kinetic modeling study of decalin oxidation and pyrolysis over a wide range of conditions

Decalin is the simplest polycyclic alkane (polynaphtenic hydrocarbon) found in liquid fuels (jet fuels, Diesel). In order to better understand the combustion characteristics of decalin, this study provides new experimental data for its oxidation in a jet-stirred reactor. For the first time, stable species concentration profiles were measured in a jet-stirred reactor at a constant mean residence time of 0.1 s and 0.5 s at respectively 1 and 10 atm, over a range of equivalence ratios (? = 0.5–1.5) and temperatures (750–1350 K). The oxidation of decalin under these experimental conditions was modeled using a semi-detailed chemical kinetic reaction mechanism (11,000 reactions involving 360 species) derived from a previously proposed scheme for the ignition of the same fuel in a shock-tube. The proposed mechanism that includes both low- and high-temperature chemistry shows reasonably good agreement with the present experimental data set. It can also represent well decalin pyrolysis and oxidation data available in the literature. Reaction path analyses and sensitivity analyses were conducted to interpret the results.     

Oxidation of H2/CO2 mixtures and effect of hydrogen initial concentration on the combustion of CH4 and CH4/CO2 mixtures: Experiments and modeling

An experimental investigation of the oxidation of hydrogen diluted by nitrogen in presence of CO₂ was performed in a fused silica jet-stirred reactor (JSR) over the temperature range 800-1050 K, from fuel-lean to fuel-rich conditions and at atmospheric pressure. The mean residence time was kept constant in the experiments: 120 ms at 1 atm and 250 ms at 10 atm. The effect of variable initial concentrations of hydrogen on the combustion of methane and methane/carbon dioxide mixtures diluted by nitrogen was also experimentally studied. Concentration profiles for O₂, H₂, H₂O, CO, CO₂, CH₂O, CH₄, C₂H₆, C₂H₄, and C₂H₂ were measured by sonic probe sampling followed by chemical analyses (FT-IR, gas chromatography). A detailed chemical kinetic modeling of the present experiments and of the literature data (flame speed and ignition delays) was performed using a recently proposed kinetic scheme showing good agreement between the data and this modeling, and providing further validation of the kinetic model (128 species and 924 reversible reactions). Sensitivity and reaction paths analyses were used to delineate the important reactions influencing the kinetic of oxidation of the fuels in absence and in presence of additives (CO₂ and H₂). The kinetic reaction scheme proposed helps understanding the inhibiting effect of CO₂ on the oxidation of hydrogen and methane and should be useful for gas turbine modeling.     

A wide-range experimental and kinetic modeling study of the pyrolysis and oxidation of 1-butyne

Allyl and propargyl radicals are involved in the production of the first aromatic ring, which is considered a crucial process in forming polycyclic-aromatic hydrocarbons and soot. 1-Butyne decomposes to propargyl radicals during its pyrolysis and oxidation. To improve our knowledge of the kinetics of 1-butyne, its pyrolysis, oxidation, and laminar flame speed properties have been measured. Pyrolysis experiments were performed in a single-pulse shock tube at 2 bar in the temperature range 1000 – 1600 K. Ignition delay times for 1-butyne/‘air’ mixtures were measured at pressures of 1, 10, 30, and 50 bar, in the temperature range 680 – 1580 K, at equivalence ratios of 0.5, 1.0, and 2.0 using rapid compression machines and shock tubes. Furthermore, laminar flame speeds were measured at ambient temperature, at p = 1, 2, 3 atm, over an equivalence ratio range of 0.6 – 1.9. A new detailed mechanism for 1-butyne has been developed and widely validated using the new experimental data and those available in the literature. Important reactions of 1-butyne pyrolysis and oxidation are determined through flux and sensitivity analyses.