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.
Catalytic Site Requirements for N2O Decomposition on Cu-, Co-, and Fe-SSZ-13 Zeolites
N2O decomposition is investigated on Cu, Co and Fe-exchanged SSZ-13 zeolite catalysts at relatively low metal loadings. The catalysts are synthesized by solution ion exchange, and subjected to X-ray diffraction (XRD), temperature-programed-reduction by H2 (H2-TPR), temperature-programed-reaction of N2O (N2O-TPR) coupled with in-situ transmission FTIR, and finally steady-state flow reaction tests. At low N2O pressures (< 0.05 kPa), all catalysts display pseudo first-order kinetics. From Arrhenius analysis, Cu and Fe-SSZ-13 display very different apparent activation energies but similar pre-exponential factors, suggesting their similar reaction mechanisms. N2O decomposition follows a dual-site mechanism, occurring on dimeric M-O-M sites in these catalysts, and O2 is formed by the combination of two O ad-atoms from two vicinal metal sites. Under low N2O pressure (0.05 kPa) and first-order kinetic regime, the reaction is limited by N-O cleavage on bare metal active sites. In comparison to Cu-SSZ-13, the much higher N2O decomposition rate over Fe-SSZ-13 is attributed to the much lower activation barriers for the N-O cleavage step. N2O decomposition occurs on isolated Co2+ ions in Co-SSZ-13. The rate-limiting step is N-O cleavage on an O-occupied Co site in the low-pressure first order kinetic regime. This single-site mechanism leads to much higher pre-exponential factors as compared to the dual-site mechanism. This beneficial factor for reaction rate enhancement, however, is compromised by the much higher activation barriers over this catalyst.
Mean field model parameterization to recover coverage-dependent kinetics
Lateral interactions between adsorbates introduce coverage dependence into adsorption energies and activation barriers of surface reactions. Lattice-based kinetic Monte Carlo (kMC) simulations can capture these interactions quantitatively but are laborious to parameterize and solve. Mean field models are more tractable, but protocols to construct and parameterize them are unclear. Here we explore the ability of a coverage-aware mean-field model to map to a lattice-kMC model of a generic two-step reaction network, including quasi-equilibrated adsorption and rate-limiting dissociation steps. We derive expressions for mean-field and coverage-dependent adsorption energies and dissociation barriers and parameterize against lattice-kMC predictions. We show that the parameterized mean-field rates correlate with ground truth lattice-kMC results across a wide range of reaction conditions and identify regions where the mean field fails. The mean field model similarly captures kMC-derived rate-order, Arrhenius and Sabatier plots at a greatly reduced computational cost. Further, the results provide guidance for parameterizing mean-field models, benchmarked against explicit lattice-based approaches for incorporating the influence of coverage effects.
Bioinspired nonwetting surfaces for corrosion inhibition over a range of temperature and corrosivity
Applications of superhydrophobic (SHS) and lubricant infused surfaces (LIS) involve exposure to corrosive environments from the acidic to the basic, at a range of temperatures, that are not fully characterized. Here, we present for the first time a multifactorial study of the effects of surface fabrication method, surface modification, surface functionalization time, temperature and pH of the immersion medium on the corrosion performance of nonwetting copper surfaces. Bioinspired SHS and LIS fabricated using facile methods of etching and electrodeposition are systematically assessed using potentiodynamic polarization measurements for their corrosion resistance in saline solution (pH ≈7) over a temperature range 23–85°C. SHS and LIS are shown to exhibit diminished corrosion rate, by up to two orders of magnitude, compared to bare copper surface. An Arrhenius model is developed for the first time, describing the temperature- dependent corrosion rate of SHS and LIS. Electrochemical impedance spectroscopy is used to show that corrosion resistance of LIS is larger by three orders of magnitude in extremely acidic (pH = 1) and by an order magnitude in extremely alkaline (pH = 14) media compared to bare copper surface. Etched LIS are generally more resistant to corrosion compared to SHS at all temperatures with excellent microstructural durability.
Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows
This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.
Constitutive behavior and processing maps of a new wrought magnesium alloy ZE20 (Mg-2Zn-0.2Ce)
ZE20 (Mg-2Zn-0.2Ce) 2 is a new wrought magnesium alloy with improved extrudability and mechanical properties [1]. To understand the constitutive behavior and workability of this new alloy, Gleeble thermomechanical testing has been carried out in this study. The flow stress behavior of ZE20 was investigated between 250 °C–450 °C and 10 –3 s –1 –1.0 s –1 in isothermal compression. Constitutive descriptions of the flow stress are provided. A new general approach at application of the extended Ludwik equation is demonstrated and was found to be more accurate than the hyperbolic sine Arrhenius model while having a similar number of model constants. Processing maps were developed based on the experimental results and are verified with microstructural investigation. A region of safe processing with non-basal texture and high activity of dynamic recrystallization (DRX) was found between 375 °C and 450 °C, from 10 –1 s –1 to 10–2.5 s –1 . A region of potentially safe processing with annealing that is associated with shear band nucleation of non-basal grains was identified for temperatures as low as 300 °C and rates as high as 10 –1 s –1 .
Multinuclear PFGSTE NMR description of 39 K, 23 Na, 7 Li, and 1 H specific activation energies governing diffusion in alkali nitrite solutions
While pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has found widespread use in the quantification of self-diffusivity for many NMR-active nuclei, extending this technique to uncommon nuclei with unfavorable NMR properties remains an active area of research. Potassium-39 ( 39 K) is an archetypical NMR nucleus exhibiting an unfavorable gyromagnetic ratio combined with a very low Larmor frequency. Despite these unfavorable properties, this work demonstrates that 39 K PFGSTE NMR experiments are possible in aqueous solutions of concentrated potassium nitrite. Further, analysis of the results indicates that 39 K NMR diffusometry is feasible when the nuclei exhibit spin–lattice and spin–spin relaxation coefficients on the order of 60–100 ms and 50–100 ms, respectively. The diffusivity of 39 K followed Arrhenius behavior, and comparative 23 Na, 7 Li, and 1 H PFGSTE NMR studies of equimolal sodium nitrite and lithium nitrite solutions led to correlations between the enthalpy of hydration with the activation energy governing self-diffusion of the cations and also of water. Realizing the feasibility of 39 K PFGSTE NMR spectroscopy has a widespread impact across energy sciences because potassium is a common alkali element in energy storage materials and other applications.
Cluster dynamics simulation of xenon diffusion during irradiation in UO 2
Diffusion of fission gas in UO 2 nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO 2 under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomic scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO 2 is adequately described by the Xe + U 2 O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U 4 O y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U 4 O y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U 4 O y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.
Determination of neutron irradiation temperatures of SiC using electrical resistivity method
The irradiation temperature of silicon carbide (SiC) was determined post-irradiation by examination of the recovery of the electrical resistivity due to thermal annealing in a rapid heating/cooling optical furnace. High-purity, high-resistivity grade SiC is routinely used as a passive temperature monitor in neutron irradiation studies at the High Flux Isotope Reactor (HFIR), and this paper presents an alternative automated technique for determination of the irradiation temperature the SiC experienced. Neutron collisions with the atoms results in displaced lattice atoms (interstitials) that act as electron donors yielding a significant decrease in electrical resistivity. The irradiation defects become thermodynamically unstable and start to recombine, when annealed above the irradiation temperature, resulting in a recovery of the electrical resistivity. The resistivity is measured at a fixed elevated temperature above ambient, which is below the target irradiation temperature. When the resistivity is plotted as a function of annealing temperatures, a clear increase is observed due to the recovery of irradiation defects. We have demonstrated that this electrical resistivity measurement of SiC is effective to determine irradiation temperature of SiC. In conclusion, energy levels of various defects in SiC were calculated from Arrhenius plots of electrical conductivity versus inverse temperature.
Fission gas diffusion and release for Cr 2 O 3 -doped UO 2 : From the atomic to the engineering scale
Here, the anticipated benefits of large grains in Cr 2 O 3 -doped UO 2 pellets include improved mechanical and fission gas retention properties. To support the assessment of fission gas release (FGR) from doped pellets, the impact of doping on fission gas diffusivity for in-reactor conditions must be understood. In this work, we tackle this issue by informing the fission gas model within the BISON fuel performance code using material models developed at the atomic scale. The investigation of intra-granular fission gas diffusivity in Cr 2 O 3 -doped UO 2 is carried out by adapting a cluster dynamics model that, accounting for UO 2 thermochemistry, is capable of describing Xe diffusion under irradiation in undoped UO 2 as the starting point. Using a thermodynamic analysis, it is shown that in stoichiometric UO 2 with additions of Cr 2 O 3 the oxygen potential is defined by the Cr-Cr 2 O 3 two-phase equilibrium. Using the cluster dynamics model, the predicted Xe diffusivity in doped UO 2 was significantly increased in both the intrinsic and irradiation-enhanced regimes compared to undoped UO 2 as a result of higher concentrations of uranium and oxygen vacancies, respectively. This is a consequence of the more oxidizing conditions at high temperature, and more reducing conditions at low temperature, as a result of doping. Arrhenius functions have been fitted to the cluster dynamics results to enable implementation of the new diffusivities in the BISON fission gas behavior model. BISON simulations were carried out, showing the competing effects of the enlarged grains and the new fission gas diffusivity model, which act to suppress and enhance fission gas release, respectively. The new physics-informed model was validated against in-reactor experimental measurements under normal operation. Additionally, benchmarking was carried out for power ramp conditions. The predicted fission gas release agreed well with the experimental data, showing noticeable improvements over the standard UO 2 model.
Determining the oxidation behavior of matrix graphite
This work presents the oxidation behavior of matrix graphite in air. Matrix graphite, graphite powder/flakes bonded by a small amount of non-graphitic carbon, surrounds coated fuel particles in order to form cylindrical fuel compacts (in prismatic core designs) or spheres (in pebble-bed reactor designs). This work focuses on oxidation tests conducted on two matrix graphite materials, one provided by Kairos Power and the other A3 matrix graphite. Some of the tests followed American Society for Testing and Materials (ASTM) oxidation testing standards using a vertical furnace system and others were performed in a thermogravimetric analyzer (TGA). It was determined that, at temperatures of 450 °C–700 °C, the oxidation rate of the Kairos matrix graphite follows the Arrhenius equation. In comparison with A3 matrix graphite, the Kairos matrix graphite shows better oxidation resistance at high temperatures (≥550 °C), but also a higher oxidation rate at low temperatures. Both the A3 matrix graphite and the Kairos matrix graphite materials may experience preferential oxidation of the partially graphitized binder. An oxygen penetration gradient was also observed when using the three characterization methods (i.e., optical microscope, x-ray tomography [XCT], and density profile by the lathe) enlisted in this research. In conclusion, the oxygen penetration depth increases with decreasing isothermal oxidation temperature, while the center of the oxidized samples (10% weight loss) remains almost untouched even at 500 °C.
In-situ determination of strain during transient burst testing and the temperature dependence of Zircaloy-4 claddings
Understanding fuel system behavior during postulated loss-of-coolant accidents is pertinent for continued safe and efficient operation of light water reactors, particularly as higher burnups are being pursued and safety margins re-evaluated. Conventional mechanical models for the incumbent Zr alloys typically rely on the assumption that steady-state creep is the dominant fuel cladding response during transient accident conditions. To investigate this assumption, simulated accident burst testing was performed on Zircaloy-4 claddings with balloon behavior measured in-situ. Here, two distinct loading conditions were utilized during burst testing: (1) constant-gas-inventory where pressure was allowed to increase with temperature and (2) constant pressure. In-situ strains and strain rates were measured via 2-dimensional digital image correlation techniques and synchronized with temperature to determine deformation dependencies. The temperature dependence of strain rate was characterized by a two segment Arrhenius relationship, with a distinct transition between the high and low temperature/strain regimes. The average activation energy of the lower temperature/strain regime was 328 ± 25 kJ/mol, in agreement with the ~320 kJ/mol used for conventional LOCA models. However, the higher temperature/strain segment, which encompassed most of ballooning, showed increased activation energies as well as a dependence on whether the burst region was in view. For tests that burst away from the camera view, the average high temperature/strain segment activation energy was 635 ± 150 kJ/mol. For samples where the rupture opening formed in view, the average activation energy was 1015 ± 179 kJ/mol. This observed shift in temperature dependence indicates a transition in deformation mechanism at the end of life, possibly to time independent failure mechanisms, which has not yet been visualized in the literature for Zr alloys. Parameters at the transition points were analyzed to determine thresholds for this change in behavior, which occurred at an average hoop strain of 6.9 ± 2.1 %.
Impact of temperature variations on BISON predictions of Ag release in AGR-1 and AGR-2 experiments
Understanding and quantifying the release of fission products like silver (Ag) from TRistructural ISOtropic (TRISO) fuel particles is important to assess the safe operation of advanced high temperature reactors. Although the silicon carbide (SiC) layer of TRISO particles is effective as the main fission product barrier, Ag can be released from intact TRISO particles. A mechanistic model for the effective Ag diffusivity, D eff , was previously developed as a function of temperature and microstructure variables informed by atomistic modeling of Ag diffusivity on the mesoscale. Here in this study, we use this model to explore how experimental temperature uncertainties impact the overall predicted Ag release. This analysis shows that temperature uncertainties have a significant impact on the overall Ag release predictions. Furthermore, we show that the time average volume average temperature (TAVA) temperature is not an appropriate proxy for temperature histories to predict fission product release. We attribute this to the Arrhenius dependence of Ag diffusivity with respect to temperature. The detailed temperature histories, therefore, provide the most accurate results are are of most importance for modeling efforts. Overall, this work shows the importance of considering the experimental uncertainty of the temperature on computational predictions of fission product transport and release and the need for more accurate temperature histories from future experiments.
First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport
The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.
Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries
Lithium plating is one of the main concerns for cell durability, as cyclable lithium and cell capacity are decreased during the process. Thermal and mechanical loads can influence cell behavior, including lithium plating driving forces, as suggested by the Arrhenius equation and the Bruggeman relation, respectively. Hence, a 1D and a 3D model are employed to study the thermal and mechanical effects on lithium plating in a lithium-ion battery with a high areal capacity (4 mAh/cm 2 ). The 1D model assumes uniform mechanical and temperature profile, while the 3D model is capable of capturing the nonuniformity of stress and temperature in the cell. Further, the effects of mechanical deformation and different temperatures on the process of lithium plating are investigated in this study. With the help of the 1D charging simulations, it is observed that lower temperatures increase the risk of lithium plating and the presence of mechanical deformations can accelerate its initiation. Furthermore, 3D charging simulations indicated that folds and boundaries of jellyrolls are more susceptible to lithium plating due to the higher stress concentration and lower temperature, respectively.
The influence of temperature on area-specific impedance and capacity of Li-ion cells with nickel-containing positive electrodes
The temperature-dependent behaviors of five nickel-containing positive electrodes (NCA, NMC811, NMC622, NMC532, and NMC111) in lithium-ion batteries are investigated using an electrochemical protocol involving rate studies, mild aging (~100 cycles), and hybrid pulse power characterization (HPPC). Tests are conducted using coin-cells with graphite negative electrodes at -20 °C, 0 °C, 20 °C, and 40 °C. Three techniques are compared for determining the area-specific impedance (ASI): i) fits to the rate study average voltages, ii) fitting to the entire voltage curves using a regularization scheme, and iii) HPPC. When fit to an Arrhenius-type equation, all methods yield similar apparent activation energies (±2 kJ/mol) for the impedance, which range from -20 to -31 kJ/mol for the electrodes. Impedance growth increases with temperature but remains at less than 0.2% per cycle for most electrodes and temperatures. NCA and NMC811 are the exceptions, which yield 0.5% and 1.5% increases in ASI per cycle, respectively, at 40 °C. For cells with the same electrodes, the capacities are similar at 20 and 40 °C but reduce at lower temperatures, with up to a 52% reduction at -20 °C and 2C. The fade in energy of the cells during C/3 cycling is attributed to decreasing capacity as opposed to increasing ASI.
Uncovering metallic glasses hidden vacancy-like motifs using machine learning
Vacancies are a ubiquitous type of defect in crystalline materials that can affect their properties and behavior, including atomic diffusivity, deformation processes, thermal conductivity, and electrical conductivity. While well-defined in crystals, it is challenging to identify such topological footprints in the amorphous structure of metallic glasses (MGs). Here, we uncover an unforeseen vacancy-like structural motif named Q7, which refers to an atomic Voronoi polyhedron with seven quadrangular faces, by investigating MG local atomic configurations with the help of machine learning. The Q7 motif exhibits characteristics similar to crystalline vacancies and plays a significant role in short-range structural disorder in MGs. The atoms at the center of the Q7 motif display larger local entropy, atomic volume, and local tension. Additionally, its concentration follows an Arrhenius-like relationship with temperature, accurately indicating the glass transition temperature, and is strongly correlated with the yield and failure of MGs during mechanical deformation. The discovery of the Q7 motif provides new insights into understanding the intricate relationship between local disorder and structural relaxation in MGs, offering a remarkable ability to predict their thermal and mechanical behavior. This potentially paves the way for enhancements in the design and manufacturing processes of MGs and other amorphous materials.
Prediction of nonlinear specific heat during single crystal HMX phase transition
Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.