A Combined Experimental and Modeling Study on Isopropyl Nitrate Pyrolysis.

Alkyl nitrates thermally decompose by homolytic cleavage of the weak nitrate bond at very low temperatures (e.g., around 500 K at reaction times of a few seconds). This provides the opportunity to study the subsequent chemistry of the initially formed radical (or its subsequent pyrolysis products, if unstable) and nitrogen dioxide at such mild conditions. In this work this idea is applied to isopropyl nitrate (iPN) pyrolysis, which is studied in a tubular reactor at atmospheric pressure, temperatures ranging from 373 to 773 K, and residence times of around 2 s. At the experimental conditions, iPN decomposition starts at 473 K with O-N bond fission producing isopropoxy radical (i-C3H7O) and NO2. i-C3H7O is rapidly converted to acetaldehyde (CH3CHO), which is the most abundant product detected, and methyl radicals. Other major products detected are formaldehyde (CH2O), methanol (CH3OH), nitromethane (CH3NO2), NO, methane, formamide (CHONH2), and methyl nitrite (CH3ONO). Four literature nitrogen chemistry models─three of those augmented with iPN specific reactions─have been tested for their ability to predict the iPN decomposition and product profiles. The mechanism by the Curran group performs best, but it still underpredicts the observed high formaldehyde and methanol yields. A rate analysis indicates that the branching ratio of the reaction between methyl radicals and nitrogen dioxide is of significant importance. Based on recent theoretical and experimental data, new rate expressions for the two reactions CH3 + NO2 → CH3O + NO and CH3 + NO2 + He → CH3ONO2 + He are calculated and incorporated in the kinetic models. It is shown that this change clearly improves the predictions, although additional work is needed to achieve good agreement between calculated and measured species profiles.

Product Identification in the Low-Temperature Oxidation of Cyclohexane Using a Jet-Stirred Reactor in Combination with SVUV-PEPICO Analysis and Theoretical Quantum Calculations.

Cyclohexane oxidation chemistry was investigated using a near-atmospheric pressure jet-stirred reactor at T = 570 K and equivalence ratio ϕ = 0.8. Numerous intermediates including hydroperoxides and highly oxygenated molecules were detected using synchrotron vacuum ultraviolet photoelectron photoion coincidence spectroscopy. Supported by high-level quantum calculations, the analysis of photoelectron spectra allowed the firm identification of molecular species formed during the oxidation of cyclohexane. Besides, this work validates recently published gas chromatography and synchrotron vacuum ultraviolet photoionization mass spectrometry data. Unambiguous detection of characteristic hydroperoxides (e.g., γ-ketohydroperoxides) and their respective decomposition products provides support for the conventional O2 addition channels up to the third addition and their relative contribution to the cyclohexane oxidation. The results were also compared with the predictions of a recently proposed new detailed kinetic model of cyclohexane oxidation. Most of the predictions are in line with the current experimental findings, highlighting the robustness of the kinetic model. However, the analysis of the recorded slow photoelectron spectra indicating the possible presence of C5 species in the kinetic model provides hints that the substituted cyclopentyl radicals from cyclohexyl ring opening might play a minor role in cyclohexane oxidation. Potentially important missing reactions are also discussed.

CO2 adsorption on pyrolysis char from protein-containing livestock waste: How do proteins affect?

Biogas generation through anaerobic digestion provides an interesting opportunity to valorize some types of animal waste materials whose management is increasingly complicated by legal and environmental restrictions. To successfully expand anaerobic digestion in livestock areas, operational issues such as digestate management must be addressed in an economical and environmentally sustainable way. Biogas upgrading is another necessary stage before intending it to add-value applications. The high concentration of CO2 in biogas results in a reduced caloric value, so the removal of CO2 would be beneficial for most end-users. The current work evaluates the CO2 uptake properties (thermogravimetry study) of low-cost adsorbent materials produced from the animal wastes generated in the livestock area itself, specifically via pyrolysis of poorly biodegradable materials, such as meat and bone meal, and the digestate from manure anaerobic digestion. Therefore, the new element in this study with respect to other studies found in the literature related to biochar-based CO2 adsorption performance is the presence of high content of pyrolyzed proteins in the adsorbent material. In this work, pyrolyzed chars from both meat and bone meal and co-digested manure have been proven to adsorb CO2 reversibly, and also the chars produced from their representative pure proteins (collagen and soybean protein), which were evaluated as model compounds for a better understanding of the individual performance of proteins. The ultra-microporosity developed in the protein chars during pyrolysis seems to be the main explanation for such CO2 uptake capacities, while neither the BET surface area nor N-functionalities on the char surface can properly explain the observed results. Although the CO2 adsorption capacities of these pristine chars (6-41.0 mg CO2/g char) are far away from data of commercially activated carbons (~80 mg CO2/g char), this application opens a new via to integrate and valorize these wastes in the circular economy of the primary sector.

Accounting for molecular flexibility in photoionization: case of tert-butyl hydroperoxide.

tert-Butyl hydroperoxide (tBuOOH) is a common intermediate in the oxidation of organic compounds that needs to be accurately quantified in complex gas mixtures for the development of chemical kinetic models of low temperature combustion. This work presents a combined theoretical and experimental investigation on the synchrotron-based VUV single photon ionization of gas-phase tBuOOH in the 9.0 - 11.0 eV energy range, including dissociative ionization processes. Computations consist of the determination of the structures, vibrational frequencies and the energetics of neutral and ionic tBuOOH. The Franck-Condon spectrum for the tBuOOH+ (X+) + e- ← tBuOOH (X) + hν transition is computed, where special treatment is undertaken because of the flexibility of tBuOOH, in particular regarding the OOH group. Through comparison of the experimental mass-selected threshold photoelectron spectra with explicitly correlated coupled cluster calculations and Franck-Condon simulations that account for the flexibility of the molecule, an estimation of the ionization energy is given. The appearance energy of the only fragment observed within the above-mentioned energy range, identified as the tert-butyl C4H9+, is also reported. Finally, the signal branching ratio between the parent and the fragment ions is provided as a function of photon energy, essential to quantify tBuOOH in gas-phase oxidation/combustion experiments via advanced mass spectrometry techniques.

Pyrolysis and Combustion Chemistry of Pyrrole, a Reference Component for Bio-oil Surrogates: Jet-Stirred Reactor Experiments and Kinetic Modeling.

Fast-pyrolysis bio-oils (FPBOs) obtained from lignocellulosic biomass are gaining attention as sustainable fuels for various applications, including the transport sector and power production. A significant fraction of bio-oils is constituted by nitrogen-containing compounds (N fuels) that should be considered when developing surrogate models for FPBOs. Moreover, the content of N fuels in FPBOs is expected to strongly contribute to the production of nitrogen oxides (NO x ) directly from fuel-bound nitrogen (fuel NO x ), in addition to the thermal NO x formation pathways typical of high-temperature combustion conditions. This work investigates the pyrolysis and combustion chemistry of pyrrole (C4H5N), a candidate reference fuel component for FPBO surrogate models. Speciation measurements in an atmospheric pressure jet-stirred reactor have been performed for both pyrolysis and oxidation conditions. Pyrolysis experiments have been performed for 1% pyrrole/helium mixtures over the temperature range T = 925-1200 K. Oxidation experiments were carried out for 1% pyrrole/oxygen/helium mixtures at three equivalence ratios (φ = 0.5, 1.0, and 2.0) over the temperature range T = 700-1200 K. These new data significantly extend the number of experimental targets for kinetic model validation available at present for pyrrole combustion. After a thorough revision of previous theoretical and kinetic modeling studies, a preliminary kinetic model is developed and validated by means of comparison to new experimental data and those previously reported in the literature. The rate of production and sensitivity analyses highlight important pathways deserving further investigations for a better understanding of pyrrole and, more in general, N fuel combustion chemistry. A critical discussion on experimental challenges to be faced when dealing with pyrrole is also reported, encouraging further experimental investigation with advanced diagnostics.

Biomass pyrolysis: thermal decomposition mechanisms of furfural and benzaldehyde.

The thermal decompositions of furfural and benzaldehyde have been studied in a heated microtubular flow reactor. The pyrolysis experiments were carried out by passing a dilute mixture of the aromatic aldehydes (roughly 0.1%-1%) entrained in a stream of buffer gas (either He or Ar) through a pulsed, heated SiC reactor that is 2-3 cm long and 1 mm in diameter. Typical pressures in the reactor are 75-150 Torr with the SiC tube wall temperature in the range of 1200-1800 K. Characteristic residence times in the reactor are 100-200 µsec after which the gas mixture emerges as a skimmed molecular beam at a pressure of approximately 10 µTorr. Products were detected using matrix infrared absorption spectroscopy, 118.2 nm (10.487 eV) photoionization mass spectroscopy and resonance enhanced multiphoton ionization. The initial steps in the thermal decomposition of furfural and benzaldehyde have been identified. Furfural undergoes unimolecular decomposition to furan + CO: C4H3O-CHO (+ M) → CO + C4H4O. Sequential decomposition of furan leads to the production of HC≡CH, CH2CO, CH3C≡CH, CO, HCCCH2, and H atoms. In contrast, benzaldehyde resists decomposition until higher temperatures when it fragments to phenyl radical plus H atoms and CO: C6H5CHO (+ M) → C6H5CO + H → C6H5 + CO + H. The H atoms trigger a chain reaction by attacking C6H5CHO: H + C6H5CHO → [C6H6CHO]* → C6H6 + CO + H. The net result is the decomposition of benzaldehyde to produce benzene and CO.

Rate rules, branching ratios, and pressure dependence of the HO2 + olefin addition channels.

In this work, we present high-pressure rate rules and branching ratios for the addition of HO2 to olefins through the concerted addition channel to form an alkyl peroxy radical (HO2 + olefin → RO2) and through the radical addition channel to form a ß-hydroperoxy alkyl radical (HO2 + olefin → ß-QOOH). These rate rules were developed by calculating rate constants for a series of addition reactions involving olefins with varying degrees of branching. The individual rate expressions were determined from electronic structure calculations performed at the CBS-QB3 level of theory combined with TST calculations. The calculated rate constants were found to be in good agreement with those reported in the literature. Next, we calculated apparent pressure- and temperature-dependent rate constants for HO2 addition to the terminal site of 1-butene using an energy-grained master equation (ME) approach and QRRK calculations with a modified strong collision (MSC) approximation. The two methods gave similar results for both reaction classes. We found that, for the radical addition reaction, the high-pressure limit for the stabilization channel is not reached until unusually high pressures (>1000 atm). Instead, this reaction leads to the direct formation of an oxirane + OH. In general, the results for the major channels are in reasonable agreement with prior theoretical and experimental data. Finally, to explicitly examine the effect of pressure, we compared concentration-time profiles for the reactions of HO2 plus butene in air that were obtained using both high-pressure and pressure-dependent mechanisms at 10 and 100 atm. These simulations showed that, contrary to general expectations, the manifestation of pressure falloff effects in kinetic modeling studies might be more prevalent at increasing pressures. This behavior is attributed to the reaction of ß-QOOH with O2, the rate of which increases with increasing pressure of air. This bimolecular reaction competes with the unimolecular reactions of ß-QOOH under conditions where falloff effects are important for that channel.

High-pressure rate rules for alkyl + O2 reactions. 2. The isomerization, cyclic ether formation, and ß-scission reactions of hydroperoxy alkyl radicals.

The unimolecular reactions of hydroperoxy alkyl radicals (QOOH) play a central role in the low-temperature oxidation of hydrocarbons as they compete with the addition of a second O(2) molecule, which is known to provide chain-branching. In this work we present high-pressure rate estimation rules for the most important unimolecular reactions of the ß-, γ-, and δ-QOOH radicals: isomerization to RO(2), cyclic ether formation, and selected ß-scission reactions. These rate rules are derived from high-pressure rate constants for a series of reactions of a given reaction class. The individual rate expressions are determined from CBS-QB3 electronic structure calculations combined with canonical transition state theory calculations. Next we use the rate rules, along with previously published rate estimation rules for the reactions of alkyl peroxy radicals (RO(2)), to investigate the potential impact of falloff effects in combustion/ignition kinetic modeling. Pressure effects are examined for the reaction of n-butyl radical with O(2) by comparison of concentration versus time profiles that were obtained using two mechanisms at 10 atm: one that contains pressure-dependent rate constants that are obtained from a QRRK/MSC analysis and another that only contains high-pressure rate expressions. These simulations reveal that under most conditions relevant to combustion/ignition problems, the high-pressure rate rules can be used directly to describe the reactions of RO(2) and QOOH. For the same conditions, we also address whether the various isomers equilibrate during reaction. These results indicate that equilibrium is established between the alkyl, RO(2), and γ- and δ-QOOH radicals.

Unimolecular thermal decomposition of phenol and d5-phenol: direct observation of cyclopentadiene formation via cyclohexadienone.

The pyrolyses of phenol and d(5)-phenol (C(6)H(5)OH and C(6)D(5)OH) have been studied using a high temperature, microtubular (µtubular) SiC reactor. Product detection is via both photon ionization (10.487 eV) time-of-flight mass spectrometry and matrix isolation infrared spectroscopy. Gas exiting the heated reactor (375 K-1575 K) is subject to a free expansion after a residence time in the µtubular reactor of approximately 50-100 µs. The expansion from the reactor into vacuum rapidly cools the gas mixture and allows the detection of radicals and other highly reactive intermediates. We find that the initial decomposition steps at the onset of phenol pyrolysis are enol/keto tautomerization to form cyclohexadienone followed by decarbonylation to produce cyclopentadiene; C(6)H(5)OH → c-C(6)H(6) = O → c-C(5)H(6) + CO. The cyclopentadiene loses a H atom to generate the cyclopentadienyl radical which further decomposes to acetylene and propargyl radical; c-C(5)H(6) → c-C(5)H(5) + H → HC≡CH + HCCCH(2). At higher temperatures, hydrogen loss from the PhO-H group to form phenoxy radical followed by CO ejection to generate the cyclopentadienyl radical likely contributes to the product distribution; C(6)H(5)O-H → C(6)H(5)O + H → c-C(5)H(5) + CO. The direct decarbonylation reaction remains an important channel in the thermal decomposition mechanisms of the dihydroxybenzenes. Both catechol (o-HO-C(6)H(4)-OH) and hydroquinone (p-HO-C(6)H(4)-OH) are shown to undergo decarbonylation at the onset of pyrolysis to form hydroxycyclopentadiene. In the case of catechol, we observe that water loss is also an important decomposition channel at the onset of pyrolysis.

High-pressure rate rules for alkyl + O2 reactions. 1. The dissociation, concerted elimination, and isomerization channels of the alkyl peroxy radical.

The reactions of alkyl peroxy radicals (RO(2)) play a central role in the low-temperature oxidation of hydrocarbons. In this work, we present high-pressure rate estimation rules for the dissociation, concerted elimination, and isomerization reactions of RO(2). These rate rules are derived from a systematic investigation of sets of reactions within a given reaction class using electronic structure calculations performed at the CBS-QB3 level of theory. The rate constants for the dissociation reactions are obtained from calculated equilibrium constants and a literature review of experimental rate constants for the reverse association reactions. For the concerted elimination and isomerization channels, rate constants are calculated using canonical transition state theory. To determine if the high-pressure rate expressions from this work can directly be used in ignition models, we use the QRRK/MSC method to calculate apparent pressure and temperature dependent rate constants for representative reactions of small, medium, and large alkyl radicals with O(2). A comparison of concentration versus time profiles obtained using either the pressure dependent rate constants or the corresponding high-pressure values reveals that under most conditions relevant to combustion/ignition problems, the high-pressure rate rules can be used directly to describe the reactions of RO(2).

Kinetic modeling of ethane pyrolysis at high conversion.

The primary objective of this study is to develop an improved first-principle-based mechanism that describes the molecular weight growth kinetics observed during ethane pyrolysis. A proper characterization of the kinetics of ethane pyrolysis is a prerequisite for any analysis of hydrocarbon pyrolysis and oxidation. Flow reactor experiments were performed with ~50/50 ethane/nitrogen mixtures with temperatures ranging from 550 to 850 °C at an absolute pressure of ~0.8 atm and a residence time of ~5 s. These conditions result in ethane conversions ranging from virtually no reaction to ~90%. Comparisons of predictions using our original mechanism to these data yielded very satisfactory results in terms of the temperature dependence of ethane conversion and prediction of the major products ethylene and hydrogen. However, there were discrepancies in some of the minor species concentrations that are involved in the molecular weight growth kinetics. We performed a series of CBS-QB3 analyses for the C(3)H(7), C(4)H(7), and C(4)H(9) potential energy surfaces to better characterize the radical addition reactions that lead to molecular weight growth. We also extended a published C(6)H(9) PES to include addition of vinyl to butadiene. The results were then used to calculate pressure-dependent rate constants for the multiple reaction pathways of these addition reactions. Inclusion of the unadjusted rate constants resulting from these analyses in the mechanism significantly improved the description of several of the species involved in molecular weight growth kinetics. We compare the predictions of this improved model to those obtained with a consensus model recently published as well as to ethane steam cracking data. We find that a particularly important reaction is that of vinyl addition to butadiene. Another important observation is that several radical addition reactions are partially equilibrated. Not only does this mean that reliable thermodynamic parameters are essential for an accurate model, but also that the reaction set describing molecular weight growth chemistry must include a final product that is sufficiently stable to shift the equilibrium toward this product despite the decrease in entropy that accompanies molecular weight growth. Another reaction, H addition to olefins, was found to inhibit molecular weight growth by leading to the production of a lower olefin plus methyl radicals.

Direct detection of products from the pyrolysis of 2-phenethyl phenyl ether.

The pyrolysis of 2-phenethyl phenyl ether (PPE, C(6)H(5)C(2)H(4)OC(6)H(5)) in a hyperthermal nozzle (300-1350 °C) was studied to determine the importance of concerted and homolytic unimolecular decomposition pathways. Short residence times (<100 µs) and low concentrations in this reactor allowed the direct detection of the initial reaction products from thermolysis. Reactants, radicals, and most products were detected with photoionization (10.5 eV) time-of-flight mass spectrometry (PIMS). Detection of phenoxy radical, cyclopentadienyl radical, benzyl radical, and benzene suggest the formation of product by the homolytic scission of the C(6)H(5)C(2)H(4)-OC(6)H(5) and C(6)H(5)CH(2)-CH(2)OC(6)H(5) bonds. The detection of phenol and styrene suggests decomposition by a concerted reaction mechanism. Phenyl ethyl ether (PEE, C(6)H(5)OC(2)H(5)) pyrolysis was also studied using PIMS and using cryogenic matrix-isolated infrared spectroscopy (matrix-IR). The results for PEE also indicate the presence of both homolytic bond breaking and concerted decomposition reactions. Quantum mechanical calculations using CBS-QB3 were conducted, and the results were used with transition state theory (TST) to estimate the rate constants for the different reaction pathways. The results are consistent with the experimental measurements and suggest that the concerted retro-ene and Maccoll reactions are dominant at low temperatures (below 1000 °C), whereas the contribution of the C(6)H(5)C(2)H(4)-OC(6)H(5) homolytic bond scission reaction increases at higher temperatures (above 1000 °C).

Detailed modeling of low-temperature propane oxidation: 1. The role of the propyl + O(2) reaction.

Accurate description of reactions between propyl radicals and molecular oxygen is an essential prerequisite for modeling of low-temperature propane oxidation because their multiple reaction pathways either accelerate the oxidation process via chain branching or inhibit it by forming relatively stable products. The CBS-QB3 level of theory was used to construct potential energy surfaces for n-C(3)H(7) + O(2) and i-C(3)H(7) + O(2). High-pressure rate constants were calculated using transition state theory with corrections for tunneling and hindered rotations. These results were used to derive pressure- and temperature-dependent rate constants for the various channels of these reactions under the framework of the Quantum Rice-Ramsperger-Kassel (QRRK) and the modified strong collision (MSC) theories. This procedure resulted in a thermodynamically consistent C(3)H(7) + O(2) submechanism, which was either used directly or as part of a larger extended detailed kinetic mechanism to predict the loss of propyl and the product yields of propylene and HO(2) over a wide range of temperatures, pressures, and residence times. The overall good agreement between predicted and experimental data suggests that this reaction subset is reliable and should be able to properly account for the reactions of propyl radicals with O(2) in propane oxidation. It is also demonstrated that for most conditions of practical interest only a small subset of reactions (e.g., isomerization, concerted elimination of HO(2), and stabilization) controls the oxidation kinetics, which makes it possible to considerably simplify the mechanism. Moreover, we observed strong similarities in the rate coefficients within each reaction class, suggesting the potential for development of relatively simple rate constant estimation rules that could be applied to analogous reactions involving hydrocarbon radicals that are too large to allow accurate detailed electronic structure calculations.

Selective removal of ethylene, a deposit precursor, from a "dirty" synthesis gas stream via gas-phase partial oxidation.

A fundamental issue in the gasification of biomass is that in addition to the desired synthesis gas product (a mixture of H(2) and CO), the gasifier effluent contains other undesirable products that need to be removed before any further downstream processing can occur. This work assesses the potential to selectively remove hydrocarbons from a synthesis gas stream via gas-phase partial oxidation. Specifically, the partial oxidation of methane-doped, ethylene-doped, and methane/ethylene-doped model synthesis gas mixtures has been investigated at ambient pressures over a temperature range of 760-910 degrees C and at residence times ranging from 0.4 to 2.4 s using a tubular flow reactor. For the synthesis gas mixtures that contain either methane or ethylene, the addition of oxygen substantially reduces the hydrocarbon concentration while only a small reduction in the hydrogen concentration is observed. For the synthesis gas mixtures doped with both methane and ethylene, the addition of oxygen preferentially removes ethylene while the concentrations of methane and hydrogen remain relatively unaffected. These results are compared to the predictions of a plug flow model using a reaction mechanism that is designed to describe the pyrolysis and partial oxidation of small hydrocarbon species. The agreement between the experimental observations and the model predictions is quite good, allowing us to explore the underlying chemistry that leads to the hydrocarbon selective oxidation. The implications of these results are briefly discussed in terms of using synthesis gas to produce liquid fuels and electrical power via a solid oxide fuel cell.

Rate constant rules for the automated generation of gas-phase reaction mechanisms.

The capability of kinetic models to predict complex chemical systems has enormously improved in the last decades, making them an increasingly important tool for process development and optimization. Extension of these approaches to even more complex systems is hindered not only by the geometrically increasing number of reactions and species to be considered but also by the necessity of assigning accurate rate constants to all of the reactions. The recent developments in automated mechanism generators can address the tedious bookkeeping issues. The requirement for development of accurate rate constant estimates remains the job of the kineticist. This task has been aided immeasurably by the combined advances in electronic structure methods and computer performance. This article describes two areas of rate estimation. First, we discuss the development of H abstraction rate estimates from C-H bonds in alkanes, cycloalkanes, and allylic systems by H atoms and point to a surprising result found for cyclopentane. Second, we briefly review our investigation of the ethyl + O(2) reaction and demonstrate the suitability of the QRRK/MSC approach for automated mechanism generation. We conclude with some suggestions for future work in this area.

Detailed modeling of the reaction of C2H5 + O2.

Modeling of low-temperature ethane oxidation requires an accurate description of the reaction of C(2)H(5) + O(2), because its multiple reaction channels either accelerate the oxidation process via chain branching, or inhibit it by forming stable, less reactive products. We have used a steady-state chemical-activation analysis to generate pressure and temperature dependent rate coefficients for the various channels of this system. Input parameters for this analysis were obtained from ab initio calculations at the CBS-QB3 level of theory with bond-additivity corrections, followed by transition state theory calculations with Wigner tunneling corrections. The chemical-activation analysis used QRRK theory to determine k(E) and the modified strong collision (MSC) model to account for collisional deactivation. This procedure resulted in a C(2)H(5) + O(2) submechanism which was either used directly (possibly augmented with a few C(2)H(5) generating and consuming reactions) or as part of a larger extended mechanism to predict the temperature and pressure dependencies of the overall loss of ethyl and of the yields of ethylene, ethylene oxide, HO(2), and OH. A comparison of the predictions using both mechanisms allowed an assessment of the sensitivity of the experimental data to secondary reactions. Except for the time dependent OH profiles, the predictions using the extended mechanism were in good agreement with the observations. By replacing the MSC model with master equation approaches, both steady-state and time dependent, it was confirmed that the MSC assumption is adequate for the analysis of the C(2)H(5) + O(2) reaction. The good overall performance of the C(2)H(5) + O(2) submechanism developed in this study suggests that it provides a good building block for an ethane oxidation mechanism.