Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics and dynamics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

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

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2025

Kinetic Deep Learning v0.1

Here, we present a method that uses protein levels to predict times series of metabolite concentrations. Understanding this type of pathway dynamics is important in order to predict the behavior of the pathway and, more pragmatically, to be able to design biological systems (such as strains bioengineered to produce chemical products) reliably. Typically, for this purpose, kinetic models consisting of differential equations based on the Michaelis-Menten dynamics have been used in the past. However, these methods can rarely produce good fits to measured data time series. Possibly, this happens because the kinetic constants are unknown or are different from the ones measured in vivo, or perhaps because Michaelis-Menten dynamics is not a satisfactory description. In order to improve the predictive nature of these kinetic models we have eliminated the Michaelis-Menten description of pathway dynamics and we have substituted it by algorithms that automatically learn these dynamics from previously obtained metabolomics and proteomics data using machine learning approaches. Specifically, kinetic deep learning uses deep learning to map proteomics time series to metabolite concentration time series, instead of learning the first metabolite derivative and integrating in (as in the first version of kinetic learning). This approach is shown to provide good to excellent results with a data set specifically collected for this purpose.

Garcia Martin, Hector [Joint BioEnergy Institute (↗
DOE OSTI2021

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

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

Atomic dynamics of gas-dependent oxide reducibility

Understanding oxide reduction is critical for advancing metal production, catalysis and energy technologies. Although carbon monoxide (CO) and hydrogen (H 2 ) are widely used reductants, the mechanisms by which they work are often presumed to be similar, both involving lattice oxygen removal. However, because of growing interest in replacing CO with H 2 to lower CO 2 emissions, distinguishing gas-specific reduction pathways is critical. Yet, capturing these atomic-scale processes under reactive gas and high-temperature conditions remains challenging. Here we use environmental transmission electron microscopy, which is capable of real-time, atomic-resolution imaging of gas–solid redox reactions to directly visualize the gas-dependent oxide reduction dynamics in NiO. We show that CO drives surface nucleation and the growth of metallic Ni islands, leading to self-limiting surface metallization. Conversely, H 2 activates a coupled surface-to-bulk transformation, where protons from dissociated H 2 infiltrate the oxide lattice to promote the inward migration of surface-generated oxygen vacancies and enabling bulk metallization. By contrast, oxygen vacancies formed by CO remain confined near the surface, where they rapidly form a metallic Ni layer that inhibits further reduction. Furthermore, these results reveal distinct atomistic pathways for CO and H 2 and provide insights that may guide metallurgical processes and catalyst design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solvated Electrons Have Multiple Personalities in Molten Salts

P-Crosscut: Solvated Electrons Have Multiple Personalities in Molten Salts [EFRC – MSEE] Alejandro Ramos-Ballesteros;1 Hung H. Nguyen;2 Kazuhiro Iwamatsu;3 Santanu Roy;4 Vyacheslav Bryantsev;4 Michael E. Woods;1 Ruchi Gakhar;1 Phillip Halstenberg;4 Bobby Layne;5 Jay A. LaVerne;6 Claudio J. Margulis:2* and James F. Wishart5* 1Idaho National Laboratory; 2The University of Iowa; 3Hunter College; 4Oak Ridge National Laboratory; 5Brookhaven National Laboratory; 6University of Notre Dame Abstract: The solvated (eS–) is a powerful reducing agent and one of the primary products of molten salt radiolysis. In these extreme high-temperature environments, the eS– can initiate cascades of redox processes that significantly alter the physical and chemical properties of a molten salt, posing challenges for molten salt reactor (MSR) performance and longevity. Consequently, mastering the fundamental behavior of the eS– could enable the design of specific molten salt mixtures with tailored Lewis acidities for controlling the speciation and chemical reactivity of the eS–, thereby mitigating its overall impact on MSR technologies. Here, we use time-resolved electron pulse radiolysis techniques for determining chemical kinetics and transient absorption spectra, combined with ab initio molecular dynamics simulations to explore the influence of multivalent metal cations on the fundamental speciation and distribution of eS– coordination environments.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
NASA NTRS1992

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.

Rock, Stacey G.↗
DOE OSTI2021

An out-of-distribution-aware autoencoder model for reduced chemical kinetics

While detailed chemical kinetic models have been successful in representing rates of chemical reactions in continuum scale computational fluid dynamics (CFD) simulations, applying the models in simulations for engineering device conditions is computationally prohibitive. To reduce the cost, data-driven methods, e.g., autoencoders, have been used to construct reduced chemical kinetic models for CFD simulations. Despite their success, data-driven methods rely heavily on training data sets and can be unreliable when used in out-of-distribution (OOD) regions (i.e., when extrapolating outside of the training set). In this paper, we present an enhanced autoencoder model for combustion chemical kinetics with uncertainty quantification to enable the detection of model usage in OOD regions, and thereby creating an OOD-aware autoencoder model that contributes to more robust CFD simulations of reacting flows.Here, we first demonstrate the effectiveness of the method in OOD detection in two well-known datasets, MNIST and Fashion-MNIST, in comparison with the deep ensemble method, and then present the OOD-aware autoencoder for reduced chemistry model in syngas combustion.

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

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2024

Coherence mapping to identify the intermediates of multi-channel dissociative ionization

Identifying the short-lived intermediates and reaction mechanisms of multi-channel radical cation fragmentation processes remains a current and important challenge to understanding and predicting mass spectra. We find that coherent oscillations in the femtosecond time-dependent yields of several product ions following ultrafast strong-field ionization represent spectroscopic signatures that elucidate their mechanism of formation and identify the intermediate(s) they originate from. Experiments on endo-dicyclopentadiene show that vibrational frequencies from various intermediates are mapped onto their resulting products. Aided by ab initio methods, we identify the vibrational modes of both the cleaved and intact molecular ion intermediates. These results confirm stepwise and concerted fragmentation pathways of the dicyclopentadiene ion. This study highlights the power of tracking the femtosecond dynamics of all product ions simultaneously and sheds further light onto one of the fundamental reaction mechanisms in mass spectrometry, the retro-Diels Alder reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Radiative properties of molecular nitrogen ions produced by helium Penning ionization and argon effects

The development of hypersonic aerospace vehicles requires a better understanding on the thermal and chemical nonequilibrium kinetics of participating species in shock layers. The computational fluid dynamic (CFD) codes developed for such flowfields overestimate the radiation in the spectral region of 300 - 600 nm. A speculation for this overestimation is that inclusion of Ar, CO2, and H2O at the upper atmosphere flight region makes a significant impact on radiative kinetics of molecular nitrogen ions. To define the effects of minority species on the radiative kinetics of N2(+), an experimental setup was made by using the helium Penning ionization. The vibrational and rotational temperature were measured by mapping the vibrational and rotational distributions of N2(+) emission with high spectroscopic resolution and absolute intensity measurements. Measured vibrational temperatures were in the range from 18,000 to 36,000 K, and rotational temperatures were in the range from 300 to 370 K. The irradiance of 391.44 nm line and rotational and vibrational temperatures were analyzed to define argon and CO2 effects on the N2(+) emission. When Ar or CO2 is injected with N2, the rotational temperature did not change. The irradiances were reduced by 34 percent and 78 percent for the 50 percent of mixture of Ar and CO2, respectively. The vibrational temperatures were increased by 24.1 percent and 82.9 percent for the 50 percent of mixture of Ar and CO2, respectively. It appears that there are no significant effects from small concentrations of Ar and CO2 at the upper atmosphere flight region.

Miller, George, III↗
DOE OSTI2023

Time-Gated Raman Spectroscopy (Intern Final Report)

Time-gated Raman spectroscopy is a powerful tool for rapid non-destructive analysis of chemical compounds, particularly when combined with other investigative methods. Time-gating the signal obtained from a sample during Raman spectroscopy allows for the removal of background emission and black body radiation, improving clarity and accuracy of scans. The main objectives of this project were to prove the lab’s time-gated Raman spectroscopy system works as proof of concept and to make select improvements to the in-house LabVIEW software used to run the time-gated system. This time-gated system will be added to the CTK’s TAP reactors in the future, but this is out of the scope of this internship project. This research is important, as understanding the mechanisms behind catalysts leads to designing better systems/materials and the perfection of catalysts, and energy-efficient chemical manufacturing is key when designing a greener future. The research for this project was done with samples at both ambient conditions in sample holders and at operando conditions (i.e., high temperature, unsteady-state) inside a small-volume chemical reactor. Numerous spectro-kinetic scans of the catalysts were obtained, but the advanced dynamics of these systems are still being interpreted. There is data that demonstrates that the time-gated system is correctly rejecting emission from samples, and one can observe real-time changes in the signal from the sample, due to coking or changes induced by redox chemistry. Chemical manufacturers and the planet are the main beneficiaries of catalysis research, as better catalysts will reduce carbon emissions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Stratosphere chemistry in a 2-D model with residual circulation

The objective of this research was to examine the effects of chemical perturbations on the stratosphere using models which can incorporate fully interactive radiative, chemical, and dynamical responses, in the context of a zonally averaged model. Model runs for the unperturbed, chlorine-perturbed and simultaneously chlorine-and CO2-perturbed cases were completed using the JPL-87 chemical kinetics data. The base case was analyzed and submitted for publication. The perturbed cases show substantial sensitivity of the predicted column ozone depletion to the perturbations affecting lower stratosphere temperature, but less to far dynamical perturbations. The column ozone distribution changed substantially when the kinetics data was changed. This implies a greater-than-expected uncertainty in predicted latitude distributions of ozone depletion, due to uncertainty about the accuracy and completeness of the chemical kinetics data set.

Guthrie, Paul D.↗
DOE OSTI2026

Universal progression of structure and dynamics in colloidal nanocrystal gels during salt-accelerated aging

Controlling the structure and function of colloidal gels requires a detailed understanding of how the various components govern network formation and aging. In particular, molecular additives like salts are widely used to tune interparticle interactions, yet their influence on gelation pathways in complex systems such as colloidal nanocrystal gels remains inadequately understood. Here, we investigate how noncoordinating salts modulate the evolution of gels formed using chemically linked tin-doped indium oxide nanocrystals. Through combined structural, dynamic, and kinetic analyses, we demonstrate that increasing salt concentration accelerates gelation. When rescaled by salt-dependent characteristic times, the evolution collapses onto universal trajectories, revealing a time-salt superposition principle. The universality extends across length scales, suggesting a consistent salt-dependent mechanism that controls both local structuring and macroscopic network formation. This observed salt modulation of structure and dynamics provides a predictive basis for controlling the kinetics of nonequilibrium nanocrystal gel assembly, enhancing the rational design of functional nanomaterials with tunable properties.

36 MATERIALS SCIENCE↗