SEARCH · Search NASA

Results for “kinetics and energetics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Kinetic versus Thermodynamic Stability of LLZO in Contact with Lithium Metal

Li 7 La 3 Zr 2 O 12 (LLZO) garnet-based oxides are a promising class of solid electrolytes used as the separator in all-solid-state batteries (ASSBs). While LLZO is considered to have a wide electrochemical stability window, its intrinsic stability in contact with lithium metal is not sufficiently well understood, and there is still a debate on the key question of whether LLZO does or does not form passivation layers before and during cycling. Utilizing both in situ and operando X-ray photoelectron spectroscopy techniques, in this work we reveal the presence of a kinetic barrier to the reduction of LLZO by Li metal, with the extent of oxygen-deficient interphase (ODI) formation depending sensitively on the energetics of Li metal arriving at the Li vertical bar LLZO interface. Despite the clear presence of a kinetic barrier to reduction, the electrochemical response of the Li vertical bar LLZO interface is unchanged by the presence of the ODI, indicating that ODI formation during electrochemical cycling does not hinder charge transfer across the Li vertical bar LLZO interface. Overall, these results reveal that the reactivity of LLZO with Li metal depends not only on the material properties of the adjoining phases (i.e., surface purity and active contact) and their resulting thermodynamic stability but also on the energy input at the interface and the resulting reaction kinetics. Furthermore, the presence of a kinetic barrier to reduction highlights the additional complexities governing the reactivity of solid-state interfaces in ASSBs and underscores the importance of operando characterization of interfacial stability to design more robust, high-performance protection strategies for solid electrolytes in contact with reactive electrodes.

25 ENERGY STORAGE↗

DOE OSTI2023

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗

DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Enhancing the connection between computation and experiments in electrocatalysis

Combining computational and experimental methods is a powerful approach to understand the variables that govern catalyst performance and ultimately design improved materials. However, the effectiveness of this approach rests on the strength of the relationships between calculated parameters and experimental measurements. These relationships are complicated by the intricacy and dynamic behaviour of catalytic active sites, and by the non-trivial relationship between calculated reaction energetics and observed rates. In this Perspective, we highlight opportunities to enhance the connection between computation and experiment in electrocatalysis. These include measuring the intrinsic kinetic behaviour of catalysts, creating precise models for the active site and its environment, and forming clear relationships between calculated reaction energetics and observed rates. As experimental and computational methods continue to become more powerful, clear connections between the two will maximize their utility to guide the design of efficient and selective electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

FLAMINGO: combining kinetic SZ effect and galaxy–galaxy lensing measurements to gauge the impact of feedback on large-scale structure

Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on $\sim$Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here, we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev–Zel’dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy–galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called $S_8$ tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect – weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Topological kinetic crossover in a nanomagnet array

Ergodic kinetics, which are critical to equilibrium thermodynamics, can be constrained by a system’s topology. Here, we studied a model nanomagnetic array in which such constraints visibly affect the behavior of the magnetic moments. In this system, magnetic excitations connect into thermally active one-dimensional strings whose motion can be imaged in real time. At high temperatures, our data showed the merging, breaking, and reconnecting of strings, resulting in the system transitioning between topologically distinct configurations. Below a crossover temperature, the string motion is dominated by simple changes in length and shape. In this low-temperature regime, the system is energetically stable because of its inability to explore all possible topological configurations. This kinetic crossover suggests a generalizable conception of topologically broken ergodicity and limited equilibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Data for Selective Oxidation of 3-Hydroxypropionic Acid to Malonic Acid over Pd/C: Mechanistic and Kinetic Study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for the sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O2 and H2O2), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R2 > 0.95). Kinetic observations also enabled a direct comparison between O2 and H2O2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

Bioproducts↗

DOE OSTI2026

Data and code from: Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

3-Hydroxypropionic Acid↗

DOE OSTI2023

How does your gyroid grow? A mesoatomic perspective on supramolecular, soft matter network crystals

We propose and describe a framework to understand the structure of supramolecular network crystals formed in soft matter in terms of mesoatomic building blocks, collective groupings of amphiphilic molecules that play a role analogous to atomic or molecular subunits of hard matter crystals. While the concept of mesoatoms is intuitive and widely invoked in crystalline arrangements of spherelike or cylinderlike (micellelike) domains, analogous notions of natural and physically meaningful building blocks of triply periodic network (TPN) crystals, like the double-gyroid or double-diamond structures are obscured by the complex, bicontinuous domain shapes and intercatenated topologies of the double networks. Focusing on the specific example of diblock copolymer melts, we propose generic rules for decomposing TPN crystals into a unique set of mesoatomic building blocks. Based on physically motivated principles, the combination of symmetries and topologies of these structures point to mesoatomic elements associated with the nodal connections, leading to mesoatomic volumes that are nonconvex and bound by smoothly curved faces, unlike the more familiar Voronoi polyhedral shapes associated with spherelike and cylinderlike mesoatoms. We analyze the shapes of these mesoatoms, their internal structure, and importantly, their local packing with neighbor mesoatomic units. Importantly, we hypothesize that mesoatoms are kinetically favored intermediate structures whose local shapes and packing template network crystal assembly on long time scales. We propose and study a minimal energetic model of mesoatom assembly for three different cubic double-network crystals, based on local shape packing, which predicts a detailed picture for kinetics of intercatenation and surface growth. Based on these analyses, we discuss several possible extensions and elaborations of the mesoatomic description of supramolecular soft matter network crystals, most notably the implications of mesoatomic malleability, a feature that distinguishes soft matter from hard matter crystals. We describe experimental observations of malleable mesoatomic units in the precursor sponge phase as well as in ordered cubic networks and suggest possibilities for observing mesoatoms in primordial, precrystalline states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O 2 and H 2 O 2 ), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R 2 > 0.95). Kinetic observations also enabled a direct comparison between O 2 and H 2 O 2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-hydroxypropionic acid↗

DOE OSTI2022

Pathways for the Formation of C 2+ Products under Alkaline Conditions during the Electrochemical Reduction of CO 2

Establishing how Cu facilitates the electrochemical CO 2 reduction reaction (CO2RR) to C 2+ products remains a critical challenge. Under typical reaction conditions, the pH near the electrode is considerably more alkaline than that in the bulk due to mass transport limitations. Challenges with probing alkaline pathways using computational methods have limited understanding of the CO2RR under experimentally relevant conditions. In this work, using the Volmer reaction on Cu (100), we demonstrate that predicted activation barriers can substantially differ between acidic and alkaline pathways. We compute reaction energetics for alkaline *CO protonation and find that, while the formation of *CHO is preferred thermodynamically, the formation of *COH is favored kinetically at high overpotential. However, we find formation of *CHO via reaction of *H and *CO feasible at room temperature. Further, we report potential-dependent energetics for forming the first C-C bond in CO2RR and find that CO dimerization likely dominates. Finally, we investigate how long-range van der Waals interactions impact our results by comparing to the meta-GGA B97M-rV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Studies of High Energy Density Discharge and Laser-Driven Deflagrating Plasma Stagnations

The goal of this proposal is to develop an understanding of the thermalization and mixing of high energy density (HED) interpenetrating and stagnating plasma deflagrations on two different platforms. The first is opposing gas-fed Z-pinch plasma expansions in a facility at Stanford University. The second is the convergence of multiple deflagration fronts generated by inverted corona laser ablations through experiments carried out by the P.I. as part of a larger team at the University of Rochester Laboratory for Laser Energetics. The link between the two platforms is the role played by kinetic effects that impact mixing and the conversion of directed to thermal energy. Such kinetic effects are important because of the relatively high temperatures despite high plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Towards engineering the growth morphology of anisotropic bicrystals

Eutectic alloys, such as Pb-free Bi-Sn-Ag or Sn-Ag-Cu, are essential for their low melting points in solder joints but suffer from brittle intermetallic phases, reducing reliability. Traditional solder design relies on equilibrium phase diagrams, often neglecting the kinetic influence of crystalline anisotropy during solidification. This study hypothesizes that the energetics at the solid-solid-liquid trijunction dictate crystalline facet selection. Using the Al-Al 3 Ni alloy as a model, this project combines atomistic simulations and experimental methods to understand interfacial energies and develop predictive models. This approach aims to mitigate the negative effects of intermetallics, enhancing component reliability in electronic applications.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Kinetics vs Thermodynamics: Engineering Photoredox Reactivity from an Upper Excited State of Fe II

Ultrafast deactivation of metal-to-ligand charge-transfer (MLCT) states into low-lying metal-centered states has long limited the utility of first-row transition-metal complexes in a broad range of applications, including photoredox catalysis. Here we bypass such limitations using a Fe II -pyridinium carbene complex to enable electron transfer reactivity directly from a higher-lying MLCT state. Time-resolved spectroscopic data revealed that the higher-energy MLCT manifold relaxes to lower-lying metal-centered states with a time constant of ca. 3 ps. The nature of the metal-centered (MC) excited state as a 3 MC (i.e., S = 1) was inferred from variable-temperature transient absorption studies. These experiments also revealed that the ground-state recovery process occurs at or very near the “Marcus barrierless regime”, thereby allowing for an estimate of its excited-state redox potential. The detailed picture of the energetics of this compound that emerged enabled fine-tuning of competing thermodynamic and kinetic pathways to effect electron transfer from the higher-energy MLCT excited state prior to relaxation to the ligand-field manifold. We believe these results open an unexplored landscape for the use of earth-abundant first-row transition-metal-based chromophores for applications in light-to-chemical energy conversion.

charge transfer↗

DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗