SEARCH · Engineering Papers

Results for “Nuclear fuel performance analysis”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Methodology for Producing a Uniform Distribution of UO2 in a Tungsten Matrix

Current work at NASA's Marshall Space Flight Center (MSFC) is focused on the development CERMET fuel materials for Nuclear Thermal Propulsion (NTP). The CERMETs consist of uranium dioxide (UO2) fuel particles embedded in a tungsten (W) metal matrix. Initial testing of W-UO2 samples fabricated from fine angular powders performed reasonably well, but suffered from significant fuel loss during repeated thermal cycling due to agglomeration of the UO2 (1). The blended powder mixtures resulted in a non-uniform dispersion of the UO2 particles in the tungsten matrix, which allows rapid vaporization of the interconnected UO2 from the sample edges into the bulk material. Also, the angular powders create areas of stress concentrations due to thermal expansion mismatch, which eventually cracks the tungsten matrix. Evenly coating spherical UO2 particles with chemical vapor deposited (CVD) tungsten prior to consolidation was previously demonstrated to provide improved performance. However, the CVD processing technology is expensive and not currently available. In order to reduce cost and enhance performance, a powder coating process has been developed at MSFC to produce a uniform distribution of the spherical UO2 particles in a tungsten matrix. The method involves utilization of a polyethylene binder during mixing which leads to fine tungsten powders clinging to the larger UO2 spherical particles. This process was developed using HfO2 as a surrogate for UO2. Enough powder was mixed to make 8 discs (2cm diameter x 8mm thickness) using spark plasma sintering. A uniaxial pressure of 50 MPa was used at four different temperatures (2 samples at each temperature). The first two samples were heated to 1400C and 1500C respectively for 5 minutes. Densities for these samples were less than 85% of theoretical, so the time at temperature was increased to 20 minutes for the remaining samples. The highest densities were achieved for the two samples sintered at 1700C (approx. 92% of theoretical). Scanning electron microscopy (SEM) of the mixed powders and the sintered samples along with energy dispersive x-ray analysis was obtained. The SEM of the powders clearly show the fine W powder adhered to the larger HfO2 particles and a uniform distribution of HfO2 particles in a tungsten matrix upon densification. Vicker's Microhardness testing was also performed on all samples using 0.5, 1.0 and 2.0 kg loads. Five indents were made at each load level. All indents were placed in the tungsten matrix to assist as a proxy in measuring densification. The highest hardness value was obtained for the 1700C specimens. The hardness average for these samples was 312.14 MPa. This powder processing method has been applied to W/UO2 powders with the SEM of the powders appearing similar to the W/HfO2 powder images.

Tucker, Dennis S.↗

DOE OSTI2025

Deployment of neural-network-based neutron microscopic cross sections in the Griffin reactor physics application

The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Involute Working Group – Validation of CFD Turbulence Models for Steady-State Safety Analysis

There are three research reactors in the world having fuel plates curved as circle-involute (a spiral generated around a circle): The Oak Ridge National Laboratory (ORNL) High Flux Isotope Reactor (HFIR) located in Tennessee, U.S.A. [ORNL, 2018]; The Laue-Langevin Institute (ILL) High Flux Reactor (RHF) located in Grenoble, France [ILL, 2018]; and, The Technical University of Munich (TUM) Research Neutron Source Heinz Maier-Leibnitz (FRM II) located in Garching, Germany [TUM, 2018]. All three reactors are currently using Highly Enriched Uranium (HEU, 235U/U ≥ 20 wt. %) as fuel and all three are actively engaged in activities to convert to Low Enriched Uranium (LEU, 235U/U < 20 wt. %) fuel. More information on these reactors can be found in Section two. For various reasons, these reactors have expressed interest in using Computational Fluid Dynamics (CFD) tools to perform their Steady-State Thermal-Hydraulic (SSTH) safety calculations. Using CFD tools for SSTH analysis represents generally a significant departure from traditional methods and requires thorough verification and validation (V&V) to be accepted by regulators. This is why Argonne National Laboratory (ANL) and the involute reactors formed an informal group to help each other in this endeavor. This so-called Involute Working Group (IWG) aims at supporting the qualification of CFD tools for SSTH safety calculations. Activities include benchmarking, code-to-code comparison, Verification and Validation (V&V) as well as technical support. More information on the IWG and the reasons why involute-plate reactors are interested in CFD can be found in Section three. The present study aims at testing three turbulence models commonly implemented in commercial software (namely K-Epsilon, K-Omega and Reynolds-Stress) and assessing their relevance for involute-plate reactor application. An experiment performed by Gambill and Bundy in 1961 was performed to measure heat transfer coefficients in conditions prototypic of the HFIR reactor [Gambill, 1961]. More details on this experiment are provided in Section four. Based on the information available, CFD models of this experiment have been created and the experimental results have been compared to the ones obtained with these CFD models. CFD models and the results of this comparison are provided in Section five.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Source term analysis of FeCrAl accident tolerant fuel using MELCOR

It has been established that incremental improvements in beyond design basis accident performance can be achieved through accident tolerant fuel (ATF). However, they have the potential to recover margin with respect to conventional fuel and therefore enhance plant economics through uprate or cycle length increase. To realize this potential, it is necessary to quantify the reduction in source term due to use of ATF, and correspondingly how this is affected by increasing the cycle length and/or burnup. This requires development of a risk informed analysis methodology for ATF under high burnup conditions, which is being developed within LWRS. To this end, a MELCOR simulation of FeCrAl ATF in a recovered Large Break LOCA (LBLOCA) scenario has been developed, using a model based on the Zion Pressurized Water Reactor (PWR). The new user defined material capability, along with the inclusion of detailed neutronics- and depletion-derived parameters such as core power distribution, decay heat behavior, and fission product inventories, allows a more detailed simulation of FeCrAl clad material properties and behavior than has previously been possible using MELCOR. These detailed FeCrAl results were compared with a zircaloy clad model to investigate the differences between the two clad types and quantify the benefits of FeCrAl ATF with an 18-month cycle. Next, the fuel cycle with ATF was extended to 24 months, to determine whether any additional safety margin provided by FeCrAl ATF could be leveraged to implement high burnup FeCrAl-clad fuel while retaining the same operating and safety limits as current zircaloy-clad fuel. For the particular scenario analyzed, the delay to fuel failure from using ATF was of the same order as the LOCA recovery time, and hence significant in reducing fission product release, with the 24-month FeCrAl core performing better than the 18-month zircaloy core. It is noted that for other transients, the reduction in release due to using FeCrAl may be less significant. Furthermore, the material model developed here can be used in such further studies in support of determining the overall source term.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Synthesis and Characterization of Pd-based Nanomaterials

Bimetallic nanoparticles (BNPs) consist of two different types of metals or alloys that are bonded together. Unique properties such as optical, electronic, thermal, and catalytic effects differ for each type of BNP. Important BNPs range from Au-Pd, Ag- Pt, Au-Pt, and Ag-Ni. Pd bimetallic nanoparticles are of interest due to their many applications such as catalysis and sensing. Bimetallic catalysts have increase reaction rates and have improved catalyst stability through the geometry and ligand distribution. Pd nanoparticles are considered to be a strong catalyst due to their high activity at low temperatures and high tolerance to moisture. The catalytic properties of bimetallic nanoparticles depend on the structural properties such as size and shape. Core-shell, hollow structure, and multi-shell alloy are three possible structures nanoparticles can form as bimetallic catalysts. BNPs can be synthesized through different methods to control the size, shape, and structure. To obtain different morphologies, a variety of methods can be performed. Different methods can range from the usage of the glancing angle deposition (GLAD) to the galvanic replacement reaction, but the methods all depend on the properties of the metals. The galvanic displacement reaction was the method used to obtain Pd-based nanoparticles. This reaction is best know for obtaining hollow shaped NPs. To determine what redox process was preformed, the activity series of metals was used. From the activity series of metals, silver (Ag) was selected to preform Pd-based nanoparticles. Objectives: Synthesize Ag nanoparticles and Ag-Pd nanoparticles to understand the morphology. Characterize the synthesized nanoparticles using scanning electron microscopy (SEM), phase analysis light scattering (PALS), dynamic light scattering (DLS), energy dispersive X-ray spectroscopy (EDS), and UV-Vis spectroscopy. Results: In the UV-Vis spectrum, the Ag-Pd bimetallic NP's plasmon band decreased as the volume of palladium increased. The surface charge increases as the concentration of palladium increases. The Pd{sup 2+} ions interact with the sodium citrate surface, and decrease the negative charge. Conclusion: Ag-Pd nanoparticles were successfully created and stabilized with sodium citrate. The addition of Pd decreased the prominent plasmon band of the Ag nanoparticles. The SEM analysis showed that Ag nanoparticles had a well-defined structure, while the Ag-Pd nanoparticles showed hollow and rough structure. The EDX analysis confirmed the presence of silver and palladium. This material can be used in many industrial and research fields such as organic synthesis, fuel cells, and environmental sensing and remediation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 2: 2D and 3D core simulations

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. The benchmark consisted of a progression in complexity from modeling single 2D pin cells, through 2D lattices and finally to whole core 2D and then 3D scenarios. The accompanying paper (Part 1) presented the results for a variety of pin cell and lattices representative of high energy cores and established the ability of the MPACT code to accurately match the predictions of the Monte-Carlo code SERPENT2 such that the former may be used as a surrogate reference code in spatially complex whole core modeling applications where using Monte-Carlo methods is somewhat impractical. This paper focuses on comparisons at the whole core level using a representative high power density 3-loop Westinghouse NSSS core operating for 24 month cycles employing high enriched fuel with use of combined burnable absorbers to hold down excess reactivity. The evolution of the core reactivity and power distribution characteristics throughout the cycle are calculated using the Westinghouse NEXUS-based code suite and compared to the results using MPACT. The results of this whole core benchmark confirm the conclusions that the Westinghouse code suite accurately models HB/HE fuel designs, with very good prediction of core reactivity and power distribution throughout the cycle. The industry can be confident that the NEXUS code system can be deployed to support PWR reload safety and operational analyses for cores featuring HB/HE fuel and 24-month cycle operation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 MARK-I Station Blackout: Part I – Single Case Analysis

Here in this analysis, the two material interaction models available in the MELCOR code are benchmarked for a severe accident at a BWR under representative Fukushima Daiichi boundary conditions. This part of the benchmark investigates the impact of each material interaction model on accident progression through a detailed single case analysis. It is found that the eutectics model simulation exhibits more rapid accident progression for the duration of the accident. The slower accident progression exhibited by the interactive materials model simulation, however, allows for a greater degree of core material oxidation and hydrogen generation to occur, as well as elevated core temperatures during the ex-vessel accident phase. The eutectics model simulation exhibits more significant degradation of core components during the late in-vessel accident phase – more debris forms and relocates to the lower plenum before lower head failure. The larger debris bed observed in the eutectics model simulation also reaches higher temperatures, presenting a more significant thermal challenge to the lower head until its failure. At the end of the simulated accident scenario, however, core damage is comparable between both simulations due to significant core degradation that occurs during the ex-vessel phase in the interactive materials model simulation. A key difference between the two models’ performance is the maximum temperatures that can be reached in the core and therefore the maximum ΔT between any two components. When implementing the interactive materials model, users have the option to modify the liquefaction temperature of the ZrO 2 -interactive and UO 2 -interactive materials as a way to mimic early fuel rod failure due to material interactions. Through modification of the liquefaction of high melting point materials with significant mass, users may inadvertently limit maximum core temperatures for fuel, cladding, and debris components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Multi-cycle reload analysis of a long cycle gas-cooled fast modular reactor

There is currently significant interest in deploying HALEU-fueled fast reactors, including the General Atomics (GA) Fast Modular Reactor (FMR). Such reactors can achieve very long fuel cycles, but with multi-batch loading will take decades to reach equilibrium. This motivates design and analysis of both the initial core and multi-cycle reload, which is typically performed using fast-running, deterministic fast reactor codes such as the Argonne Reactor Computation (ARC) codes. In this paper, multicycle reload of the GA FMR is analyzed using the ARC codes. The GA FMR utilizes 19.75 % enriched fuel in a 16 year cycle with a three-batch strategy, with twice-burned fuel placed on the core periphery. The GA FMR has a softened neutron spectrum due to reflecting elements in the core, so the neutronic solution is first benchmarked against the OpenMC Monte Carlo code. Discrepancy on k eff is 400–600 pcm, likely due to the softened neutron spectrum, heterogeneous fuel assembly design and central reflector. However, the rms discrepancy on the assembly power distribution is only 0.6 %, despite the presence of the central reflector. A reload strategy is devised for the first three cycles of such a reactor, ultimately spanning the first 45–48 years of its operation. The fresh core uses 19.75 %, 19.25 % and 16.75 % enriched fuel in place of fresh, once-burned and twice-burned and is then subsequently refueled with only 19.75 % enriched fuel. The cycle length is varied over 3 cycles of operation to balance fuel utilization and reactor availability, specifically with use of an extended 18-year Cycle 1, followed by a shortened 11-year Cycle 2. Cycle 3 is close to the target 16-year length. Finally, placing twice burned assemblies next to the GA FMR central reflector can reduce power peaking by 3 %, at the expense of slightly reducing the cycle length.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of MCNP, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) confirm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results analysis methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exemplified through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Sulfur Pellets Responses to a Bare and Steel Reflected Pulse of the Oak Ridge National Laboratory Health Physics Research Reactor

The experiments analyzed in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the $\textit{Fast Burst Reactor}$. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. The HPRR was an unmoderated, unshielded fast reactor that used highly enriched uranium and molybdenum alloy as fuel. The reactor was initially sent to the Nevada Test Site in 1962, where it was used to evaluate radiation doses received as a result of the Hiroshima and Nagasaki bombings during World War II. A few years later, the reactor was sent back to ORNL to be part of the Dosimetry Application Research (DOSAR) facility shown in Figure 1, which included a reactor building shown on the left (west) of the picture and a control and laboratory building in the upper right corner (northeast). The critical assembly was used for numerous technical studies, including systems calibration, dosimetry, radiobiology of plants and animals, testing of radiation alarms, as well as teaching and training in radiation dosimetry and nuclear engineering. Between 1963 and 1987, the HPRR was operated for thousands of hours, achieving criticality close to 10,000 times and motivating many publications. The HPRR was decommissioned in 1987. The goal of this effort was to use historical data from operation of the HPRR to create a criticality accident alarm system (CAAS) benchmark to be included in the $\textit{International Handbook of Evaluated Criticality Safety Benchmark Experiments}$ (ICSBEP Handbook). A thorough inspection was performed of all available documentation and information available. The most promising experiments that were selected for evaluation were those described in the 1987 ORNL report entitled $\textit{Health Physics Research Reactor Reference Dosimetry}$, ORNL-6240. The report includes reference dosimetry results of the shielded and unshielded configurations of the HPRR after burst operations. Because of changes to the reactor positioning and storage systems that were made in 1985, the previous dosimetry reports became obsolete, and the newly designed experiments were needed to create the HPRR’s adjusted dosimetry data. The various results reported in ORNL-6240 include reference doses and dose equivalents from different conventions at different distances and elevations as determined using the detected neutron fluence and conversion factors. The HPRR neutron fluence was obtained through different methods, including sulfur pellet analysis and threshold detector unit data. Information about the HPRR spectrum was also obtained through Bonner sphere measurements. This benchmark is focused on a part of the measured sulfur fluences reported in Appendix H of ORNL-6240. Standard commercial sulfur pellets were placed at different distances from the HPRR centerline during burst operation and were activated due to the 32 S(n,p) 32 P reaction. The resulting 32 P activity was then measured and the information about the corresponding sulfur fluence and/or neutron dose could be extracted. Many of those measurements have 7 been performed with the HPRR in its bare configuration or with different shields (combinations of Lucite, concrete, steel). All the necessary, precise information about material and/or dimensions of the different shields was not found, so it was decided to focus only on the unshielded and steel-shielded configurations to minimize the benchmark uncertainty. A total of 31 cases (24 unshielded and 7 shielded cases at different positions) of sulfur fluence were selected before evaluation to develop the benchmark.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

High-temperature steam oxidation study of irradiated FeCrAl defueled specimens

Post irradiation examinations (PIE) were performed on irradiated iron-chromium-aluminum (FeCrAl) specimens. These FeCrAl specimens were fabricated at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The experimental setup involved subjecting FeCrAl cladding, along with UO 2 pellets, to irradiation in the Idaho National Laboratory Advanced Test Reactor (ATR). In parallel, the FeCrAl alloy tubing without UO 2 pellets was irradiated at ORNL's High Flux Isotope Reactor (HFIR). After irradiation, the ATR-irradiated rodlet was transported to an ORNL hot cell, where it was sectioned into multiple samples for the PIE and severe-accident testing. The sectioning process revealed that the fuel was not bonded to the cladding and could be easily detached from sectioned cladding slices. Microstructural analysis of the fuel cross sections demonstrated no significant interaction between the fuel and the cladding. Additionally, high-temperature steam oxidation tests on defueled cladding segments showed minimal oxygen uptake even at 1200 °C. Here, the ATR-irradiated specimens began to exhibit signs of enhanced oxidation upon reaching a temperature of 1300 °C. Furthermore, enhanced oxidation was observed on the inner surface of the ATR-irradiated FeCrAl specimen, which had been subjected to 1300 °C for a duration of 1 min. By contrast, high-temperature steam oxidation experiments indicated that the HFIR-irradiated FeCrAl cladding provided good thermal stability when exposed to 1300 °C for up to 4 h. Comparative analysis encompassing the oxidation behavior of the ATR-irradiated fueled FeCrAl, HFIR-irradiated unfueled FeCrAl, and unirradiated FeCrAl suggests that the fuel–cladding interaction, although not visible via standard microscale electron microscopy measurements, may accelerate the deterioration of FeCrAl cladding in beyond-design-basis accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

AGC-4 Experiment Irradiation Monitoring Data Qualification Final Report

The Graphite Technology Development Program ran a series of six experiments to quantify the effects of irradiation on nuclear grade graphite. This report focuses on the fourth experiment, Advanced Graphite Creep 4 (AGC 4). The Advanced Reactor Development (ARD) Technology Development Office (TDO) Program for research and development activities require documentation of qualified monitoring data to design and license the first high-temperature reactor nuclear plant. Qualified data meets the requirements for use as described in the experiment planning and quality assurance documents. Failed data do not meet the requirements and provide no useable information. Trend data may not meet all requirements but still provide some useable information. Use of Trend data requires assessment of how any deficiencies affect a particular use of the data. AGC-4 began with Advanced Test Reactor (ATR) Cycle 157D on May 30, 2015. After irradiating the graphite for two cycles, the capsule was removed from the reactor after ATR Cycle 158A, which ended on January 2, 2016, due to interference with another experiment. Irradiation was resumed with Cycle 162A on October 7, 2017 after the interfering experiment was removed from the reactor. AGC-4 irradiation completed by the end of Cycle 166B on January 10, 2020. Between Cycles 162A and 166B, AGC-4 capsule was removed from the reactor core during two Powered Axial Locator Mechanism (PALM) Cycles 163A and 165A to avoid overheating. All thermocouples (TCs) have functioned throughout the AGC-4 experiment. A total of 9,256,184 out of all 11,167,213 temperature records (or 82.9%) are Qualified for use by the ARD Program and 1,911,029 missing values are Failed. Argon, helium, and total gas flow data were within expected ranges, except only eight out-of-range values occurred during outages. A total of 21,240,627 out of all 22,288,5103 gas flow rates (or 92.8%) are Qualified for use by the ARD Program and 1,644,476 are Failed records mostly due to missing values. Discharge gas line moisture values were consistently low during full ATR power, except for 162B, when moisture content was increasing to more than 200 ppmv by the end of this cycle. During outages, moisture content reached as high as 700 ppmv. 1,089,941 out of a total of 1,107,331 moisture values (or 98.4%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 17,390 missing moisture values are Failed. Graphite creep specimens were subjected to one of three loads: 393, 491, or 589 lbf. For a brief period during Cycle 157D between 12:19 on June 2, 2015, and 08:23 on June 11, 2015, the load cells were wired incorrectly, resulting in missing stack load data. Missing stack loads were estimated from measured ram pressures using regression equations developed from the existing data from Cycle 157D. Estimated stack loads during this period are considered to be an accurate representation of actual load applied to the stacks. These loads deviate slightly from the planned loads. This deviation does not prevent the data from being Qualified for use but must be taken into account when analyzing the effect of load on creep. 6,095,728 out of a total of 6,403,985 moisture values (or 95.2%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 308,257 missing load values are Failed. Stack displacement increased consistently throughout the eight cycles, with total displacement reached highest value of 2.4 in by the end of irradiation. During ATR outages, a set of pneumatic rams raised the stacks of graphite creep specimens to ensure the specimens were not stuck within the test train. This stack raising was performed seven times throughout irradiation. All stacks were raised successfully each time. 4,744,974 out of a total of 6,094,513 displacement values (or 77.9%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 1,349,539 displacement values are Failed mostly due to missing values. Analyses were conducted on correlations between TCs to look for trends and step changes that might indicate instrument degradation or failure. Correlation analysis was used to identify instances when TCs form short circuits, referred to as virtual junctions, which result in TCs reporting temperatures from some location in the capsule other than the location where they were intended to read. No evidence of virtual junctions was found. Analyses were also conducted on control charts of temperature differences of between two TCs, which are expected to behave consistently throughout the entire irradiation period. Upward or downward trend over time indicates at least one TC in the pair was drifted. Examining control charts for all 66 possible pairs out of twelve TCs installed in the AGC-4 capsules reveal no clear drift failures occurred, except unstable behavior of two TCs in Zone 3, TC-7 and TC-8, over irradiation time. In conclus

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Analysis of the MUSIC 3 He Multiplicity Data

A measurement campaign called the Measurement of Uranium Subcritical and Critical (MUSiC) was performed on a range of configurations of highly-enriched uranium (HEU) from December 2020 through April of 2021. While part of the focus was to measure reactor kinetics parameters on delayed supercritical systems, an additional focus was performing neutron noise measurements on subcritical configurations from deeply subcritical to nearly delayed critical. Multiple detector systems were used to perform these measurements, such as a 3 He multiplicity detector called the Neutron Multiplicity Array Detector (NoMAD) and a liquid scintillator system called the Rossi-α Measurement Rapid Organic Discriminating Detector (RAM-RODD). Also included were a scintillator system from the University of Michigan and a set of small 3He tubes that have previously been used to measure Rossi-α values on near-critical systems. The focus of this paper will be a comparison of prospective analysis methods for the NoMAD measurements. Previous subcritical measurements at the National Criticality Experiments Research Center (NCERC) submitted to the International Criticality Safety Benchmark Evaluation Project (ICSBEP) used the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This relies on the time correlations of neutron detections to infer the spontaneous fission rate and neutron multiplication of a system through binning the time tagged detections and analyzing resulting histograms of the numbers of counts. However, there are other neutron noise methods that rely on similar processes, such as the Hansen-Dowdy formalism which uses a slightly different methodology to extract multiplication from the neutron multiplicity counting moments. Comparisons can be made between these experimental results and those obtained through simulations to validate or identify deficiencies in analysis, detection methods, or the underlying nuclear data. Different time gating strategies and their effects on count rate uncertainties are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2023

Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer

With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.

Microgravity↗

DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Application of BISON to UO 2 MiniFuel fission gas release analysis

There has been a recent push to accelerate fuel qualification by developing revolutionary capabilities to reduce irradiation periods, and thereby, reduce the time required to qualify a new fuel system. One such capability is the MiniFuel irradiation capsule designed to miniaturize fuel samples and irradiate “mini” fuel samples under isothermal temperature conditions. MiniFuel allows steady-state irradiations to decouple the traditionally coupled fission rate (i.e., power) and temperature parameters to understand and generate microstructures observed in fuel operated in a commercial reactor. Furthermore, this process offers the possibility to gather in situ data as well as postirradiation or transient data such as thermal conductivity, specific heat, fission gas diffusion and release, etc. However, accelerating fuel qualification is not solely reliant on generating large amounts of data but also on developing an informed test matrix designed to rapidly generate impactful data. Additionally, this process is reliant on fuel performance codes, such as BISON, to evaluate MiniFuel irradiations using existing material models. This process pinpoints model/data gaps, identifies desired irradiation conditions, and subsequently supports model validation and development. This work describes the use of BISON to perform a number of sensitivity studies designed to understand conditions that lead to fission gas release (FGR) under steady-state isothermal irradiation conditions and temperature transient conditions. The model is applied to a UO 2 MiniFuel example and shows an overall good qualitative agreement with experimental FGR annealing tests under different temperature conditions. It also accounts well for microstructural effects on FGR. When quantitatively compared with FGR data from previously irradiated 103 MWd/kgU UO 2 discs under thermal annealing, the model shows a less satisfactory agreement with the experimental data. Finally, a UO 2 MiniFuel test matrix is proposed to help to extend the model's operational range and validate the new FGR model capabilities to higher burnups and transient conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗