Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Graphite”

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

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

AGR-3/4 TRISO Fuel Compact Ceramography

The combined third and fourth irradiation in the Advanced Gas Reactor (AGR) program (AGR-3/4) contained tristructural isotropic (TRISO)-coated particle fuel and designed-to-fail (DTF) fuel particles. The DTF particles were only coated with highly anisotropic pyrocarbon (PyC) so they would purposely fail during the AGR-3/4 irradiation and provide a source of fission products for measurement. To observe the post-irradiation morphology of these DTF particles and the TRISO-coated particles, three AGR-3/4 compacts were mounted in epoxy, sectioned above their centerlines, and ground/polished. Three rounds of grinding/polishing and optical microscopy were performed so that the particles could be observed at multiple planes. Each compact contained approximately 1,872 TRISO particles and exactly 20 DTF particles. The three AGR-3/4 compacts examined covered a wide temperature range and featured both the hottest average irradiation temperature (1375°C) and the coldest/lowest burnup (872°C and 5.5% fissions per initial metal atom) of any compacts to undergo post-irradiation examination (PIE) in the AGR program to date. A total of 29 DTF particles were located and observed via microscopy across the three compacts. All the DTF particles observed via microscopy had completely failed PyC coatings. Some different DTF kernel and PyC morphologies were observed that have been attributed to the differences in irradiation temperature and/or burnup. At low and medium irradiation temperatures (roughly 850 to 1050°C), it appears irradiation-induced dimensional changes in the DTF PyC caused it to fracture and fold in on itself, and the kernel deformed to accommodate this PyC deformation. The extent of DTF kernel deformation to accommodate DTF PyC buckling tended to be greater for the medium-temperature fuel compared to the cold fuel. In the high-irradiation-temperature compact (1375°C), the DTF PyC layer appears to have completely reacted chemically with the kernel material. The only material surrounding the DTF particles in this hot compact resembles the compact matrix material. Observations of the TRISO-coated driver fuel particles and fuel compact graphitic matrix were also made. The TRISO particle kernels in the medium-temperature compact (1047°C), and many in the high-temperature compact (1375°C), showed morphologies consistent with what was seen in AGR-1 and AGR-2. In other TRISO particles from the high-temperature compact, a spatial gradient in the kernel appearance was evident. Here the cool side of the kernel (away from the center of the compact) had a much darker appearance (like that of the buffer and pyrocarbon). This is the first observation of this kind of kernel spatial gradient in AGR fuel. It is believed that the high irradiation temperature of that compact (an average of 1375°C) activates or accelerates an unidentified chemical reaction, and the temperature gradient within the fuel causes a sharp spatial gradient. The low-temperature/low-burnup compact also showed unique kernel morphologies where the TRISO-coated particles had clearly distinguishable oxide and carbide phases at the center of the kernel, an oxide rind surrounding this, and remnants of a carbide skin surrounding the oxide rind. These are features observed in the as-fabricated fuel that are no longer present in higher burnup fuel. Occasional gaps between the outer pyrolytic carbon (OPyC) and graphitic matrix material were observed in all compacts. Gaps were most often found in the small, matrix-filled spaces between adjacent particles. Finally, buffer morphologies in AGR-3/4 TRISO particles were like those observed in AGR-1 and AGR-2, and the buffer fracture frequencies in AGR-3/4 and AGR-2 TRISO particles were plotted versus fast neutron fluence and irradiation temperature. In AGR-3/4 the buffer fracture rate was highest (23%) for an irradiation temperature of 1047°C and a fast fluence of 5.18E25 n/m 2 . Increasing the irradiation temperature to 1375°C reduced the buffer fracture rate to 14%. Reducing the irradiation temperature to 872°C and the fast fluence by a factor of 3 reduced the buffer fracture rate to 7%. The temperature and fluence dependencies observed for AGR-3/4 buffer fracture were consistent with those observed in AGR-2. Lower fluences and/or higher temperatures significantly reduced the buffer fracture frequencies. Lower buffer fracture frequencies were also observed at lower temperatures as long as the fluence was also significantly reduced. This is believed to be due to lower fluence resulting in less irradiation-induced dimensional change (shrinkage) of the buffer, and a higher temperature promoting enhanced creep relaxation of stresses within the buffer, leading to less buffer fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

2023 Critical Materials Strategy

The global effort to curb carbon emissions is accelerating demand for clean energy technologies and the materials they rely on. Demand for these materials will only continue to grow, especially as some nations aim to achieve net zero emissions by 2050. While some major materials like steel, copper, and aluminum are already powering the fossil fuel economy, others are more minor materials with potential supply risks. These risks could jeopardize the ability to reduce greenhouse gas emissions within the desirable timeframe to avoid significant climate change. In some cases, it may be necessary to take action to improve the resilience of material supply chains and mitigate supply risks. Understanding the importance of individual materials to clean energy and the supply risks associated with them is necessary to identify which materials may serve as potential roadblocks to a clean energy future. The U.S. Department of Energy (DOE) issued a series of 13 supply chain deep dive assessment reports on various energy technologies in 2022 in response to President Biden’s Executive Order on America’s Supply Chains (E.O. 14017). These reports emphasized that supply chain bottlenecks can occur at any stage of the value chain from mining and refining to component and even sub-system manufacturing. The bottlenecks are a combination of factors such as material availability, equipment availability, work force availability and quality, logistics, regulatory framework, and market conditions. These bottlenecks were worsened during the global Covid-19 pandemic. Its lingering impacts have hindered capacity expansion for material supply chains and prevented product lead-time recovery. One approach to reduce supply chain risks for the United States is to have a strong domestic manufacturing sector with a diverse set of producers. Boosting responsible domestic production would require leveraging the latest science not only in material extraction but also in developing substitutes, recycling, reuse, and remanufacturing. This report is an updated analysis of previous Critical Materials Strategy (CMS) reports published by the DOE in 2010, 2011, and 2019 based on national and global priorities, technology advancement, and technology adoption trends. Like the CMS reports, this analysis presents the results of a formal material criticality assessment to identify which materials are critical to the continued deployment of clean energy technologies globally. The analysis in this report leveraged the DOE supply chain deep dive assessments to develop the initial list of materials to evaluate. This DOE Critical Materials Assessment (CMA) is conducted independently of criticality assessments performed by other U.S. government agencies, such as that conducted by the U.S. Geological Survey (USGS). This analysis complements the USGS critical minerals determination in three aspects. First, the DOE assessment is performed from a global perspective, while the USGS analysis focusses on the importance of minerals to the U.S. economy. Second, this report focuses on the importance of materials to clean energy technologies, rather than to the economy in general. Lastly, this study is forward looking to 2035 based on clean energy deployment scenarios, whereas the USGS assessment is retrospective. Materials evaluated in this report that do not appear in the USGS Critical Minerals List include copper, uranium, electrical steel, and SiC. A draft version of this report received ~80 public comments related to supporting data and methodological improvement. Those comments have been incorporated as much as possible where appropriate. Highlights of findings from this 2023 CMA include: Rare earth materials (neodymium, praseodymium, dysprosium, and terbium) used in magnets in electric vehicle (EV) motors and wind turbine generators continue to be critical. While dysprosium (Dy) and terbium (Tb) are both heavy rare earth elements that serve the same function in magnets, the criticality of Tb is slightly lower than that for Dy in the short term due to the widespread use of Dy in high-grade magnets and Tb’s present role as a substitute. Similarly, praseodymium (Pr) is critical in the medium term but only near critical in the short term because it is more substitutable in magnets than neodymium (Nd); Materials used in batteries for EVs and stationary storage are now considered to be critical. While cobalt (Co) was found to be critical in this and previous reports, lithium (Li) becomes critical in the medium term due to its broader use in various battery chemistries and the rampant growth of the EV industry. Natural graphite is a new addition in this assessment and is also found to be critical; Platinum group metals used in hydrogen electrolyzers, such as platinum (Pr) and iridium (Ir), are critical due to an increased focus on hydrogen technologies to achieve net zero carbon emissions, while those used in catalytic converters, such as rhodium (Rh) and palladium (Pd), were screened out due to the decreased importance of catalytic converters in the medium term; Gallium (Ga) continues to be critical due to its use in light-emitting diodes (LEDs). In addition, the use of Ga has increased in magnet manufacturing and in semiconductor in forms such as gallium arsenide (GaAs) or gallium nitride (GaN); Major materials like Aluminum (Al), copper (Cu), nickel (Ni), and silicon (Si) move from noncritical in the short term to near critical in the medium term due to their importance in electrification; Electrical steel is near critical due to its use in transformers for the grid and electric motors in EVs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Investigating the AGN-201M research reactor's unique dominance ratio

The AGN-201M reactor is a small, thermal spectrum reactor located at the University of New Mexico (UNM). It is moderated by polyethylene, reflected by graphite, with fuel comprised of uranium microspheres embedded in polyethylene plates that are separated by an aluminum baffle. The dominance ratio of a multiplying system is the ratio of the first higher-mode eigenvalue to the fundamental, k{sub 1}/k{sub 0}. It is a key parameter in characterizing the rate of convergence, analyzing system stability, and determining how tightly coupled the system's neutron fission is. The purpose of this paper is to investigate the system's unusually low dominance ratio of 0.632. This investigation varied the reactor's physical parameters, such as the fuel geometry, fuel density, and reflector thickness, and examined their impact on the reactor's dominance ratio. Additionally, this study examined neutronically similar systems in hopes of identifying common causes for systems with low dominance ratios. The reason why the dominance ratio of the reactor is small compared to larger thermal reactors was determined to be because of its significantly smaller geometry and the composition of the fuel plates. The reflector's effect on the dominance ratio is small in comparison to the other factors, but was found to have a non-zero effect. Furthermore, the AGN-201M was found to have a significantly lower dominance ratio than systems in which it shares a very high (c{sub k} >95%) degree of neutronic similarity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

The stellar origins of 96 Zr excesses in presolar graphites from the Murchison meteorite

Context. Zirconium-96 is a stable isotope that can be synthesized under different neutron-rich nucleosynthetic conditions. Astrophysical models predict its production to occur in various stellar environments: from low-to-intermediate-mass asymptotic giant branch (AGB) stars to massive stars and core-collapse supernovae. Aims. Detections of 96Zr excesses, in combination with other isotopic measurements from presolar grains can provide unique constraints on its stellar origin. Presolar grains are microscopic particles found in primitive Solar System materials, which formed in stellar winds and supernova ejecta. The isotopic composition of each grain can provide us a snapshot of the nucleosynthetic processes that took place during the parent star’s lifetime. Methods. In this study, we measured the stable isotopes of C, N, O, Mo, Zr, and Ru in high-density presolar graphite grains from the Murchison meteorite and found four grains that contain positive isotopic anomalies in 96 Zr carried by their internal subgrains. We analyzed multi-element isotopic datasets from each grain to explore the source of the observed 96 Zr excesses. Results. Comparisons with stellar models indicate that two grains likely condensed in an intermediate-mass AGB star with initial metallicity of Z ≤ 0.014. Their 96 Zr/ 94 Zr ratios also match those predicted for born-again AGB stars undergoing a very late thermal pulse and rapidly accreting white dwarfs. After comparing the relative populations of the aforementioned dust-producing stars, we propose rapidly accreting white dwarfs as a new, and more likely, stellar source for one of the presolar grains. The remaining two grains could have originated in the supernova ejecta of massive stars, due to correlated excesses in the p-nuclides, 92, 94 Mo. Thus, grains with 96 Zr anomalies can have a variety of stellar origins, in agreement with theoretical studies. Conclusions. Our study highlights the importance of multi-element analysis in constraining the types of stars where presolar grains have condensed. These data will help improve our understanding of various nucleosynthesis processes in different stellar phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Thermomechanical Processing of Uranium Alloys with 10 and 50 Weight Percent Zirconium

This work focuses on casting and thermomechanical processing for moderate sized (500-1000g) 10% by weight and 50% by weight uranium-zirconium alloys. Three castings of uranium – 10wt% zirconium and a fourth casting of uranium – 50wt% zirconium were completed using different castings processes. The first casting process relied on only induction mixing within the VIM system thereby eliminating the need for an arc-melting feedstock preparation step. The second casting process used the arc-melting process prior for feedstock to the same VIM system. In this study, the induction mixing process alone with a graphite crucible appears insufficient for the fabrication of UZr alloys. Graphite is suitable with the use of master alloy UZr feedstock produced by arc-melting. There was a moderate increase of approximately 100 ppm of carbon within the U10Zr alloys. However, there was a near 1000 ppm increase of carbon in the U50Zr alloy. This carbon content is considered too high for forming processes, and therefore, warrants the use of a ceramic crucible for UZr castings near 50wt% zirconium. Zirconia does not appear suitable for UZr castings near 50wt% zirconium. The use of Beryllia or Thoria may be a more suitable option. However, due to the significantly lower liquidus and lower zirconium content, Zirconia crucibles do appear viable for carbon reduction in the U10Zr system. Two thermomechanical processes are also studied in this work - hot extrusion of the U50Zr system and heat treatment in an argon atmosphere. This work demonstrated that extrusion of the U50Zr alloy near 680°C is a viable choice even with +1000 ppm carbon. Extrusion of the U50Zr system near 630°C was unsuccessful and does not appear viable unless significant force and tooling strength is available. There appears to be a moderate increase in extrudate hardness when compared to as-cast. Heat-treatment reduced the hardness of as-cast alloys in the U10Zr castings as well as resolve some variance in the microstructures. However, for the U50Zr billet that did not successfully extrude at 630°C, heat treatment at 900°C for 96 hours did not have a significant impact on hardness or microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Experimental validation of a high fidelity Monte Carlo neutron transport model of the MIT graphite exponential pile

High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Proton branching ratios of 23 Mg levels

Background: The anomalous 22 Ne abundance measured in certain presolar graphite grains is thought to arise from the decay of 22 Na that was synthesized at high temperatures in core-collapse supernovae. To better interpret this abundance anomaly, the primary destruction mechanism of 22 Na, the 22 Na(p,γ) 23 Mg reaction, must be better understood. Purpose: Determine proton branching ratios of several 23 Mg excited states that play a role in the high-temperature 22 Na(p,γ) 23 Mg reaction rate. Methods: Particle decays of 23 Mg excited states populated with the previously reported 24 Mg(p,d) 23 Mg transfer reaction measurement [Kwag et al., Eur. Phys. J. A 56, 108 (2020)] were analyzed to extract proton branching ratios. The reaction was studied using a 31-MeV proton beam from the Holifield Radioactive Ion Beam Facility of Oak Ridge National Laboratory and 24 Mg solid targets. Results: Proton branching ratios of several 23 Mg excited states in the energy range of Ex = 8.044 - 9.642 MeV were experimentally determined for the first time for the p0 and p1'(p1+p2+p3) decay channels. Conclusions: These new branching ratios for 23 Mg levels can provide an experimental foundation for an improved high-temperature rate of the 22 Na(p,γ) 23 Mg reaction needed to understand production of anomalously high 22 Ne abundance in core-collapse supernovae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

ENDF/B-VIII.1: Thermal Neutron Scattering Sublibrary

The thermal neutron scattering law (TSL) sublibrary aims to describe the interaction of incident neutrons at thermal or sub-thermal energies with different compound materials such fuels, moderators and special-purpose materials. In ENDF/B-VIII.1 there was a large number of new and updated TSL evaluations, including traditional moderators (light water, Beryllium metal, Beryllium Oxide, Calcium Hydride, plastics (Polystyrene and Lucite), graphite (reactor-grade and crystalline), anhydrous Hydrogen Fluoride, and heavy paraffinic oil); exotic moderators (Beryllium Carbide, Zirconium Hydride, Yttrium Hydride, Lithium-7 Hydride and Deuteride), FLiBe molten salt, structural materials and cladding (Silicon Carbide, Silicon Dioxide, Zirconium Carbide), fuels (Plutonium Dioxide, Uranium Carbide, Uranium metal, Uranium Nitride, Uranium Dioxide, Uranium Hydride), and special purpose materials. In ENDF/B-VIII.1 we also distribute alongside the evaluated files, a comma-separated file (CSV), named TSL_MAT_numbers.csv, which lists all evaluated files in the current release and their corresponding unique MAT number.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Fission Product Quantification and Burnup Analysis via Gamma Spectrometry of the AGR-5/6/7 TRISO Fuel Experiment

Five capsules containing varying configurations of TRISO fuel compacts were irradiated in the most recent Advanced Gas Reactor (AGR) TRISO test campaign, AGR-5/6/7. As part of the post-irradiation examination (PIE) of this work all the fuel compacts and capsule graphite holders were measured with the Precision Gamma Scanner (PGS) to determine isotopic activities and locations using gamma spectrometry. Measuring the compacts allowed for the comparison of physics model depletion calculations with the measured compact activity to validate models. These measurements also are used to find compacts of interest with low fission product retention that may be tested in subsequent PIE. Measurements of the holders can confirm fission product distribution and migration. Of the isotopes measured, Cs-137 and Cs-134 were used to estimate location of failed particles as well as determine a measured burnup estimation. Results from Capsule 1 compacts and the holder indicate many fuel particle failures occurred. There are high amounts of cesium in the holder near Stack 9 of Capsule 1 that corresponds to low fractions?measured-to-calculated (M/C) ratio?of cesium remaining in those compacts. All other capsules had cesium retention ratios near 1. The isotope 110mAg had less consistent results. Capsule 2 compacts exhibited a M/C of 0.20, much lower than anticipated for the temperatures these compacts underwent. Capsule 3 experienced much higher temperatures, >1300C and had an average M/C of 0.75. These are the opposite results than what was predicted by past AGR experiments for the relationship between temperature and isotope retention. Capsule 4 and 5 showed good fission product retention and agreement with the physics model predictions. The burnup results from the PGS measurements showed strong agreement between the physics models and the measured values. Using a ratio of Cs-134/Cs-137 to calculate the burnup match closely with the physics models. Additionally, the Cs-137 burnup results were less accurate but extremely close to the model results.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Electrochemical stability and corrosion behavior of Ni-Cr alloy in molten LiCl-KCl salt

This study investigated the electrochemical properties and potential-dependent stability of Cr(II) and Ni(II) in a eutectic LiCl-KCl molten salt at 500 °C under an argon atmosphere using a three-electrode cell configuration with a glassy carbon working electrode, Ag/Ag + reference electrode, and graphite counter electrode. Cyclic voltammetry was employed to characterize the Cr(0)/Cr(II)/Cr(III) and Ni(0)/Ni(II) redox transitions, yielding diffusion coefficients of 0.50 (±0.01) × 10 −5 cm 2 s −1 for Cr(II) and 0.91 (±0.16) × 10 −5 cm 2 s −1 for Ni(II) at 500 °C. The redox stability domains of these species provided insight into the corrosion mechanisms of a Ni-Cr alloy (80–20 wt %) in LiCl-KCl salt. The measured open circuit potential of the alloy corresponded to the range of potentials where the Ni(0) and Cr(II) species are stable, indicating a preferential dissolution of Cr. Immersion testing of the alloy for 336 h confirmed significant Cr depletion (or Ni-enrichment) at the surface layer, corroborating Cr dealloying as the dominant corrosion reaction under open circuit conditions. Additionally, the application of a small anodic overpotential (0.11 V) resulted in selective Cr dissolution while a large anodic overpotential (0.25 V) resulted in co-dissolution of both Cr and Ni. Furthermore, these results demonstrate that the corrosion mechanisms of a Ni-Cr alloy in LiCl-KCl depend upon the applied potential, and thus the redox conditions of the molten salt should be carefully controlled in the design of a redox buffer to mitigate corrosion in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Interfacial Pressure Improves Calendar Aging of Lithium Metal Anodes

Lithium metal is a very attractive anode material because its theoretical specific capacity is approximately 10 times higher than conventional graphite anodes. Despite great promise, Li anodes suffer from capacity fade due to instabilities with the electrolyte as well as stranding of active Li. We have previously shown that applied interfacial pressure improves Li anode cycling because the pressure reduces the propensity for Li isolation and enables easier reconnection. Many researchers have also shown that calendar aging can lead to Li capacity loss and this has been attributed to either electrolyte decomposition with concurrent Li corrosion or to the formation of stranded Li. Our prior research focused on calendar aging during cycling suggests the mechanism for calendar aging is largely related to stranding of Li during rest and reconnection of the stranded Li upon further cycling, evidenced by similar average Coulombic efficiencies and Li loss in cells with and without rest. Because our calendar aging studies suggest Li stranding as a major cause of Coulombic efficiency drops and our Li cycling studies suggest this can be mitigated partially through applied interfacial pressure, we hypothesized that applied pressure would improve calendar aging by reducing stranded Li and enabling reconnection. We systematically varied applied pressure (0-1000 kPa) on Li metal anodes during cycling tests with and without intermittent calendar aging periods. Though the Coulombic efficiency decreases during aging periods, the lost capacity is recovered during subsequent cycles, as shown though average Coulombic efficiency and cumulative Li capacity loss analysis. We find that application of pressure partially mitigates calendar aging, in accordance with our hypothesis that calendar aging is caused by Li standing and can be mitigated to some degree with interfacial pressure. This is further supported by our results showing that the average Coulombic efficiency and cumulative Li capacity losses are similar over 50 cycles for cells that were continuously cycled and cells with periodic calendar aging periods. This result indicates that the losses during aging are reversible, which is consistent with Li stranding and reconnection. We show that pressure is one mitigation technique that helps reduce Li calendar aging in this study, but our finding that calendar aging is primarily governed by the stranding and reconnection of dead Li has wider implications. This research suggests that other mitigations which have been shown to prevent dead Li formation or encourage reconnection during cycling would also likely be successful for the purpose of improving calendar aging. The authors were supported by a Laboratory Directed Research and Development (LDRD) program. This work was performed, in part, at the Center for Integrated Nanotechnologies, an Office of Science User Facility operated for the U.S. Department of Energy (DOE) Office of Science. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC (NTESS), a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration (DOE/NNSA) under contract DE-NA0003525. This written work is authored by an employee of NTESS. The employee, not NTESS, owns the right, title and interest in and to the written work and is responsible for its contents. This work was authored in part by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Any subjective views or opinions that might be expressed in the written work do not necessarily represent the views of the U.S. Government.

applied pressure↗
DOE OSTI2021

Effects of 18 O Exchange on Neutron Emission Rates of Aging 238 PuO 2 LWRHUs

The Lightweight Radioisotope Heater Unit (LWRHU) provides about 1 Watt of thermal power to critical electrical and mechanical systems in the extreme temperatures of space, maintaining the operating temperature of the spacecraft’s scientific and general operating equipment. Los Alamos National Laboratory (LANL) designed and tested the LWRHU in the 1970s through early 1980s. The first production campaign at LANL was in the early 1980s in support of the Galileo mission in which over 100 LWRHUs were manufactured. In the mid-1990s, almost 200 LWRHUs were manufactured in support of the Cassini mission. After the primary Cassini manufacturing campaign, the Department of Energy (DOE) approved the fabrication of more than 10 additional LWRHUs as spares in 1998. The LWRHU assembly contains the following components: a fuel pellet, a vented capsule, a pyrolytic graphite insulator, and a fine-weave pierced fabric (FWPF) graphite aeroshell. The fuel pellet is a hot-pressed 238 PuO 2 cylindrical pellet that is sintered to create a ceramic pellet. The pellet is encapsulated in a platinum-rhodium encapsulation. The welded fueled clad is placed inside a pyrolytic graphite thermal insulation to protect the fuel from the heat produced in reentry events. The exterior layer of the LWRHU assembly is the FWPF aeroshell. Before sintering, PuO 2 granules are heated in a furnace at temperatures greater than or equal to 700°C under an enriched 16 O atmosphere in order to reduce trace quantities of 18 O. Because 18 O undergoes an (α,n) reaction, 18 O needs to be removed as much as possible to meet radiation emission requirements set by the Jet Propulsion Laboratory (JPL) and the National Aeronautics and Space Administration (NASA), as well as other agencies. The relative abundances of the oxygen isotopes are indicated by neutron emission rate (NER) measurements. LWRHUs are vented with a sintered platinum frit vent (George 1986). During extraterrestrial operation, the vent in the cladding allows helium from the decay of plutonium to escape and relieves pressure from the system while preventing the solids from escaping (Rinehart 1996, Tate 1982 & 1985). This vent is protected during manufacture by a platinum-30 rhodium cover. The capsule vent is activated by milling a 0.025-inch diameter hole through the protective cover to a nominal depth of 0.015 inches, generally immediately before the LWRHU is loaded into the aeroshell (George 1986). Prior to long-term storage, LWRHUs are vented and thus potentially allow infiltration of 18 O from the ambient atmosphere. Although NER decreases with age due to the reduction in emitted alpha particles, it might be expected that NER reduction is mitigated by the exchange of 18 O. However, this is not necessarily the case and no trend was observed relating NER to age other than the natural decay of 238 Pu (Mulford 2021). The most probable and dominant factor for determining specific NER was manufacturing discrepancies between individual units. It is important to understand the lack of observed 16 O- 18 O exchange and to compare LWRHU pellets as much as possible.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗