Mesoscale modeling of Microstructural Evolution in High Burnup UO2 Fuel and Transient Response
Overview of work done to model UO2 HBS modeling at an upcoming DOE-CEA meeting.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Engineering topics
Matthews, Christopher
Publications and source records attributed to Matthews, Christopher.
Empirical and mechanistic transient fission gas release model for high-burnup LOCA conditions
In response to the nuclear industry desire to extend burnup beyond current licensing practices, the US Nuclear Regulatory Commission (NRC) released a research information letter (RIL) that provides a basis for analyzing fuel fragmentation, relocation, and dispersal (FFRD) in light-water reactors. Of the five elements discussed, the most ambiguous is the significance of transient fission gas (FG) release (FGR) (tFGR) and its effects on fuel performance under loss-of-coolant accident conditions. In fresh fuel, FG migration and eventual release is primarily governed by diffusion-based mechanisms at higher temperatures (>1,000 °C). However, the mechanisms governing FGR changes as burnup increases. More recent research indicates that FGR increases as burnup increases, specifically under temperature transient conditions, and this release occurs at lower temperatures with a new release mechanism. This behavior has been attributed to microcracking in the fuel and is likely related to microstructure embrittlement with the presence of over pressurized FG bubbles. The NRC RIL outlines the complexity of the phenomenon and a need for a deeper understanding to adequately address FFRD for regulatory application. Here, therefore, this manuscript intends to summarize the publicly available tFGR data and discuss the observed dependencies (e.g., burnup, heating rate, sample geometry, terminal temperature). An empirical model has been developed and benchmarked against recently published experimental data. However, this empirical model is limited to conditions for which fitting data exist and, therefore, a high-level discussion is included to provide a roadmap for atomistically-informed multiscale modeling in conjunction with experimental data collection to develop a mechanistic tFGR model widely applicable to a broad range of nuclear fuel conditions at high burnup.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Development of bubble evolution model for new mechanistic transient fission gas release capability in BISON
This report summarizes efforts within NEAMS to investigate the mechanisms that govern fission gas behavior in UO 2 . In particular, the focus is on understanding how fission gas behavior causes transient fission gas release and fragmentation/pulverization of high burnup structure (HBS) in UO 2 . HBS forms in the periphery of the pellet where temperatures are relatively low. Previously, MD simulations were performed to determine the reaction energies for various Xe and U defects with bubbles, as a function of Xe to vacancy ratio or, equivalently, pressure. As had been shown in FY22, it was found that the unmodified version of the Simple Integrated Fission Gas Release and Swelling (SIFGRS) model within BISON greatly over-predicted the number gas atoms per vacancy in the bubbles in the outer rim of the pellet (a ratio of > 1 million). This was due to slow grain boundary vacancy diffusivity and not accounting for the pressure-dependent reaction energy for Xe interstitials with bubbles. The application of the pressure dependent reaction energies was able to restrict Xe to vacancy ratios to 2:1, which is far more realistic than those originally obtained from SIFGRS.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Development and formulation of physics based metallic fuel models and comparison to integral irradiation data
Metallic fuel has an important historical significance in the development of nuclear reactors and continues to be relevant to the progression of advanced test and power reactors. A number of models, ranging from empirical to mechanistic, have been developed and implemented in various fuel performance codes to describe U-Zr and U-Pu-Zr fuel and typical fast reactor cladding materials. One challenge of using these models to simulate fuel performance is the inevitable tangling of coupled phenomena that can cloud proper implementation, calibration, and eventual utilization of new models. Here in an effort to provide a baseline capability that will facilitate the use of advanced models, new capabilities have been implemented into the fuel performance code BISON specific to metallic fuel simulations, ranging from materials properties, fission gas release and swelling calculations, coolant channel models, and cladding correlations. These models have been applied to the X441/X441A EBR-II experimental assembly data, a set of irradiated metallic UPuZr fuel rods of varying pin designs. The models implemented in BISON are able to capture the general trend of the expected response of the fuel and cladding to irradiation in EBR-II, especially when considering the spread in experimental measurements and the uncertainties inherited from the historical material models. Ultimately, the models outlined here provide the baseline capabilities on which new models can build upon in order to improve the prediction of metallic fuel performance simulations in off-normal designs or operations.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Surrogate-assisted Data-Free Inference with summary statistics for predicting xenon diffusivity in uranium oxide nuclear fuel
Explore the source record for details and available documents.
Robbe, Pieterjan Magda M↗
The National Criticality Experiments Research Center and its role in support of advanced reactor design
The National Criticality Experiments Research Center (NCERC) located at the Nevada National Security Site (NNSS) in the Device Assembly Facility (DAF) and operated by Los Alamos National Laboratory (LANL) is the only general purpose critical experiments facility in the United States. Experiments from subcritical to critical and above prompt critical are carried out at NCERC on a regular basis. In recent years, NCERC has become more involved in experiments related to nuclear energy, including the Kilopower/KRUSTY demonstration and the recent Hypatia experiment. Multiple nuclear energy related projects are currently ongoing at NCERC. This paper discusses NCERC’s role in advanced reactor design and how that role may change in the future.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
2.0 - MOOSE: Enabling massively parallel multiphysics simulation
The last 2 years have been a period of unprecedented growth for the MOOSE community and the software itself. The number of monthly visitors to the website has grown from just over 3,000 to now averaging 5,000. In addition, over 1,800 pull requests have been merged since the beginning of 2020, and the new discussions forum has averaged 600 unique visitors per month. The previous publication has been cited over 200 times since it was published 2 years ago. This paper serves as an update on some of the key additions and changes to the code and ecosystem over the last 2 years, as well as recognizing contributions from the community.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
CALPHAD modeling of uranium nitride (UN) fabrication routes enabled by first-principles calculations
The thermochemical details of fabricating uranium nitride (UN) by ammonolysis of uranium tetraflouride (UF 4 ) were determined using density functional theory (DFT) and CALculation of PHAse Diagrams (CALPHAD) computational methods. The thermochemical data of all binary, ternary, and quaternary U-H-N-F phases were computed using DFT, and the data for the phases that have not been measured experimentally, including UN 2 and NH 4 F(g), were combined with existing experimentally-determined data for CALPHAD modeling. The DFT data were benchmarked using experimental Gibbs energy of reaction and experimental thermochemical data for individual species. Phase diagrams relevant to the ammonolysis reaction are depicted, showing regions of stability for solid U-N, U-F and U-N-F phases. An unidentified phase produced in a previous experiment was identified as UN 0.95 F 1.2 (UNF) by comparing its X-ray diffraction spectrum to the experimental spectrum, and its formation during the fabrication of UN from UF 4 is supported by the simulated phase diagram. Here, it is calculated that UN 2 can be produced by the ammonolysis of UF 4 , but requires elevated temperatures, high NH 3 (g) partial pressure, and large amounts of flowing NH 3 (g) to avoid solid flu oride impurities in the uranium nitride. Likewise, U 2 N 3 can be produced instead at temperatures greater than 980 K. The use of silane (SiH 4 ) gas was investigated as a potential additive in the ammonolysis fabrication route to speed removal of fluorine. The addition of SiH 4 (g) offers little advantage to the removal of fluorine, and adds the complication of Si 3 N 4 formation. The use of DFT to fill in missing data to perform CALPHAD calculations demonstrated here allows for the determination of more comprehensive and trustworthy phase diagrams than the use of existing experimental data alone.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Bayesian calibration for summary statistics with applications to a cluster dynamics model.
Abstract not provided.
Robbe, Pieterjan↗
Marine Energy to Hydrogen Working Meeting (Workshop Report)
The Marine Energy to Hydrogen Working Meeting was held virtually on February 17, 2021. The meeting was sponsored by the U.S. Department of Energy Hydrogen and Fuel Cell Technologies Office (HFTO).
08 HYDROGEN↗
Marmot V2
MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.
Aagesen, LarryK.↗
Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel
Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.
Aagesen, Larry K. (ORCID:000000034936676X)↗
Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels
TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Phase-field simulations of fission gas bubbles in high burnup UO2 during steady-state and LOCA transient conditions
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Metallic fuel cladding degradation model development and evaluation for BISON
Fuel cladding chemical interaction (FCCI) and coolant cladding chemical interaction (CCCI), both consume stainless-steel cladding in sodium-cooled fast reactors (SFRs), and effectively degrade cladding mechanical performance through thickness reduction. Therefore, FCCI/CCCI and corresponding cladding degradation models are essential for advanced fuel performance codes, such as BISON, to reliably predict cladding behavior in SFRs. Here, we report the development and implementation efforts for BISON's FCCI/CCCI correlations and cladding degradation models based on U.S. legacy metallic fuel data. The models were evaluated using Integral Fast Reactor (IFR) program X447/A experiment data supported by the ongoing integration project enabling coordinated application of BISON and EBR-II fuel irradiation and physics database (FIPD). Furthermore, the implemented models were demonstrated to improve BISON's capabilities of predicting cladding damage and degradation behavior that is more consistent with post-irradiation examination observations. Additionally, some limitations of current BISON modules are identified, which are to be overcome through the BISON-FIPD integration efforts.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Metal Hydride Simulations Using SWIFT
Transition metal dihydrides are considered to be good candidate materials for nuclear reactor moderator applications due to their ability to maintain high hydrogen densities up to high temperatures. Owing to the recent interest in metal hydrides for use in microreactor applications where the size requirements make advanced moderators exceedingly attractive, a clear need for capturing the evolving thermochemistry of these materials is evident. In an effort to provide a tool for reactor designers and experimentalists, the MOOSE-based SWIFT code (Stoichiometry With Internally Fluctuating Temperature) has been adapted to tackle the thermochemical behavior of yttrium hydride at elevated temperatures.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
FY2021 Progress Report on BISON Metallic Fuel Model Development and V&V Using EBR-II Legacy Data
In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification project in FY2021 are summarized. The cladding degradation model based on the FCCI/CCCI wastage calculations has been developed and implemented into BISON. A comprehensive evaluation of the cladding degradation model was performed based on FIPD data of the IFR experiment X447. BISON objects were also developed to enable direct use of time-varying cladding outer surface temperature profile as temperature boundary conditions, which proved to provide more accurate temperature predictions for the metallic fuel pins irradiated in EBR-II. Additionally, a new BISON object was implemented to enable direct comparison between BISON predicted data and FIPD-based post-irradiation examination (PIE) results, which would significantly facilitate BISON metallic fuel verification and validation (V&V) activities. These new BISON-FIPD integration features were used in the establishment of a low-burnup fuel swelling evaluation framework as demonstration. The framework was successfully used to evaluate current fuel swelling models based on the IFR experiment X423.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗