SEARCH · Search NASA

Results for “Fuel Rod Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Development of a high-fidelity multi-cycle model of the NuScale small modular reactor using VERA

With growing renewables penetration, there is increased interest in flexible power operation for nuclear reactors. For multi-unit SMRs, in particular the NuScale SMR, which is an integral pressurized water reactor, there are opportunities to optimize flexible power operation across multiple units to limit the degradation of structural and control components. Here, we focus on degradation of in-core components, specifically the control rods and reactor pressure vessel. To perform these studies a high-fidelity, multi-cycle representation of the NuScale SMR is required, with a detailed representation of the structural and control components. To this end, the NuScale SMR has been modelled using the Virtual Environment for Reactor Applications (VERA) software. The entire transition to equilibrium is simulated, from Cycle 1 through to the equilibrium cycle. The equilibrium cycle model shows a good agreement with the NuScale design certification application (DCA) results, with differences attributable to a combination of using public domain data for the present study, and methodological differences. K-effective, power distributions, reactivity coefficients, and boron letdown curves are compared and all found to closely match. This shows that the VERA model is suitable for further studies. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

ATF material model development and validation for priority fuel concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the Bison fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on the focus was on U 3 Si 2 fuel and FeCrAl cladding as part of a high impact problem through the Nuclear Energy and Advanced Modeling Simulation (NEAMS) program. Then, developments for Cr 2 O 3 -doped UO 2 fuel, and SiC-SiC and Cr-coated zirconium-based claddings began based upon industry interests. In late fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors (CASL) took over further ATF work in Bison in support of the Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC identified their list of priority fuel and cladding concepts, which included Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and Cr-coated zirconium-based and FeCrAl claddings. In particular, the NRC suggested that reports similar in form to NUREG/CR-7024 [1] that was developed for traditional LWR materials UO 2 and zirconium-based claddings (i.e., Zircaloy-4, M5®, ZIRLO) be created for the priority ATF concepts. The approach to ATF capability development in Bison since the beginning has been two-fold: (1) empirical correlations and (2) multiscale model development. Both approaches have uncertainty inherent to them. Uncertainty in empirical correlations is bounded by the experimental data upon which with the correlation was developed. Models developed through a multiscale approach have uncertainty associated with the lower length scale calculations and input parameters that must be propagated to the engineering scale model in Bison. In this report, the recommended models, their range of applicability (e.g., temperature, burnup), and associated uncertainty for the NRC priority fuel concepts Cr 2 O 3 -doped UO 2 and U 3 Si 2 are presented in a manner similar to the aforementioned NUREG. In addition, the Cr 2 O 3 -doped UO 2 models are validated to the Halden Reactor tess IFA-677.1 rods 1 and 5 and IFA-716.1 rod 1. For U 3 Si 2 the models are validated to the recent post-irradiation examination (PIE) data ATF-13 and ATF-15 rods that were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory. Finally, an uncertainty quantification (UQ) and sensitivity analysis (SA) is completed on a single rod for each fuel concept that takes into account the defined uncertainty in select models. The results of the UQ and SA analyses indicate that given the large uncertainty in some of the input ATF models that the uncertainty on the Bison predictions of fuel performance metrics of interest bound the available experimental data. Correlation coefficients are also reported that identified the uncertain inputs with the strongest correlation (positive or negative) with the outputs of interest. Fuel performance metrics investigated include fuel elongation, rod internal pressure, fuel centerline temperature, and fission gas release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

TRISO SiC Failure Probability for Reactivity Initiated Accidents in High-Temperature Gas-Cooled Reactors

This work analyzes the failure process of the silicon carbide (SiC) layer in tristructural isotropic (TRISO) during reactivity-initiated accident scenarios for a high-temperature gas-cooled reactor (HTGR) with BISON. Two cases are considered—a group control rod withdrawal (CRW) and a control rod ejection (CRE)—reproduced from a previous study. Failure probability is modeled using Weibull statistics, and worst-case scenario Weibull parameters are adopted to simulate the envelopes in BISON with a one-dimensional TRISO model. CRW scenario results are characterized by higher values of maximum energy deposition and final temperature and volumetric strain with respect to the CRE ones, but the latter have remarkably higher SiC failure probability, mainly due to the offset in strain rates between the two cases. This work also confirms the validity and conservatism of the performance envelopes produced in a previous work by replicating the envelope formulation using RELAP5-3D and RAVEN with a different sampling technique and obtaining consistent results. A sensitivity analysis using the Sobol variance decomposition method on SiC failure probability is then performed involving a set of inputs on both CRW and CRE. The two most important parameters are Weibull modulus and characteristic stress, and their relative importance depends on the specific case. The proposed interpretation of the results is that both energy deposition and strain rate influence the relative degree of importance of the failure parameters. Computation of 95% confidence intervals around worst-case scenario SiC failure probability values is also carried out for four different sets of Weibull parameters. Heren a new criterion for SiC TRISO quality classification built upon safety-based ranges of Weibull parameters is proposed to be integrated in future Fuel-Production Quality Assurance Plans.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

SATS Enhanced Capabilities and Demonstration of Improved Ramp Rates for In-Cell Testing

This report presents the development of the second-generation Severe Accident Test Station, referred to as SATS 2.0. SATS 2.0 represents a major advancement in addressing critical research needs for the United States nuclear industry, particularly in relation to fuel fragmentation, relocation, and dispersal (FFRD) concerns associated with the nuclear industry’s goal to extend fuel burnup beyond a peak rod average of 62 GWd/tU. The central objective of SATS 2.0 was the deployment a new furnace capable of achieving heating rates up to 100°C/s. This achievement addresses critical capability gaps, enabling assessment of prioritized research objectives established by the Electric Power Research Institute’s (EPRI) Collaborative Research on Advanced Fuel Technologies (CRAFT) working group. These research objectives are directly related to high-burnup loss-of-coolant accident (LOCA) conditions and the effects of prolonged exposure to elevated temperatures, both of which are crucial to enhancing nuclear safety and efficiency. The SATS 2.0 system builds upon prior experience with the original SATS system, which focused on evaluating accident-tolerant fuel (ATF) cladding concepts during accident conditions. However, the new system not only expands upon prior core capabilities by achieving higher heating rates with a better furnace but also plans to incorporate novel auxiliary systems to create a versatile platform for addressing current research needs. For example, a system was developed for the quantification and characterization of fission gas released during high-temperature transients. Additionally, in situ measurement capabilities were developed, such as digital image correlation which enabled the real-time capture of strain related to balloon and burst events and fiber-optic sensors that allowed for high-fidelity characterization of temperature gradients. Demonstration tests of the new 12-lamp furnace achieved heating rates up to 120°C/s. Notably, SATS 2.0 demonstrated its capacity to simulate LOCA burst tests at different pressures, yielding burst data that align with historical empirical models. Moreover, the system exhibited its capability to simulate complex conditions observed in anticipated operational occurrences, while effectively mitigating temperature overshoots. These accomplishments mark significant progress toward overcoming FFRD challenges and advancing the United States nuclear industry's safety basis and technical capabilities for extended burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Employing MACS/ViBRANT as a Surrogate MARVEL Reactor for Startup Reactivity Tuning and Supervisory Control Processes

Advanced nuclear reactors are a key part of the future of nuclear energy both in the United States and globally. They offer unique benefits for various energy-demanding applications, including use in remote locations, compact size, modular manufacturing, remote monitoring, low and/or variable power rating operation, and reliance on novel technologies to enhance operational safety. To achieve economic feasibility, advanced reactors must significantly reduce their workforces in comparison with the current fleet. Achieving this reduction will occur through reducing staff workloads using technology to achieve autonomous or semi-autonomous operations, demonstrated by comprehensive testing and validation activities. These operations will require both software and hardware platforms during the design and testing phases. While simulations are useful during the design phase, their performance can significantly deviate during actual deployment on hardware. This report presents the outcomes of a collaborative technical initiative between the U.S. Department of Energy (DOE) Microreactor Program (MRP) and Advanced Sensors and Instrumentation (ASI) Program. The collaboration utilized the Microreactor Automated Control System (MACS) hardware platform to bridge the gap between theoretical reactor design and actual startup and control operations. Two key use cases were investigated: facilitating the startup testing period and demonstrating supervisory control. The first use case details the key Microreactor Applications Research Validation and Evaluation (MARVEL) reactor startup physics testing activities conducted using the MACS platform. These activities included drum worth measurements, shutdown margin assessment, temperature feedback analysis, and scram time evaluation, as well as unique testing that would apply to the MARVEL reactor to demonstrate the testing methodologies in a low-risk environment. The MACS platform, serving as a surrogate representation of the MARVEL reactor, proved instrumental in performing these tests. The exercise revealed aspects that led to optimized processes, refined hardware design, and enhanced base software capabilities. By maturing methods and technologies in this manner, the initiative promises to reduce wasted time in the actual on-site reactor deployment effort, thereby saving significant time and resources. The second use case focuses on the development and implementation of supervisory control methods aimed at managing core tilt, which can result from asymmetrical operations or manufacturing imperfections in fuel rods or reactivity control devices. A key objective was to assess and compare the use of artificial intelligence (AI) for supervisory control. The effort aimed to define the role of supervisory control to enhance performance without risking control instability. This effort explored three distinct approaches: rules-based (RB) methods, optimization techniques, and reinforcement learning (RL) algorithms. Each approach was evaluated for its ease of implementation, its usability, and its effectiveness in responding to asymmetries in neutron flux. Comparative analysis of these approaches provided valuable insights into their applicability and effectiveness, offering a robust framework for advanced reactor operations. Together, these two use cases highlight the potential of hardware test beds to help streamline the design, operation, and control of advanced nuclear reactors. This collaborative effort underscores the importance of continued innovation and experimentation in achieving the next generation of safe, reliable, and economically viable nuclear energy solutions.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

CTF: A modernized, production-level, thermal hydraulic solver for the solution of industry-relevant challenge problems in pressurized water reactors

CTF is a thermal hydraulic (T/H) subchannel tool that has been extensively developed over the past ten years as part of the Consortium for Advanced Simulation of Light Water Reactors (CASL) program. The code was selected early in the CASL program for support of high-impact challenge problems that were found to be relevant to the nuclear industry and its currently operating fleet of pressurized water reactors (PWRs), including issues such as departure from nucleate boiling (DNB), crud-induced power shifts (CIPSs), and reactivity-insertion accidents (RIAs). By incorporating CTF into the multiphysics Virtual Environment for Reactor Application (VERA) core simulator software developed by CASL, CTF has become the primary means of providing fluid and fuel thermal feedback, as well as T/H figure-of-merits (FOMs) in large-scale reactor simulations. With the goal of solving industry challenge problems, CASL placed great emphasis on developing high-quality, high-performance, validated software tools that offer higher fidelity than what is currently possible with current industry methods. In support of this effort, CTF was developed from a research tool into an nuclear quality assurance (NQA-1)–compliant, production-level software tool that is capable of addressing the stated challenge problems and goals of CASL. This work presents a review of the major technological achievements that were realized in developing CTF over the past decade of the CASL program and presents an overview of the code solution approach and closure models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2015

Development of Adaptive Model Refinement (AMoR) for Multiphysics and Multifidelity Problems

This project investigated the development and utilization of Adaptive Model Refinement (AMoR) for nuclear systems simulation applications. AMoR refers to utilization of several models of physical phenomena which differ in prediction fidelity. If the highest fidelity model is judged to always provide or exceeded the desired fidelity, than if one can determine the difference in a Quantity of Interest (QoI) between the highest fidelity model and lower fidelity models, one could utilize the fidelity model that would just provide the magnitude of the QoI desired. Assuming lower fidelity models require less computational resources, in this manner computational efficiency can be realized provided the QoI value can be accurately and efficiently evaluated. This work utilized Generalized Perturbation Theory (GPT) to evaluate the QoI, by convoluting the GPT solution with the residual of the highest fidelity model determined using the solution from lower fidelity models. Specifically, a reactor core neutronics problem and thermal-hydraulics problem were studied to develop and utilize AMoR. The highest fidelity neutronics model was based upon the 3D space-time, two-group, nodal diffusion equations as solved in the NESTLE computer code. Added to the NESTLE code was the ability to determine the time-dependent GPT neutron flux. The lower fidelity neutronics model was based upon the point kinetics equations along with utilization of a prolongation operator to determine the 3D space-time, two-group flux. The highest fidelity thermal-hydraulics model was based upon the space-time equations governing fluid flow in a closed channel around a heat generating fuel rod. The Homogenous Equilibrium Mixture (HEM) model was used for the fluid and Finite Difference Method was applied to both the coolant and fuel pin energy conservation equations. The lower fidelity thermal-hydraulic model was based upon the same equations as used for the highest fidelity model but now with coarse spatial meshing, corrected somewhat by employing effective fuel heat conduction values. The effectiveness of switching between the highest fidelity model and lower fidelity model as a function of time was assessed using the neutronics problem. Based upon work completed to date, one concludes that the time switching is effective in annealing out differences between the highest and lower fidelity solutions. The effectiveness of using a lower fidelity GPT solution, along with a prolongation operator, to estimate the QoI was also assessed. The utilization of a lower fidelity GPT solution was done in an attempt to avoid the high computational burden associated with solving for the highest fidelity GPT solution. Based upon work completed to date, one concludes that the lower fidelity adjoint solution is not sufficiently accurate with regard to estimating the QoI; however, a formulation has been revealed that may provide a path for addressing this shortcoming.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2018

Final Scientific/Technical Report, Multi-Sensor Inspection and Robotic Systems for Dry Storage Casks

Once removed from service, fuel rod assemblies are cooled in a pool and then placed in dry storage until they can be transported to a repository for disposal. The extended duration that used nuclear fuel spends in dry storage makes it prudent to be able to nondestructively inspect the dry storage casks. The degradation mechanisms of most concern are chloride-induced stress corrosion cracking of the stainless steel canister that confines the spent fuel and distributed cracking of the concrete in the overpack structure. Multi-sensor inspection and robotic systems were researched, developed, built, and tested. One system operates completely outside the cask to inspect the concrete. The other system, known as the PRINSE system, for Proactive Robotic Inspection of Nuclear Storage Enclosures, enters the cask through an outlet vent and navigates around the outer surface of the stainless steel canister. The PRINSE system delivers sensors for: (i) real-time characterization of the canister surface in order to assess whether chlorides are present, (ii) detecting the presence of cracks using shear-horizontal ultrasonic guided waves, (iii) measuring the temperature of the canister surface, and (iv) a gamma radiation counter. The four sensing modalities are applicable to any cask containing a stainless steel canister, however the robotic delivery system is designed specifically for Holtec’s HI-STORM family of casks. Horizontal axis casks have exposed concrete surfaces. A fully contactless air-coupled scanning ultrasonic test platform and a frequency-wavenumber domain data analysis method are used for these concrete structures. Analysis, numerical simulation and experimental results demonstrate that distributed cracks in concrete can be detected and cracked regions localized within a scanned region of concrete. Vertical axis casks (like the HI-STORM family) have steel-clad concrete overpacks. In a separate effort, a semi-coupled test platform and a guided wave attenuation-based data analysis method are used for disbond detection. Analysis, numerical simulation and experimental results show that the bonding conditions of steel-concrete interfaces can be quantitatively characterized using the guided wave signal analysis scheme. The PRINSE system is designed for vertical axis dry storage casks, specifically the HI-STORM family. The double pulse laser-induced breakdown spectroscopy (LIBS) system takes advantage of an optical fiber to deliver intense light to the canister surface that enables surface composition characterization in-situ in real time. Chlorine concentrations as low as 10 mg/m2 were detected on stainless steel. Additionally, a simpler laser setup was created to detect sodium as a surrogate for chlorine, and sodium was detectable when the chlorine concentration was as low as 0.5 mg/m2. The electromagnetic acoustic transducers (EMATs) enable detection of stress corrosion cracks. The EMATs generate shear-horizontal waves that propagate around the canister and are reflected by cracks enabling 100% of the canister weld lines to be inspected. The EMAT-based crack detection system was tested at the Pacific Northwest National Laboratory on a canister containing blind defects and was shown to detect the crack sizes it was designed to detect. The PRINSE system is a novel marsupial robot. The command center and winch outside the cask control the delivery arm and train of cars that traverse the outside of the canister from the cask ventilation system. The sensors are housed in the cargo bays of three cars that comprise the train, which is lowered to the bottom of the canister by the winch. Encoders on the winch enable the sensors to be delivered to a prescribed location, e.g., for re-inspection of a problem area. The PRINSE system does not mar the canister surface, can be extracted manually in the case of a power failure, and is designed for an elevated temperature and gamma radiation environment, although not all components of the prototype are made from the final materials. The PRINSE system was tested on a full canister that was placed inside a HI-STORM cask at the Holtec Manufacturing Division. It was also demonstrated on a full-height mockup at Penn State to a live audience.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

ZPPR-15 Small Reactivity Worth Experiments

In addition to measuring the worth of sodium voiding and the worth of simulated control rods, the ZPPR staff measured the Doppler worth of four standard samples, the worth of axial expansion and the worth of radial bowing in ZPPR-15. Doppler sample worths were measured in ZPPR-15A, ZPPR-15B and ZPPR-15D. The Doppler samples were natural UO 2 , depleted uranium metal, depleted U-10Zr alloy and 33% enriched UO 2 . These samples were placed in a sample capsule in the inner core of ZPPR-15 and heated to various temperatures between 300 K and 1100 K. The resulting changes in reactivity relative to the reference configuration were measured to determine the Doppler worths of the samples. Special segmented drawers were constructed to measure the worth of axial expansion in ZPPR15. The segments were connected to each other, and the back segment was connected to a spring-loaded cable. Pulling the cable introduced small gaps between drawer segments to simulate axial expansion, and releasing the cable eliminated those gaps. The worth of this simulated axial expansion was measured in ZPPR-15A, ZPPR-15B and ZPPR-15D. The worth of radial bowing was measured in ZPPR-15, but the experimenters regarded the result as unsatisfactory. Radial bowing was measured by two methods in ZPPR-17A, so the ZPPR-17A bowing measurements were analyzed in place of the unpublished ZPPR-15 radial bowing measurement. The first ZPPR-17A measurement consisted of rearranging the plates in 96 core drawers at the core-radial blanket boundary to simulate outward motion of the fuel during bowing. The second ZPPR-17A measurement used a new bowing oscillator to move the core plates a small distance vertically in a special drawer. The published uncertainties for the Doppler worth, axial expansion worth and radial bowing worth measurements are the statistical uncertainties in the measurements. In reality, the published uncertainty for each of these measurements is just one component of the total uncertainty. There are additional uncertainties specific to each measurement and a common uncertainty related to conversion from the natural measurement units, cents, to pcm for calculations. A full uncertainty analysis was performed for each measurement. The significant uncertainties were quantified, and a total uncertainty was determined for each measurement. The ZPPR-15 and ZPPR-17A experimental records were used to create detailed as-built Monte Carlo models of the Doppler worth, axial expansion and radial bowing measurement configurations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

Fundamental Studies of Tritium Diffusivity and Formation on the Surfaces of Defective LiAlO 2 Pellets

One of the important components of tritium-producing burnable absorber rods (TPBAR) is the γLiAlO 2 layer used in the form of an annular ceramic pellet enriched with the 6 Li isotope. Within TPBAR architecture, this oxide is located between the zircaloy-4 liner and nickel-plated zircaloy-4 tritium getter. When irradiated in a pressurized water reactor (PWR), the 6 Li pellets absorb neutrons, simulating the nuclear characteristics of a burnable absorber rod, and produce tritium ( 3 H, T). The 3 H specie chemically reacts with the metal getter, leading to the formation of a metal hydride. Accurate analysis of the 3 H transport through the ceramic pellets and the barrier/cladding system is hampered by the lack of fundamental data about the hydrogen isotope solubility and diffusivity. In TPBAR architecture, between the γ-LiAlO 2 pellet and the Ni-plated Zircaloy-4 getter (or liner), there is a gas gap. After 3 H diffuses from bulk γ-LiAlO 2 to its surface, it is possible that different 3 H species (such as T, T 2 , T 2 O, and OT, etc.) are formed on this surface which can then get into the gas gap region. Furthermore, upon recombination and desorption molecular tritium can be formed on the γ-LiAlO 2 surface. Therefore, exploring the formation of 3 H-related species on γ-LiAlO 2 surface becomes very important as an integral part of the study of 3H solubility and diffusivity in the pellets.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

OECD/NEA MPCMIV Benchmark - Preliminary fuel performance results

The on-going OECD/NEA Multi-physics Pellet Cladding Mechanical Interaction Validation (MPCMIV) benchmark aims to provide guidance on multi-physics validation through the modelling of two cold ramps. In this paper, the first results for the base irradiation of the father rod and fuel rodlet (refabricated from the father rod), and the first cold ramp are presented. The base irradiation consists of 3 years of irradiation in the Forsmark-2 reactor. The cold ramp encompasses a steady-state pre-ramp period of less than one hour at a low constant linear heat rate (LHR) followed by a ramp test (< 1 min) with a much higher maximal LHR. The base irradiation is modelled using the fuel performance codes FRAPCON and FAST, while the cold ramp modeling is using the fuel performance code FRAPTRAN. Several missing parameters for FRAPCON/FAST base irradiation models have been selected using multiple references such as the OECD/NEA light water reactor Uncertainty Analysis in Modelling (UAM) benchmark specifications and the FRAPCON Integral Assessment report. The obtained results for the base irradiation such as the cladding outer diameter show reasonable agreement with the experimental measurements. The results for the cold ramp show larger differences with the experimental measurements (e.g., cladding axial elongation). However, such differences have been observed as well in other studies involving pellet cladding mechanical interaction analyses and are attributed to the fuel performance modelling assumptions and to the tuned modelling parameters that were not covered by the specifications. Uncertainty and sensitivity analysis might allow a better quantification of these missing parameters. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Sodium Void Worth and Control Rod Worth Measurements in ZPPR 15

Sodium void worth and control rod worth were measured in Phases A, B, C and D of the ZPPR15 experimental program. Sodium voiding was simulated by replacing stainless steel cans containing sodium by empty stainless steel cans of the same sizes. The control rod worth measurements were measurements of simulated control rods (CRs) and simulated empty control rod positions (CRPs). A CRP consisted of a 2 x 2 array of drawers filled with sodium plates. A CR consisted of a 2 x 2 array of drawers containing alternating columns of boron carbide plates and sodium plates in the front 18 inches and just sodium plates in the back 18 inches. The ZPPR staff published the measured sodium void worths and control rod worths and the statistical uncertainties in the measured values. In reality, the published uncertainty for each of these measurements is only one component of the total uncertainty. There are additional contributions from uncertainties related to measurement technique, plate compositions and reproducibility. However, the dominant uncertainty for the sodium void worth and control rod worth measurements is the uncertainty in βeff which is required to convert from the natural measurement units, cents, to the calculation units, pcm or Δ/k 1 k 2 . A full uncertainty analysis was performed for the sodium void worth and control rod worth measurements. The significant uncertainties were quantified, and a total uncertainty was derived for each sodium void worth or control rod worth measurement. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models of the sodium void worth and control rod worth configurations. These models were used to compute the sodium void worths and control rod worths. (C/E – 1) and the related uncertainty in (C/E – 1) were derived for the measured worths.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2015

Monte-Carlo Simulations of the Nuclear Energy Deposition Inside the CARMEN-1P Differential Calorimeter Irradiated into OSIRIS Reactor

The nuclear heating measurements in Material Testing Reactors (MTRs) are crucial for the study of nuclear materials and fuels under irradiation. The reference measurements of this nuclear heating are especially performed by a differential calorimeter including a graphite sample material. These measurements are then used for other experimental conditions in order to predict the nuclear heating and thermal conditions induced in the irradiation devices. Nuclear heating is a great deal of interest at the moment as the measurement of such heating is an important issue for MTRs reactors. This need is especially generated by the new Jules Horowitz Reactor (JHR), under construction at CEA/Cadarache 'French Alternative Energies and Atomic Energy Commission'. This new reactor, that will be operational in late 2019, is a new facility for the nuclear research on materials and fuels. Indeed the expected nuclear heating rate is about 20 W/g for nominal capacity of 100 MW. The present Monte Carlo calculation works belong to the IN-CORE (Instrumentation for Nuclear radiation and Calorimetry On line in Reactor): a joint research program between the CEA and Aix- Marseille University in 2009. One scientific aim of this program is to design and develop a multi-sensors device, called CARMEN, dedicated to the measurements of main physical parameters simultaneously encountered inside JHR's experimental channels (core and reflector) such as neutron fluxes, photon fluxes, temperature, and nuclear heating. A first prototype was already developed. This prototype includes two mock-ups dedicated respectively to neutronic measurements (CARMEN-1N) and to photonic measurements (CARMEN-1P) with in particular a specific differential calorimeter. Two irradiation campaigns were performed successfully in the periphery of OSIRIS reactor (a MTR located at Saclay, France) in 2012 for nuclear heating levels up to 2 W/g. First Monte Carlo calculations reduced to the graphite sample of the calorimeter were carried out. A preliminary analysis shows that the numerical results overestimate the measurements by about 20 %. A new approach has been developed in order to estimate the nuclear heating by two methods (energy deposition or KERMA) by considering the whole complete geometry of the sensor. This new approach will contribute to the interpretation of the irradiation campaign and will be useful to improve the out-of-pile calibration procedure of the sensor and its thermal response during irradiations. The aim of this paper is to present simulations made by using MCNP5 Monte-Carlo transport code (using ENDF/B-VI nuclear data library) for the nuclear heating inside the different parts of the calorimeter (head, rod and base). Calculations into two steps will be realized. We will use as an input source in the model new spectra (neutrons, prompt-photons and delayed-photons) calculated with the Monte Carlo code TRIPOLI-4{sup R} inside different experimental channels (water) located into the OSIRIS periphery and used during the CARMEN-1P irradiation campaign. We will consider Neutrons- Photons-Electrons and Photons-Electrons modes. We will begin by a brief description of the differential-calorimeter device geometry. Then the MCNP5 model used for the calculations of nuclear heating inside the calorimeter elements will be introduced. The energy deposition due to the prompt-gamma, delayed-gamma and neutrons, the neutron-activation of the device will be considered. The different components of the nuclear heating inside the different parts of the calorimeter will be detailed. Moreover, a comparison between KERMA and nuclear energy deposition estimations will be given. Finally, a comparison between this total nuclear heating Calculation and Experiment in graphite sample will be determined. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2017

Fuel pin zoning to improve the accuracy of the palisades spent fuel pool criticality analysis

Most PWR spent fuel pool criticality analyses assume that the burnup is uniformly distributed among the fuel pins in the assembly. This assumption is accurate for PWR cores that have the typical small inter-assembly gap. The Palisades reactor core, however, is different in that it uses cruciform control rods in between fuel assemblies. Every interior assembly has two sides having a wide water gap and two sides having a narrow water gap. After depleting under these conditions, the pins near the wide water gap will have a higher burnup than the pins in the center. If the burnup is distributed uniformly among the pins, the resulting k calculation will be conservative but excessively so. The excess conservatism was found to be on the order of 1% in k. For this analysis, SCALE 6.1 depletion with TRITON is used followed by the CSAS5 module for the Palisades pool analysis. Several burnup zoning arrangements (uniform enrichment) were investigated ranging from 3 zones to as many as 19 zones. The 19 zone analysis agreed well with a CASMO analysis where every pin was its own depletion zone. It was found that 8 fuel zones gave essentially the same result as 19 zones with less complexity. The depletion calculation using 8 fuel zones instead of a single zone is more complicated but the resulting k calculation no longer has the excessive conservatism seen in a one zone model. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

Dose Rate Evaluation for the ES-3100 Package with HEU Content Using MCNP, ADVANTG, Monaco, and MAVRIC

This paper presents a comparative study of dose rate calculations for the ES-3100 package with highly enriched uranium (HEU) content for different source configurations using computer codes: (1) MCNP; (2) Automated Variance Reduction Generator (ADVANTG)/MCNP; (3) Monaco; and (4) Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC). The Model ES-3100 package was developed at the Y-12 National Security Complex for domestic and international transportation of Type B fissile radioactive material. In this study, six different source configurations (i.e., solid cylinder, cylindrical hemi-shell, cylindrical shell, rectangular plate, cylindrical rod, or cylindrical segment form) having 36 kg of HEU metal inside the package containment vessel (based on configurations in the ES-3100/HEU safety analysis report for packaging) are evaluated. Dose rates at the package surfaces and 1 m from the package surfaces are calculated for these different source configurations. MCNP and Monaco cases are run without any biasing options to accelerate the convergence. Consistent Adjoint Driven Importance Sampling (CADIS) and Forward-Weighted CADIS (FW-CADIS) methods developed at the Oak Ridge National Laboratory are implemented in ADVANTG/MCNP and MAVRIC codes to accelerate the convergence. ADVANTG generates variance reduction parameters using Denovo code, and MCNP is used with the variance reduction parameters to accelerate the convergence. MAVRIC uses Denovo code to construct an importance map and a biased source distribution that are supplied to Monaco to accelerate the Monte Carlo simulation. The FW-CADIS option in ADVANTG and MAVRIC is used to accelerate the convergence in this study. The accelerated convergence cases (ADVANTG/MCNP and MAVRIC) are about 100 times faster with 100 times less particle simulation than those cases run without biasing options (analog MCNP and analog Monaco). The MCNP, ADVANTG/MCNP, Monaco, and MAVRIC calculated dose rates at the package surfaces and at 1 m from the package surfaces for the different source configurations are compared and are found to be in general agreement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗