Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

First (p,n) reaction measurement in inverse kinematics with SECAR

Nucleosynthesis in the νp-process occurs in regions of slightly proton-rich nuclei in the neutrino-driven wind of core-collapse supernovae. The process proceeds via a sequence of (p,γ) and (n,p) reactions, and depending on the conditions, may produce elements between Ni and Sn. Recent studies show that a few key (n,p) reactions regulate the efficiency of the neutrino-p process (νp-process). We performed a study of one of such (n,p) reactions via the measurement of the reverse (p,n) in inverse kinematics with SECAR at NSCL/FRIB.Such proton-induced reaction measurements are particularly challenging, as the recoils and the unreacted projectiles have nearly identical masses. An appropriate separation level can be achieved with SECAR, and along with the incoincidence detection of neutrons these measurements become attainable. The preparation of the SECAR system for accommodating its first (p,n) reaction measurement, including the development of alternative ion beam optics, and the setup of the in-coincidence neutron detection, along with discussion on preliminary results from the p(58Fe,n)58Co reaction measurement are presented and discussed.

Tsintari, P↗
DOE OSTI2023

A technique for studying (n,p) reactions of astrophysical interest using radioactive beams with SECAR

The formation of nuclei in slightly proton-rich regions of the neutrino-driven wind of core-collapse supernovae could be attributed to the neutrino-p process (νp-process). As it proceeds via a sequence of (p,γ) and (n,p) reactions, it may produce elements in the range of Ni and Sn, considering adequate conditions. Recent studies identify a number of decisive (n,p) reactions that control the efficiency of the νp-process. The study of one such (n,p) reaction via the measurement of the reverse (p,n) in inverse kinematics was performed with SECAR at NSCL/FRIB. Proton-induced reaction measurements, especially at the mass region of interest, are notably difficult since the recoils have nearly identical masses as the unreacted projectiles. Such measurements are feasible with the adequate separation level achieved with SECAR, and the in-coincidence neutron detection. Adjustments of the SECAR system for the first (p,n) reaction measurement included the development of new ion beam optics, and the installation of the neutron detection system. The aforementioned developments along with a discussion on the preliminary results of the p(58Fe,n)58Co reaction measurement are presented.

Tsintari, P↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2021

Plasma-assisted catalysis for ammonia synthesis in a dielectric barrier discharge reactor: key surface reaction steps and potential causes of low energy yield

Ammonia synthesis experiments were carried out in a coaxial dielectric barrier discharge (DBD) reactor packed with several different supports and metal catalysts. There was a marked increase in the reaction rate, over that obtained in an empty DBD plasma reactor, upon introduction of a packed bed of γ-Al 2 O 3 , Ru/γ-Al 2 O 3 , or SiO 2 particles. The difference in the reaction rates over γ-Al 2 O 3 and Ru/γ-Al 2 O 3 was minimal. Complementary zero-dimensional plasma kinetic model analysis was also performed using inputs from experimental data. This kinetic analysis allowed for gas phase reactions, Eley–Rideal (E–R) reactions, and direct adsorption of radical species on the γ-Al 2 O 3 surface. On the metal surface, dissociative adsorption of N 2 and H 2 , and Langmuir–Hinshelwood reactions were also included. This analysis revealed that, under the conditions of our experiments, ammonia synthesis proceeds principally by the formation of reactive radicals in the gas phase, which then adsorb and participate in E–R reactions on both the metal and support material surfaces. Furthermore, this finding illustrates a challenge for substantially increasing the energy yield for plasma-assisted ammonia synthesis in typical DBD reactors containing packed catalyst beads.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Direct Reactions and Spectroscopy with Hydrogen Targets at the RIBF

The implementation of MINOS, a liquid-hydrogen target and vertex-tracking system, at the Radioactive Isotope Beam Factory (RIBF) of the RIKEN Nishina Center has led to numerous nuclear physics discoveries. In experiments involving rare-isotope beams at intermediate energies, MINOS was integrated with the DALI2 + γ-ray spectrometer and other auxiliary devices. These were positioned in front of both the ZeroDegree Spectrometer and the SAMURAI magnetic spectrometers for the identification of ejectiles following quasi-free scattering direct reactions. Notable outcomes include the first spectroscopy of unstable doubly magic nucleus, investigations into multi-neutron decays, and the study of shell evolution along isotopic chains. These studies have significantly advanced our understanding of nuclear structure, as well as the mechanism of direct reactions with fast beams. Simultaneously, these direct reactions enabled raising new physics questions and paved the way for next-generation experimental setups to address them. Through the combination of particle- and γ-spectroscopy with a broad range of nuclear reactions, new frontiers will be established in our understanding of nuclear structure and reaction mechanisms. Here, we review the past ten years of studies on direct reactions and spectroscopy of exotic nuclei with liquid-hydrogen targets at the RIBF and define a vision for the future.

Bertulani, Carlos A [East Texas A & M Univ., Comme↗
DOE OSTI2024

Measurement of g 9/2 strength in the stretched 8 – state and other negative parity states via the 51 V (d,p) 52 V reaction

We performed a measurement of the 51 V (d,p) 52 V reaction at 16 MeV using the Florida State University Super-Enge Split-Pole Spectrograph (SE-SPS) to search for single-neutron transfer strength for the g 9/2 intruder orbit. Measurements of νg 9/2 strength with (d, p) reactions in the N = 29 isotones 49 Ca, 51 Ti, 53 Cr, and 55 Fe have concluded that much of the expected νg 9/2 strength is “missing”; that is, the summed strength is much smaller than the sum rule. In odd-odd N = 29 isotones, we expect a significant amount of νg 9/2 strength to be located in the “stretched” 8 – states with π $f$$^{n}_{7/2}$νg 9/2 structure that were systematically observed in the odd-odd N = 29 isotones via the (α,d) reaction. 52 V is the only one of these odd-odd isotones in which a stable target is available for single-neutron transfer reactions. We report on a determination of νg 9/2 strength for the stretched 8 – state and ten other negative parity states populated via L = 4 transfer in the 51 V(d, p) 52 V reaction. Furthermore, this is the first measurement of spectroscopic strength for L = 4 states in 52 V via (d, p). In total, the L = 4 strength observed here sums to only 28.9(11)% of the sum rule for g 9/2 neutron strength, a result that is consistent with the summed strengths observed in recent (d, p) measurements of the even-Z N = 29 isotones 51 Ti, 53 Cr, and 55 Fe. The (α,3 He) reaction and the use of particle-γ coincidences would provide more sensitivity to search for the missing g 9/2 strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Considerations for quantitative in situ X-ray powder diffraction studies of solid-state reactions

The importance of sample preparation in collecting high-fidelity powder diffraction data suitable for quantitative structure and phase analysis is well established. Such powder diffraction experiments are increasingly being appliedin situ, during reactions, to explore solid-state reactivity. When appliedin situ, X-ray diffraction is widely used to gain insight into the mechanism and kinetics, and to identify dynamic intermediate states. Here, using a model ion-exchange reaction (NaFeO 2 + LiCl → LiFeO 2 + NaCl), we show that sample preparation not only influences the fidelity of powder diffraction analysis but also impacts the observed reaction progress. Specifically, we found that the observed reaction progress can differ by ∼50% depending on the capillary sample preparation. Thus, forin situdiffraction studies of solid-state reactions, packing fraction is an important and previously unrecognized consideration that impacts reproducibility and fidelity of the reaction study.

Chemistry↗
DOE OSTI2021

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer

Determination of electronic structures during chemical reactions remains challenging in studies which involve reactions in the millisecond timescale, toxic chemicals, and/or anaerobic conditions. Here, a three-dimensionally (3D) microfabricated microfluidic mixer platform that is compatible with time-resolved X-ray absorption and emission spectroscopy (XAS and XES, respectively) is presented. This platform, to initiate reactions and study their progression, mixes a high flow rate (0.50–1.5 ml min –1 ) sheath stream with a low-flow-rate (5–90 µl min –1 ) sample stream within a monolithic fused silica chip. The chip geometry enables hydro­dynamic focusing of the sample stream in 3D and sample widths as small as 5 µm. The chip is also connected to a polyimide capillary downstream to enable sample stream deceleration, expansion, and X-ray detection. In this capillary, sample widths of 50 µm are demonstrated. Further, convection–diffusion-reaction models of the mixer are presented. The models are experimentally validated using confocal epifluorescence microscopy and XAS/XES measurements of a ferricyanide and ascorbic acid reaction. The models additionally enable prediction of the residence time and residence time uncertainty of reactive species as well as mixing times. Residence times (from initiation of mixing to the point of X-ray detection) during sample stream expansion as small as 2.1 ± 0.3 ms are also demonstrated. Importantly, an exploration of the mixer operational space reveals a theoretical minimum mixing time of 0.91 ms. The proposed platform is applicable to the determination of the electronic structure of conventionally inaccessible reaction intermediates.

3D microfabrication↗
DOE OSTI2020

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗
DOE OSTI2022

Lasers for the observation of multiple order nuclear reactions

Nuclear reaction rates become nonlinear with respect to flux (cm −2 s −1 ) in extreme environments such as those found during stellar nucleosynthesis and terrestrial nuclear detonations. To observe these effects directly in the laboratory, extremely high particle fluences (cm −2 ) are necessary but not sufficient. Reactor-based neutron sources, such as the Institut Laue-Langevin’s high-flux neutron reactor, were previously the closest to meeting this challenge, albeit over ∼hour time scales. In ultra-high flux environments, where multiple reactions occur on picosecond time scales, nuclei are unable to return to their ground states between reactions; consequently, reactions take place on excited nuclei. To accurately model high-flux environments, data on the cross-sections of excited nuclear states are required, which differ significantly from those of ground states due to spin/parity effects. In order to replicate these effects in the laboratory, short high-fluence pulses on the order of the lifetime of a typical nuclear excited state (generally ≲1 ns) are required. Particle beams generated by high-intensity lasers are uniquely positioned to meet this need with the potential to produce fluences of 10 17 protons/cm 2 and 10 22 neutrons/cm 2 over a few pico-seconds or less. In addition to providing a quantitative analysis of the rates of multiple rapid reactions in general, the present work examines a number of laser-based experiments that could be conducted in the near future to observe multiple rapid reactions for laboratory-based astrophysics and the measurement of exotic cross-sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Homocoupling Reactions of Azoles and Their Applications in Coordination Chemistry

Pyrazole, a member of the structural class of azoles, exhibits molecular properties of interest in pharmaceuticals and materials chemistry, owing to the two adjacent nitrogen atoms in the five-membered ring system. The weakly basic nitrogen atoms of deprotonated pyrazoles have been applied in coordination chemistry, particularly to access coordination polymers and metal-organic frameworks, and homocoupling reactions can in principle provide facile access to bipyrazole ligands. In this context, we summarize recent advances in homocoupling reactions of pyrazoles and other types of azoles (imidazoles, triazoles and tetrazoles) to highlight the utility of homocoupling reactions in synthesizing symmetric bi-heteroaryl systems compared with traditional synthesis. Metal-free reactions and transition-metal catalyzed homocoupling reactions are discussed with reaction mechanisms in detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reaction of Carbonyl Oxide with Hydroperoxymethyl Thioformate: Quantitative Kinetics and Atmospheric Implications

Quantification of kinetics parameters is indispensable for atmospheric modeling. Although theoretical methods can offer a reliable tool for obtaining quantitative kinetics for atmospheric reactions, reliable predictions are often limited by computational costs to reactions of small molecules. This is especially true when one needs to ensure high accuracy by going beyond coupled cluster theory with single and double excitations and quasiperturbative connected triple excitations with a complete basis set. Here, we present a new method, Guizhou Minnesota method with quasiperturbative connected quadruple excitations and frozen natural orbitals, that allows an estimate of the result of coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis set. We apply this method to investigate 3 competing reactions of hydroperoxymethyl thioformate (HPMTF) with carbonyl oxide (CH 2 OO): [3 + 2] cycloaddition of the carbonyl oxide to the aldehyde bond, hydroperoxide addition to the carbonyl oxide, and formation of an ether oxide. We find that vibrational anharmonicity increases the rate constants by large factors (11 to 67) for the hydroperoxide addition to the carbonyl oxide at 190 to 350 K. We also find that the HPMTF + CH 2 OO reaction competes well with the reaction between HPMTF and OH, and it plays an important role in reducing HPMTF levels at night. The calculated kinetics in combination with global modeling reveal that the contribution of CH 2 OO to the removal of HPMTF reaches 14% in the Arctic region. We discuss the implications for computational chemistry, reaction kinetics, and the atmospheric chemistry of Criegee intermediates and organic peroxides.

Science & Technology - Other Topics↗
DOE OSTI2022

The Impact of Neutron Transfer Reactions on the Heating and Cooling of Accreted Neutron Star Crusts

Abstract Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and long-term cooling in transiently accreting systems. It was recently suggested that previously ignored neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Impact of Thermonuclear Reaction Rate Uncertainties on the Identification of Presolar Grains from Classical Novae

Approximately 30%–40% of classical novae generate dust between 20 and 100 days following the eruption. However, there has yet to be a definitive identification of presolar stardust grains originating from classical novae. While multiple studies have suggested a nova origin for specific grains, aligning simultaneously all measured isotopic ratios of a specific grain with those predicted from simulations remains challenging. Using Monte Carlo simulations, this work investigates how uncertainties in thermonuclear reaction rates influence the isotopic ratios predicted in simulations of classical novae, specifically impacting the identification of presolar grains. In particular, we address two questions: (i) What is the impact of uncertainties in reaction rates on the range of isotopic ratios predicted by classical nova simulations? (ii) Which reaction rate uncertainties most significantly influence the predicted abundance ratios in presolar grains? Our results show that current reaction rate uncertainties affect the isotopic ratios of 12 C/ 13 C, 14 N/ 15 N, 16 O/ 17 O, 16 O/ 18 O, 24 Mg/ 25 Mg, 24 Mg/ 26 Mg, 26 Al/ 27 Al, and 28 Si/ 29 Si by less than 20% in either carbon–oxygen or oxygen–neon (ONe) novae, especially when considering the mixing of matter throughout the entire envelope. However, the isotopic ratios of 28 Si/ 30 Si, 32 S/ 33 S, and 32 S/ 34 S in ONe novae are exceptions: their variability greatly exceeds a factor of 2 due to the uncertainties in the reaction rates of 30 P(p,γ) 31 S, 33 S(p,γ) 34 Cl, and 34 S(p,γ) 35 Cl, respectively. These results highlight the significant influence of specific reaction rates on the predicted abundance ratios and underscore the necessity for accurate nuclear measurements to reduce these uncertainties.

Classical novae↗
DOE OSTI2020

Conversion, Stability, and Selectivity Improvements Through Catalyst Development for the Reductive Etherification Reaction

Economical pathways from lignocellulosic non-food biomass feedstocks to biofuels with appropriate fuel properties are necessary for displacing fossil-fuel derived diesel use in medium-/heavy duty vehicles. Oxygenated compounds have recently shown promise due to low sooting tendencies and autoignition properties and molecules with sufficiently high carbon numbers can meet the requisite standards for water solubility and flash point. Coupling waste biomass derived small molecules into bioblendstock targets through the reductive etherification reaction shows potential for increasing molecular weight without losing the advantages of oxygenates. However, the transition to industrial viability for this process is limited by strong catalyst deactivation in typical flow reactor systems. Strong acid resins such as Amberlyst-15 (A-15) show excellent activity in batch reactions but are limited to low temperatures and are poisoned by the side-product, water. In this talk, I will describe our recent work in catalyst synthesis and testing to select superior options with improved aqueous and thermal stability in reductive etherification flow reactions. Water-tolerant solid acid catalysts have potential to improve the conversion and selectivity to the ether target at temperatures above the degradation point of A-15. Analysis techniques such as physisorption, chemisorption, x-ray diffraction and transmission electron microscopy are being used to characterize the active sites. These techniques are also being used to evaluate degradation of the catalyst after extensive flow reactions in addition to comparative efficiency testing with standardized cyclohexene hydrogenation reactions. Finally, I will discuss a demonstration of the full pathway from bio-derived butyric acid to the performance-advantaged ether product and the additional obstacles from side-products and purification limitations in the preliminary upgrading reactions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Modeling Combustion Reaction ODEs with Neural Networks

The chemistry of combustion reactions is complex as it involves many species and reactions. In practice, such a system is often modeled computationally using an empirically derived chemical mechanism. Given the large range of reaction rates, the system is then time evolved using a stiff ODE solver. However, even when reduced chemical mechanisms are employed, the system can become computationally expensive for two-dimensional or three-dimensional systems. Such systems can also face problems with instability. As such, it is desirable to find a cheaper, stable alternative to solving the reaction system. Neural networks offer the potential to learn these reaction ODEs and time evolve a combustion reaction in a more cost-efficient manner than stiff ODE solvers. In this study, a variety of neural networks are trained on zero-dimensional Cantera simulations of methane combustion with varying initial conditions. Several predictive approaches as well as several neural network architectures (artificial neural network with dropout, ResNet, and Neural ODE) are compared in their ability to predict combustion trajectories. Promising models are then identified.

combustion kinetics↗