Results for “differential kinetics” — Page 4

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↗

NASA NTRS1999

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods. Dilatometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions.

Bayuzick, R. J.↗

NASA NTRS1998

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods. Dilitometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post-processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions. Levitation processing of refractory materials in ultra high vacuum provides an avenue to conduct these kinetic studies on single samples. Two experimental methods have been identified where ultra high vacuum experiments are possible; electrostatic levitation in ground-based experiments and electromagnetic processing in low earth orbit on TEMPUS. Such experiments, reported here, were conducted on zirconium. Liquid zirconium is an excellent solvent and has a high solubility for contaminants contained in the bulk material as well as those contaminants found in the vacuum environment. Oxides, nitrides, and carbides do not exist in the melt, and do not form on the surface of molten zirconium, for the materials and vacuum levels used in this study. Ground-based experiments with electrostatic levitation have shown that the statistical nucleation kinetic experiments are viable and yield results which are consistent with classical nucleation theory. The advantage of low earth orbit experiments is the ability to vary the flow conditions in the liquid prior to nucleation. The put-pose of nucleation experiments in TEMPUS was to examine.

Bayuzick, R. J.↗

NASA NTRS1977

Kinetics of imidization and crosslinking in PMR-polyimide resin

Infrared spectroscopy and differential scanning calorimetry were employed to study the imidization and crosslinking kinetics of norbornenyl-capped, addition-type polyimide resins (designated PMR for polymerization of monomer reactants). The spectral and thermal analyses were performed on resin specimens which had been isothermally aged at temperatures appropriate for imidization (120 to 204 C) and crosslinking (275 to 325 C). Imidization occurs rapidly (approximately 0.01/min) at short times, while at times longer than approximately 0.5 hour, the rate decreases significantly (approximately 0.0001/min). The crosslinking reaction exhibits first order kinetics during the initial portion of the reaction and its rate appears to be limited by the reversion of the norbornenyl Diels-Alder adduct. The total heat evolved per mole of endcap during crosslinking shows an inverse dependence on the molecular weight of the imide prepolymers. This reflects the effect of endcap dilution and decreased mobility of the larger oligomers.

Lauver, R. W.↗

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↗

NASA NTRS1967

On the accuracy of determining rate constants in enzymatic reactions.

Rate constant estimates accuracy assessed for differential equations in connection with kinetics of enzymatic reaction

Aris, R.↗

NASA NTRS1977

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗

NASA NTRS1979

Crystallization kinetics in a multiply saturated basalt magma - An experimental study of Luna 24 ferrobasalt

The paper deals with the effects of kinetics on mineral/melt partitioning and on fractional crystallization for a Luna-24 ferrobasalt. The composition is nearly multiply saturated under lunar surface conditions, making it possible to study the response of several mineral phases to kinetic factors during cooling. The differential suppression of the temperature of appearance of olivine, clinopyroxene, and plagioclase causes changes in the liquid line of descent. The course of liquid line of descent is mapped as a function of the cooling rate, and the partitioning of elements between pyroxene/liquid and olivine/liquid is examined.

Grove, T. L.↗

NASA NTRS2011

Stellar Mixing: I. Formalism

In this paper we use the Reynolds stress models (RSM) to derive algebraic expressions for the following variables: a) heat fluxes; b) J.l fluxes; and c) momentum fluxes. These relations, which are fully 3D, include: 1) stable and unstable stratification, represented by the Brunt-Vaislila frequency, N(exp 2) =-g/H(sub p_(del - del(sub ad))(1 - RI(sub mu)); 2) double diffusion, salt-fingers, and semi-convection, represented by the density ratio R(sub mu) = del(sub mu)/(del - del(sub ad)); 3) shear (differential rotation), represented by the mean squared shear Sigma(exp 2) or by the Richardson number, Ri =N(exp 2)Sigma(exp -2); 4) radiative losses represented by a Peclet number, Pe; 5) a complete analytical solution of the ID version of the model. In general, the model requires the solution of two differential equations for the eddy kinetic energy K and its rate of dissipation, epsilon. In the local and stationary cases, when production equals dissipation, the model equations are all algebraic.

Canuto, V .M.↗

NASA NTRS1999

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods [1-6]. Dilitometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics [7]. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov [8]. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions.

Bayuzick, Robert J.↗

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↗

NASA NTRS1984

Kinetics of crystal nucleation and growth in Pd(40)Ni(40)P(20) glass

Samples of Pd(40)Ni(40)P(20) glass, produced by cooling the melt at 1 or 800 K/s, are heated in a differential scanning calorimeter to determine the crystallization kinetics. Optical microscopy shows that eutectic crystallization proceeds both by growth from the surface of the samples and by the growth of spherical regions around preexisting nuclei in the interior. A modified Kissinger (1957) analysis is used to obtain the activation energy for crystal growth (3.49 eV). The steady state homogeneous nucleation frequency at 590 K is about 10 million/cu m per sec. This is estimated to be the maximum nucleation frequency: it is too low to account for the observed population of quenched-in nuclei, which are therefore presumed to be heterogeneous. The major practical obstacle to glass formation in this system is heterogeneous nucleation.

Drehman, A. J.↗