Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential kinetics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Parameter Identification for Battery Abuse Reaction Kinetics Models Using DSC Analysis [SWR-23-87]

This code estimates the parameters as inputs for NREL’s abuse reaction kinetics model using Differential scanning calorimetry (DSC) experimental data. The parameters for calculating the rates of abuse reactions among battery cell components at elevated temperatures include frequency factor, activation energy, and reaction orders. The least-squared method is employed to fit the reaction model to DSC data. The code is implemented in MATLAB.

Yang, Chuanbo↗
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

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

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

Surface Equilibration Mechanism Controls the Stability of a Model Codeposited Glass Mixture of Organic Semiconductors

While previous work has identified the conditions for preparing ultrastable single-component organic glasses by physical vapor deposition (PVD), little is known about the stability of codeposited mixtures. Here, we prepared binary PVD glasses of organic semiconductors, TPD (N,N'-Bis(3-methylphenyl)-N,N'-diphenylbenzidine) and m-MTDATA (4,4',4"-Tris[phenyl(m-tolyl)amino]triphenylamine), with a 50:50 mass concentration over a wide range of substrate temperatures (T sub ). The enthalpy and kinetic stability are evaluated with differential scanning calorimetry and spectroscopic ellipsometry. Binary organic semiconductor glasses with exceptional thermodynamic and kinetic stability comparable to the most stable single-component organic glasses are obtained when deposited at T sub = 0.78–0.90T g (where T g is the conventional glass transition temperature). When deposited at 0.94T g , the enthalpy of the m-MTDATA/TPD glass equals that expected for the equilibrium liquid at that temperature. Thus, the surface equilibration mechanism previously advanced for single-component PVD glasses is also applicable for these codeposited glasses. Furthermore, these results provide an avenue for designing high-performance organic electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2021

CSPlib - A Software Toolkit for the Analysis of Dynamical Systems and Chemical Kinetic Models

CSPlib is an open source software library for analyzing general ordinary differential equation (ODE) systems and detailed chemical kinetic ODE systems. It relies on the computational singular perturbation (CSP) method for the analysis of these systems. The software provides support for: General ODE models (gODE model class) for computing source terms and Jacobians for a generic ODE system; TChem model (ChemElemODETChem model class) for computing source term, Jacobian, other necessary chemical reaction data, as well as the rates of progress for a homogenous batch reactor using an elementary step detailed chemical kinetic reaction mechanism. This class relies on the TChem [2] library; A set of functions to compute essential elements of CSP analysis (Kernel class). This includes computations of the eigensolution of the Jacobian matrix, CSP basis vectors and co-vectors, time scales (reciprocals of the magnitudes of the Jacobian eigenvalues), mode amplitudes, CSP pointers, and the number of exhausted modes. This class relies on the Tines library; A set of functions to compute the eigensolution of the Jacobian matrix using Tines library GPU eigensolver; A set of functions to compute CSP indices (Index Class). This includes participation indices and both slow and fast importance indices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ReactionMechanismSimulator.jl: A modern approach to chemical kinetic mechanism simulation and analysis

Abstract We present ReactionMechanismSimulator.jl (RMS), a modern differentiable software for the simulation and analysis of chemical kinetic mechanisms, including multiphase systems. RMS has already been applied to problems in combustion, pyrolysis, polymers, pharmaceuticals, catalysis, and electrocatalysis. RMS is written in Julia, making it easy to develop and allowing it to take advantage of Julia's extensive numerical computing ecosystem. In addition to its extensive library of optimized analytic Jacobians, RMS can generate and use Jacobians computed using automatic differentiation and symbolically generated analytic Jacobians. RMS is demonstrated to be faster than Cantera and Chemkin in several benchmarks. RMS also implements an extensive set of features for analyzing chemical mechanisms, including a library of easy‐to‐call plotting functions, molecular structure resolved flux diagram generation, crash analysis, traditional sensitivity analysis, transitory sensitivity analysis, and an automatic mechanism analysis toolkit. RMS implements efficient adjoint and parallel forward sensitivity analyses. We also demonstrate the ease of adding new features to RMS.

Johnson, Matthew S.↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Curl up with a good B : detecting ultralight dark matter with differential magnetometry

Ultralight dark matter (such as kinetically mixed dark-photon dark matter or axionlike dark matter) can source an oscillating magnetic-field signal at the Earth’s surface, which can be measured by a synchronized array of ground-based magnetometers. The global signal of ultralight dark matter can be robustly predicted for low masses, when the wavelength of the dark matter is larger than the radius of the Earth, $λ_{DM}$ $\gg$ R. However, at higher masses, environmental effects, such as the Schumann resonances, can become relevant, making the global magnetic-field signal B difficult to reliably model. In this work, we show that ∇ × B is robust to global environmental details, and instead only depends on the local dark matter amplitude. We therefore propose to measure the local curl of the magnetic field at the Earth’s surface, as a means for detecting ultralight dark matter with $λ_{DM}$ ≲ R. As this measurement requires vertical gradients, it can be done near a hill/mountain. Our measurement scheme not only allows for a robust prediction, but also acts as a background rejection scheme for external noise sources. We show that our technique can be the most sensitive terrestrial probe of dark-photon dark matter for frequencies 10 Hz ≤ $f_{A'}$ ≤ 1 kHz (corresponding to masses 4 × 10 –14 eV ≤ $m_{A'}$ ≤ 4 × 10 –12 eV). It can also achieve sensitivities to axionlike dark matter comparable to the CAST helioscope, in the same frequency range.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Pulse Radiolysis Investigation of Radicals Derived from Water-Soluble Cyanine Dyes: Implications for Super-resolution Microscopy

Light-induced blinking, an inherent feature of many forms of super-resolution microscopy, has been linked to transient reduction of the fluorescent cyanine dye used as an imaging agent. There is, however, only scant literature information related to one-electron reduced cyanine dyes, especially in an aqueous environment. Here, we examine a small series of cyanine dyes, possessing disparate π-conjugation lengths, under selective reducing or oxidizing conditions. The experiment allows recording of both differential absorption spectra and decay kinetics of the resultant one-electron reduced or oxidized transient species in water. Relative to the ground state, absorption transitions for the various radicals are weak and somewhat broadened but do allow correlation with the π-conjugation length. In all cases, absorption maxima lie to the blue of the main ground-state transition. Under anaerobic conditions, the transient species decay on the microsecond to millisecond time scale, with the mean lifetime depending on molecular structure, radiation dose, and dye concentration. The experimental absorption spectra recorded for the one-electron reduced radicals and the presumed dimer cation radical compare well to spectra obtained from time-dependent density functional theory calculations. Furthermore, the results allow conclusions to be drawn regarding the plausibility of the reduced species being responsible for light-induced blinking in direct stochastic optical reconstruction microscopy.

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

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

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

Compositional effects on cure kinetics, mechanical properties and printability of dual-cure epoxy/acrylate resins for DIW additive manufacturing

Interest in 3D printing of thermoset resins has increased significantly in recent years. One approach to additive manufacturing of thermoset resins is printing dual-cure resins with direct ink write (DIW). Dual-cure resins are multi-component resins which employ an in situ curable constituent to enable net-shape fabrication while a second constituent and cure mechanism contribute to the final mechanical properties of the printed materials. In this work, the cure kinetics, green strength, printability, and print fidelity of dual-cure epoxy/acrylate thermoset resins are investigated. Resin properties are evaluated as a function of acrylate concentration and in situ UV exposure conditions. The acrylate cure kinetics are probed using photo-differential scanning calorimetry and the impacts of resin composition and UV cure profile on the acrylate extent of conversion are presented. Continuous and pulsed UV cure profiles are shown to affect total conversion due to variances in radical efficiency at different UV intensities and acrylate concentrations. The effects of acrylate concentration on the kinetics of the epoxy thermal cure and the final mechanical properties are also investigated using dynamic mechanical analysis and three-point bend measurements. The glass transition temperature is dependent on formulation, with increasing acrylate content decreasing the T g . However, the room temperature shear moduli, flexural moduli, strength, strain-to-failure, and toughness values are relatively independent of resin composition. The similarity of the final properties allows for greater flexibility in resin formulation and in situ cure parameters, which can enable the printing of complex parts that require high green strength. We found that the in situ UV print intensities and exposure profiles that are necessary to achieve the best print quality are not, in most cases, the conditions that maximize conversion of the acrylate network. Lastly, this highlights the importance of developing optimized resin compositions which enable complete cure of the acrylate network by promoting acrylate dark cure or thermal cure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An SECM-Based Spot Analysis for Redoxmer-Electrode Kinetics: Identifying Redox Asymmetries on Model Graphitic Carbon Interfaces

The fundamental process in non-aqueous redox flow battery (NRFB) operation revolves around electron transfer (ET) between a current collector electrode and redox-active organic molecules (redoxmers) in solution. Here, we present an approach utilizing scanning electrochemical microscopy (SECM) to evaluate interfacial ET kinetics between redoxmers and various electrode materials of interest at desired locations. This spot-analysis method relies on the measurement of heterogeneous electron transfer rate constants (k f or k b ) as a function of applied potential (E-E 0 '). As demonstrated by COMSOL simulations, this method enables the quantification of Butler-Volmer kinetic parameters, the standard heterogeneous rate constant, k 0 , and the transfer coefficient, α. Our method enabled the identification of inherent asymmetries in the ET kinetics arising during the reduction of ferrocene-based redoxmers, compared to their oxidation which displayed faster rate constants. Similar behavior was observed on a wide variety of carbon electrodes such as multi-layer graphene, highly ordered pyrolytic graphite, glassy carbon, and chemical vapor deposition-grown graphite films. However, aqueous systems and Pt do not exhibit such kinetic effects. Our analysis suggests that differential adsorption of the redoxmers is insufficient to account for our observations. Displaying a greater versatility than conventional electroanalytical methods, we demonstrate the operation of our spot analysis at concentrations up to 100 mM of redoxmer over graphite films. Looking forward, our method can be used to assess non-idealities in a variety of redoxmer/electrode/solvent systems with quantitative evaluation of kinetics for applications in redox-flow battery research.

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

Stochastic hydrodynamics and hydro-kinetics: Similarities and differences

The hydro-kinetic formalism has been used as a complementary approach to solving the Stochastic Differential Equations (SDE) corresponding to noisy hydrodynamics. The hydro-kinetic formalism consists of a deterministic set of relaxation type equations that tracks the evolution of 2-point correlation functions of stochastic hydrodynamic quantities. Hence they are comparatively easier to solve than the SDEs, which are computationally intensive and need to deal with arbitrarily large gradients. This work compares the two approaches for the propagation and diffusion of conserved charge fluctuations in the Bjorken hydrodynamic model. For white noise, the two approaches agree. For colored Catteneo noise, which is causal, the two approaches diverge. This is because white noise only induces two-point correlations, while Catteneo noise also induces higher-order correlations. Furthermore, this difference is quantified from the effects of causal evolution and influence from higher-order correlations induced by the Catteneo noise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗