Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Fuel Performance Modeling”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Credible Criticality Safety Margin in the 30B Package with LEU+ UF 6 and Hypothetical Water Ingress

The commercial nuclear industry is pursuing advancements in fuel and reactor design that increase the uranium enrichment above 5 wt. % 235 U. These advancements will necessitate the ability to transport bulk quantities of UF 6 at increased enrichments. Currently, the 30B cylinder is the primary container used by the industry for UF 6 storage and transportation and has a long history of successful shipments. This container can support up to 2,277 kg of UF 6 at a maximum enrichment of 5 wt. % 235 U. Previous evaluations have assessed the potential impact of criticality safety for 30B transport at higher enrichments but assumed moderator intrusion would not require evaluation. Although current regulations allow for the exception of moderator intrusion for UF 6 packages through the design and quality control of the package content, this exception is limited to enrichments up to 5 wt. % 235 U. Thus, an investigation of moderator intrusion into a 30B cylinder should be performed. Moderator intrusion into a 30B cylinder is a unique condition for criticality safety evaluation in transportation because of the violent chemical reactions that occur between UF 6 and H 2 O. The resulting intrusion is strongly dependent on the breach size, breach location, breach interface solid/ullage of the UF 6 content, the UF 6 distribution (which is temperature dependent), and temperature/pressure conditions which are dynamic in accident conditions. Additionally, the complexity of the HF-UO 2 F 2 -H 2 O interface that occurs during the event can influence the potential solubility of uranium in the system, as well as influence the amount of UF 6 reaction with H 2 O. With the operating experience from Orano Federal Services LLC, the UF 6 chemical expertise from Oak Ridge National Laboratory (ORNL), and the criticality safety expertise from ORNL, this paper evaluates the neutronic conditions (i.e., k eff ) that apply the understood chemistry and experimental conditions that occur during moderator intrusion of a breached 30B cylinder under postulated accidents. This report examines the historical evaluations of UF 6 transport and expands these evaluations for the enrichments expected for nuclear industry advancement. These simulations primarily feature a homogeneous mixture of UF 6 and H 2 O as an infinite media system and as a sphere with water reflection given an impurity limit of 0.5 wt. % UF 6 . This report demonstrates that this H/U limit is valid up to 8 wt. % enrichment for both H 2 O and HF as the moderating mixture. The water-reflected homogeneous mixed spheres evaluated demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. Additionally, this report evaluates moderator ingress scenarios of a 30B cylinder. The 30B cylinder simulations vary the cylinder orientation, the mass of UF 6 in the system, the mass of H 2 O in the system, and how much H 2 O has reacted with UF 6 . The moderator ingress of a 30B cylinder was evaluated for the water-reflected homogeneous mixed spheres to demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. These simulations are standard practice and are independent of the UF 6 -H 2 O reaction. This work also explored more complex simulations that layer the H 2 O over the UF 6 , which incorporates the amount of reacted UF 6 as a separate layer and considers solubility limits of UO 2 F 2 in H 2 O as it is produced. This approach is intended to simulate the postulated event of a large hairline crack occurring on a submerged cylinder. Finally, this work considers the potential of H 2 O mixing into the UF 6 solid heterogeneously using a sponge-like model. The water ingress into the UF 6 is treated as random size spheres, and the reaction products form into a layer over the spheres.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of the Monte Carlo N-Particle radiation transport simulation tool, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) confirm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results assessment methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exemplified through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Criteria for Retention of 3013 S1 Containers Based on Relative Risk and Expert Judgment

An evaluation was performed to assess the suitability of thirty-three 3013 containers proposed for retention. These containers have moisture levels greater than 0.08 wt.% – the S1 population. The remainder of the S1 population stored at SRS will be down blended and disposed of by the end of 2028. Based on field surveillance and shelf-life data available to date as well as informed technical judgment, no container is currently expected to fail in its 50-year storage period. However, corrosion risk varies across the S1 population. Relative risks were evaluated using predicted Consensus Scores and their 95% Upper Prediction Limits (UPLs). The predicted values are based on a statistical model of Consensus Score as a function of moisture, chloride, and whether the packaged material was electrorefining scrap packaged at Hanford. Consensus Score has been shown to be a useful indicator of corrosion potential, and the UPL captures uncertainty in the model predictions, providing a conservative indicator of corrosion potential. Using UPLs to determine relative risks, together with expert review, three containers were identified as not suitable for retention, and the remainder were determined to be suitable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of the Monte Carlo N-Particle radiation transport simulation tool, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) con firm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results assessment methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exempli fied through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Design Status and Experimental Strategy for Initial Molten Salt Irradiation Experiments in HFIR

The development and deployment of molten salt reactor (MSR) technologies require experimental capabilities that can evaluate molten fuel salt behavior and structural material performance under representative irradiation conditions. Although modeling and separate effects testing provide important insight, there remains a critical lack of in-pile data that capture the coupled effects of neutron irradiation, temperature, salt chemistry, and time. Informed by lessons learned from historical MSR programs and recent international irradiation efforts, this report presents a structured approach to addressing existing gaps in molten salt irradiation testing using the High Flux Isotope Reactor (HFIR). A phased irradiation strategy is presented that incrementally increases experimental complexity while managing cost, risk, and facility constraints. The framework progresses from passive, nonfueled static capsule experiments to fueled, instrumented, and ultimately circulating molten salt systems, providing a pathway for materials screening, mechanistic understanding, and qualification-relevant testing. The report defines a near-term Phase 1 passive capsule concept, associated irradiation conditions, and a conceptual post-irradiation examination strategy. Ongoing modeling, out-of-pile testing, and safety analysis activities are described to support continued capability maturation and to enable future phases of molten salt irradiation testing in support of MSR technology deployment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Neutronic Consideration of TREAT Facility Fuel SiC Recladding

The Transient Reactor Test (TREAT) Facility is an air-cooled reactor composed of graphite-urania fuel encapsulated in zirconium alloy canisters and surrounded by ~2 ft (~61 cm) of nuclear grade graphite. The primary role of TREAT is to provide in-pile nuclear-heated safety research via transient testing [1]. Benchmark efforts have demonstrated sensitivity of TREAT neutronics models to graphite cross section data and key neutron poisons found in the fuel assemblies [2]. During fabrication of TREAT fuel assemblies, a small concentration of boron was inadvertently introduced with variable concentration into the aggregate of fuel blocks prepared for TREAT use. The concentration of boron, and its respective uncertainty, represent the most sensitive component in a neutronics study; unfortunately, nothing can be done to reduce or otherwise mitigate its impact. Hydrogen in the fuel blocks reduces the mean free path of the core and its true quantity and uncertainty is approximated. The hydrogen remains from pitch binder used in the fuel fabrication process and was not completely released due to incomplete fuel graphitization. Destructive fuel testing is necessary to quantify the actual hydrogen content and to improve benchmark model quality. A third neutron poison of some significance is the hafnium impurity in the Zircaloy-3 (Zry-3) assembly cladding. The quantity of hafnium is expected to be within nominal manufacturing limits for nuclear grade zircaloy. Recladding the fuel assemblies with material that is less neutron absorbing is expected to benefit reactor neutron economy. Moreover, TREAT’s total transient energy is restricted based on the rate of Zry-3 oxidation in TREAT’s air coolant environment. Hence, cladding with more benign air reaction kinetics could dramatically uprate the facility’s experimental capability. Silicon carbide (SiC) has been proposed as a candidate material for recladding TREAT assemblies, as SiC/SiC demonstrate favorable properties within nuclear environments [3]. The current operational limits of TREAT are restricted by the oxidation of Zry-3. Change of the cladding to SiC could feasibly double the current temperature limitation. Recladding opportunities allow for previously damaged fuel assemblies to be returned to service, as Zry-3 is not commercially available. The initial effort required to consider proceeding with any recladding effort is to neutronically characterize key reactor parameters while replacing Zry-3 with SiC. A characterization summary, with respective results, is provided herein.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

U.S. Efforts in Support of Examinations at Fukushima Daiichi- November 2020 Meeting Notes with Updated Information Requests

Much is still not known about the end-state of core materials in each unit that was operating on March 11, 2011 at the Fukushima Daiichi Nuclear Power Station (Daiichi). Information obtained from Daiichi is required to inform Decontamination and Decommissioning (D&D) activities, improving the ability of the Tokyo Electric Power Company Holdings, Incorporated (TEPCO Holdings) to characterize potential hazards and to ensure the safety of workers involved with cleanup activities. This information also has important implications for the safety and operation of U.S. commercial nuclear power plants. This document summarizes results from the Fiscal Year 2021 (FY2021) U.S. effort to review Daiichi information and extract insights to enhance the safety of existing and future nuclear power plant designs. This U.S. effort, which was initiated in 2014 by the Department of Energy Office of Nuclear Energy (DOE-NE), is completed by a group of experts in reactor safety and plant operations that identify examination needs and evaluate recent Daiichi examination data to address these needs. Since its inception, annual reports were issued that document significant safety insights being obtained in areas of special emphasis: system and component performance, radionuclide surveys and sampling, debris end-state location, combustible gas effects, and plant operations and maintenance. In addition to reducing uncertainties related to severe accident modeling progression, these insights have and continue to be used to update guidance for severe accident prevention, mitigation, and emergency planning. Reduced uncertainties in modeling the events at Daiichi improve the realism of reactor safety evaluations that inform future D&D activities. A key aspect of prior U.S. efforts, the updated list of information requests, is included in this FY2021 report to ensure that they are transmitted to organizations within Japan. This report also continues to emphasize how information obtained from the affected reactors at Daiichi has been and will continue to be used to update severe accident management strategies and reduce uncertainties in systems analysis code models. In addition, recommendations are included that would expand the use of this information to provide insights regarding maintenance, radiation protection, design, and siting activities for existing and new reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

AGR-5/6/7 Irradiation As Run Predictions Using PARFUME

The PARticle FUel ModEl (PARFUME), a fuel performance modeling code used for high-temperature gas-cooled reactors, was used to model the Advanced Gas Reactor (AGR)-5/6/7 irradiation test using as-run physics and thermal data. The AGR-5/6/7 irradiation test consists of the combined fifth, sixth, and seventh planned irradiations of the AGR Fuel Development and Qualification Program. The AGR-5/6/7 test train is a multi-capsule, instrumented experiment that is designed for irradiation in the 133.4-mm diameter northeast flux trap position of the Advanced Test Reactor (ATR) at Idaho National Laboratory. Each capsule contains compacts filled with uranium oxycarbide unaltered fuel particles. This report documents the calculations performed to predict the failure probability of tristructural isotropic (TRISO)-coated fuel particles during the AGR-5/6/7 experiment. In addition, this report documents the calculated fission product release fraction from the fuel. The calculations include modeling of the AGR 5/6/7 irradiation that occurred from February 2018 to July 2020 over nine ATR cycles, including six normal cycles and three power axial locator mechanism cycles, for a total of approximately 376 effective full power days (EFPD). The irradiation conditions and material properties of the AGR-5/6/7 test predicted zero fuel particle failures in Capsules 1, 3, and 4. Fuel particle failures were predicted in two of the compacts in Capsule 2 and one particle failure is predicted in each one of the compacts in Capsule 5. All compacts that exhibited fuel particle failures predicted by PARFUME were caused by localized stress concentrations in the silicon carbide (SiC) layer caused by cracking in the inner pyrolytic (IPyC) layer. In addition, shrinkage of the buffer and IPyC layer during irradiation resulted in formation of a buffer-IPyC gap. Compacts with a lower irradiation temperature and fluence experienced the smallest buffer-IPyC gap formation. Conversely, higher irradiated temperature compacts with a high fluence experienced the largest buffer-IPyC gap formation. Compact 3-6-3 experienced the largest buffer IPyC gap formation of just under 21.7 µm. The release fraction of fission products silver (Ag), cesium (Cs), and strontium (Sr) vary depending on capsule location and irradiation temperature. The maximum release fraction of Ag occurs in Capsule 3, reaching up to 59.5% for the TRISO fuel particles (compact 3-6-3). The release fraction of the other two fission products, Cs and Sr, are much smaller. A maximum Cs release fraction of 1.1% occurred in compact 3-4-3 and 4.4% for Sr in compact 3-6-3.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The impact of anisotropic thermal expansion on the isothermal annealing of polycrystalline α -uranium

Although grain growth impacts microstructural evolution in a wide variety of materials systems, the effect of anisotropic thermal expansion on grain boundary mobility and texture evolution has not been widely studied. Anisotropic thermal expansion occurs in multiple non-cubic metals, and the thermomechanical processing behavior of these materials can be better understood with further study into the impact of thermal expansion on grain boundary mobility and texture evolution. In this work, we develop a mesoscale phase field model of grain growth that includes the effect of anisotropic thermal expansion, which is applied to study polycrystalline α-uranium, a highly anisotropic metal. Three-dimensional simulations on polycrystalline a-uranium with and without thermal expansion eigenstrains are performed to study the grain boundary mobility and texture evolution as a function of temperature. A strain-free temperature of 933 K is selected, and the system is studied at within the range of 873 – 933 K at intervals of ten degrees, resulting in increasing thermal eigenstrain with decreasing temperature. We also estimate a grain boundary mobility prefactor and activation energy based on existing experimental data of isothermal annealing of α-uranium. The grain boundary mobility is found to display significant deviation from Arrhenius behavior with the inclusion of thermal expansion eigenstrain as the amount of thermal eigenstrain (and thus elastic strain energy within the system) increases. This result explains an experimentally observed grain boundary mobility deviation from Arrhenius behavior. Furthermore, the texture evolution is affected, such that the grain orientations become less random with increasing thermal eigenstrain, which could explain experimentally observed texture behavior. These results indicate that the effect of thermal expansion should be considered when predicting the thermomechanical processing behavior of α-uranium and other materials with anisotropic thermal expansion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Techno-economic analysis of advanced small modular nuclear reactors

Here, small modular nuclear reactors (SMRs) represent a robust opportunity to develop low-carbon and reliable power with the potential to meet cost parity with conventional power systems. This study presents a detailed, bottom-up economic evaluation of a 12 × 77 MW e (924 MW e total) light-water SMR (LW-SMR) plant, a 4 × 262 MW e (1,048 MW e ) gas-cooled SMR (GC-SMR) plant, and a 5 × 200 MW e (1,000 MW e total) molten salt SMR (MS-SMR) plant. Cost estimates are derived from equipment costs, labor hours, material inputs, and process-engineering models. The advanced SMRs are compared to natural gas combined cycle plants with and without post-combustion carbon capture and a conventional large nuclear reactor. Overnight capital cost (OCC) and levelized cost of energy (LCOE) estimates are developed. The OCC of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$4,844$/kW, $\$4,355$/kW, and $\$3,985$/kW respectively. The LCOE of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$89.6$/MWh, $\$81.5$/MWh, and $\$80.6$/MWh respectively. A Monte Carlo analysis is performed, for which the OCC and construction time of the LW-SMR is found to have a lower mean and standard deviation than a conventional large reactor. The LW-SMR OCC is found to have a mean of $\$5,233$/kW with a standard deviation of $\$658$/kW and a 90 % probability of remaining between $\$4,254$/kW and $\$6,399$/kW, while the construction duration is found to have a mean of 4.5 years with a standard deviation of 0.8 years and a 90 % probability of remaining between 3.4 and 6.0 years. The economic impact of economies of scale, simplification, modularization, and construction time for SMRs are discussed. Additionally, policy implications for direct SMR capital subsidies and the impact of a carbon tax on natural gas emissions are explored.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Towards an Agent-Based Blackboard System for Reactor Design Optimization

The field of reactor design is rich with opportunities for applications of computational optimization algorithms; these applications can range from preliminary core design to reactor shuffling patterns. Many of these schemes rely on sets of previously generated solutions (sometimes referred to as “generations”) to inform future decisions. While it is important to build upon prior knowledge, this process requires a full generation of solutions to be formed before future solutions can be examined. Rather than relying on a generational scheme to perform an optimization, we propose using an agent-based approach in conjunction with a blackboard framework for performing reactor design optimizations. Utilizing an agent-based approach allows agents to perform tasks independently, while retaining the ability to build off of previous solutions. We develop an agent-based blackboard system (ABBS) for determining the Pareto front (PF) in sodium fast reactor design optimization problems and compared this with the Non-Dominated Sorting Genetic Algorithm II (NSGA-II). Our goal is to evaluate the viability of the ABBS in producing a PF that is comparable with the NSGA-II algorithm. The design space consists of the fuel height, fuel smear, and plutonium fraction in the core, and we seek to minimize the reactivity swing and plutonium mass, while maximizing the burnup. The diversity, coverage, and spread of the PFs generated by the two methods are examined, and the ABBS is able to converge to the same PF as the NSGA-II algorithm. These results show that the ABBS is able to find optimal designs that are similar to those found by the NSGA-II algorithm. We conclude our study by applying the ABBS to the design of a sodium-cooled fast reactor to dispose of weapons-grade plutonium. The ABBS finds a core design that can burn upwards of 17.5 kg of weapons-grade plutonium per year and degrade an additional 195 kg of weapons-grade plutonium per year into non-weapons-grade material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Engineering defect energy landscape of CoCrFeNi high-entropy alloys by the introduction of additional dopants

The concept of high-entropy alloys (HEAs) focusing on tuning the overall chemical complexity represents a novel alloy design strategy. In contrast, alloying of a metallic matrix with minor doping elements with limited and localized tunability has been a common practice to improve material performance. Combining the idea of globally engineering defect energy landscape in HEAs and the localized doping strategy in dilute alloys, in this work, we explore doping effects of minor elements in a HEA matrix to further enhance the overall and localized chemical tunability, aiming to improve its irradiation resistance. Specifically, we study the influence of minor Al, Cu, Ti, and Pd substitutional doping elements on defect energetics in a CoCrFeNi model HEA based on density-functional theory (DFT) calculations. The DFT results indicate that the formation and migration energies of vacancies can be strongly influenced when a dopant is introduced at the first nearest neighbor shells around a vacancy. On the other hand, interstitial energetics are only slightly affected. Among the four elements, Ti and Pd generally decrease vacancy formation energies and increase vacancy migration energies more significantly than Al and Cu. The doping effects become more pronounced when the concentration of the substitutional dopants increases. Based on the energy distributions obtained from DFT, we build a kinetic Monte Carlo (kMC) model to assess the impact of dopants on vacancy-mediated diffusivity in the doped HEAs. Our results suggest that Ti and Pd can lower the tracer diffusivity in the considered HEAs and act as trapping sites, whereas Cu may enhance the atomic transport. This work indicates that substitutional doping in HEAs is an effective strategy in metallurgy to further tune the defect and transport properties of complex alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing irradiation performance of the Tristructural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in uncertainty of calculated fuel temperature is the gas gap uncertainty due to the nub-to-shell clearance because of fabrication error. A revision to the thermal model was made to examine the most probable offset position for four different dates of interest. The offset position options varied in magnitude and azimuthal direction of offset of the holder top and bottom position. The best-fit offset of the holder is estimated based on the minimum root mean square error of residuals for all operational thermocouples (TCs). From these results, the following conclusions were made: 1) During earlier cycles (162A – 164A) when numerous TCs were still operational, the best-fit offset distance varied in range [0.002 – 0.0035 in.] for both the top and bottom holder while offset azimuthal direction varied widely, especially for the holder bottom. 2) Holder offset led to slightly lower Capsule 1 average temperature, but wider temperature variation (lower minimum and higher peak fuel temperatures). 3) The maximum offset of 0.006 in. to the northwest direction for both the top and bottom holder yielded a much higher temperature than the original peak fuel temperature, which increased from 1422 to 1557 ? (or 135 ? increase). This could be one of the reasons for massive particle failure near the end of Cycle 165A. 4) Even though, the highest temperature at the TC-1-7 tip reached slightly higher than 1000?, the temperature along the TC wire reached as high as 1335? assuming a holder offset of 0.006 in. to the northwest. Consequently, at this temperature nickel from this TC wire could cause particle failures in Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Technical Evaluation of Accelerated Basin De-Inventory Material Addition to Sludge Batch 11 (Rev. 1)

The Accelerated Basin De-inventory (ABD) program involves discarding spent nuclear fuel that is currently stored in L-Basin to the Defense Waste Processing Facility (DWPF) for vitrification. The first ABD discards will occur during the preparation of Sludge Batch (SB) 11. Savannah River Mission Completion has requested that the Savannah River National Laboratory assess the technical gaps related to the increased gadolinium poisoning requirement and the impacts of performing the Low Temperature Aluminum Dissolution (LTAD) process in Tank 51 with H-Canyon discards present. The following summarizes the evaluation of the impacts of increasing the quantity of gadolinium (and related topics) from what was previously evaluated in the SRNL studies of gadolinium-poisoned ABD material solubility, the overall ABD flowsheet review, and increasing the fissile mass loading in glass: 1) Based on literature surveys, there is no indication that organic interactions with gadolinium will be significant at the high pH (typically >13) conditions of the Concentration, Storage, and Transfer Facilities. Any interactions of gadolinium with organics in DWPF are not expected to adversely impact DWPF or downstream facilities. Thus, there is little-to-no residual risk from organic interactions with gadolinium [Gap closed]; 2) Adding depleted uranium to ABD material, targeting 235 U enrichment of 4.90% within each transfer window, will mitigate potential impacts from an increase in soluble 235 U enrichment during sludge washing and LTAD. The plan to take advantage of previous transfers and allow 235 U enrichment of >5% during the final transfer window carries a risk that Tank 51 supernate will have a 235 U enrichment of >5%, which should be evaluated for acceptance; 3) Increasing the gadolinium mass ratio to 3.0:1 Gd: 235 U(eq SLU ) should lead to the same or higher partitioning of gadolinium into the solid phase within the DWPF Chemical Process Cell, resulting in both liquid and solid phases with expected partitioning of Gd consistent with the prior solubility study [Gap closed for SB11]; 4) There are no expected impacts on DWPF melt temperature and melter operations due to the minimal ~0.2 weight percent (wt%) increase in Gd concentration relative to previous sludge batches [Gap closed for SB11]; 5) As observed previously, Gd is expected to enter the off-gas system via physical entrainment, but at a slightly higher concentration than what was observed for SB9 melter off-gas pluggage deposits (0.07 wt%) [Gap closed for SB11] ; 6) There are no expected impacts on DWPF recycle or the Recycle Collection Tank glycolate destruction process. [Gap closed for SB11]; 7) Gd is projected to be a trace component in the SB11 glass (<0.5 wt%) and can be ignored for process control. Trace components do not significantly impact glass durability, thus the conclusions of the previous Product Consistency Test evaluation at a fissile mass loading of 2,500 g fissile/m3 glass still applies to SB11. The ~0.1 wt% increase in Gd2O3 concentration relative to the previous study will not impact the predictability of SB11 glass with the DWPF Product Composition Control System (PCCS) models for durability or the acceptability of glass according to the Waste Acceptance Product Specifications (WAPS) criterion for product consistency [Gap closed for SB11]; 8) No additional Toxicity Characteristic Leaching Procedure testing is necessary for SB11 and the hazardous waste specification of the SB11 DWPF waste form is unchanged after the addition of the ABD stream [Gap closed for SB11]. The following summarizes the evaluation of the impacts of adding two-thirds of the ABD material to Tank 51 prior to LTAD: 1) The addition of two-thirds of the ABD increases overall aluminum mass from 1.39×10 4 kg to 1.64×10 4 kg (15.5% ABD Al). The form of the insoluble portion of the Al resulting from ABD addition should be the more readily dissolved Al(OH) 3 and amorphous forms. The portion of the ABD aluminum that is processed by LTAD is expected to be completely soluble, thus requiring that less of the boehmite in the sludge be dissolved to reach the same Al target in the SB. [Gap closed for SB11]; The expected LTAD impact on other components, as related primarily to the components in ABD, are discussed. Gd is expected to remain insoluble during LTAD and not impact the solubility of other components. [Gap closed for SB11]; The addition of two-thirds of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 13,100 kg (63% ABD U) and the projected plutonium mass from 86.0 kg to 89.5 kg (3.9% ABD Pu). The addition of all of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 16,100 kg (70% ABD U) and the projected plutonium mass from 86.0 kg to 90.4 kg (5.3% ABD Pu). The 235 U enrichment will be ≤5%. The fissile uranium will be adequately poisoned by Gd and the fissile Pu will be adequately poisoned by Fe from the sludge. [Gap closed for SB11]; There is a low risk that ABD addition will impact the rheology or pumpability of the slurry. There is a low but higher risk of ABD addition prior to LTAD impacting the settling rate; Based on the evaluation of adding two-thirds of the ABD material and all of the ABD material prior to the LTAD process, there is no volume or mass limit that would need to be imposed on ABD additions prior to LTAD. [Gap closed for SB11]. Revision 1 of this report addresses a variation on the ABD additions and LTAD strategy where sodium hydroxide additions for LTAD may be performed intermittently or concurrently with an ABD addition window. The proposed change does not alter the conclusions of this evaluation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Equilibrium Core Model for Micro Pebble Bed Reactors Using OpenMC

Estimating the equilibrium state for pebble bed reactors (PBRs) presents complex challenges as it requires simultaneous consideration of changes in the pebbles’ movement as well as their fuel compositions. Whereas traditional approaches use multigroup diffusion codes for neutronics calculations of PBRs’ equilibrium state, the double-heterogeneity of PBRs complicates neutron cross-section generation. Continuous-energy Monte Carlo (MC) methods are better suited for detailed PBR analysis because of their natural handling of double-heterogeneity, but they demand substantially more computational resources. Here, this study introduces a novel method for efficiently estimating the equilibrium state in small and micro PBRs with reduced computational cost. The method is anticipated to accelerate the processes of core design and performing parametric studies for utilizing advanced fuel and structural materials. The HTR-10 reactor design was used for validating the method’s predictions and evaluating its computational efficiency. When compared to reference calculation values from the literature, criticality (k-effective) was predicted to be approximately within the margin of error of the MC transport calculation, average core power density (in megawatts per cubic meter) was predicted within 2.5% relative error, and maximum thermal flux (10 13 n/cm 2 .s −1 ) was predicted within 1.8% relative error. The calculated inventory of fission products and fuel composition in the equilibrium core were within 15% and 16.6%, respectively, when compared to reported values from the literature. The difference is attributed to variance in the considered values of the core temperature, which was found to significantly affect the depletion analyses.

Equilibrium core↗