Origin of anomalous high-rate Na-ion electrochemistry in layered bismuth telluride anodes
Not Available
SEARCH · Search NASA
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.
Not Available
Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.
The front propagation formulation (FPF) for large eddy simulation (LES) of turbulent premixed flames is presented and discussed. The FPF method adopts a regularized Dirac delta function to minimize the spurious propagation of under-resolved filtered flame fronts in LES, while reproducing the reaction characteristics of premixed flame fronts. Here, the extension of the FPF method to LES and the relation to existing sub-filter combustion models are discussed. An algebraic model for the regularized Dirac delta function is presented and validated using direct numerical simulation (DNS) data for statistically planar premixed flames propagating in homogeneous isotropic turbulence. A posteriori validation for a filtered 1-D planar laminar flame is performed to demonstrate the characteristics of FPF in reproducing the propagation speed and the thickness of the under-resolved flame front. LES of a turbulent Bunsen flame is also performed to demonstrate the performance of the proposed method.
Reactive burn models for propagating detonation waves in an explosive are based on the ZND (Zel'dovich-Von Neumann-Doring) theory. A key property of the ZND theory, which leads to self-sustaining detonation waves, is a sonic point relative to the front in the reaction zone. Moreover, burn models assume the explosive is a homogeneous material. This assumption is used for plastic-bonded explosives, even though they are heterogeneous materials. The heterogeneities are accounted for only with an empirical burn rate. Furthermore, the reaction-zone width can be less than the length scale of heterogeneities. This raises the question of the effect of heterogeneities on the reaction zone. In particular, whether the reaction zone can be steady and have a well defined sonic point. To examine this question we focus on PBX 9501 since there is data on the grain-binder heterogeneities and on the reaction-zone velocity time history of a steady planar detonation wave. Using pseudo-mesoscale 1-D simulations, we examine the effect of a detonation wave transversing explosive and binder segments of random length. The simulations show that the reaction zone is quasi-steady with only short wavelength small amplitude variations in the detonation wave pressure and speed. The main effect of the heterogeneities is to introduce noise in the pressure and velocity fields that propagates into the Taylor wave behind the detonation. This is consistent with the scatter in the reaction-zone measurements for PBX 9501.
Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.
Cast Stone is a cementitious waste form being considered for the solidification of low activity waste at the DOE Hanford Site. Under near-surface disposal conditions in an arid environment, Cast Stone is subject to drying and carbonation which may impact retention of waste constituents. Here this study investigates the effects of environmental CO 2 concentration and relative humidity (RH) on the carbonation of Cast Stone. The rate of carbonation front ingress and the extent of carbonation reaction were characterized for samples aged up to 48 weeks at three RH levels and two CO 2 concentrations. While the 68% RH environment allowed the greatest reaction extent, the 15% RH environment yielded the deepest carbonation front ingress. At 68% RH, there was a linear relationship between the ingress rate and CO 2 concentration. Carbonation reactions increased the drying rate of Cast Stone. Redistribution of sodium toward the wetting front was observed under drying and carbonation conditions.
The evolution of single-molecule and -particle fluorescence microscopy imaging techniques for the investigation of chemical reactions yielded achievements at all fronts, from synthetic organic chemistry to materials science. The removal of ensemble averaging effects uncovers unique, otherwise hidden reaction intermediates, reactivity, and kinetics. Thus, new insights can be obtained even for well-explored chemical reactions and processes. The aim of this article is to serve as a guide and a source of inspiration for scientists who wish to develop single-molecule and -particle fluorescence microscopy to answer questions in chemistry. Here, we describe the current status of technical accomplishments in the field and point out technical and conceptual obstacles for overcoming in the near future, with focus on improvement of bench-scale synthetic reactions.
Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.
Achieving the desired combination of superior detonation performance and insensitivity to shock initiation has been a long-standing goal in high explosive (HE) science and technology. Having previously established the shock insensitivity of 1,1-diamino-2,2-dinitroethene (also known as DADNE or FOX-7) single crystals to 20 GPa (extended to 25 GPa in this work), the FOX-7 detonation response was determined through wave profile measurements in ~250 μm thick single crystals shock compressed to 64 GPa. Quite unexpectedly, FOX-7 demonstrated the classic Chapman–Jouguet (C–J) detonation response—reaction completion in the detonation front (<0.7 ns) at pressures of 44 GPa and higher—not observed in other insensitive high explosives. Here, the experimentally determined C–J pressure (35 GPa), detonation wave velocities and the detonation products equation of state—together with shock insensitivity to 25 GPa—demonstrate that FOX-7 single crystals display a near-optimal combination of high performance and shock insensitivity, not observed in another HE crystal.
The international thermonuclear experimental reactor (ITER) is a worldwide project currently being built in France for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. ITER is of the tokamak based design using strong magnetic fields to confine the very hot plasma needed to induce the fusion reaction. Tokamak devices are currently the front leading designs. Building a successful magnetic fusion device for energy production is of great challenge. A key obstacle to such design is the performance during abnormal events including plasma disruptions and so-called edge-localized modes (ELMs). In these events, a massive and sudden release of energy occurs quickly, due to loss of full or partial plasma confinement, leading to very high transient power loads on the reactor surface boundaries. A successful reactor design should tolerate several of these transient events without serious damages such as melting and vaporization of the structure. This paper highlights, through comprehensive state-of-the-art computer simulation of the entire ITER interior design during such transient events, e.g., ELMs occurring at normal operation and disruptions during abnormal operation, potential serious problems with current plasma facing components (PFCs) design. The HEIGHTS computer package is used in these simulations. The ITER reactor design was simulated in full and exact 3D geometry including all known relevant physical processes involved during these transient events. The current ITER divertor design may not work properly and may requires significant modifications or new innovative design to prevent serious damage and to ensure successful operation.
This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.
Here, we demonstrate that a paradigm shift from considering the deuteron as a system of a bound proton and neutron to that of a pseudovector system in which we observe a proton and neutron results in the possibility of probing a new “incomplete” P-statelike structure on the light front (LF). This occurs at extremely large internal momenta, which can be achieved in a high energy transfer electrodisintegration of the deuteron. Investigating the deuteron on the light front, where the vacuum fluctuations are suppressed, we found that this new structure, together with the conventional S and D states, is leading order in transferred energy of the reaction and thus not suppressed on the light front. The incompleteness of the observed P state results in a violation of the angular condition that can happen only if the deuteron contains non-nucleonic structures, such as Δ Δ, N*N or hidden color components. We demonstrate that experimentally verifiable signatures of incomplete P states are angular anisotropy of the light front momentum distribution of the nucleon in the deuteron, as well as an enhancement of the tensor polarization strength beyond the S- and D-wave predictions at large internal momenta in the deuteron.
Heterolayered structures consist of two or more different types of layer and can exhibit exceptional physical properties. Rational routes to synthesize new members of such compounds are required because most of these compounds have been discovered unintentionally. So far there is no generic method to vertically stack chemically different layers to form two-dimensional compounds owing to a lack of understanding of the synthesis of these materials. Here we report the use of molten hydroxides as unconventional solutions for the rapid stacking of oxide and chalcogenide layers with precise composition control. In addition, the crystal growth of heterolayered phases can be achieved by the reaction of different components at their diffusion front in molten hydroxides. This approach creates conditions in which the building blocks for each heterolayer can coexist, enabling heterolayered structures and bypassing the challenges of traditional solid-state chemistry methods where short reactant diffusion lengths predominate. This crystal growth methodology for heterolayers is also applicable to systems that do not form congruent melts at high temperatures.
Research into the early-stage uranium-hydrogen degradation mechanism lacks unified interpretation and understanding due to inadequate temporal and spatial characterization resolution required to unify reaction observations. In this study, white-light interferometry was adapted to characterize the uranium-hydrogen reaction at a temperature of 50 °C and hydrogen pressure of 13.8 kPa. The uranium surface remained unchanged for 62 min during the hydride induction period before the first hydride blister appeared. This blister was monitored from induction through critical spallation and UH₃ powder release during rapid surface hydriding. The critical spallation dimensions measured for height, diameter, area, and perimeter of the first blister were 1.79 µm, 17.58 µm, 187.7 µm 2 and 59.79 µm. Hydriding growth kinetics were assessed by tracking the first hydride site over the reaction time. Analysis shows that the spall-front velocity of the hydride was changing at a rate of 0.91 µm/min which compares favorably with legacy Sievert’s experiments in literature. Total percentage of the surface area hydrided after 242 min was 42.7%, and post-characterization of the UH 3 powders shows both α-UH 3 and β-UH 3 . This work highlights the strength of characterizing the early-stage uranium-hydrogen reaction using white-light interferometry, but more importantly, unifying the understanding and mechanism of the uranium-hydrogen reaction kinetics.
Due to downsizing trends, flame-wall interaction (FWI) is increasingly prominent in gas turbines (GTs). FWI has direct consequences on flame stabilization and pollutant emissions, but it is not well understood in turbulent flows representative of GTs. We present results from a direct numerical simulation (DNS) of a turbulent CH4/H2 model GT low-swirl laboratory-scale flame interacting with an inclined wall. The results from the laboratory flame include simultaneous measurements of velocity using stereo particle imaging velocimetry and OHxCH2O planar laser induced fluorescence. The adaptive-mesh refinement solver PeleLMeX is used, with 24-species, 105-reaction reduced Aramco chemical kinetics mechanism. The premixed fuel-air mixture consists of hydrogen-enriched methane with 70% hydrogen volume fraction and 0.4 equivalence ratio. The inflow is prescribed to match experimental measurements at the burner exit. Karlovitz and turbulent Reynolds numbers are 300 and 400, respectively. The simulation and experimental results show excellent agreement. The flame features a bowl-shape stabilization, with a corrugated, continuous flame front at the leading edge, followed by fragmented reaction zones downstream. A large diffuse cloud of CH2O is formed downstream of the quenching point. The simulation results indicate that the cloud of CH2O is the result of incomplete methane combustion, with CH2O "leaking" from the locally quenched reaction zones.The DNS provides fine-grain resolution of turbulence-flame-wall interaction that cannot be captured with experimental measurements. With access to the entire solution vector at each cell of the computational domain, the local quenching.
Inspired by nature's morphogenesis, a new 3D printing process –growth printing (GP)– takes advantage of a self‐propagating curing front to produce 3D polymeric parts following a growth‐like development plan. The propagation of the curing front is driven by the exothermic polymerization of dicyclopentadiene (DCPD), which transforms the liquid resin into a stiff polymer as it propagates at 1 mm s −1 . GP is triggered when a heated initiator contacts the uncured liquid resin in an open container. The initiator nucleates the frontal polymerization reaction and the isotropic radial propagation of the growth front. Simultaneously, the initiator is moved up across the free surface of the resin, pulling the cured object out of the uncured resin. The motion trajectory of the initiator with respect to the free resin surface controls the growth morphology of the 3D part. An inverse design algorithm is developed to produce 3D parts by modeling the reaction‐diffusion‐driven solidification process. This process has substantial energy savings and high printing speeds.
Reactive burn models for a high explosive (HE) have three components: an equation of state (EOS) for the products, an EOS for the reactants and a reaction rate. For the products, we discuss the uncertainty in data for the planar detonation state and its release isentrope. For the reactants, the uncertainty in the principal shock locus data is discussed. For the reaction rate, we discuss uncertainties in shock initiation data (distance-of-run to detonation as a function of shock pressure), and curvature effect data (detonation speed as a function of front curvature). An important source of uncertainty for calibrating the reaction rate for a plastic-bonded explosive (PBX) results from variations in heterogeneities; that is the mesoscale structure of HE grains, binder and porosity. This shows up in the scatter of both shock initiation data and curvature effect data with small changes in initial density and initial temperature, and systematic changes with PBX lot. As illustrative examples, we use data for the conventional HMX based PBX 9501 and the insensitive TATB based PBX 9502 since these explosives have been extensively studied and high quality data are readily available in the literature.
In combustion simulations, it is usually assumed that highly-energized collision complexes (either radical or stable species) formed in exothermic reactions are always under thermal equilibrium as they evolve through competing networks of reactions. However, in practical flames, non-trivial amounts of reactive radicals such as H, O and OH are present apart from O2. As a result, collisions of these reactive species with the energized collision complexes have been shown in recent studies to induce non-thermal reactivity. These studies have also demonstrated that such non-thermal reactions can be suitably represented in macroscopic kinetics models as chemically termolecular reactions. The present work was focused on identifying and quantifying the effects of including such reactions on the evolution of an initial deflagration front to a developing detonation in H2/CH4-air mixtures under boosted internal combustion (IC) engine conditions. Specifically, fully resolved simulations, with and without non-thermal reactivity, were performed for a constant volume reactor containing stoichiometric H2/CH4-air mixture. It was found that inclusion of termolecular chemistry resulted in a delayed onset of end-gas auto-ignition. Concurrently, the developing detonation intensity was observed to be significantly higher. Chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions responsible for the observed variation in the results.