A detailed chemical kinetic model for high temperature ethanol oxidation

A detailed chemical kinetic model for ethanol oxidation has been developed and validated against a variety of experimental data sets. Laminar flame speed data (obtained from a constant volume bomb and counterflow twin‐flame), ignition delay data behind a reflected shock wave, and ethanol oxidation product profiles from a jet‐stirred and turbulent flow reactor were used in this computational study. Good agreement was found in modeling of the data sets obtained from the five different experimental systems. The computational results show that high temperature ethanol oxidation exhibits strong sensitivity to the fall‐off kinetics of ethanol decomposition, branching ratio selection for C₂H₅OH + OH ↔ Products, and reactions involving the hydroperoxyl (HO₂) radical. The multichanneled ethanol decomposition process is analyzed by RRKM/Master Equation theory, and the results are compared with those obtained from earlier studies. The ten‐parameter Troe form is used to define the C₂H₅OH(+M) ↔ CH₃ + CH₂OH(+M) rate expression as k^∞ = 5.94E23 T^−1.68 exp(−45880 K/T) (s⁻¹) k^o = 2.88E85 T^−18.9 exp(−55317 K/T) (cm³/mol/sec) F_cent = 0.5 exp(−T/200 K) + 0.5 exp(−T/890 K) + exp(−4600 K/T) and the C₂H₅OH(+M) ↔ C₂H₄ + H₂O(+M) rate expression as k^∞ = 2.79E13 T^0.09 exp(−33284 K/T) (s⁻¹) k^o = 2.57E83 T^−18.85 exp(−43509 K/T) (cm³/mol/sec) F _cent = 0.3 exp(−T/350 K) + 0.7 exp(−T/800 K) + exp(−3800 K/T) with an applied energy transfer per collision value of <ΔE_down> = 500 cm⁻¹. An empirical branching ratio estimation procedure is presented which determines the temperature dependent branching ratios of the three distinct sites of hydrogen abstraction from ethanol. The calculated branching ratios for C₂H₅OH + OH, C₂H₅OH + O, C₂H₅OH + H, and C₂H₅OH + CH₃ are compared to experimental data. © 1999 John Wiley & Sons, Inc. Int J Chem Kinet 31: 183–220, 1999     

A comprehensive mechanism for methanol oxidation

A comprehensive detailed chemical kinetic mechanism for methanol oxidation has been developed and validated against multiple experimental data sets. The data are from static-reactor, flow-reactor, shock-tube, and laminar-flame experiments, and cover conditions of temperature from 633–2050 K, pressure from 0.26–20 atm, and equivalence ratio from 0.05–2.6. Methanol oxidation is found to be highly sensitive to the kinetics of the hydroperoxyl radical through a chain-branching reaction sequence involving hydrogen peroxide at low temperatures, and a chain-terminating path at high temperatures. The sensitivity persists at unusually high temperatures due to the fast reaction of CH₂OH+O₂=CH₂O+HO₂ compared to CH₂OH+M=CH₂O+H+M. The branching ratio of CH₃OH+OH=CH₂OH/CH₃O+H₂O was found to be a more important parameter under the higher temperature conditions, due to the rate-controlling nature of the branching reaction of the H-atom formed through CH₃O thermal decomposition. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet 30: 805–830, 1998     

Methanol oxidation in a flow reactor: Implications for the branching ratio of the CH3OH+OH reaction

The oxidation of methanol in a flow reactor has been studied experimentally under diluted, fuel-lean conditions at 650–1350 K, over a wide range of O₂ concentrations (1%–16%), and with and without the presence of nitric oxide. The reaction is initiated above 900 K, with the oxidation rate decreasing slightly with the increasing O₂ concentration. Addition of NO results in a mutually promoted oxidation of CH₃OH and NO in the 750–1100 K range. The experimental results are interpreted in terms of a revised chemical kinetic model. Owing to the high sensitivity of the mutual sensitization of CH₃OH and NO oxidation to the partitioning of CH₃O and CH₂OH, the CH₃OH + OH branching fraction could be estimated as α = 0.10 ± 0.05 at 990 K. Combined with low-temperature measurements, this value implies a branching fraction that is largely independent of temperature. It is in good agreement with recent theoretical estimates, but considerably lower than values employed in previous modeling studies. Modeling predictions with the present chemical kinetic model is in quantitative agreement with experimental results below 1100 K, but at higher temperatures and high O₂ concentration the model underpredicts the oxidation rate. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 40: 423–441, 2008     

Theoretical study on the gas‐phase reaction mechanism between rhodium monoxide and methane for methanol production

The gas‐phase reaction mechanism between methane and rhodium monoxide for the formation of methanol, syngas, formaldehyde, water, and methyl radical have been studied in detail on the doublet and quartet state potential energy surfaces at the CCSD(T)/6‐311+G(2d, 2p), SDD//B3LYP/6‐311+G(2d, 2p), SDD level. Over the 300–1100 K temperature range, the branching ratio for the Rh(⁴F) + CH₃OH channel is 97.5–100%, whereas the branching ratio for the D‐CH₂ORh + H₂ channel is 0.0–2.5%, and the branching ratio for the D‐CH₂ORh + H₂ channel is so small to be ruled out. The minimum energy reaction pathway for the main product methanol formation involving two spin inversions prefers to both start and terminate on the ground quartet state, where the ground doublet intermediate CH₃RhOH is energetically preferred, and its formation rate constant over the 300–1100 K temperature range is fitted by k_CH3RhOH = 7.03 × 10⁶ exp(?69.484/RT) dm³ mol^?1 s^?1. On the other hand, the main products shall be Rh + CH₃OH in the reactions of RhO + CH₄, CH₂ORh + H₂, Rh + CO +2H₂, and RhCH₂ + H₂O, whereas the main products shall be CH₂ORh + H₂ in the reaction of Rh + CH₃OH. Meanwhile, the doublet intermediates H₂RhOCH₂ and CH₃RhOH are predicted to be energetically favored in the reactions of Rh + CH₃OH and CH₂ORh + H₂ and in the reaction of RhCH₂ + H₂O, respectively. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

A Mass-spectrometric study of propionaldehyde oxidation in the negative temperature coefficient region

The oxidation of propionaldehyde has been investigated in a 1-L Pyrex reactor at total pressures of 50–120 torr and temperatures 553–713 K. Detection of reactants and products was principally by molecular beam mass spectrometry, although certain species could only be measured by gas-chromatographic analysis. At 553 K the yield of water was ～83% of the propionaldehyde consumed, leading to the conclusion that OH is the principal chain carrier near the beginning of the negative temperature coefficient region. Many oxygenated organics (CH₂O, CH₃CHO, C₂H₅OH, C₂H₅O₂H, CH₃O₂H) and C₂H₄ are formed during the oxidation process. These oxidation products are consistent with the important role of O₂ addition to C₂H₅ radicals at 553 K followed by subsequent reactions of the C₂H₅O₂ radical. As the temperature is increased, the product concentrations smoothly change to a much simpler distribution in which C₂H₄, H₂O₂, and CO are the dominant products.     

Ethylene pyrolysis and oxidation: A kinetic modeling study

Ethylene oxidation and pyrolysis was modeled using a comprehensive kinetic reaction mechanism. This mechanism is an updated version of one developed earlier. It includes the most recent findings concerning the kinetics of the reactions involved in the oxidation of ethylene. The proposed mechanism was tested against ethylene oxidation experimental data (molecular species concentration profiles) obtained in jet stirred reactors (1–10 atm, 880–1253 K), ignition delay times measured in shock tubes (0.2–12 atm, 1058–2200 K) and ethylene pyrolysis data in shock tube (2–6 atm, 1700–2200 K). The general prediction of concentration profiles of minor species formed during ethylene oxidation is improved in the present model by using more accurate kinetic data for several reactions (principally: HO₂ + HO₂ → H₂O₂ + O₂, C₂H₄ + OH → C₂H₃ + H₂O, C₂H₂ + OH → Products, C₂H₃ → C₂H₂ + H).     

A Computational Study of the Kinetics and Mechanism for the C2H3 + CH3OH Reaction

The mechanism for the C₂H₃ + CH₃OH reaction has been investigated by the Gaussian‐4 (G4) method based on the geometric parameters of the stationary points optimized at the B3LYP/6–31G(2df, p) level of theory. Four transition states have been identified for the production of C₂H₄ + CH₃O (TSR/P1), C₂H₄ + CH₂OH (TSR/P2), C₂H₃OH + CH₃ (TSR/P3), and C₂H₃OCH₃ + H (TSR/P4) with the corresponding barriers 8.48, 9.25, 37.62, and 34.95 kcal/mol at the G4 level of theory, respectively. The rate constants and branching ratios for the two lower energy H‐abstraction reactions were calculated using canonical variational transition state theory with the Eckart tunneling correction at the temperature range 300–2500 K. The predicted rate constants have been compared with existing literature data, and the uncertainty has been discussed. The branching ratio calculation suggests that the channel producing CH₃O is dominant up to about 1070 K, above which the channel producing CH₂OH becomes very competitive.     

Deuterium Isotope Effect on Excess Enthalpies of Methanol or Ethanol and Their Deuterium Derivatives at 298.15 K

Excess enthalpies of six binary mixtures of CH₃ OD+CH₃ OH, CH₃ OD+CD₃ OD, CD₃ OD+CH₃ OH, C₂ D₅ OD+C₂ H₅ OH, C₂ D₅ OD+C₂ H₅ OD, C₂ H₅ OD+C₂ H₅ OH have been determined over the whole range of mole fractions at 298.15 K in order to know the isotopic effect on hydrogen-bonding accurately, although there are many reports on the differences in the strength of hydrogen-bonding between OH and OD. All excess enthalpies measured are very small and endothermic. The mixtures of CH₃ OD+ CH₃ OH, and C₂ D₅ OD+C₂ H₅ OH showed the largest excess enthalpies among each methanol and ethanol mixtures. The difference of intermolecular interaction between OH and OD in methanol and ethanol was almost same value of (1.82±0.04) J mol^-1 Excess enthalpies of 1,4-dimethylbenzene+1,3-dimethylbenzene and 1,4-dimethylbenzene+1,2-methylbenzene were measured by three different principle calorimeters at 298.15 K in order to know the precision of calorimetry for a small enthalpy change. This revised version was published online in July 2006 with corrections to the Cover Date.     

Kinetics and mechanism associated with the reactions of hydroxyl radicals and of chlorine atoms with 1‐propanol under near‐tropospheric conditions between 273 and 343 K

Rate constants for the reactions of OH radicals and Cl atoms with 1‐propanol (1‐C₃H₇OH) have been determined over the temperature range 273–343 K by the use of a relative rate technique. The value of k(Cl + 1‐C₃H₇OH) = (1.69 ± 0.19) × 10^?12 cm³ molecule^?1 s^?1 at 298 K and shows a small increase of 10% between 273 and 342 K. The value of k(OH + 1‐C₃H₇OH) increases by 14% between 273 and 343 K with a value of (5.50 ± 0.55) × 10^?12 cm³ molecule^?1 s^?1 at 298 K, and further when combined with a single independent experimentally determined value at 753 K gives k(OH + 1‐C₃H₇OH) = 4.69 × 10^?17T^1.8 exp(422/T) cm³ molecule^?1 s^?1, which fits each data point to better than 2%. Two well‐established structure–activity relationships for H abstraction by OH radicals give accurate predictions of the rate constant for OH + 1‐C₃H₇OH, provided the β‐CH₂ group is given an increased reactivity of a factor of about 2 over that for the structurally equivalent CH₂ group in alkanes at 298 K. A quantitative product analysis was carried out at 298 K for the Cl‐initiated photooxidation of 1‐C₃H₇OH, using both FTIR and gas chromatography. HCHO, CH₃CHO, and C₂H₅CHO were the only major organic primary products observed, although HCOOH was found in much smaller amounts as a secondary product. A key characteristic of the analysis was that the initial values of the product ratio [CH₃CHO]/[C₂H₅CHO] were effectively constant for NO pressures between 0.15 and 0.3 Torr, but fell by about 35% as the pressure fell to 0.0375 Torr. From a detailed consideration of the mechanism for the oxidation, it is suggested that C₂H₅CHO, CH₃CHO (+HCHO), and 3 molecules of HCHO are formed uniquely from CH₃CH₂CHOH, CH₃CHCH₂OH, and CH₂CH₂CH₂OH radicals, respectively. On this basis, use of the product yields gives the branching ratios of 56, 30, and 14% for Cl atom reaction at the α‐, β‐, and γ‐C? H positions in 1‐C₃H₇OH at 298 K. Given the very low temperature coefficients involved, little change will occur over tropospheric temperature ranges. © 2001 John Wiley & Sons, Inc. Int J Chem Kinet 34: 110–121, 2002     

Bimolecular QRRK analysis of methyl radical reactions

Reactions which proceed through energized adducts, including radical recombinations, insertions, and addition to unsaturates, frequently exhibit unusual kinetic behavior. The branching ratios among various product channels are often complex functions of both temperature and pressure. Four such reactions involving methyl radicals are analyzed by combining chemical activation distribution functions with QRRK methods to predict rate constants for each channel. These include three oxidation paths, CH₃ + O, CH₃ + O₂, CH₃ + OH, and the addition reaction CH₃ + C₂H₂. These predictions are compared to experiments wherever possible; generally, the agreement is quite satisfactory. Analysis of the energetics of the various reaction channels, using parameters which are readily available, provides a convenient framework for prediction. Suggested rate constants for the various channels for the four reactions are given at three pressures, 20, 760, and 7600 Torr, for the temperature range 300–2500 K. The approach used here can easily be applied to other reactions.     

Master Equation Modeling of the Unimolecular Decompositions of α‐Hydroxyethyl (CH3CHOH) and Ethoxy (CH3CH2O) Radicals

The unimolecular decomposition of two radical isomers of C₂H₅O (CH₃CH₂O/ethoxy, CH₃CHOH/α‐hydroxyethyl) are investigated by means of Rice–Ramsperger–Kassel–Marcus/master equation simulations in helium and nitrogen bath gases on an accurate one‐dimensional potential energy surface. For ethoxy, simulations are carried out between temperatures of 406 and 1200 K and pressures of 0.001 and 100 atm. For CH₃CHOH, simulations are carried out between temperatures of 800 and 1500 K and pressures of 0.001 and 100 atm. Results are compared with available experimental data, with good agreement. The dominant product of α‐hydroxyethyl decomposition is CH₃CHO + H, with C₂H₃OH + H and CH₃ + CH₂O, being minor channels. Rate coefficients are strongly dependent on temperature and pressure and are recommended with attendant uncertainty factor estimates. The relative roles of vinyl alcohol and acetaldehyde in the context of combustion chemistry are also discussed.     

Vapor — liquid and vapor — liquid — liquid phase equilibria of binary and ternary systems of nitrogen,ethene and methanol: Experiment and data evaluation

Zeck, S. and Knapp, H., 1986. Vapor—liquid and vapor—liquid—liquid phase equilibria of binary and ternary systems of nitrogen, ethene and methanol: experiment and data evaluation. Fluid Phase Equilibria, 26: 37–58.VLE and VLLE of three binary and one ternary system containing the components N₂, C₂H₄ and CH₃OH are investigated in a high-pressure phase equilibrium apparatus with vapor recirculation at temperatures 240 < T < 298 K and pressures 4 < p < 100 bar. Immiscibilities in the liquid phase are observed in the binary system C₂H₄CH₃OH with a lower critical end point and in the ternary system N₂C₂H₄CH₃OH.The experimental results are reported and compared with the results of other investigators and of available correlations.     

Role of CH, CH3, and OH Radicals in Organic Compound Decomposition by Water Plasmas

Decomposition of acetone, methanol, ethanol, and glycerine by water plasmas at atmospheric pressure has been investigated using a direct current discharge. At torch powers of 910–1,050 W and organic compound concentrations of 1–10 mol%, the decomposition rate of methanol and glycerine was over 99%, while those of acetone and ethanol was 95.4–99%. The concentrations of H₂ obtained were 60–80% in the effluent gas for any compounds by pyrolysis. Based on the experimental results, the decomposition mechanism of organic compounds in water plasmas was proposed and the roles of intermediate species such as CH, CH₃, and OH have been investigated; CH radical generated from organic compounds decomposition was easily oxidized to form CO; incomplete oxidation of CH₃ leads to C₂H₂ generation as well as soot formation; and negligible amount of soot observed from glycerine decomposition even at high concentration indicated that oxidation of CH_×(×:1–3) was enhanced by OH radical.     

Study of the ethane oxidation reaction by the kinetic tracer method

The kinetics of ethane oxidation was studied at 320, 340, 353 and 380°C, mixture composition 2 C₂H₆ + 1 O₂, and total pressure 609 torr. It was found that at 320°C CH₂O and CH₃CHO were branching agents. A series of experiments was conducted on 2C₂H₆ + O₂ oxidation in the presence of 0.7% ¹⁴C-labeled ethylene. The ethylene oxide was found to form only from C₂H₄, formaldehyde formed from C₂H₄ and C₂H₆; and CH₃CHO, C₂H₅OH, and CH₃OH formed only from ethane. The formation rates of C₂H₄, C₂H₄O, and CH₂O were calculated by the kinetic tracer method. At 320°C the fraction of oxygen-containing products formed from C₂H₄ was 16–18%, and at 353 and 380°C it was 30–40%.     

Toward a comprehensive mechanism for methanol pyrolysis

A single kinetic mechanism for methanol pyrolysis is tested against multiple sets of experimental data for the first time. Data are considered from static, flow, and shock tube reactors, covering temperatures of 973 to 2000 K and pressures of 0.3 to 1 atmosphere. The model results are highly sensitive to the rates of unimolecular fuel decomposition and of various chain termination reactions that remove CH₂OH and H radicals, as well as to experimental temperature uncertainties. The secondary fuel decomposition reaction CH₃OH = CH₂OH + H, which has previously been included only in mechanisms for high temperature conditions, is found to have a significant effect at low temperatures as well, through radical recombination. The reaction CH₃O + C = CH₃ + CO₂, rather than CH₃OH + H = CH₃ + H₂O, is found to be the dominant source of CH₃ at low temperatures. The reverse of CH₃ + OH = CH₂OH + H is important to CH₃ production at high temperatures.     

Experimental determination of the equilibrium constant of the reaction CH3 + O2 ⇄ CH3O2 during the gas-phase oxidation of methane

For the experimental determination of the equilibrium constant of the reaction CH₃ + O₂ ? CH₃O₂ (1), the process of methane oxidation has been studied over the temperature range of 706–786 K. The concentration of CH₃O₂ has been measured by the radical freezing method, and that of CH₃ from the rate of accumulation of ethane, assuming that C₂H₆ is produced by the reaction CH₃ + CH₃ → C₂H₆ (2). The equilibrium constant of reaction (1) has been obtained at four temperatures. For the heat of the reaction the value Δ?H₂₉₈ = -32.2 ± 1.5 kcal/mol is recommended.     

The decomposition of ethanol vapour by infrared radiation from a pulsed HF-laser

The decomposition of ethanol vapour induced by infrared radiation from a pulsed HF-laser has been studied as a function of pressure. At high pressures, above 10 torr, the main primary processes appear to be:C₂H₅OH → H₂ + CH₃CHO,C₂H₅OH → C₂H₄ + H₂O,C₂H₅OH → CH₃ + CH₂OHin a ratio of 3:2:1 which is independent of pressure. At low pressures the process yielding C₂H₄ and H₂O becomes dominant. The results suggest that the high pressure behaviour involves a “thermal” decomposition with collisional processes dominating, whereas at low pressures the decomposition is due to multiple photon absorption which at the lowest pressures approaches a collision-free unimolecular decomposition.     

Mechanism and Rate Constants of the CH3+ CH2CO Reaction: A Theoretical Study

The mechanism of the reaction of ketene with methyl radical has been studied by ab initio CCSD(T)‐F12/cc‐pVQZ‐f12//B2PLYPD3/6‐311G** calculations of the potential energy surface. Temperature‐ and pressure‐dependent reaction rate constants have been computed using the Rice–Ramsperger–Kassel–Marcus (RRKM)–Master Equation and transition state theory methods. Three main channels have been shown to dominate the reaction; the formation of the collisionally stabilized CH₃COCH₂ radical and the production of the C₂H₅ + CO and HCCO + CH₄ bimolecular products. Relative contributions of the CH₃COCH₂, C₂H₅ + CO, and HCCO + CH₄ channels strongly depend on the reaction conditions; the formation of thermalized CH₃COCH₂ is favored at low temperatures and high pressures, HCCO + CH₄ is dominant at high temperatures, whereas the yield of C₂H₅ + CO peaks at intermediate temperatures around 1000 K. The C₂H₅ + CO channel is favored by a decrease in pressure but remains the second most important reaction pathway after HCCO + CH₄ under typical flame conditions. The calculated rate constants at different pressures are proposed for kinetic modeling of ketene reactions in combustion in the form of modified Arrhenius expressions. Only rate constant to form CH₃COCH₂ depends on pressure, whereas those to produce C₂H₅ + CO and HCCO + CH₄ appeared to be pressure independent.     

Validation of a thermal decomposition mechanism of formaldehyde by detection of CH2O and HCO behind shock waves

The thermal decomposition of formaldehyde was investigated behind shock waves at temperatures between 1675 and 2080 K. Quantitative concentration time profiles of formaldehyde and formyl radicals were measured by means of sensitive 174 nm VUV absorption (CH₂O) and 614 nm FM spectroscopy (HCO), respectively. The rate constant of the radical forming channel (1a), CH₂O + M → HCO + H + M, of the unimolecular decomposition of formaldehyde in argon was measured at temperatures from 1675 to 2080 K at an average total pressure of 1.2 bar, k_1a = 5.0 × 10¹⁵ exp(‐308 kJ mol^?1/RT) cm³ mol^?1 s^?1. The pressure dependence, the rate of the competing molecular channel (1b), CH₂O + M → H₂ + CO + M, and the branching fraction β = k_1a/(kA_1a + k_1b) was characterized by a two‐channel RRKM/master equation analysis. With channel (1b) being the main channel at low pressures, the branching fraction was found to switch from channel (1b) to channel (1a) at moderate pressures of 1–50 bar. Taking advantage of the results of two preceding publications, a decomposition mechanism with six reactions is recommended, which was validated by measured formyl radical profiles and numerous literature experimental observations. The mechanism is capable of a reliable prediction of almost all formaldehyde pyrolysis literature data, including CH₂O, CO, and H atom measurements at temperatures of 1200–3200 K, with mixtures of 7 ppm to 5% formaldehyde, and pressures up to 15 bar. Some evidence was found for a self‐reaction of two CH₂O molecules. At high initial CH₂O mole fractions the reverse of reaction (6), CH₂OH + HCO ? CH₂O + CH₂O becomes noticeable. The rate of the forward reaction was roughly measured to be k₆ = 1.5 × 10¹³ cm³ mol^?1 s^?1. © 2004 Wiley Periodicals, Inc. Int J Chem Kinet 36: 157–169 2004     

Kinetics of ethane oxidation

Ethane oxidation in jet-stirred reactor has recently been investigated at high temperature (800–1200 K) in the pressure range 1–10 atm and molecular species (H₂, CO, CO₂, CH₄, C₂H₂, C₂H₄, C₂H₆) concentration profiles were obtained by probe sampling and GC analysis. Ethane oxidation was modeled using a comprehensive kinetic reaction mechanism including the most recent findings concerning the kinetics of the reactions involved in the oxidation of C₁? C₄ hydrocarbons. The proposed mechanism is able to reproduce experimental data obtained in our high-pressure jet stirred reactor and ignition delay times measured in shock tube in the pressure range 1–13 atm, for temperatures extending from 800 to 2000 K and equivalence ratios of 0.1 to 2. It is also able to reproduce atoms concentrations (H,O) measured in shock tube at ≈2 atm. The same detailed kinetic mechanism can also be used to model the oxidation of methane, ethylene, propyne, and allene in similar conditions.