Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arrhenius”

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

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2023

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

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

Structure–transport correlations in Na 11 Sn 2 SbSe 12 and its sulfide solid solutions

We report a new ion conductor—Na 11 Sn 2 SbSe 12 —as a possible candidate for a solid catholyte in composite cathodes for all-solid state Na-batteries, which exhibits a room temperature ionic conductivity of 0.15 ± 0.03 mS cm −1 and an activation energy of 0.39 ± 0.02 eV. The sulfide solid solutions of Na 11 Sn 2 SbSe 12 , namely, Na 11 Sn 2 SbS 12−x Se x (x = 1 and 6), were also investigated through a combination of Rietveld refinement against powder x-ray diffraction data and electrochemical impedance spectroscopy to reveal the complex structure–property relationships governing ion transport in this class of materials. Meanwhile, broadening of the Na-ion diffusion pathways in Na 11 Sn 2 SbSe 12 is expected to facilitate Na-ion transport compared to the sulfide-rich member of the solid solution, the opposite holds: increasing the Se fraction in Na 11 Sn 2 SbS 12−x Se x leads to a little change in the activation energy but a reduction in the ionic conductivity. We ascribe this to the lowering of the prefactor, σ 0 , in the Arrhenius relationship with increasing lattice “softening” as a function of higher Se content.

Ramos, Erika P.↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

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

Self-diffusion of liquid deuterium hydride and liquid tritium

Here, we present a quasi-elastic neutron scattering study of liquid deuterium hydride carried out using the Disk Chopper Spectrometer at the National Institute of Standards and Technology. Under saturated vapor pressure, the self-diffusion constant of deuterium hydride obeys an Arrhenius law D = D 0 exp(-E A /k B T), where the prefactor D 0 is given by D 0 = 9.5 ± 1.2 Å 2 /ps and the activation energy is given by E A = 58 ± 2 K. We apply the quantum law of corresponding states to the known diffusion constants of the hydrogen isotopologues. From this application, we estimate that D 0 ≈ 9.1 Å 2 /ps and E A ≈ 75 K in liquid tritium. Young’s theory of quantum-mechanical effects in van der Waals fluids is shown to apply to the diffusion constants of the liquid hydrogens. Our results underscore the importance of nuclear quantum effects in shaping the properties and behavior of the hydrogen isotopologues.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Effect of a collapsing gas bubble on the shock-to-detonation transition in liquid nitromethane

We studied the shock-induced collapse of butane gas bubbles in the homogeneous explosive nitromethane (NM) to investigate the effects of hot spot formation on the detonation process. A butane bubble was injected into a sample of NM, and a shock wave from a flat plate impactor compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and bubble collapse processes. A multiband optical fiber pyrometer measured the time-resolved thermal radiance, and we used the results and emissivity values extracted from spectral fits to estimate temperatures. We measured the characteristics of the shock-to-detonation transition in NM with and without a bubble. All experiments were performed at shock pressures near 8 GPa, where neat NM can detonate. A single bubble in this system was shown to sensitize NM, leading to a reduced run-to-detonation time. We used hydrodynamic modeling to predict shock wave propagation, the extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. We used a temperature-dependent Arrhenius burn model for simulations, and it yielded much better results than reactive burn models that depend only on pressure and density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Translational diffusion in supercooled water at and near the glass transition temperature—136 K

The properties of amorphous solid water at and near the calorimetric glass transition temperature, T g , of 136 K have been debated for years. One hypothesis is that water turns into a “true” liquid at T g (i.e., it becomes ergodic) and exhibits all the characteristics of an ergodic liquid, including translational diffusion. A competing hypothesis is that only rotational motion becomes active at T g , while the “real” glass transition in water is at a considerably higher temperature. To address this dispute, we have investigated the diffusive mixing in nanoscale water films, with thicknesses up to ∼100 nm, using infrared (IR) spectroscopy. The experiments used films that were composed of at least 90% H 2 O with D 2 O making up the balance and were conducted under conditions where H/D exchange was essentially eliminated. Because the IR spectra of multilayer D 2 O films (e.g., thicknesses of ∼3–6 nm) embedded within thick H 2 O films are distinct from the spectrum of isolated D 2 O molecules within H 2 O, the diffusive mixing of (initially) isotopically layered water films could be followed as a function of annealing time and temperature. The results show that water films with total thicknesses ranging from ∼20 to 100 nm diffusively mixed prior to crystallization for temperatures between 120 and 144 K. The translational diffusion had an Arrhenius temperature dependence with an activation energy of 40.8 ± 3.5 kJ/mol, which indicates that water at and near T g is a strong liquid. The measured diffusion coefficient at 136 K is 6.25 ± 1.4 × 10 −21 m 2 /s.

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

Synthesis and composition control of epitaxial FeWO4 thin films using pulsed laser deposition

We report the growth of epitaxial FeWO4 thin films on c-plane sapphire via pulsed laser deposition using a Fe 2 O 3 /WO 3 target ablated in an O 2 atmosphere. At a constant O 2 pressure of 1 mTorr, x-ray diffraction (XRD) confirms FeWO 4 as the major crystalline phase for substrate temperatures from 500 to 800 °C. As temperature increases, x-ray fluorescence (XRF) shows the Fe/W ratio remains nearly constant at 0.90 ± 0.02, while x-ray photoelectron spectroscopy (XPS) shows Fe 3+ to Fe 2+ conversion increases from 20% to 35%. Morphological analysis reveals phase separation, likely of amorphous Fe 3+ oxide from crystalline FeWO 4 , increasing with Fe 3+ conversion. This is attributed to an O-rich laser ablated flux, where conversion is driven by the Arrhenius temperature dependence of Fe 3+ → Fe 2+ reduction on the film surface. At a constant substrate temperature of 750 °C, XRD confirms FeWO 4 formation for O 2 pressures from 0.5 to 10 mTorr. As pressure increases, XRF shows the Fe/W ratio decreases from 0.98 to 0.70, while XPS shows Fe 3+ conversion rises from 15% to 70%. Morphology shows phase separation decreasing with increasing Fe 3+ conversion. This is attributed to scattering, where higher O 2 pressure makes the laser ablated flux O-deficient relative to Fe and W, facilitating Fe 2+ formation. Films with Fe 3+ conversion above ∼30% and Fe/W ratios from 0.86 to 0.96 exhibit FeWO 4 optical and electronic properties suitable for photoanode applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

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

Sulfidation kinetics of titanium and Ti-6Al-4V with elemental sulfur

The sulfidation kinetics of titanium and Ti-6Al-4V (Ti64) using elemental sulfur were studied by measuring the amount of metal consumed instead of the change in mass or scale thickness. Here this metal consumption approach combines the parabolic oxidation law, Arrhenius equation, and Pilling-Bedworth ratio to determine the apparent activation energy and the pre-exponential factor of consumption of Ti64 during sulfidation. The estimations of apparent activation energy and pre-exponential factor were 110.4 kJ·mol -1 and 1.3×10 -9 m 2 s -1 for the titanium-sulfur, and 116.2 kJ·mol -1 and 1.2×10 -9 m2s -1 for the Ti64-sulfur system. Finally, a material removal predictor was developed to predict the amount of Ti64 consumed for a given sulfidation temperature and time.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quasi-equilibrium predictions of water desorption kinetics from rapidly-heated metal oxide surfaces

Controlling sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for pulsed power applications. Despite the short time scales involved, quasi-equilibrium ideas based on transition state theory (TST) and Arrhenius temperature dependence have been widely applied to fit desorption activation free energies. Here, we apply molecular dynamics (MD) simulations in conjunction with equilibrium potential-of-mean-force (PMF) techniques to directly compute the activation free energies (ΔG*) associated with desorption of intact water molecules from Fe 2 O 3 and Cr 2 O 3 (0001) surfaces. The desorption free energy profiles are diffuse, without maxima, and have substantial dependences on temperature and surface water coverage. Incorporating the predicted ΔG* into an analytical form gives rate equations that are in reasonable agreement with non-equilibrium molecular dynamics desorption simulations. We also show that different ΔG* analytical functional forms which give similar predictions at a particular heating rate can yield desorption times that differ by up to a factor of four or more when the ramp rate is extrapolated by 8 orders of magnitude. This highlights the importance of constructing a physically-motivated ΔG* functional form to predict fast desorption kinetics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Low-temperature high-frequency dynamic magnetic susceptibility of classical spin-ice Dy 2 Ti 2 O 7

In this study, radio-frequency (14.6 MHz) AC magnetic susceptibility, , of Dy 2 Ti 2 O 7 was measured using self-oscillating tunnel-diode resonator. Measurements were made with the excitation AC field parallel to the superimposed DC magnetic field up to 5 T in a wide temperature range from 50 mK to 100 K. At 14.6 MHz, a known broad peak of from kHz—range audio-frequency measurements around 15 K for both [111] and [110] directions shifts to 45 K, continuing the Arrhenius activated behavior with the same activation energy barrier of E a ≈ 230 K. Magnetic field dependence of along [111] reproduces previously reported low-temperature two-in-two-out to three-in-one-out spin configuration transition at about 1 T, and an intermediate phase between 1 and 1.5 T. The boundaries of the intermediate phase show reasonable overlap with the literature data and connect at a critical endpoint of the first order transition line, suggesting that these features are frequency independent. An unusual upturn of the magnetic susceptibility at T → 0 was observed in magnetic fields between 1.5 T and 2 T for both magnetic field directions, before fully polarized configuration sets in above 2 T.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Non-equilibrium dynamics in geometrically frustrated spin glass Bi 2 Fe 3 GaO 9 with a Cairo lattice

We have explored the magnetic relaxation of the spin glass phase in Bi 2 Fe 3 GaO 9 , a geometrically frustrated magnet with a unique Cairo lattice made of pentagonal building blocks. Using dc magnetization measurements, we find the relaxation behavior in single crystals follows the Arrhenius law as a function of temperature but remains relatively constant under different magnetic fields. Through a carefully designed protocol, we observe significant rejuvenation and memory effects, reminiscent of classical dilute spin glasses rather than spin jamming systems. Our finding suggests that Bi 2 Fe 3 GaO 9 is probably located in the crossover between conventional spin glass and perfect spin jam and thus provides an excellent platform for investigating the evolution of underlying magnetic structure, paving the way for further neutron scattering research.

aging↗
DOE OSTI2023

Helium bubble size effects on the surface morphological response of plasma-facing tungsten

Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of an atomic layer deposition system for deposition of alumina as a hydrogen permeation barrier

Tritium permeation into and through materials poses a critical challenge for the development of nuclear fusion reactors. Minimizing tritium permeation is essential for the safe and efficient use of available fuel supplies. In this work, we present the design, construction, and validation of custom atomic layer deposition (ALD) and deuterium permeation measurement systems aimed at developing thin-film hydrogen permeation barriers. Using the ALD system, we deposited conformal $\mathrm{Al}_{2}\mathrm{O}_{3}$ films on copper foil substrates and characterized their growth behavior, morphology, and composition. ALD growth rates of ∼1.1 Å/cycle were achieved for temperatures between 100 ∘ C and 210 ∘ C. Permeation measurements on bare and alumina coated copper foils revealed a significant reduction in deuterium flux with the addition of a ∼10 nm $\mathrm{Al}_{2}\mathrm{O}_{3}$ layer. While bare copper followed diffusion-limited transport consistent with Sievert’s law, the alumina-coated samples exhibited surface-limited, pore-mediated transport with linear pressure dependence. Arrhenius analysis showed distinct differences in activation energy for the two transport regimes, and permeation reduction factors exceeding an order of magnitude were observed. These results demonstrate the potential of ALD-grown $\mathrm{Al}_{2}\mathrm{O}_{3}$ films as effective hydrogen isotope barriers and provide a foundation for future studies on film optimization and integration into fusion-relevant components.

atomic layer deposition↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗