Multitarget Global Sensitivity Analysis of n-Butanol Combustion

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single-molecule kinetics. Using a stochastic Markov-state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated for model problems and for a physically realistic mechanism for an alkene hydrogenation reaction on a single-atom catalyst. The spectral analysis permits a hierarchy of timescales to be identified from the single-molecule signal, which correspond to specific relaxation modes in the network.

The role of the three body photodissociation channel of water in the evolution of dioxygen in astrophysical applications

A recent experiment at the Dalian Coherent Light Source (DCLS) has provided measurements of the partial cross sections for the photodissociation of water vapor over an unprecedented range of wavelengths in the vacuum ultraviolet (VUV) region. It was found that the three body dissociation channel, H + H + O(3P/1D), becomes prominent at wavelengths shorter than the Lyman α-line at 121.6 nm. The present work explores the kinetic consequences of this discovery for several astrophysically motivated examples. The irradiation of a dilute low-temperature gas by unscreened solar radiation, similar to early stage photochemical processing in a comet coma, shows significant increase in the production of O2-molecules at shorter times, <1 day, that might physically correspond to the photochemical reaction zone of the coma. Several examples of planetary atmospheres show increased O-atom production at high altitudes but relatively little modification of the equilibrium O2 concentrations predicted by conventional models.

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH₃C(O)CH₃) at a nominal temperature of 1800 K. The major product observed is ketene (CH₂CO). Minor products identified include acetaldehyde (CH₃CHO), propyne (CH₃CCH), propene (CH₂CHCH₃), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. The pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

Sparsity Facilitates Chemical-Reaction Selection for Engine Simulations

Analysis of large-scale, realistic models incorporating detailed chemistry can be challenging because each simulation is computationally expensive, and a complete analysis may require many simulations. This paper addresses one such problem of this type, chemical-reaction selection in engine simulations. In this computationally challenging case, it is demonstrated how the important concept of sparsity can facilitate chemical-reaction selection, which is the process of finding the most important chemical reactions for modeling a chemical process. It is difficult to perform accurate reaction selection for engine simulations using realistic models of the chemistry, as each simulation takes processor weeks to complete. We developed a procedure to efficiently accomplish this selection process with a relatively small number of simulations using a form of global sensitivity analysis based on sparse regression. The chemical-reaction selection leads to an analysis of the ignition chemistry as it evolves within the compression-ignition engine simulations and allows for the spatial development of the selected chemical reactions to be studied in detail.

Global Sensitivity Analysis of Large Reaction Mechanisms Using Fourier Amplitude Sensitivity Test

Global sensitivity analysis (GSA) of large chemical reaction mechanisms remains a challenge since the model with uncertainties in the large number of input parameters provides large dimension of input parameter space and tends to be difficult to evaluate the effect of input parameters on model outputs. In this paper, a criterion for frequency selection to input parameter is proposed so that Fourier amplitude sensitivity test (FAST) method can evaluate the complex model with a low sample size. This developed FAST method can establish the relationship between the number of input parameters and sample size needed to measure sensitivity indices with high accuracy. The performance of this FAST method which can allow both the qualitative and quantitative analysis of complex systems is validated by a H2/air combustion model and a CH4/air combustion model. This FAST method is also compared with other GSA methods to illustrate the features of this FAST method. The results show that FAST method can evaluate the reaction systems with low sample size, and the sensitivity indices obtained from the FAST method can provide more important information which the variance-based GSA methods cannot obtain. FAST method can be a remarkably effective tool for the modelling and diagnosis of large chemical reaction.

Global Sensitivity Analysis with Small Sample Sizes: Ordinary Least Squares Approach

A new version of global sensitivity analysis is developed in this paper. This new version coupled with tools from statistics, machine learning, and optimization can devise small sample sizes that allow for the accurate ordering of sensitivity coefficients for the first 10-30 most sensitive chemical reactions in complex chemical-kinetic mechanisms, and is particularly useful for studying the chemistry in realistic devices. A key part of the paper is calibration of these small samples. Because these small sample sizes are developed for use in realistic combustion devices, the calibration is done over the ranges of conditions in such devices, with a test case being the operating conditions of a compression ignition engine studied earlier. Compression-ignition engines operate under low-temperature combustion conditions with quite complicated chemistry making this calibration difficult, leading to the possibility of false positives and false negatives in the ordering of the reactions. So an important aspect of the paper is showing how to handle the trade-off between false positives and false negatives using ideas from the multiobjective optimization literature. The combination of the new global sensitivity method and the calibration are sample sizes a factor of approximately 10 times smaller than were available with our previous algorithm.

Effects of new Ab initio rate coefficients on predictions of species formed during n-butanol ignition and pyrolysis

Experimental, time-resolved species profiles provide critical tests in developing accurate combustion models for biofuels such as n-butanol. A number of such species profiles measured by Karwat et al. [ Karwat, D. M. A.; et al. J. Phys. Chem. A 2011 , 115 , 4909 ] were discordant with predictions from a well-tested chemical kinetic mechanism developed by Black et al. [ Black, G.; et al. Combust. Flame 2010 , 157 , 363 ]. Since then, significant theoretical and experimental efforts have focused on determining the rate coefficients of primary n-butanol consumption pathways in combustion environments, including H atom abstraction reactions from n-butanol by key radicals such as HO2 and OH, as well as the decomposition of the radicals formed by these H atom abstractions. These reactions not only determine the overall reactivity of n-butanol, but also significantly affect the concentrations of intermediate species formed during n-butanol ignition. In this paper we explore the effect of incorporating new ab initio predictions into the Black et al. mechanism on predictions of ignition delay time and species time histories for the experimental conditions studied by Karwat et al. The revised predictions for the intermediate species time histories are in much improved agreement with the measurements, but some discrepancies persist. A rate of production analysis comparing the effects of various modifications to the Black et al. mechanism shows significant changes in the predicted consumption pathways of n-butanol, and of the hydroxybutyl and butoxy radicals formed by H atom abstraction from n-butanol. The predictions from the newly revised mechanism are in very good agreement with the low-pressure n-butanol pyrolysis product species measurements of Stranic et al. [ Stranic, I.; et al. Combust. Flame 2012 , 159 , 3242 ] for all but one species. Importantly, the changes to the Black et al. mechanism show that concentrations of small products from n-butanol pyrolysis are sensitive to different reactions than those presented by Stranic et al.

First analysis of the 1-v″ progression of the Ångström (B1Σ+-A1Π) band system in the rare 13C17O isotopologue

The 1-v″ progression of the Ångström band system, so far unobserved in the rare (13)C(17)O isotopologue, was obtained under high resolution as an emission spectrum using a high accuracy dispersive optical spectroscopy. In the studied region 22,700-24,500 cm(-1), 146 spectral lines were observed, among which 118 were interpreted as belonging to the 1-0 and 1-1 bands of B-A system, and the next 28 were interpreted as extra lines belonging to the 1-1 band of B(1)Σ(+)-e(3)Σ(-) intercombination system, also unobserved in the (13)C(17)O molecule so far. All those lines were precisely measured with an estimated accuracy better than 0.0025 cm(-1), and rotationally analyzed. As a result the following in the (13)C(17)O molecule were calculated for the first time: the merged rotational constants B1 = 1.790227(23) cm(-1), D1 = 6.233(47) × 10(-6) cm(-1), and ΔG1/2 = 2010.9622 (69) cm(-1) and the equilibrium constants, ωe = 2076.04(57) cm(-1), ωexe = 32.54(28) cm(-1), Be = 1.824678(15) cm(-1), αe = 2.2967(24) × 10(-2) cm(-1), De = 5.226(25) × 10(-6) cm(-1), and βe = 6.71(48) × 10(-7) cm(-1) for the B(1)Σ(+) Rydberg state, as well as the individual rotational constant B0 = 1.50485(78) cm(-1), and the equilibrium constants ωe = 1463.340(21) cm(-1), Be = 1.49902(12) cm(-1), αe = 1.7782(49) × 10(-2) cm(-1), De = 7.36(56) × 10(-6) cm(-1) for the A(1)Π state, and σe = 21,854.015(51) cm(-1), RKR turning points, Franck-Condon factors (FCF), relative intensities, and r centroids for the Ångström band system. With the help of the strong and vast A(1)Π (v = 0) ∼ e(3)Σ(-) (v = 1) interaction, the experimental parameters of the e(3)Σ(-) (v = 1) perturbing state were established in the (13)C(17)O molecule for the first time.

Quantum Tunneling Affects Engine Performance

We study the role of individual reaction rates on engine performance, with an emphasis on the contribution of quantum tunneling. It is demonstrated that the effect of quantum tunneling corrections for the reaction HO2 + HO2 = H2O2 + O2 can have a noticeable impact on the performance of a high-fidelity model of a compression-ignition (e.g., diesel) engine, and that an accurate prediction of ignition delay time for the engine model requires an accurate estimation of the tunneling correction for this reaction. The three-dimensional model includes detailed descriptions of the chemistry of a surrogate for a biodiesel fuel, as well as all the features of the engine, such as the liquid fuel spray and turbulence. This study is part of a larger investigation of how the features of the dynamics and potential energy surfaces of key reactions, as well as their reaction rate uncertainties, affect engine performance, and results in these directions are also presented here.