Search NASA⌕ Search

SEARCH · Search NASA

Results for “initialization”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Shadow of the Future: Developing Trust and Software within the Exascale Computing Project

Collaboration and team science are emerging areas of interest in software production. Historically, multi-institutional research collaborations are difficult to initiate and maintain, negatively impacting communication, negotiation, and dialogue between industry, government, and academic researchers. The Exascale Computing Project (ECP), a massive, multi-team, high-stakes initiative, facilitated broader research collaboration under a shared funding structure and extended timeline to support scientific discovery. Here, we conducted interviews with ECP teams, representing a variety of domain specialties, research institutions, and programming backgrounds. Using thematic analysis, we assessed how ECP’s structure created an environment of increased trust among projects and how software shared between teams facilitated sustained collaboration. We found that the expectation of future collaboration, i.e., the shadow of the future, greatly enhanced trust among teams and the quality of scientific software produced. Based on our findings within ECP projects, we connect to the existing literature on trust in software engineering and share recommendations for sustainable multi-institutional collaboration and shared best software practices.

Exascale computing project↗
DOE OSTI2025

Temperature dependence of uranium and thorium partitioning in igneous zircons

Zircon is a key mineral in geochronology because of its chemical and physical durability and tendency to incorporate radioactive trace elements such as U and Th. Quantifying the partitioning of the actinide elements is critical to constrain initial non-secular equilibrium amounts of 230 Th in zircon. An excess or deficit of 206 Pb will be produced from such an initial excess/deficit of 230 Th from the secular equilibrium condition, which influences the calculated 206 Pb/ 238 U age (Schärer, 1984, Mattinson, 1973). However, there is no standard way to calculate Th/U partitioning ratios when applying age corrections to young igneous zircon, making uncertainties hard to estimate. To better understand U and Th partitioning of zircon/melt that is the cornerstone of the secular disequilibrium correction, zircon was synthesized in one-atmosphere experiments using basaltic andesite, andesite, and rhyolite starting materials, doped with Zr, U, and Th. Different experimental temperatures and oxygen fugacity conditions (ΔQFM−4 to ΔQFM+4) were explored to examine their effects on U and Th partitioning. In addition, we specifically quantify the effects of sector zoning, fractional crystallization, and melt composition on U and Th partitioning. By combining experimental and natural zircon data, we find that temperature has the primary control on the partitioning of U and Th in the zircon and calibrate an inverse relationship between these partition coefficients and crystallization temperatures. The calibrated equation can be applied to the 230 Th correction for an improvement in the accuracy of Th-corrected 206 Pb/ 238 U dates when the zircon crystallization temperature is known.

206Pb/238U dating↗
DOE OSTI2025

Experimental platforms for investigating feature-driven jets for HED mix model validation

High-energy-density (HED) systems, such as inertial confinement fusion (ICF), are susceptible to hydrodynamic instabilities that can significantly affect both experimental results and modeling predictions. Isolated features, such as fill tubes or divots in the capsule, can cause material to jet as a result of the compressive shock exciting the Richtmyer–Meshkov instability, and serve as one of the primary degradation mechanisms in ICF yield. Simulations of feature-driven jets and how they mix require extensive experimental validation, particularly for understanding to what degree the initial size and shape of a feature influence jet dynamics, and how much instability feeds through downstream layers. A better understanding of feature-driven jetting can improve our mix modeling capabilities and increase hydrodynamic simulation accuracy. This manuscript describes a series of experimental platforms fielded by Los Alamos National Laboratory as a part of the Mshock Omega 60 and ModCons Omega EP campaigns to explore feature-driven jetting. These platforms are designed to benchmark jet evolution and growth as a function of initial feature size and shape, investigate jet-layer interactions leading to instability feedthrough, and will be used to characterize jet-jet interactions resulting from clusters of features. In conclusion, preliminary results for both platforms are shown. The ModCons experiments are on-going, and a discussion of future work directions is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Extraction of lithium from battery recycling wastewater using synergistic D2EHPA and TBP

The recycling of spent lithium-ion batteries (LIBs) poses significant challenges, including the generation of large volumes of chemically complex wastewater. The composition of this wastewater is influenced by both the intrinsic chemistry of the batteries and the specific recycling processes employed. Notably, this wastewater contains economically valuable components, such as lithium, which can be recovered. Here, in this study, a solvent extraction (SX) process was investigated as a method to recover lithium from battery recycling wastewater, especially from anode-washing stage. Initially, various commercial extractants were evaluated, including di(2-ethylhexyl)phosphoric acid (D2EHPA), mono-2-ethylhexyl (2-ethylhexyl)phosphonate (PC88A), bis(2,4,4-trimethylpentyl)phosphinic acid (Cyanex 272), 2-hydroxy-5-nonylacetophenone oxime (LIX 84-I), tri-butyl phosphate (TBP), and their combinations. Among these, D2EHPA + TBP demonstrated a synergism to advance and maximize the lithium extraction. Subsequently, the effects of key parameters, including D2EHPA concentration, TBP concentration, contact time, initial pH, and aqueous-to-organic (A/O) phase ratio, were systematically investigated and optimized. A two-stage SX approach was employed to enhance lithium recovery. Under the optimized conditions of 30.0 vol% D2EHPA, 10.0 vol% TBP, 10 min of contact time, and a 1:1 A/O phase ratio, a lithium extraction efficiency of more than 88% in a two-stage solvent extraction was achieved. Lithium was subsequently stripped from the loaded organic solution using sulfuric acid (H 2 SO 4 ). Using 2.0 M H 2 SO 4 , lithium stripping was achieved after two counter-current stripping stages at an organic-to-aqueous (O/A) phase ratio of 6:1. This stripping process enriched the lithium concentration by a factor of four compared to the original lithium concentration in the anode-washing wastewater. The recyclability of the synergistic D2EHPA + TBP system was also evaluated over four extraction-stripping cycles. The results demonstrated that the system maintained high extraction and stripping efficiencies.

D2EHPA↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗
DOE OSTI2025

High and Ultra-High Temperature Reaction Kinetics by Single Nanoparticle Mass Spectrometry

Methodology is presented for non-destructive, optically-detected single nanoparticle (NP) mass spectrometry, with the goal of extracting surface reaction kinetics for single NPs at high temperatures. Methods for determining the NP charge, mass, and temperature as a function of time are discussed, and the data are used to extract both the absolute kinetics for mass change, as well as the efficiencies of the surface processes that cause them. Factors that contribute to the uncertainties in absolute and relative mass determination, and in the resulting kinetic parameters, are discussed. The method allows the NP-to-NP variations in initial reactivity to be measured directly, along with the time evolution of reactivity resulting from NP structural/compositional changes that occur under reaction conditions. The strengths and limitations of single nanoparticle mass spectrometry as a high temperature surface kinetics tool are discussed in the context of sublimation and O2 oxidation kinetics experiments for single hafnium (Hf) NPs at temperatures ranging above 2400 K. The Hf oxidation kinetics are compared to analogous oxidation experiments for silicon, graphite, and carbon black NPs. In all four cases, the oxidation chemistry was dominated by processes that result in net mass loss, and the distinct mechanisms responsible are discussed. All four NPs also eventually passivated, i.e., the efficiencies for oxidative etching decreased by at least two orders of magnitude, relative to the initial efficiencies. Furthermore, the passivation mechanisms, which are quite different for carbon, compared to silicon or hafnium, are discussed. Carbon NP passivation is attributed to structural isomerization leading to fully coordinated, fullerene-like NP surfaces, while for silicon and hafnium, passivation results from delayed formation of an oxide layer, triggered by accumulation of oxygen in the NP sub-surface region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The role of annealing and grain boundary controls on the mechanical properties of limestones and marbles

Chemical and mechanical processes are coupled in many geological and geochemical environments. Reactive processes anneal defects and restructure grain boundaries, modifying their elastic properties, levels of internal friction, wave propagation rates and fracture behaviors. The nature of these changes is, however, contingent on the initial state of the rock. Here, in this study, impulse excitation was used to measure changes in mechanical properties as a function of dry and steam heating time at 300 °C for three carbonate rocks: Carrara marble, Carthage marble (Burlington Limestone), and Texas Cream limestone (Edwards limestone), with initial porosities of about 1 %, 3 %, and 27 %, respectively. Frequency-dependent phenomena along with mineral recrystallization were observed. This was coupled with small-angle X-ray and neutron scattering analysis to determine the relationship between changes in bulk mechanical and microstructural properties. An observed decrease in both the Young's and shear moduli in the experimental limestones as a function of heating time reflects and quantifies a reduction in the stiffness of the rock due to annealing. The internal friction of the samples first increases then decreases with reaction time, reflecting defect annealing, but this was balanced against grain boundary dissolution apparently driven by condensation of steam in the confined grain-boundary environment. This suggests an increase in the boiling point under confinement leading to dissolution, increased porosity and widening of the grain boundaries. In addition, comparison of small-angle X-ray and neutron scattering results suggests that it is inappropriate to assume that pores are empty for quantitative analysis of typical small-angle scattering samples.

58 GEOSCIENCES↗
DOE OSTI2026

Cleavage at the nsp5–nsp6 site of SARS-CoV-2 main protease intermediate precursor is faster from a monomer than a dimer form

Our previous studies of severe acute respiratory syndrome coronavirus 2 main protease (MPro) precursor monomer indicate that the initial N-terminal nonstructural protein (nsp)4/nsp5 cleavage occurs intramolecularly, with a small fraction of the active site loop equilibrium being in the active state. To understand the influence of dimer formation of MPro upon N-terminal cleavage on the subsequent C-terminal nsp5/nsp6 intermolecular cleavage kinetics, the stepwise processing of a monomeric, inactive precursor containing the native terminal cleavage sites of MPro (MBP- (−6) MPro C145A(+3) -GB1-6H, 86.2 kDa) by mature WT MPro (MPro WT ) was investigated. Differential scanning fluorimetry and analytical ultracentrifugation measurements of various MPro constructs suggest that the C145A mutation decreases the dimer dissociation constant (K dimer ) by ∼26-fold, relative to WT C145 and H41A. The monomeric precursor’s nsp4–nsp5 site appears to saturate MProWT’s active sites and cleave faster, followed by a slower first-order cleavage at the C-terminal site. No detectable product resulting from the C-terminal cleavage is observed until most of the N-terminal cleavage is complete. The initial intermediate product (termed MPro C145A-IP ) is a homodimer with an estimated K dimer of <0.05 μM. In contrast, the first-order kinetics observed for the cleavage of the monomeric form of the intermediate product is at least 300 times faster than that of the dimer form. Room-temperature X-ray structure of the MPro C145A-IP –ensitrelvir complex is like that of the MPro WT –ensitrelvir complex and reveals a dynamic C-terminal region including MPro residues 302 to 306. These results are interpreted from the point of view of a mechanism in which nsp5–nsp6 cleavage may occur from a monomeric intermediate, and dimer formation restricts this cleavage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗
DOE OSTI2024

Effect of die design on microstructure and mechanical joint strength in friction self-piercing riveted AA7055-T76 and AA7055-T76

In the present study, a unique friction self-piercing riveting (F-SPR) was applied to join high-strength, low-ductility aluminum alloy (AA) 7055-T76 and AA7055-T76 to mitigate a cracking issue. Also, the effect of die design on both joint formation and the mechanical joint performance were investigated. A crack-free joint was achieved when joining AA7055-T76 and AA7055-T76 by F-SPR. As the bottom diameter of the die increased, the mechanical interlocking distance initially increased then became almost saturated. At the same time, the solid-state bonding gradually disappeared, and the lap-shear fracture changed from bottom aluminum fracture to rivet pull-out with partial bottom aluminum fracture. The average maximum lap-shear joint strength of 11.57 kN and cross tension of 5.65 kN were achieved from the D2 die design owing to the combined contributions from high mechanical interlocking distance and good solid-state bonding at the joint interface. Ultrafine grain refinement was observed at the region next to the rivet shank outer surface and along the materials flow line between the top and bottom aluminum sheet in the rivet cavity. Further, because of the same joint-process parameters, the hardness profile at the joints made by different dies did not show obvious differences. As the distance from the base metal to the rivet outer surface decreased, frictional heat caused the hardness to initially decrease from 190 to 150 HV in the heat-affected zone. Hardness increased up to 175 HV in the thermomechanical affected zone (TMAZ-I) owing to large plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of post-weld heat treatment on microstructure and mechanical properties of rotary inertia friction welded dissimilar 422-4140 martensitic steel piston joints

The study investigated how post‐weld heat treatment (PWHT) temperature affects the microstructure and localized deformation/fracture during bend testing of rotary inertia friction welds (RIFW) between AISI 422 stainless steel and AISI 4140 steel. RIFW produced a fully martensitic interface with approximately 550 HV hardness in both the thermo-mechanically affected (TMAZ) and heat‐affected zones (HAZ). Due to differences in temper resistance, the 4140 TMAZ/HAZ softened progressively under PWHT temperatures from 525 °C to 700 °C, while the 422 TMAZ unexpectedly maintained about 550 HV up to 600 °C before significantly softening at temperatures ≥625 °C. This asymmetric softening generated steep hardness gradients across the interface at temperatures ≤600 °C. Furthermore, carbon migration across the interface was minimal up to 600 °C, moderate at 625 °C, and by 700 °C produced a carbide‐rich eutectoid layer in the 422 TMAZ alongside a carbon‐depleted soft ferrite layer in the 4140 TMAZ. Strain during bending was PWHT‐dependent, concentrating on the 4140 side; in as-welded joints, the high hardness led to deformation and crack initiation in the base metals, whereas in PWHT samples, cracking initiated in the softened 4140 TMAZ near the interface. The intermediate PWHT temperature of 625 °C offered the best balance of limited carbon diffusion across the interface, relatively low peak weld hardness and minimized hardness gradients across the interface, more homogenous deformation, and good bend test performance.

4140 low-alloy steel↗
DOE OSTI2024

Phase stability and microstructure of neutron-irradiated substoichiometric yttrium dihydrides

The impact of the neutron-displacement damage on phase stability and microstructure of substoichiometric yttrium dihydrides (YH x , x <2) were investigated to assess their use as solid moderator in high-temperature nuclear reactors. YH x specimens were, thus, subjected to neutron irradiations in the range of 0.1–2 displacements per yttrium atom (dpa-Y) in the temperature range of 536–878°C at the Oak Ridge National Laboratory's (ORNL's) High Flux Isotope Reactor (HFIR). YH x specimens were initially prepared at stoichiometry (H/Y) ratios of 1.69 and 1.83. HFIR-irradiated specimens were characterized by variety of techniques to investigate H retention characteristics including dimensional analysis, optical microscopy, scanning electron microscopy electron back scatter diffraction (EBSD), transmission electron microscopy, thermal desorption spectroscopy (TDS), and high-energy x-ray diffraction (HE-XRD) characterizations. Overall, YH x exhibited notable structural and phase stability under short-term neutron-irradiation, except for the samples with significant silicon carbide (SiC) interaction at high doses and temperatures. Basic dimensional and mass measurements were misleading for accurate assessment of H retention, as confirmed by EBSD phase maps, XRD line profiles, and TDS signals. Thus, it was discussed that a robust H retention metric is needed to assess irradiated hydrides. Further, nanoscale cavities were observed as a result of the neutron irradiation in all samples. Although no clear impact of dose and irradiation temperature was determined, the initial H/Y ratio had an impact on the cavity number density where low H/Y specimens had high-resistance to cavity formation. The Y-vacancy cluster formation at the collision stage of the displacement cascade and their stabilization by H were considered to be the likely underlying mechanisms for the observed cavity microstructure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Application of pH-dependent environmental protection agency test method 1313 to low-activity nuclear waste glass

Here, in this study, we present results of a room temperature, pH-dependent leach test established by United States Environmental Protection Agency (EPA) Leaching Environmental Assessment Framework (LEAF) (EPA Method 1313 and further abbreviated EPA 1313) to study the alteration behavior of low-activity nuclear waste (LAW) glasses. Implementation of EPA 1313, which is closer to LAW disposal temperatures, could allow for more flexibility in the formulation of Hanford LAW glasses, specifically in terms of increased waste loading. A set of 48 LAW glasses were selected to test with EPA 1313. Briefly, the test method consists of adding glass to a number of solutions at different initial pH values and at a fixed mass-to-solution-volume ratio. The glass is then allowed to corrode for two days and the final solution pH and concentration of different glass elements in solution is measured. The test was modified from the published procedure to constrain particle size to provide more reproducible results, reduce the studied pH range to a more applicable range, and to reduce solution volume to reduce costs associated with material preparation. Compositional modeling was used to predict the measured elemental releases at various pH. The model predictions improved in tests where there was a relatively small difference between initial and final solution pH indicating that pH controls would result in better predictive models.

Accelerated aging↗
DOE OSTI2025

Enabling BWR fuel rod analysis in the BISON fuel performance code

Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

NOx mediated formation of ε – UO 3 during thermal decomposition of uranyl nitrate hexahydrate under static air conditions

The thermal decomposition of uranyl nitrate hexahydrate in air at temperatures of 650 – 800 °C is expected to yield α – U 3 O 8 through intermediate UO 3 phases. Here, n this study, Raman spectroscopy complimented by powder X-ray diffraction revealed the unexpected formation of ε – UO 3 as an accompanying phase during thermal decomposition of UNH at 700 °C under static air conditions. Experiments demonstrated that ε – UO 3 does not form during the initial ramp up stage nor during the high temperature heating period; but instead forms during the cooling stage between 250 – 400 °C. Additional experiments revealed that UNH initially decomposes through an amorphous UO 3 intermediate prior to α – U 3 O 8 formation. Attempts to bypass the U 3 O 8 precursor in the formation of ε – UO 3 were unsuccessful, supporting the necessity of U 3 O 8 (α – U 3 O 8 in this study) in the formation of ε – UO 3 . The observed phase evolution is proposed to result from the NOx species generated during UNH thermal decomposition, which create localized oxidizing conditions within the furnace under static air conditions. Furthermore, a predominantly phase pure ε – UO 3 was synthesized under static air conditions through the addition of an extra plateau at 250 °C during the cooling step. These findings demonstrate the importance of gas-phase chemistry and cooling conditions in uranium oxide phase evolution and provide additional insight into ε – UO 3 formation pathways.

Epsilon uranium trioxide↗