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 451 records · Page 25

Quantifying deformation during Zry-4 burst testing: a comparison of BISON and a combined in-situ digital image correlation and infrared thermography method

The development and verification of nuclear fuel-cladding performance codes require extensive testing to establish empirical correlations, necessitating more accelerated testing methods that can replace large validation studies. Here, in the present work, a technique was developed and implemented to experimentally quantify the relationship between temperature, cladding internal pressure, and strain in-situ during a burst test, with the specific aim of generating data for code validation and model refinement. Digital image correlation was used to measure hoop, axial, and radial strain, and infrared thermography to quantify axial temperature gradients. Several key experimental modifications to traditional LOCA burst testing were necessary, but are shown to effectively have little impact on burst temperatures. The time/temperature dependent pressure and time/spatially dependent thermal gradient data were used with BISON to simulate cladding burst and strain rates. Although the measured DIC strains reach appreciable amounts at slightly lower temperatures than BISON simulated strains, there is good agreement between the measured and simulated strains in terms of magnitude and rates leading up to burst. However, in the few seconds prior to and during burst, the BISON simulated strains are significantly lower than the measured strains, indicating the potential for model improvement. The combined experiment and simulation technique may be applied to any new cladding concept to accelerate fuel qualification.

36 MATERIALS SCIENCE↗

Interface Problem Formulation Improvements with Application to Nuclear Fuel Performance Analysis

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation 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 documents new developments and robustness improvements in mechanical and thermal (gap heat transfer) contact formulations. The improvements range from the migration of industrial level (“assessment”) nuclear fuel model setups to the usage of mortar formulations, the addition of frictional contact to one-dimensional layered representations of fuel and cladding components, and the addition of the Petrov-Galerkin approach to dual mortar, which improves performance on curved, relatively coarse meshes. In addition, the Lagrange-multiplier enforcement of mechanical mortar contact constraints has been extended to two additional types of enforcement: penalty and augmented Lagrange-Uzawa. We show that the latter approach yields the same interface results as dual mortar in the Multiphysics Object-Oriented Simulation Environment, with the advantage of not worsening the condition number of the system matrix—thereby enabling the use of some general implementations of iterative preconditioners, at the expense of additional system evaluations (i.e., augmentations).

42 ENGINEERING↗

Fuel Performance Analysis of Chromium-Coated Cladding under Burst Conditions

To reduce the oxidation of zirconium-based alloy cladding at high temperatures, accident tolerant fuel systems have been proposed. Of the concepts identified, chromium-coated cladding has been shown to slow oxidation without greatly impacting the fuel system geometry or neutronic performance. To determine how coated-cladding tubes will perform under high-temperature accident conditions, pressurized-tube burst tests have been performed using the Severe Accident Test Station at Oak Ridge National Laboratory. To begin modeling these tubes to better understand how the coating will impact cladding behavior, these burst tests were simulated with the BISON fuel performance code. Cladding tube surface temperatures for the burst test were developed by fitting thermocouple data into axial and azimuthal profiles, while pressure data were compared until cladding failure. The temperatures at failure and the pressure evolution show relatively good agreement between the simulation and experiment results. This is the first step of a larger effort to simulate the cladding deformation process under high-temperature transient conditions and assess the cladding margin to failure more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Distant Mirror: Solar Oscillations Observed on Neptune by the Kepler K2 Mission

Starting in 2014 December, Kepler (K2) observed Neptune continuously for 49 days at a 1-minute cadence. The goals consisted of studying its atmospheric dynamics, detecting its global acoustic oscillations, and those of the Sun, which we report on here. We present the first indirect detection of solar oscillations in intensity measurements. Beyond the remarkable technical performance, it indicates how Kepler would see a star like the Sun. The result from the global asteroseismic approach, which consists of measuring the oscillation frequency at maximum amplitude max velocity and the mean frequency separation between mode overtones delta velocity, is surprising as the max velocity measured from Neptune photometry is larger than the accepted value. Compared to the usual reference max velocity of the sun equal to 3100 microhertz, the asteroseismic scaling relations therefore make the solar mass and radius appear larger by 13.8 plus or minus 5.8 percent and 4.3 plus or minus 1.9 percent, respectively. The higher max velocity is caused by a combination of the value of max velocity of the sun, being larger at the time of observations than the usual reference from SOHO/VIRGO/SPM (Variability of solar IRradiance and Gravity Oscillations / on board SOHO (Solar and Heliospheric Observatory) / Sun PhotoMeters) data (3160 plus or minus 10 microhertz), and the noise level of the K2 time series, being 10 times larger than VIRGO's. The peak-bagging method provides more consistent results: despite a low signal-to-noise ratio (S/N), we model 10 overtones for degrees iota equal 0, 1, 2. We compare the K2 data with simultaneous SOHO/VIRGO/SPM photometry and Bison (Birmingham Solar-Oscillations Network) velocity measurements. The individual frequencies, widths, and amplitudes mostly match those from VIRGO and BiSON within 1 sigma, except for the few peaks with the lowest S/N.

planets and satellites: individual (Neptune) – s↗

RELAP-7 Development Updates

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory, and is the next generation tool in the RELAP reactor safety/systems analysis application series. RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). During Fiscal Year (FY) 2015, the RELAP-7 code has been further improved with expanded capability to support boiling water reactor (BWR) and pressurized water reactor NPPs analysis. The accumulator model has been developed. The code has also been coupled with other MOOSE-based applications such as neutronics code RattleSnake and fuel performance code BISON to perform multiphysics analysis. A major design requirement for the implicit algorithm in RELAP-7 is that it is capable of second-order discretization accuracy in both space and time, which eliminates the traditional first-order approximation errors. The second-order temporal is achieved by a second-order backward temporal difference, and the one-dimensional second-order accurate spatial discretization is achieved with the Galerkin approximation of Lagrange finite elements. During FY-2015, we have done numerical verification work to verify that the RELAP-7 code indeed achieves 2nd-order accuracy in both time and space for single phase models at the system level.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a gas-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory high performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON MOOSE-based applications. Griffin provides the reactor physics capabilities including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient solutions. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions. BISON solves thermal conduction problems for the average pebble and TRISO in the pebble-bed core thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel and moderator temperatures. We present results for the coupled steady state equilibrium core and a protected loss of flow event. Although this model is to be considered prototypical with regard to the capabilities in BlueCRAB, its results are consistent with published work by the Institute of Nuclear and New Energy Technology (INET) in China and other research entities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamical System Scaling Application to Zircaloy Cladding Thermal Response During Reactivity-Initiated Accident Experiment

New fuel design and development currently requires 20 to 25 years to be qualified for use by the nuclear power industry. The thermal-hydraulics community has taken advantage of scaling theory to design reduced scale experiments that correctly preserve dominant key phenomena while quantifying distorted phenomena. These techniques can be leveraged in the design and analysis of fuel performance experiments to help reduce the timeline associated with fuel design and development. This study uses the Dynamical System Scaling (DSS) method to analyze cladding temperature data from the recent SETH-C experiment in the TREAT facility and accompanying BISON simulations to assess dynamic distortions occurring throughout the fast power excursion transient. The DSS analysis revealed that on the cool down from peak cladding temperature that the fuel radial power profile is the most sensitive modeling parameter with a heterogenous radial peaking factor corresponding to the lowest distortion compared to a uniform energy deposition. For the heat up to peak cladding temperature the heterogeneous radial power profile corresponded to the shortest process action. Finally, for the heat up to peak cladding temperature, the gap conductance model sensitivity was quantified using process action and shows that the default Light Water Reactor gap conductance model corresponded to the longest process action.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multi-Physics Surrogate Modeling with Dimensionality Reduction

This summary provides a new algorithm to construct surrogate models for multi-physics, nonlinear, high dimensional models which typically suffer from the curse of dimensionality. To deal with the high computational cost and difficult error estimation, we employ recent advances in reduced order modeling (ROM) techniques to do an error-preserving reduction for the number of influential model parameters, followed by a conventional surrogate construction process. To implement ROM for a multiphysics model, a range finding algorithm (RFA) is used to identify a linear transformation in which the number of influential parameters can be minimized and the errors resulting from the non-influential parameters can be upper-bounded. The proposed algorithm is applied to the MAMMOTH coupled code system developed at the Idaho National Laboratory. In particular, we investigate the reducibility of Rattlesnake and BISON codes and their multi-physics coupling. The most important result of this work is that the reduction obtained with coupled physics is more significant than that with individual physics models. This summary explores the application of ROM techniques on multi-physics coupling codes system. Concise methodology is provided and selected numerical results are provided to show that coupling system provides additional reduction at interfaces. Future work could be invested to determine the optimum surrogate model format. (authors)

97 MATHEMATICS AND COMPUTING↗

Gaseous swelling of U 3 Si 2 during steady-state LWR operation: A rate theory investigation

Rate theory simulations of fission gas behavior in U 3 Si 2 are reported for light water reactor (LWR) steady-state operation scenarios. We developed a model of U 3 Si 2 and implemented into the GRASS-SST code based on available research reactor post-irradiation examination (PIE) data, and density functional theory (DFT) calculations of key material properties. Simplified peripheral models were also introduced to capture the fuel-cladding interaction. The simulations identified three regimes of U 3 Si 2 swelling behavior between 390 K and 1190 K. Under typical steady-state LWR operating conditions where U 3 Si 2 temperature is expected to be below 1000 K, intragranular bubbles are dominant and fission gas is retained in those bubbles. The consequent gaseous swelling is low and associated degradation in the fuel thermal conductivity is also limited. Those predictions of U 3 Si 2 performance during steady-state operations in LWRs suggest that this fuel material is an appropriate LWR candidate fuel material. Fission gas behavior models established based on this work are being coupled to the thermo-mechanical simulation of the fuel behavior using the BISON fuel performance multi-dimensional finite element code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A DISTANT MIRROR: SOLAR OSCILLATIONS OBSERVED ON NEPTUNE BY THE KEPLER K 2 MISSION

Starting in 2014 December, Kepler K 2 observed Neptune continuously for 49 days at a 1 minute cadence. The goals consisted of studying its atmospheric dynamics, detecting its global acoustic oscillations, and those of the Sun, which we report on here. We present the first indirect detection of solar oscillations in intensity measurements. Beyond the remarkable technical performance, it indicates how Kepler would see a star like the Sun. The result from the global asteroseismic approach, which consists of measuring the oscillation frequency at maximum amplitude ν {sub max} and the mean frequency separation between mode overtones Δ ν , is surprising as the ν {sub max} measured from Neptune photometry is larger than the accepted value. Compared to the usual reference ν {sub max,⊙} = 3100 μ Hz, the asteroseismic scaling relations therefore make the solar mass and radius appear larger by 13.8 ± 5.8% and 4.3 ± 1.9%, respectively. The higher ν {sub max} is caused by a combination of the value of ν {sub max,⊙}, being larger at the time of observations than the usual reference from SOHO /VIRGO/SPM data (3160 ± 10 μ Hz), and the noise level of the K 2 time series, being 10 times larger than VIRGO’s. The peak-bagging method provides more consistent results: despite a low signal-to-noise ratio (S/N), we model 10 overtones for degrees ℓ = 0, 1, 2. We compare the K 2 data with simultaneous SOHO /VIRGO/SPM photometry and BiSON velocity measurements. The individual frequencies, widths, and amplitudes mostly match those from VIRGO and BiSON within 1 σ, except for the few peaks with the lowest S/N.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Development of Accelerated Steady-state Test Capsule Experiments to Replicate EBR-II Fuel Behavior Using BISON Fuel Performance Analysis

Here in this work, BISON fuel performance calculations were performed to predict the fuel behavior of accelerated burnup U-Pu-Zr fuel, with temperature operation conditions of the fuel and the cladding mirroring conditions within EBR-II fuel pins. The temperature operating conditions within the FAST accelerated burnup rods were aimed at replicating EBR-II X447/X447A fuel surface and inner cladding surface temperatures. Due to the FAST capsule design, these temperatures can be replicated with fission rate densities being significantly increased. fuel performance modeling has not been assessed for novel experiments such as accelerated burnup utilizing the FAST capsule within ATR. This is an important step in understanding accelerated irradiation methods as many performance models are empirical models conforming to the results of PIE but do not always include physical models that would represent the changes in irradiation tests.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Demonstration of Load-Follow Simulation with VERA-CS and Standalone 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.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Pellet Cladding Mechanical Interaction Modeling Using the Extended Finite Element Method

As a brittle material, the ceramic UO2 used as light water reactor fuel experiences significant fracturing throughout its life, beginning with the first rise to power of fresh fuel. This has multiple effects on the thermal and mechanical response of the fuel/cladding system. One such effect that is particularly important is that when there is mechanical contact between the fuel and cladding, cracks that extending from the outer surface of the fuel into the volume of the fuel cause elevated stresses in the adjacent cladding, which can potentially lead to cladding failure. Modeling the thermal and mechanical response of the cladding in the vicinity of these surface-breaking cracks in the fuel can provide important insights into this behavior to help avoid operating conditions that could lead to cladding failure. Such modeling has traditionally been done in the context of finite-element-based fuel performance analysis by modifying the fuel mesh to introduce discrete cracks. While this approach is effective in capturing the important behavior at the fuel/cladding interface, there are multiple drawbacks to explicitly incorporating the cracks in the finite element mesh. Because the cracks are incorporated in the original mesh, the mesh must be modified for cracks of specified location and depth, so it is difficult to account for crack propagation and the formation of new cracks at other locations. The extended finite element method (XFEM) has emerged in recent years as a powerful method to represent arbitrary, evolving, discrete discontinuities within the context of the finite element method. Development work is underway by the authors to implement XFEM in the BISON fuel performance code, and this capability has previously been demonstrated in simulations of fracture propagation in ceramic nuclear fuel. These preliminary demonstrations have included only the fuel, and excluded the cladding for simplicity. This paper presents initial results of efforts to apply XFEM to model stress concentrations induced by fuel fractures at the fuel/cladding interface during pellet cladding mechanical interaction (PCMI). This is accomplished by enhancing the thermal and mechanical contact enforcement algorithms employed by BISON to permit their use in conjunction with XFEM. The results from this methodology are demonstrated to be equivalent to those from using meshed discrete cracks. While the results of the two methods are equivalent for the case of a stationary crack, it is demonstrated that XFEM provides the additional flexibility of allowing arbitrary crack initiation and propagation during the analysis, and minimizes model setup effort for cases with stationary cracks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fuel Performance Calculations of Tagged UO 2 MiniFuel Disks

Isotopic taggants are being considered as a novel form of forensic identifiers for nuclear fuels. As part of this effort, eighteen tagged MiniFuel disks were irradiated in the High Flux Isotope Reactor at Oak Ridge National Laboratory. In this work, the BISON fuel performance code was used to recreate the pellets and predict their behavior and their post-irradiation material properties. These predictions are to serve as a baseline to help investigators identify results of interest during post-irradiation examination of actual MiniFuel disks. The simulation methodology and models are described herein, along with the individual disk predictions. Several general trends were identified in this analysis. Based on the models, it is predicted that the pellets’ internal temperatures will rise over time, fission gas release will remain constant at around 1%, regardless of the burnup, and no grain growth will occur in any of the disks. The reasons for these predictions are discussed below. Probable shortcomings in the BISON predictions are addressed, and future work is proposed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fuel performance of chromium-coated zirconium alloy and silicon carbide accident tolerant fuel claddings

The U.S. Department of Energy’s Accident Tolerant Fuel program is focused on extending the time for fuel failure during postulated severe accidents compared to the standard UO2-Zr alloy fuel system. This research investigates the feasibility of two different chromium-coated cladding concepts, one of which is zirconium-alloy based and the other is composite-SiC based. Both claddings had 50µm coatings, deducted from the base layer thicknesses. The claddings were reported on, using the multi- physics fuel performance tool BISON, under steady-state PWR operating conditions as well as two transients: power ramp and a loss-of-coolant accident (LOCA). The chromium-coated claddings showed comparable thermo-mechanical performance to the reference Zircaloy-4 cladding. As chromium is reported to provide an order of magnitude improvement in oxidation resistance, it is expected to be a better alternative in accident scenarios. Simulation results for both concepts call for further experimental investigation as well as modeling of beyond design-basis accidents.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗