Results for “fuel performance modeling” — Page 22

DOE OSTI2025

The SAS4A/SASSYS-1 Version 5.8 Safety Analysis Code System

SAS4A/SASSYS-1 is a software simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced nuclear reactors. Detailed, mechanistic models of steady-state and transient thermal, hydraulic, kinetic, and mechanical phenomena are employed to describe the response of the reactor core, the reactor primary and secondary coolant loops, the reactor control and protection systems, and the balance-of-plant to accidents caused by changes in coolant flow, loss of heat rejection, or reactivity insertion. The consequences of single and double-fault accidents can be modeled, including fuel and coolant heating, fuel and cladding mechanical behavior, core reactivity feedbacks, coolant loop performance including natural circulation, and decay heat removal. Analyses are typically terminated upon demonstration of reactor and plant shutdown to permanently coolable conditions, or upon violation of design basis margins. The objective of the analysis is to quantify accident consequences as measured by the transient behavior of system performance parameters, such as fuel and cladding temperatures, reactivity, and cladding strain. Originally developed for analysis of sodium cooled reactors with oxide fuel clad by stainless steel, the models were subsequently extended and specialized to metallic fuel clad with advanced alloys and to several other coolant options, including lead, LBE, and water.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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

Pressurized Water Reactor Control Rod Ejection Analysis Using PARCS, RELAP5-3D, and BISON for High Burnup Fuel

This paper describes safety assessment of high burnup fuel (HBF) during loss of coolant accidents (LOCA), reactivity-initiated accidents (RIA), and anticipated operational occurrences (AOOs) for Southern Nuclear’s Vogtle Electric Generating plant (VEGP). A goal of this project is to create a multiphysics modeling approach in aspects of thermal hydraulics, neutronics, and fuel performance, as well as to use it in combination with experimental studies (Figure 1). In support of the project, this paper seeks to address fuel performance of a HBF pin during steady state conditions and during a control rod ejection (CRE) scenario using a loose coupling of PARCS, RELAP5-3D, and BISON.. The primary goal of this work was to analyze high burnup fuel behavior during a CRE scenario. The results indicate fuel failure is not expected to occur under the conditions presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Understanding Carbon Dioxide Transfer in Direct Methanol Fuel Cells Using a Pore-Scale Model

Abstract The gas flow of carbon dioxide from the catalyst layer (CL) through the microporous layer (MPL) and gas diffusion layer (GDL) has great impacts on the water and fuel management in direct methanol fuel cells (DMFCs). This work has developed a liquid–vapor two-phase model considering the counter flow of carbon dioxide gas, methanol, and water liquid solution in porous electrodes of DMFC. The model simulation includes the capillary pressure as well as the pressure drop due to flow resistance through the fuel cell components. The pressure drop of carbon dioxide flow is found to be about two to three orders of magnitude higher than the pressure drop of the liquid flow. The big difference between liquid and gas pressure drops can be explained by two reasons: volume flowrate of gas is three orders of magnitude higher than that of liquid; only a small fraction of pores (<5%) in hydrophilic fuel cell components are available for gas flow. Model results indicate that the gas pressure and the mass transfer resistance of liquid and gas are more sensitive to the pore size distribution than the thickness of porous components. To buildup high gas pressure and high mass transfer resistance of liquid, the MPL and CL should avoid micro-cracks during manufacture. Distributions of pore size and wettability of the GDL and MPL have been designed to reduce the methanol crossover and improve fuel efficiency. The model results provide design guidance to obtain superior DMFC performance using highly concentrated methanol solutions or even pure methanol.

Electrochemistry↗

DOE OSTI2025

Systems Analysis of Biomass and Coal Co-firing Power Plants with Deep Carbon Capture Toward Net-zero Emissions

Achieving a net-zero emission economy in the United States requires integrating diverse low-carbon and negative-emission technologies into the existing fossil fuel-dominant power fleet. Potential technologies from the low-carbon portfolio include renewable power, fossil power with carbon capture and storage (CCS), bioenergy with CCS (BECCS), and direct air capture (DAC). Renewable power is a clean energy source but has to pair with costly battery storage to provide dispatchable electricity. Fossil power with CCS offers dispatchable electricity yet still relies on DAC to offset residual emissions, even when deploying deep CCS with more than 90% CO2 capture. Coal-biomass co-firing with CCS, a subset of BECCS, is a reliable energy production technology that can be retrofitted from existing electricity generation units (EGUs). Power plant retrofit maximizes the use of the current U.S. coal power fleet without the need for large-scale deployment of new renewable power, battery storage, or DAC. Retrofitting coal-biomass co-firing with deep CCS in EGUs is a promising option, but not a universal solution. Biomass co-firing at a power plant introduces economic challenges and indirectly poses pressure on land and water resources. Meanwhile, retrofitting deep CCS affects plant efficiency and raises electricity generation costs. Overall, the technical feasibility and economic viability of plant retrofits vary across EGUs, as they are contingent upon the regional availability of biomass, unit-specific characteristics, site-specific fuel supply costs, and adjacent CO2 storage potential. Government incentives like 45Q can improve the retrofit viability, though the impact requires further quantification. A comprehensive analysis at the unit level is essential to address the question regarding the fate of the U.S. coal-fired electricity generation fleet toward the net-zero emission goal. This study conducts a systematic techno-economic-environmental assessment of EGUs to identify the viability of biomass co-firing and deep CCS retrofits in the U.S. coal-fired power fleet. Specifically, it characterizes the techno-economic performance of deep carbon capture, estimates life cycle greenhouse gas (GHG) emissions, and conducts a fleet-level assessment on retrofit viability. The key objectives are (1) to estimate the unit-specific performance and retrofitted cost under various biomass co-firing levels and CO2 capture rates; (2) to determine the possibility of reaching net-zero emission at the fleet level; (3) to quantify the cumulative capacities that are suitable for plant retrofits under current and future biomass supply scenarios; and (4) to improve the understanding of policy impacts on such retrofits to help the power sector’s transition to a net-zero economy. Techno-economic Model of Deep Carbon Capture. This study develops the performance and economic models for Monoethanolamine-based post-combustion CO2 capture at 95–99% capture rates. The process is simulated in Aspen Plus, analyzing the performance of carbon capture technology by varying the plant sizes, solvent lean loading, CO2 concentrations, and flue gas inlet temperature. Based on the key inputs and output parameters of CO2 capture, a reduced-order performance model of deep carbon capture is formulated. In addition, an engineering-economic model integrating the performance metrics is developed to estimate the capital as well as operation and maintenance (O&M) costs. Capital cost estimations follow the framework of the Integrated Environmental Control Model (IECM) and incorporate data regressions from three technical reports by IECM, the National Energy Technology Laboratory (NETL), and the National Renewable Energy Laboratory. The O&M cost estimation utilizes the actual inventory consumption rate and labor requirements. Both performance and cost models are embedded into IECM v13.0-beta, a fossil-fuel power plant modeling tool. Life Cycle Assessment of Power Plants. This study estimates the GHG emissions of power plants through life cycle assessment (LCA). The LCA scope includes fuel supply, combustion-based power generation, and CO2 transport and storage. The fuel-based life cycle module is designed following the framework of the NETL Unit Process Library and CO2U LCA Guidance Toolkit. The module is then incorporated into IECM v13.0-beta. The process-based LCA is applied to estimate the GHG emissions of coal and biomass supply, coal- and coal-biomass co-firing power plant operation, as well as CO2 pipeline transport and geographical sequestration. An uncertainty analysis is conducted to quantify the variability and uncertainty associated with the LCA using the Latin Hypercube Sampling (LHS) method. Fleet-level Assessment. This study evaluates the technical and economic feasibility of selected coal-fired EGUs, examines the role of tax credits in retrofit viability, and assesses the competitiveness of retrofitted units against other low-carbon options. Unit screening identifies EGUs for the study, focusing on new, efficient baseload units with air pollution controls. The power plant databases are then established to organize unit-specific information on performance and operating conditions from the relevant public databases. Biomass for co-firing retrofits is selected based on home and neighboring county availability, ensuring sustained operation with at least a 5% co-firing level. The CO2 storage site is determined by state-level storage potential, with ArcGIS Pro and NETL CO2 Saline Storage Cost Model used to identify the optimal balance between the nearest transport distances and affordable storage costs. The latest IECM v13.0-beta is then employed to configure and evaluate the eligible EGUs with or without the deployment of deep CCS and biomass co-firing. A supply curve is established to illustrate the cumulative installed capacity suitable for retrofits at different cost levels. A sensitivity analysis on tax credits for carbon sequestration is performed. Finally, a unit-level cost comparison is conducted among retrofitted plants, renewable power with battery storage, and abated fossil fuels with DAC. Expected Results. This study evaluates the technical, economic, and environmental metrics of each EGU across an array of CO2 capture rates and biomass co-firing level scenarios. Unit-level comparisons will identify critical factors influencing technical performance. The supply curves with and without tax incentives will provide insights into the impact of tax credits on biomass co-firing and CCS deployment. The cost comparisons with renewables and DAC-retrofit will assess the competitiveness of the retrofitted units. Life cycle emissions from each unit will be assessed to identify the scenarios under which net-zero emissions can be achieved. These analyses are expected to determine the total coal-fired capacity suitable for serving as a low-carbon energy source with or without tax incentives. The study results are novel in identifying optimal unit-specific strategies for producing carbon-neutral power, whether through retrofitting EGUs with deep CCS, biomass co-firing, DAC, or installing renewable power with battery. The findings will provide insight into nationwide efforts to ensure reliable, affordable, and low-carbon electricity. It also will inform investment decisions and policies in the deployment of deep carbon capture and negative emission technologies for a net-zero energy future.

Biomass Co-firing↗

DOE OSTI2021

BISON As-run AGR-3/4 Irradiation Test Predictions

BISON, a nuclear fuel performance application built using the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library, was used to model the Advanced Gas Reactor (AGR)-3/4 irradiation test using as-run physics and thermal hydraulics data. The AGR-3/4 test consists of the combined third and fourth planned irradiations of the AGR Fuel Development and Qualification Program. The AGR-3/4 test train consists of twelve separate and independently controlled and monitored capsules. Each capsule contains four compacts filled with both uranium oxycarbide (UCO) unaltered “driver” fuel particles and UCO designed-to-fail (DTF) fuel particles. The DTF fraction was specified to be 1×10-2. This report documents the calculations performed to predict the failure probability of tristructural isotropic (TRISO) coated driver fuel particles during the AGR-3/4 experiment on a single compact. This report will demonstrate the capabilities of BISON to model the complex AGR-3/4 irradiation test and identify further development needed to capture the fuel particle failure probability and source term on every compact for further comparison from post-irradiation examination (PIE) data. The calculations include the modeling of the AGR-3/4 irradiation that occurred from December 2011 to April 2014 in the Advanced Test Reactor (ATR) over a total of 10 ATR cycles including seven normal cycles, one low power cycle, one unplanned outage cycle, and one Power Axial Locator Mechanism (PALM) cycle.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Implications of point defect accumulation on UO 2 thermal conductivity and fission gas release under accelerated fuel irradiation

Evaluation of thermal properties is a crucial factor for nuclear fuel performance. During reactor operation, the accumulation of fission products and irradiation-induced lattice defects are responsible for degradation in thermal conductivity. Consequently, it affects fuel temperature and fission gas release (FGR) among other Multiphysics processes important for economics and safety analysis. We analyze the implications of point defects (PD) accumulation described using a rate theory (RT) Model on lattice thermal conductivity of UO 2 . Here, we demonstrate that fission rate-dependent point defect concentrations have the largest impact on in-pile thermal conductivity in the periphery of light water reactor fuels below a temperature threshold governed by the migration barrier of defects. Our analysis provides a mechanistic description of this phenomena which current fuel performance codes treat empirically. The reduction of thermal conductivity in the low -temperature rim region acts as additional thermal resistance and leads to a temperature notably larger than suggested by Lucuta thermal conductivity correlation. These effects are anticipated to have notable impacts when fuels are exposed to accelerated radiation. The impact of such point defect-informed treatment of thermal conductivity on fuel performance is evaluated by a detailed analysis of fission gas behavior and its release. We consider several models capturing different stages of fission gas bubble evolution and fission gas release (FGR). Finally, a new fission rate-dependent correction to the Lucuta correlation is proposed. The results show a significant reduction in thermal conductivity at the fuels’ periphery and an increase in fuel centerline temperature specifically at low burnups. Ultimately a modified LC shows a higher FGR compared to the original LC, while the acceleration process results in a reduction in overall FGR.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Implementation and testing of physics-based pulverization model in BISON

This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.

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 the irradiation performance of tri-structural 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 the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

advanced gas reactors↗

DOE OSTI2023

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

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural 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 the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

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

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural 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 the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach. Vol. I: PWR Fuel

Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Mortar-based Cohesive Zone Model with Application to TRISO Particle Debonding

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report outlines work towards the development of an integrated, mortar-based cohesive zone framework with mechanical contact. We model debonding of tri-structural isotropic (TRISO) layers using a mortar approach, extending the framework used for mechanical contact, using weighted quantities. A bilinear mixed mode traction approach, already used in Jiang et al. (2021), is applied analogously herein to model debonding between TRISO particle layers. The examples include the debonding of the inner pyrolytic carbon (IPyC) from the silicon carbide (SiC) layers, as well as partial and full debonding of the buffer from the IPyC layers. Both 2D and 3D examples are utilized to demonstrate the superior numerical performance of the new mortar approach, which allows for modeling arbitrarily meshed interfaces (i.e., non-matching discretizations).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Evaluation of Mechanistic and Empirical Models against Existing FFRD and LOCA Experimental Databases

The desire of the nuclear industry to improve the economics of existing nuclear power plants has necessitated research into the potential of a phenomenon known as fuel fragmentation, relocation, and dispersal (FFRD). This phenomenon is possible during a loss-of-coolant accident (LOCA) at relatively high burnup. The Nuclear Energy Advanced Modeling and Simulation program has been developing simulation capabilities for FFRD and LOCA in the BISON fuel performance code for multiple years. This year, the effort has been focused on evaluating new lower length scale informed pulverization thresholds as well as updating and adding new empirical models for various phenomena. Models added or updated this year include a preliminary transient fission gas release model, new high-temperature Zircaloy creep models, a Zircaloy rupture opening area model, and a temperature-dependent emissivity during radiation from the fuel rod to the surrounding atmosphere. The models are verified through implementation tests to demonstrate code correctness after addition to BISON. The models are then assessed against a subset of the existing BISON validation suite for LOCA and FFRD cases. Two new cases, Studsvik Rods 192 and 193, have been added. The results indicate that the inclusion of a bubble pressure evolution model in the bubbles in the high-burnup structure has the largest impact compared to the 3D fracture criterion on reducing the calculated amount of pulverized fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Pellet cladding mechanical interaction as a potential failure mechanism during a control rod drop accident in a boiling water reactor

Boiling water reactors (BWRs) represent approximately one-third of the operating fleet in the United States, contributing significantly towards the global effort in reducing carbon emissions. Even though most of the operating fleet has been in operation for quite some time, continued advancements in new nuclear fuel (such as accident tolerant fuel) or operating regimes (such as power up-rates and higher burnup operation) necessitates similar advancements in modeling and simulation capabilities. Bison, a component of the Virtual Environment for Reactor Applications (VERA), is a high-fidelity fuel performance code able to explore the fuel performance of a wide variety of fuel types in one-, two-, and three-dimensions. Until recently, the code had not been used for analyses of BWRs. Modeling capabilities have been added for Gd-bearing UO{sub 2} and pure zirconium liners. New models have been added based upon the U.S. Nuclear Regulatory Commission (NRC) guide-lines for hydrogen pickup in Zircaloy-2 claddings and pellet-clad mechanical interaction (PCMI) failure during a reactivity insertion accident (RIA), known as a rod drop accident (CRDA) in BWRs. Implementing and/or improving Bison modeling capabilities extended its analytical reach to areas beyond its original intended purpose. This paper demonstrates one of these capabilities as a proof of concept. Recently developed models enable Bison to provide an alternate approach for cladding integrity determination in CRDA evaluations, which currently use bounding conservative estimates. As part of this demonstration, NRC guidance on hydrogen-pickup and PCMI failure models were utilized in this research. Even though more research is needed in establishing right inputs and process in this area, this paper demonstrates Bison's ability to determine cladding integrity in a CRDA evaluation. This first of a kind demonstration is a proof of concept in this area, which could potentially be extended to a number of other areas where a more accurate cladding integrity determination would be needed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

2022 Accomplishments Report Template

Chromium-doped UO 2 fuels is an industrially important Accident Tolerant Fuel (ATF) due to the ability to increase the average grain size relative to standard UO 2 . However, there is a lack of knowledge of the thermophysical properties of these doped fuels, which is necessary for accurately modeling the in-pile performance. The work accomplished this FY focused on developing large grain Cr-doped UO 2 pellets as a function of Cr content and subsequently measuring the thermophysical and mechanical properties. The microstructural variations were also investigated to evaluate chromium formation within the microstructure of the sintered pellets.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Evaluation of Radiolysis Data for Hydrogen Gas Generation During Gamma Irradiation of Pre-Corroded and Pristine Aluminum Samples - An Aluminum SNF Dry Storage Study Interim Report

Information and data from radiolysis testing to measure hydrogen (H₂) generated from hydrated oxides on aluminum exposed to ⁶⁰Co radiation were reviewed to evaluate hydrogen generation with radiation dose. Hydrogen generation rate is a primary input to the coupled thermal-chemical system model of the performance of aluminum-clad spent nuclear fuel in a dry storage canister (ASNF-in-canister) over its storage life. Hydrogen generation data and hydrogen generation rate (hydrogen generated per absorbed radiation dose, $\Delta$(H₂)/$\Delta$(dose)) were reported in two reports for radiolysis testing [1, 2]. The testing involved first-time data¹ for radiolytic yield of hydrogen from hydrated-oxides-on-aluminum substrates, and from pristine (non-corroded) aluminum substrates. The laboratory-grown hydrated-oxides-on-aluminum substrates contained a ~ 5 μm film consisting of the trihydroxide bayerite (Al(OH)₃) with amounts of the oxyhydroxide boehmite (AlOOH). These specimens were assumed to also contain physisorbed water (unquantified). The pristine substrates were flat coupons of aluminum that were ground to a 600-grit finish but were not immersed in water to grow a hydrated oxide. The pristine substrates were tested as companion specimens to provide information on hydrogen generation from material without chemisorbed water and were assumed to contain physisorbed water (unquantified).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Nuclide inventory validation: radiochemical assay data quality and modeling challenges in benchmark models development

Oak Ridge National Laboratory is conducting radiochemical assay experiments, using high-precision analytical protocols validated with a comprehensive quality assurance plan, to expand the nuclide inventory validation basis for high burnup spent nuclear fuel. Preliminary measurement data for key actinides and fission products in two pressurized water reactor spent fuel samples are being used to investigate the impact of measurement data uncertainty on the sample burnup estimation. These measurement data are also being used to examine the impact of assumptions applied when developing best-estimate computational models to simulate fuel irradiation history. The simulations are being performed using depletion capabilities in the SCALE code system. Comparison of calculated and measured nuclide concentrations shows good agreement for the considered nuclides. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗