Search NASA⌕ Search

SEARCH · Search NASA

Results for “BISON 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessment of buffer-IPyC thermomechanical debonding behavior using new experimental strength data in BISON

TRIstructural ISOtropic (TRISO) fuel is a nuclear fuel commonly used in High Temperature Gas-cooled Reactors (HTGRs). A single sub-millimeter-diameter TRISO fuel particle consists of a spherical fuel kernel surrounded by four coating layers: a low-density pyrocarbon buffer layer, an inner pyrolytic carbon (IPyC) layer, a silicon carbide (SiC) layer, and an outer pyrolytic carbon (OPyC) layer. The kernel is commonly made of UO2 or a mixture of uranium carbide and uranium oxide (UCO). During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. Not all irradiated particles will experience buffer-IPyC gap formation. The debonding may be partial, or it may be nearly total. However, from post-irradiation examination of UCO TRISO fuels irradiated as part of the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, it was concluded that partial buffer-IPyC debonding was the most common type of buffer-IPyC interaction. To predict TRISO thermomechanical performance, multi-physics models have been built that are being continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for 1D, 2D, and 3D TRISO particle simulations. This code is used to calculate fuel temperature, kernel swelling, buffer densification, thermal and irradiation creep, fracture, and fission gas production and release. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated AGR TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. That observation is also implemented in the TRISO interface debonding model. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Full-core high-burnup BWR LOCA fuel performance analysis and FFRD susceptibility

The susceptibility of the boiling water reactor (BWR) Limerick Unit 1 to fuel fragmentation, relocation, and dispersal during a postulated large-break loss-of-coolant accident (LBLOCA) was calculated using a multiphysics framework. The simulations include full-core, rod-resolved neutronic, thermal hydraulic, and fuel performance models using the VERA, TRACE, and BISON codes. This work focused on the transient BISON simulations, which include both the normal operation and LBLOCA periods in the same simulations. Cladding integrity was assessed using two correlations that are included with BISON. make page break Several new BWR-specific features were recently added to BISON. This work represents the first time these features have been included in a core-scale set of simulations. This study hence evaluates the performance of these new models for an operating reactor with realistic operating conditions. Simulation results showed that cladding integrity was maintained (i.e., no rods burst). Finally, future work to improve BWR and PWR predictions using this framework is suggested.

BISON↗
DOE OSTI2021

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Measurements of fission gas release from nuclear fuel in support of BISON fuel performance analysis

This project established a fuel irradiation and fission gas release (FGR) loop that has been designed for implementation at the PULSTAR reactor of North Carolina State University (NCSU). The FGR loop is capable of analyzing various types of fuel samples including, e.g., uranium dioxide fuel (UO 2 ) and uranium silicide (U 3 Si 2 ) samples. The fuel samples were produced and characterized at Idaho National Laboratory (INL) and shipped NCSU. The samples are ready for irradiation in the FGR loop to a given total thermal neutron fluence and under controlled temperature conditions. An innovative high resolution gamma-ray spectrometry approach and instrumentation have been established to monitor the loop’s sweep gas and perform high accuracy measurements of the gas release rates relative to predicted birth rates to establish the release-to-birth ratio (R/B) for a given radionuclide. In addition, uncertainty minimization is achieved through the formulation of the measurement observables into relative release-to-birth indicators for selected radionuclides. The experimental findings will be utilized to enhance and validate the fission gas release models that are currently implemented in the BISON fuel performance code. The experimentally deduced physical parameters will be incorporated in BISON’s fission gas release models. Atomistic simulations of fuel behavior have been performed to investigate the mechanisms of fission gas migration in the fuel. At this stage, performance of the fuel irradiation experiments awaits licensing approval by the US NRC.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

PolyPole-1: An accurate numerical algorithm for intra-granular fission gas release

This paper describes the development of a new numerical algorithm (called PolyPole-1) to efficiently solve the equation for intra-granular fission gas release in nuclear fuel. The work was carried out in collaboration with Politecnico di Milano and Institute for Transuranium Elements. The PolyPole-1 algorithms is being implemented in INL's fuels code BISON code as part of BISON's fission gas release model. The transport of fission gas from within the fuel grains to the grain boundaries (intra-granular fission gas release) is a fundamental controlling mechanism of fission gas release and gaseous swelling in nuclear fuel. Hence, accurate numerical solution of the corresponding mathematical problem needs to be included in fission gas behaviour models used in fuel performance codes. Under the assumption of equilibrium between trapping and resolution, the process can be described mathematically by a single diffusion equation for the gas atom concentration in a grain. In this work, we propose a new numerical algorithm (PolyPole-1) to efficiently solve the fission gas diffusion equation in time-varying conditions. The PolyPole-1 algorithm is based on the analytic modal solution of the diffusion equation for constant conditions, with the addition of polynomial corrective terms that embody the information on the deviation from constant conditions. The new algorithm is verified by comparing the results to a finite difference solution over a large number of randomly generated operation histories. Furthermore, comparison to state-of-the-art algorithms used in fuel performance codes demonstrates that the accuracy of the PolyPole-1 solution is superior to other algorithms, with similar computational effort. Finally, the concept of PolyPole-1 may be extended to the solution of the general problem of intra-granular fission gas diffusion during non-equilibrium trapping and resolution, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

R&D Plan for RISMC Industry Application #1: ECCS/LOCA Cladding Acceptance Criteria

The Nuclear Regulatory Commission (NRC) is finalizing a rulemaking change that would revise the requirements in 10 CFR 50.46. In the proposed new rulemaking, designated as 10 CFR 50.46c, the NRC proposes a fuel performance-based equivalent cladding reacted (ECR) criterion as a function of cladding hydrogen content before the accident (pre-transient) in order to include the effects of higher burnup on cladding performance as well as to address other technical issues. A loss of operational margin may result due to the more restrictive cladding embrittlement criteria. Initial and future compliance with the rule may significantly increase vendor workload and licensee costs as a spectrum of fuel rod initial burnup states may need to be analyzed to demonstrate compliance. The Idaho National Laboratory (INL) has initiated a project, as part of the DOE Light Water Reactor Sustainability Program (LWRS), to develop analytical capabilities to support the industry in the transition to the new rule. This project is called the Industry Application 1 (IA1) within the Risk-Informed Safety Margin Characterization (RISMC) Pathway of LWRS. The general idea behind the initiative is the development of an Integrated Evaluation Model (IEM). The motivation is to develop a multiphysics framework to analyze how uncertainties are propagated across the stream of physical disciplines and data involved, as well as how risks are evaluated in a LOCA safety analysis as regulated under 10 CFR 50.46c. This IEM is called LOTUS which stands for LOCA Toolkit for US, and it represents the LWRS Program’s response to the proposed new rule making. The focus of this report is to complete an R&D plan to describe the demonstration of the LOCA/ECCS RISMC Industry Application # 1 using the advanced RISMC Toolkit and methodologies. This report includes the description and development plan for a RISMC LOCA tool that fully couples advanced MOOSE tools already in development in order to characterize and optimize plant safety and operational margins. Advanced MOOSE tools that are needed to complete this integrated evaluation model are: RAVEN, RELAP-7, BISON, and MAMMOTH.

10CFR50.46C↗
DOE OSTI2016

FUEL PERFORMANCE MODELING OF FECRAL CLADDING FAILURE DURING LOSS OF COOLANT CONDITIONS

Accident tolerant materials are defined as those that provide significantly increased response time in the event of an accident while providing similar or improved performance as the conventional UO2/Zircaloy-4 fuel rods during normal operation. One avenue for materials to demonstrate enhanced accident tolerance is to have improved reaction (e.g., oxidation) kinetics with steam, resulting in slower hydrogen (or other combustible gases) generation rate, while maintaining acceptable thermo-mechanical properties, fuel-clad interactions, and fission product behavior. Due to their improved reaction kinetics with steam compared to Zircaloy-4, iron-chromium-aluminum (FeCrAl) alloys are being considered as an accident tolerant cladding material. In order to determine whether FeCrAl alloys are indeed accident tolerant materials, comparisons between conventional zirconium-based claddings under accident conditions is required. To this end, we have developed a burst (failure) criterion for FeCrAl using existing experimental data for first-generation alloys. This model was incorporated into the BISON fuel performance code to perform comparative studies with Zircaloy-4 cladding tubes under non-oxidizing and oxidizing LOCA conditions. The separate effects LOCA conditions considered were adopted from the PUZRY and REBEKA test series. Due to the fact that FeCrAl alloys have a larger thermal absorption cross-section and higher stiffness than Zircaloy-4, cladding thickness will be smaller for FeCrAl than for Zircaloy-4. Three different cladding thicknesses were analyzed in this work. The results indicate that for both non-oxidizing (PUZRY) and oxidizing (REBEKA) LOCA conditions, FeCrAl has similar burst behavior (i.e., time to burst and pressure at burst) as Zircaloy-4. It should be noticed that this study used data for first-generation FeCrAl alloys to develop the burst criterion. Although no burst data?currently exists for C35M (a second-generation FeCrAl alloy under development at Oak Ridge National Laboratory), improved burst behavior is anticipated. As burst data for C35M becomes available the burst criterion will be updated and the simulations revisited to determine the accident tolerance of second-generation FeCrAl alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Fuel performance modeling of FeCrAl cladding failure during loss of coolant conditions

Accident tolerant materials are defined as those that provide significantly increased response time in the event of an accident while providing similar or improved performance as the conventional UO2/Zircaloy-4 fuel rods during normal operation. One avenue for materials to demonstrate enhanced accident tolerance is to have improved reaction (e.g., oxidation) kinetics with steam, resulting in slower hydrogen (or other combustible gases) generation rate, while maintaining acceptable thermo-mechanical properties, fuel-clad interactions, and fission product behavior. Due to their improved reaction kinetics with steam compared to Zircaloy-4, iron-chromium-aluminum (FeCrAl) alloys are being considered as an accident tolerant cladding material. In order to determine whether FeCrAl alloys are indeed accident tolerant materials, comparisons between conventional zirconium-based claddings under accident conditions is required. To this end, we have developed a burst (failure) criterion for FeCrAl using existing experimental data for first-generation alloys. This model was incorporated into the BISON fuel performance code to perform comparative studies with Zircaloy-4 cladding tubes under non-oxidizing and oxidizing LOCA conditions. The separate effects LOCA conditions considered were adopted from the PUZRY and REBEKA test series. Due to the fact that FeCrAl alloys have a larger thermal absorption cross-section and higher stiffness than Zircaloy-4, cladding thickness will be smaller for FeCrAl than for Zircaloy-4. Three different cladding thicknesses were analyzed in this work. The results indicate that for both non-oxidizing (PUZRY) and oxidizing (REBEKA) LOCA conditions, FeCrAl has similar burst behavior (i.e., time to burst and pressure at burst) as Zircaloy-4. It should be noticed that this study used data for first-generation FeCrAl alloys to develop the burst criterion. Although no burst data?currently exists for C35M (a second-generation FeCrAl alloy under development at Oak Ridge National Laboratory), improved burst behavior is anticipated. As burst data for C35M becomes available the burst criterion will be updated and the simulations revisited to determine the accident tolerance of second-generation FeCrAl alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

NEAMS Workbench and Bison Fuel Performance Remote Application Configuration

The mission of the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program is to develop, apply, deploy, and support state-of-the-art predictive modeling and simulation tools for the design and analysis of current and future nuclear energy systems. This is accomplished using computing architectures that range from laptops to leadership class facilities. The NEAMS Workbench [1,2,3] initiative will facilitate the transition from conventional computer codes to high-fidelity ones by providing a common analysis environment and interface for model creation, review, execution, output review, and visualization for integrated codes. The BISON fuel performance code [4] is a finite element method (FEM)-based software being developed by multiple organizations led by the Idaho National Laboratory for the US Department of Energy Office of Nuclear Energy. BISON can simulate a variety of fuel and cladding compositions and geometries in varying dimensionalities, depending upon the problem type and complexity, including one-dimensional (1D), 1.5D, 2D axially symmetric (2D-rz), 2D radial-azimuthal (2D-rθ), or 3D models. The NEAMS Workbench allows execution of the BISON fuel performance code on local and remote compute resources. In this report, the detailed instructions of configuring the NEAMS Workbench for remote execution of the BISON fuel performance code on a portable batch system (PBS)-scheduled compute resource, as well the transfer of associated Bison files, is illustrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

BISON Methods and Contact Algorithms

This paper reviews five critical algorithms for BISON’s simulation of nuclear fuel: the burnup model, the coolant channel model, fuel cracking models, mechanical contact, and thermal contact.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Development of a BISON validation case for the TRISO transient irradiations in NSRR using effective heat capacity methods

The current tristructural isotropic (TRISO) fuel assessment and validation database in BISON primarily covers steady- state irradiation and high-temperature furnace testing. Transient assessment cases are potentially needed to support U.S. industry efforts in designing and deploying commercial reactors using TRISO fuels. Historical transient tests in- volving TRISO fuels used highly conservative conditions compared to the typical high-temperature gas-cooled reactor accident scenarios. Despite this, modeling historical transient tests is fundamental for evaluating BISON’s predictive capabilities, adapting material properties for high-temperature and high-particle-power regimes, and developing a sys- tematic validation approach for TRISO transient applications. This work developed a 1D model of transient experiments carried out at the Nuclear Safety Research Reactor using BISON. BISON’s predictions of energy deposition, UO 2 melting onset, and molten volume fractions are compared against experimental measurements. Melting was modeled using an effective specific heat capacity model for UO 2 . We found that BISON’s predictions are in reasonable agreement with experimental data for low-energy-deposition cases, and that BISON overpredicts melting at higher energy depositions. We also discuss the potential causes of discrepancies between the simulated and measured results and propose ways to further develop the model. Although these simulations used conservative conditions compared to those expected for actual TRISO-fueled reactors, they extended the range of conditions reflected in the data in the existing BISON database for TRISO fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Coupling MOOSE-Wrapped MPACT to BISON

As part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort is being made to leverage codes like Virtual Environment for Reactor Analysis (VERA) and Michigan Parallel Characteristics Transport (MPACT) by coupling them with other NEAMS codes. To facilitate coupling with other Multiphysics Objected Oriented Simulation Environment (MOOSE) applications, a MOOSE-wrapped MPACT app is created. The capability of this app, named Trogdor, is demonstrated by coupling it with another MOOSE app, BISON. Several single pin problems were run to test the Trogdor and BISON coupling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

(Project 14-6770) An Investigation to Establish Multiphysical Property Dataset of Nuclear Materials Based on in-situ Observations and Measurements

In-core nuclear materials including fuel pins and cladding materials fail due to issues including corrosion, mechanical wear, and pellet cladding interaction. In most such scenario microstructure dependent and corrosioninduced chemistry dependent property changes significantly affect performance of cladding, pellet, and housing. Emphasis of this work was on replace conventional pellet-cladding material models with a new straingradient viscoplasticity model that is informed by transmission electron microscopy (TEM) based measurements and by nanomechanical Raman spectroscopy (NMRS) based measurements. The TEM measurements are quantitative in nature and therefore reveal stress-strain relations with simultaneous insights into mechanisms of deformation at nanoscale. The NMRS measurements reveal the similar information at mesoscale along with additional information on relating local microstructural stresses with applied stresses. The resulting information is used to fit constants in the strain gradient viscoplasticity model as well as to validate one. During TEM measurements, a micro-electro-mechanical system based setup was developed with mechanical actuation, sensing, heating, and electrical loading. Contrary to post-mortem analysis or qualitative visualization, this setup combines direct visualization of the mechanisms behind deformation with measurement of stress, strain, thermal and electrical properties. The unique research philosophy of visualizing the microstructure at high resolution while measuring the properties led to fundamental understanding in grain size and temperature effects on measured mechanical properties such as fracture toughness. A key contribution is the role of mechanical loading boundary conditions to deconvolute the insitu TEM based nanoscale and NMRS based mesoscale data to bulk behavior. First the literature based pellet cladding mechanical interaction model based on the work of Retel’s and Williamson’s in literature work to predict tempurature and stress distribution in cladding and pellet at normal operating condition was analyzed. Later the data was fitted to find constants for a viscoplastic strain gradient model. The developed model still needs to be refined and calibrated using various experimental results. That remains the focus of future work. Overall, a major thrust of the work was therefore on active control of the microstructure (grain size, defect density and types) exploiting the multi-physics coupling in materials. In particular, using experiments the synergy of current density, mechanical stress and temperature were studied to annihilate defects and recrystallize grains. The developed model is being examined for implementation in BISON. Multiple invited talks, international journal publications, and conference publications were performed by students supported on this work. Another output is support multiple PhD and masters thesis students who will be an important asset for future basic nuclear research. Future Work Recommendations: A nuclear reactor operates under significant variations of thermal loads due to energy cycling and mechanical loads due to constraint effects. Significant thermal and chemical diffusion takes place at the pallet-cladding level. While the proposed work established new experimental approach and new dataset for Zircaloy-4, the irradiation level was in the range of 1-2 dpa. Samples with higher dpa need to be examined. Therefore, a continual of support of the performed work is essential. Currently, these are the only experiments that can measure the produced data. The work also needs to be extended to different fuel types and cladding types such as SiC and FeCrAl based claddings. A combination of datasets for these materials can then be used to analyze accurately predict behavior of critical pellet cladding systems in accident scenario with high heat flux and high thermal loads. This is a BIG unknown as if now.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Implementing Fuel Cladding Bonding and Assessing its Impact on Axial Gas Communication

Fuel rods irradiated in light-water reactors to burnup values above 45 MWd/kgU are subject to the formation of a chemical and mechanical bond between the fuel and cladding upon gap closure. The formation of the bond subsequently inhibits the ability of fission gases released from the fuel to flow freely to the plenum of the rod. The flowing of gases within fuel rods is referred to as axial gas communication. During transients, such as loss of coolant accidents, the bond may influence cladding deformation prior to breaking. Upon bond breakage, gases are able to more freely communicate to the lower pressure regions of the rod. The impact of bonding on gas communication and the ballooning behavior of the cladding during a loss of coolant accident is of interest to the nuclear industry in support of burnup extension for the existing light-water reactor fleet. In this report, a model to capture the effects of fuel cladding bonding in the BISON fuel performance code is presented. The theory of the model along with implementation testing is provided. A summary of a previously developed axial gas communication is given to set the stage for how the two models may be coupled together. A full-length pressurized-water reactor fuel rod demonstration is highlighted to evaluate the impact of including bonding on axial gas communication calculations using a preliminary coupling methodology. An overview of the next steps regarding the modeling of bonding, axial gas communication, and a more tightly coupled framework is also provided.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Modeling and Measurement of Axial Gas Transport in Nuclear Fuels

This document is an American Nuclear Society (ANS) transaction for the ANS annual conference. This transaction provides an overview of the Modeling and Measurement of Axial Gas Transport Laboratory Directed Research and Development project including preliminary modeling results. The goal is to provide model for axial gas transport in the Bison Fuel Performance code to support industry's desire to achieve higher burnups in existing light water reactors. A brief overview of the experimental approach and associated modeling activities is presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

The application of MAMMOTH for a detailed tightly coupled fuel pin simulation with a station blackout

The accurate calculation of desired quantities to predict fuel behavior in a nuclear reactor requires the solution of inter linked equations representing different physics phenomena. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding outer surface. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. In this work the MAMMOTH application is employed in the demonstration of an approach that more consistent with first principles, in a simulation that couples the neutronics application Rattlesnake, the thermal fluids application RELAP-7 and the fuels performance application BISON. A single fuel pin was modeled based on the dimensions of a fuel rod from a Westinghouse 17*17 optimized PWR fuel assembly. The simulation consisted of a depletion period of 1343 days, representing three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted with a total pin power of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shutdown). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a tightly coupled approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain, which includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2-D axisymmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1-D flow equations. For this work, the mesh on the Rattlesnake side is either be 3-D (for low order PN transport or diffusion) or 2-D axisymmetric (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burnup are transferred accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burnup distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh to perform a fuels performance solution. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. For this paper, sets of results from the detailed calculation are provided for both depletion and the SBO event. We demonstrate that a detailed simulation akin to first principles can be achieved using MAMMOTH and the various MOOSE application on differing domains. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗