1
|
Amiri V, Asatryan R, Swihart M. Automated Generation of a Compact Chemical Kinetic Model for n-Pentane Combustion. ACS OMEGA 2023; 8:49098-49114. [PMID: 38162756 PMCID: PMC10753700 DOI: 10.1021/acsomega.3c07079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 11/21/2023] [Accepted: 11/24/2023] [Indexed: 01/03/2024]
Abstract
We have employed automated mechanism generation tools to construct a detailed chemical kinetic model for combustion of n-pentane, as a step toward the generation of compact kinetic models for larger alkanes. Pentane is of particular interest as a prototype for combustion of alkanes and as the smallest paraffin employed as a hybrid rocket fuel, albeit at cryogenic conditions. A reaction mechanism for pentane combustion thus provides a foundation for modeling combustion of extra-large alkanes (paraffins) that are of more practical interest as hybrid rocket fuels, for which manual construction becomes infeasible. Here, an n-pentane combustion kinetic model is developed using the open-source software package Reaction Mechanism Generator (RMG). The model was generated and tested across a range of temperatures (650 to 1350 K) and equivalence ratios (0.5, 1.0, 2.0) at pressures of 1 and 10 atm. Available thermodynamic and kinetic databases were incorporated wherever possible. Predictions using the mechanism were validated against the published laminar burning velocities (Su) and ignition delay times (IDT) of n-pentane. To improve the model performance, a comprehensive analysis, including reaction path and sensitivity analyses of n-pentane oxidation, was performed, leading us to modify the thermochemistry and rate parameters for a few key species and reactions. These were combined as a separate data set, an RMG library, that was then used during mechanism generation. The resulting model predicted IDTs as accurately as the best manually constructed mechanisms, while remaining much more compact. It predicted flame speeds to within 10% of published experimental results. The degree of success of automated mechanism generation for this case suggests that it can be extended to construct reliable and compact models for combustion of larger n-alkanes, particularly when using this mechanism as a seed submodel.
Collapse
Affiliation(s)
- Venus Amiri
- Department of Chemical and Biological
Engineering, University at Buffalo, The
State University of New York, Buffalo, New York 14260, United States
| | - Rubik Asatryan
- Department of Chemical and Biological
Engineering, University at Buffalo, The
State University of New York, Buffalo, New York 14260, United States
| | - Mark Swihart
- Department of Chemical and Biological
Engineering, University at Buffalo, The
State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
2
|
Zádor J, Martí C, Van de Vijver R, Johansen SL, Yang Y, Michelsen HA, Najm HN. Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces. J Phys Chem A 2023; 127:565-588. [PMID: 36607817 DOI: 10.1021/acs.jpca.2c06558] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.
Collapse
Affiliation(s)
- Judit Zádor
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Carles Martí
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | | | - Sommer L Johansen
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Yoona Yang
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Hope A Michelsen
- Department of Mechanical Engineering, University of Colorado Boulder, Boulder80309, Colorado, United States
| | - Habib N Najm
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| |
Collapse
|
3
|
Automatic construction of transition states and on-the-fly accurate kinetic calculations for reaction classes in automated mechanism generators. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112852] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
4
|
Gao CW, Liu M, Green WH. Uncertainty analysis of correlated parameters in automated reaction mechanism generation. INT J CHEM KINET 2020. [DOI: 10.1002/kin.21348] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Connie W. Gao
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts
| | - Mengjie Liu
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts
| | - William H. Green
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts
| |
Collapse
|
5
|
Vernuccio S, Broadbelt LJ. Discerning complex reaction networks using automated generators. AIChE J 2019. [DOI: 10.1002/aic.16663] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Sergio Vernuccio
- Department of Chemical and Biological Engineering Northwestern University Evanston Illinois
| | - Linda J. Broadbelt
- Department of Chemical and Biological Engineering Northwestern University Evanston Illinois
| |
Collapse
|
6
|
Sun X, Zong W, Wang J, Li Z, Li X. Pressure-dependent rate rules for cycloaddition, intramolecular H-shift, and concerted elimination reactions of alkenyl peroxy radicals at low temperature. Phys Chem Chem Phys 2019; 21:10693-10705. [PMID: 31086861 DOI: 10.1039/c9cp01207a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reactions of cycloaddition, intramolecular H-shift and concerted elimination of alkenyl peroxy radicals are three kinds of important reactions in the low temperature combustion of alkenes. In this study, the cycloaddition reactions are divided into classes considering endo-cycloaddition, exo-cycloaddition and the size of the transition states; the intramolecular H-shift reactions are divided into classes depending upon the ring size of the transition states and the type of C-H bonds from which the hydrogen atom is transferred; the concerted elimination reactions are divided into classes according to the type of H-CβCαOO bond that is broken. All geometry optimizations are performed at the B3LYP/6-31G(2df,p) level. With the electronic structure calculations being performed using the composite Gaussian-4 (G4) method, high pressure limit rate constants and pressure-dependent rate constants at pressures varying from 0.01 to 100 atm are calculated by using canonical transition state theory and the Rice-Ramsberger-Kassel-Marcus/master equation method, respectively. All rate constants are given in the form of the modified Arrhenius expression. The high pressure limit rate rules and the pressure-dependent rate rules are derived by averaging the rate constants of a representative set of reactions in each class. The results show that the rate rules for these three classes of reactions have a large uncertainty and the impact of the pressure on the rate constants increases as temperature increases.
Collapse
Affiliation(s)
- Xiaohui Sun
- School of Chemical Engineering, Sichuan University, Chengdu, Sichuan 610065, China.
| | - Wengang Zong
- School of Chemical Engineering, Sichuan University, Chengdu, Sichuan 610065, China.
| | - Jingbo Wang
- School of Chemical Engineering, Sichuan University, Chengdu, Sichuan 610065, China.
| | - Zerong Li
- College of Chemistry, Sichuan University, Chengdu, Sichuan 610064, China.
| | - Xiangyuan Li
- School of Chemical Engineering, Sichuan University, Chengdu, Sichuan 610065, China.
| |
Collapse
|
7
|
Döntgen M, Schmalz F, Kopp WA, Kröger LC, Leonhard K. Automated Chemical Kinetic Modeling via Hybrid Reactive Molecular Dynamics and Quantum Chemistry Simulations. J Chem Inf Model 2018; 58:1343-1355. [PMID: 29898359 DOI: 10.1021/acs.jcim.8b00078] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An automated scheme for obtaining chemical kinetic models from scratch using reactive molecular dynamics and quantum chemistry simulations is presented. This methodology combines the phase space sampling of reactive molecular dynamics with the thermochemistry and kinetics prediction capabilities of quantum mechanics. This scheme provides the NASA polynomial and modified Arrhenius equation parameters for all species and reactions that are observed during the simulation and supplies them in the ChemKin format. The ab initio level of theory for predictions is easily exchangeable, and the presently used G3MP2 level of theory is found to reliably reproduce hydrogen and methane oxidation thermochemistry and kinetics data. Chemical kinetic models obtained with this approach are ready to use for, e.g., ignition delay time simulations, as shown for hydrogen combustion. The presented extension of the ChemTraYzer approach can be used as a basis for methodological advancement of chemical kinetic modeling schemes and as a black-box approach to generate chemical kinetic models.
Collapse
Affiliation(s)
- Malte Döntgen
- Chair of Technical Thermodynamics , RWTH Aachen University , 52062 Aachen , Germany.,Molecular Science, Department of Chemistry , University of Helsinki , 00560 Helsinki , Finland
| | - Felix Schmalz
- Chair of Technical Thermodynamics , RWTH Aachen University , 52062 Aachen , Germany
| | - Wassja A Kopp
- Chair of Technical Thermodynamics , RWTH Aachen University , 52062 Aachen , Germany
| | - Leif C Kröger
- Chair of Technical Thermodynamics , RWTH Aachen University , 52062 Aachen , Germany
| | - Kai Leonhard
- Chair of Technical Thermodynamics , RWTH Aachen University , 52062 Aachen , Germany
| |
Collapse
|
8
|
Vereecken L, Aumont B, Barnes I, Bozzelli J, Goldman M, Green W, Madronich S, Mcgillen M, Mellouki A, Orlando J, Picquet-Varrault B, Rickard A, Stockwell W, Wallington T, Carter W. Perspective on Mechanism Development and Structure-Activity Relationships for Gas-Phase Atmospheric Chemistry. INT J CHEM KINET 2018. [DOI: 10.1002/kin.21172] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- L. Vereecken
- Institute for Energy and Climate Research: IEK-8 Troposphere; Forschungszentrum Jülich GmbH; Jülich Germany
| | - B. Aumont
- Laboratoire Interuniversitaire des Systèmes Atmosphériques (LISA); UMR 7583 CNRS; Universités Paris-Est Créteil et Paris Diderot; Institut Pierre-Simon Laplace; Créteil Cedex France
| | - I. Barnes
- School of Mathematics and Natural Sciences; Physical & Theoretical Chemistry; University of Wuppertal; Wuppertal Germany
| | - J.W. Bozzelli
- Department of Chemistry and Environmental Science; New Jersey Institute of Technology; Newark NJ 07102
| | - M.J. Goldman
- Department of Chemical Engineering; Massachusetts Institute of Technology; Cambridge MA 02139
| | - W.H. Green
- Department of Chemical Engineering; Massachusetts Institute of Technology; Cambridge MA 02139
| | - S. Madronich
- Atmospheric Chemistry Observations and Modeling Laboratory; National Center for Atmospheric Research; Boulder CO 80307
| | - M.R. Mcgillen
- School of Chemistry; University of Bristol; Cantock's Close; Bristol BS8 1TS UK
| | - A. Mellouki
- Institut de Combustion; Aérothermique, Réactivité et Environnement (ICARE); CNRS/OSUC; 45071 Orléans Cedex 2 France
| | - J.J. Orlando
- Atmospheric Chemistry Observations and Modeling Laboratory; National Center for Atmospheric Research; Boulder CO 80307
| | - B. Picquet-Varrault
- Laboratoire Interuniversitaire des Systèmes Atmosphériques (LISA); UMR 7583 CNRS; Universités Paris-Est Créteil et Paris Diderot; Institut Pierre-Simon Laplace; Créteil Cedex France
| | - A.R. Rickard
- Wolfson Atmospheric Chemistry Laboratories; Department of Chemistry; University of York; York YO10 5DD UK
- National Centre for Atmospheric Science; University of York; York YO10 5DD UK
| | - W.R. Stockwell
- Department of Physics; University of Texas at El Paso; El Paso TX 79968 USA
| | - T.J. Wallington
- Research & Advanced Engineering; Ford Motor Company; Dearborn MI 48121-2053
| | - W.P.L. Carter
- College of Engineering; Center for Environmental Research and Technology (CE-CERT); University of California; Riverside CA 92521
| |
Collapse
|
9
|
Bhoorasingh PL, Slakman BL, Seyedzadeh Khanshan F, Cain JY, West RH. Automated Transition State Theory Calculations for High-Throughput Kinetics. J Phys Chem A 2017; 121:6896-6904. [PMID: 28820268 DOI: 10.1021/acs.jpca.7b07361] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A scarcity of known chemical kinetic parameters leads to the use of many reaction rate estimates, which are not always sufficiently accurate, in the construction of detailed kinetic models. To reduce the reliance on these estimates and improve the accuracy of predictive kinetic models, we have developed a high-throughput, fully automated, reaction rate calculation method, AutoTST. The algorithm integrates automated saddle-point geometry search methods and a canonical transition state theory kinetics calculator. The automatically calculated reaction rates compare favorably to existing estimated rates. Comparison against high level theoretical calculations show the new automated method performs better than rate estimates when the estimate is made by a poor analogy. The method will improve by accounting for internal rotor contributions and by improving methods to determine molecular symmetry.
Collapse
Affiliation(s)
- Pierre L Bhoorasingh
- Department of Chemical Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Belinda L Slakman
- Department of Chemical Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Fariba Seyedzadeh Khanshan
- Department of Chemical Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Jason Y Cain
- Department of Chemical Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Richard H West
- Department of Chemical Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| |
Collapse
|
10
|
Vereecken L, Glowacki DR, Pilling MJ. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem Rev 2015; 115:4063-114. [DOI: 10.1021/cr500488p] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Luc Vereecken
- Max Planck Institute for Chemistry, 55128 Mainz, Germany
| | - David R. Glowacki
- PULSE
Institute and Department of Chemistry, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department
of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| | | |
Collapse
|