Accident Tolerant Fuel Modeling in Bison
An overview of accident tolerant fuel modeling capabilities in Bison.
SEARCH · Engineering Papers
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.
An overview of accident tolerant fuel modeling capabilities in Bison.
The renewed interest in metallic U-Zr nuclear fuel alloy has led to a drive for deeper understanding of the mechanisms driving the phenomena observed under irradiation conditions. The Department of Energy Advanced Fuel Campaign has developed infrastructure to support metallic fuel development, including Post Irradiation Examination (PIE) of legacy Fast Flux Test Facility (FFTF) Mechanistic Fuel Failure (MFF) experiments. The PIE performed on legacy FFTF MFF experiments gives insight on metallic fuel performance and can address the lack of knowledge and scarcity of reliable data identified in several studies over recent years. Unfortunately, PIE efforts can cost significant time and resources which can impede the progress of metallic U-Zr fuel development. Metallic U-Zr fuel performance modeling can be used to inform PIE efforts on regions of interest for relevant investigations and can help understand phenomena observed in PIE. This report demonstrates the current progress of FFTF MFF fuel performance simulations using the BISON fuel performance code and discusses the support provided by simulation to PIE efforts. Progress in temperature, profilometry, fission gas release, plenum pressure, and zirconium redistribution simulation results have been demonstrated.
A Modeling and Simulation (M&S) exercise is being performed for the High burnup Experiments for Reactivity initiated Accident (HERA) project under the Nuclear Energy Agency (NEA) Framework for Irradiation Experiments (FIDES) program. The goal of the M&S exercise is to improve M&S and experiment integration, facilitate community involvement in experiment design and interpretation, facilitate community collaboration, and aid in ensuring program data meet fuel performance code needs. The M&S exercise will compile and compare results from over 20 international organizations using 14 different fuel performance codes. This paper presents the results from the BISON fuel performance code generated by the Idaho National Laboratory participants.
Fission product mass transfer within and between TRISO particle and compact layers is an important phenomenon. It directly impacts fission product release predictions, which are used as source terms for safety and licensing calculations. Modeling fission product transport within layers and across bonded interfaces is relatively straightforward under the assumptions of isotropic Fickian diffusion, concentration continuity, and flux continuity. Modeling fission product transport across gaps and debonded layers is more difficult. Gaps are known to form between the buffer and inner PyC (IPyC), and debonding may occur between the IPyC and silicon carbide (SiC). A new mass transfer model was developed to provide a more accurate and robust fission product release calculation by accounting for interlayer sorptivity. One importance of the new model is to account for temperature and material changes across the gap based on sorption isotherm. This work presents the development of the sorption behavioral model and a fission product trapping model, and analyses the subsequent fission product diffusion behavior, specifically cesium, through TRISO particles and compact materials using the BISON fuel performance code.
Presentation requested by our collaborators on the technology commercialization fund project. This presentation provides an overview of the axial gas communication capabilities in the BISON fuel performance code. A discussion on models that impact gas communication, the theory behind the implementation, verification, and validation cases are presented.
Lunchtime presentation requested by Advanced Fuels Campaign colleagues for their visit with visitors from the Japan Atomic Energy Agency. The presentation provides an overview of the BISON fuel performance code and covers a few already published examples related to accident tolerant fuel research using the software.
A refined BISON-Griffin workflow for designing targeted TRISO transient analyses in the Transient Reactor Test Facility
NEAMS/AGR BISON modeling of TRISO updates
The process for developing and qualifying nuclear fuels for commercial nuclear application requires fundamental material development, characterization, and design; out-of-pile testing on unirradiated materials; integral fuel rod irradiations, testing, and postirradiation examinations; and transient analyses. The historical approach depends on the generation of large empirical datasets and series of integral fuel rod irradiations, and this approach ultimately takes ~20 years—or sometimes longer—to acquire data through extensive sequential testing. Thus, the qualification and eventual deployment of new fuel systems constitute a long process. However, recent technological advancements have provided researchers the opportunity to perform out-of-cell, in situ measurements to assess material performance for the duration of the experiment. One such example of this capability is the use of digital image coordination and thermal imaging to assess Zircaloy cladding performance under a simulated loss-of-coolant accident (LOCA) transient condition. In situ measurements generally provide high-fidelity strain, strain rates, and temperature surface maps. This is critical for the US nuclear industry, which is actively developing a technical basis to support extending the peak rod average burnup from 62 to ~75 GWd/tU and the deployment of accident-tolerant fuel. However, the US Nuclear Regulatory Commission (NRC) outlined in its research information letter several technical issues that the industry must address before extending burnup. One topic of specific interest is understanding the cladding balloon and rupture geometry during the LOCA heat-up phase. By leveraging these advanced in situ capabilities, this work used in situ data generated from a simulated LOCA to better understand high-temperature creep and its effect on Zircaloy balloon and rupture performance. Here, this work used the BISON fuel performance code to assess the high-temperature creep model predictions with in situ data.
During the initial rise to power in light water reactors (LWR), thermal gradients within a pellet cause radial and axial cracks to form in the fuel. The effect of these cracks is to reduce the pellet-cladding gap and accelerate the interaction between the fuel and cladding. This process is known as fuel relocation and may also include contributions from pellet eccentricity and cladding ovality. Since the cladding experiences both elevated temperatures and high external pressure due to the coolant, the cladding typically creeps inward further reducing the pellet-cladding gap. Once the gap is closed and pellet cladding mechanical interaction (PCMI) begins, both the thermal and mechanical behavior of the fuel is affected. The compressive forces exerted on the fuel due to contact with the cladding cause the fractured fuel sections to move back toward their original position which is termed relocation recovery. A model for this phenomenon is implemented in the BISON fuel performance code and applied to a set of validation test cases. In conclusion, the predicted fuel rod diameter is compared to experimental measurements to evaluate the influence of relocation recovery over a range of operating conditions.
This report summarizes a NEAMS (Nuclear Energy Advanced Modeling and Simulation) project focused on sensitivity analysis of a fuels performance benchmark problem. The benchmark problem was defined by the Uncertainty Analysis in Modeling working group of the Nuclear Science Committee, part of the Nuclear Energy Agency of the Organization for Economic Cooperation and Development (OECD ). The benchmark problem involv ed steady - state behavior of a fuel pin in a Pressurized Water Reactor (PWR). The problem was created in the BISON Fuels Performance code. Dakota was used to generate and analyze 300 samples of 17 input parameters defining core boundary conditions, manuf acturing tolerances , and fuel properties. There were 24 responses of interest, including fuel centerline temperatures at a variety of locations and burnup levels, fission gas released, axial elongation of the fuel pin, etc. Pearson and Spearman correlatio n coefficients and Sobol' variance - based indices were used to perform the sensitivity analysis. This report summarizes the process and presents results from this study.
This report describes the coupling scheme involving power distribution from MCNP to Bison.
In this report, load-follow simulations using VERA-CS with one-way coupling to standalone BISON has been demonstrated including both a single rod with a full cycle of load-follow operations and a quarter-core model with a single month of load-follow.
The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program under the US Department of Energy (DOE) consists of three product lines. The Fuels Product Line (FPL), hosted at Idaho National Laboratory (INL), has produced a multiphysics framework with the name Multiphysics Object Oriented Simulation Environment (MOOSE) (Gaston et al. 2015). The FPL has also produced several high-fidelity physics codes including the fuel performance code Bison (Hales et al. 2016). The Reactors Product Line (RPL), hosted at Argonne National Laboratory (Argonne), has developed the high-order spectral-element Computational Fluid Dynamics (CFD) code Nek5000 (Fischer et al. 2015). The Integration Product Line (IPL), hosted at Oak Ridge National Laboratory (ORNL), is tasked with ensuring tools from the other product lines can be coupled together so that users may perform high-fidelity multiphysics simulations.
This report presents the results of a parametric thermo-mechanical parametric analysis of SiC/SiC cladding performed using the BISON fuel performance code. In this analysis different steady state linear heat rates (LHR) (18kW/m, 22kW/m, 26kW/m and 30kW/m) and initial fuel rod-gap thicknesses (50μm, 65μm, 80μm and 95μm) were considered along with three different model implementation of fuel creep. These include 1) no fuel creep 2) fuel creep based on the Hagrman model and 3) fuel creep based on a modified model. The fuel rod geometries used in this analysis are representative of a generic 17 x 17 fuel rod assembly used in a Westinghouse PWR.
The Virtual Environment for Reactor Applications (VERA) [1] being developed by the Consortium for Advanced Simulation of Light-Water Reactors (CASL) allows a user to runnumerous modeling codes with a single input file. This includes modeling codes such as MPACT/CTF [2,3], a neutronics/thermal hydraulics code, and Bison [4,5], a fuel performance code. This is achieved through the use of preprocessors which act upon the common VERA input and subsequent inputs.
SiC/SiC composites are considered leading candidates for accident tolerant fuel cladding in light water reactors. However, when SiC-based materials are exposed to neutron irradiation, they experience significant changes in dimensions and physical properties. In the extreme nuclear reactor environment, SiC-based fuel cladding will see neutron damage, significant heat flux, as well as a corrosive environment. In this report, we evaluate the effect of non-uniform changes in the dimensions caused by neutron irradiation with variable temperatures, along with the closing of the fuel – clad gap, on the stress development in the clad over the course of time. To ensure reliable and safe operation of accident tolerant fuel cladding concepts such as SiC-based materials, it is important to assess its thermo-mechanical performance under in-reactor conditions of irradiation and realistic temperature distributions. In this work, we have used the BISON fuel performance modeling code to perform the thermo-mechanical analysis of SiC/SiC cladding. A constitutive model is constructed and solved numerically to predict the stress distribution and variation in the cladding under normal operating conditions. The dependence of dimensions and physical properties variation with irradiation and temperature has been incorporated. These initial results from the analyses provides insights about the stress distribution and variation with time in the cladding as well as the interaction of fuel pellet with the cladding under different conditions of initial fuel pellet-cladding gap and steady state power levels.
In order to improve the accident tolerance of light water reactor (LWR) fuel, coated cladding materials have been proposed which exhibit much slower oxidation kinetics in high-temperature steam than Zr-alloys. This behavior should decrease the energy release due to oxidation and allow the cladding to remain integral longer in the presence of high-temperature steam, making accident mitigation more likely. In order to determine how these coating materials will adjust the behavior of the Zr-alloy cladding during normal reactor operation, simulations have been performed using the Bison fuel performance code. This analysis focused on modeling and comparing the integral thermo-mechanical performance of coated and uncoated zircaloy cladding with uranium dioxide fuel under normal operating conditions. This analysis demonstrates the ability to model a discrete meshed coating bonded to the cladding exterior and provides a comparison with an uncoated fuel rod. A generic range of coating material properties are implemented given the differences and unknowns regarding the proprietary coating properties used by the vendors. A main-effects analysis was performed to identify trends and sensitivities in how coating material properties affect key fuel performance metrics. Lastly, because residual stresses are expected to form in the coating, arising from the high temperature coating application process and application techniques, an additional analysis was performed to identify the impact of these stresses on the cladding behavior. Future analysis efforts will target more comprehensive documentation of this main effects analysis by sampling the results at the beginning and middle-of-life in order to capture important behavior before and during the onset of gap closure.