Results for “Nuclear Fuel Performance Modeling” — Page 63

DOE OSTI2020

Risk-Informed ATF and FLEX Analysis for an Enhanced Resilient BWR Under Design-Basis and Beyond-Design-Basis Accidents

This report documents the activities performed by the Idaho National Laboratory (INL) during fiscal year (FY) 2020 for the U.S. Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program, Risk-Informed System Analysis (RISA) Pathway, Enhanced Resilient Plant (ERP) Systems research. The purpose of the RISA Pathway research and development is to support plant owner-operator decisions with the aim to improve the economics, reliability, and maintain the high levels of safety of current nuclear power plants over periods of extended plant operations. The concept of ERP refers to the combinations of accident-tolerant fuel (ATF), optimal use of diverse and flexible coping strategy (FLEX), enhancements to plant components and systems, and the incorporation of augmented or new passive cooling systems, as well as improved fuel cycle efficiency. The objective of the ERP research effort is to use the RISA methods and toolkit in industry applications, including methods development and early demonstration of technologies, in order to enhance existing reactors’ safety features (both active and passive) and to substantially reduce operating costs through risk-informed approaches to plant design modifications to the plant and their characterization. One main focus of the FY 2020 efforts documented in this report was to extend the analyses conducted in FYs 2018 and 2019 for a pressurized water reactor (PWR) to a boiling water reactor (BWR). The same analysis process, risk analysis approaches, and analysis tools as in the previous work for PWR were used for a generic BWR with near-term ATF cladding (i.e., Iron-Chromium-Aluminum [FeCrAl] cladding and Chromium [Cr]-coated cladding) designs under the postulated station blackout (SBO) and medium loss-of-coolant (MLOCA) accident scenarios. In addition, a FLEX model was developed and incorporated into a generic BWR probabilistic risk assessment (PRA) model using the INL-developed software tool, Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE), to assess the risk impact from FLEX. The other main focus of the FY 2020 efforts was to advance analysis methods, including developing dynamic approach for FLEX human reliability analysis (HRA) using the INL-developed software tool, Event Modeling Risk Assessment using Linked Diagrams (EMRALD), as well as developing a multicriterion benefit evaluation (MCBE) method for evaluating costs and benefits of safety enhancements in nuclear power plants (NPPs). As a case study, the MCBE method was applied to evaluate the costs and benefits brought by FLEX implementation.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗

DOE OSTI2021

Thermal properties of U-Mo alloys irradiated under high fission power density

A variety of physical and thermal property measurements made on irradiated U-Mo alloy monolithic fuel samples with a Zr diffusion barrier and clad in Al alloy 6061 are reported. These properties form the basis for additional investigations in thermal behavior of U-Mo fuel under high performance research reactor irradiation conditions. Measurements were performed on samples harvested from fuel plates irradiated under high fission power density and surface heat flux and complement previous measurements on samples subjected to moderate fission power density and surface heat flux. Specific heat capacity of irradiated monolithic U-Mo increases with increasing temperature and fission density, but carries large uncertainties associated with thickness measurements used to determine the mass fraction of the layer constituents for extraction of the U-Mo specific heat capacity, specifically the very thin and nonuniform Zr diffusion barrier. Density of irradiated monolithic U-Mo decreases with increasing fission density resulting from fission gas swelling and porosity formation. Thermal diffusivity of irradiated monolithic U-Mo increased with increasing temperature and decreased with increasing fission density. Thermal conductivity of irradiated monolithic U-Mo increased with increasing temperature, and in general, appeared to decrease and become less sensitive to temperature with increasing fission density. Thermal conductivity measurements compared to a semi-empirical model agreed within ±30%. It is suspected that the main driver in deviation between model calculations and measurements is the fission gas swelling determination, which does not account for fission power density or irradiation temperature.

uranium (U), molybdenum, nuclear fuel, thermal con↗

DOE OSTI2020

Total Measurement Uncertainty in Neutron Coincidence Multiplicity Analysis

Neutron multiplicity counting is the most commonly used nondestructive assay technique for determining the plutonium mass within containers of scrap PuO 2 or mixed oxide (MOX). In multiplicity analysis, the 240 Pu eff mass, leakage multiplication, and alpha ratio (the ratio of [α, n]-to-spontaneous fission neutron production) are the three primary unknown sample properties. They must be determined simultaneously. To solve for these three unknowns in a multiplicity assay, three measured values are needed: the singles, doubles, and triples neutron count rates. While the analysis is limited to solving for three unknowns, there are many additional factors that impact the observed count rates and contribute to the measurement uncertainty. In this study we investigate the various uncertainty contributors for the multiplicity analysis through a combination of traditional uncertainty propagation techniques supplemented by Monte Carlo simulations to address the dependences not explicitly expressed by the point source model. Uncertainties arising from counting statistics, calibration parameters, calibration method, nuclear data, and various material characteristics (isotopic abundances, chemical form, density, and impurities) are considered. A Total Measurement Uncertainty (TMU) estimate is then developed from these uncertainty contributors. This study is confined to multiplicity analysis of items commonly encountered in international safeguards applications. That is, the study focused on Pu oxides and MOX materials for the masses ranging up to 4000 grams total Pu. Multiplicity measurements were simulated using MCNP V6 based on the Plutonium Scrap Multiplicity Counter (PSMC), Epithermal Multiplicity Counter (ENMC), Pyrochemical Multiplicity Counter, and Large Epithermal Multiplicity Counter (LEMC) for this study; however, this report focuses on the parameterization of the uncertainties for the PSMC. The performance differences between the PSMC and the other multiplicity counting systems are relatively small, primarily manifesting in the impact on measurement precision so that the evaluation developed for the PSMC can be applied to the other multiplicity counting systems. Finally an analysis tool, the Multiplicity TMU Estimator, was developed from this study to serve as an aid for evaluation of the total measurement uncertainty of multiplicity assay results obtained from the commonly used INCC acquisition and analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

STAT7 v2.0 User Guide

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software includes non-power research and test reactor analysis, and it has been used for the conversion to low- enriched uranium fuel of U.S. High Performance Research Reactors such as the Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat reactor safety analysis, such as during fuel reloading, STAT7 accommodates flexibility in analyzing many practical aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for reactor thermal hydraulics analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to fitted functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze every axial node in each lateral stripe of all fuel plates and coolant channels in every fuel element of an entire reactor core. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. In a single execution, STAT7 can be used to perform full core analysis by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Dual-Purpose Canister Filling Demonstration Project Progress Report

This report discusses the progress made at the Oak Ridge National Laboratory to support direct disposal of dual-purpose canisters (DPCs). Acceptable filler materials should demonstrate that the probability of criticality in DPCs during the disposal timeframe is below the probability threshold for inclusion in a repository performance assessment. This effort, which will ultimately result in a full-scale demonstration, includes computational fluid dynamics (CFD) modelling developed to gauge the filling process and to uncover any unforeseen issues. Filling simulations of the lower region (mouse holes) of a prototypic DPC show successfully simulate filling of the void space inside the canister and a smooth, even progression of the liquid level. Flow through a pipe that is similar to the drainpipe in a DPC is being investigated separately to gain valuable insight of the flow regime inside a pipe. Three physical experiments validating the computational filling model have been completed using surrogate liquids. One experiment using molten metal has been completed and the results demonstrate adequate filling of the canister void spaces. Although the scale experiment observed some void spaces related to shrinkage of the metal during cooling, the volume filled is expected to be sufficient to meet the purpose of moderator exclusion. Further experiments will scale up the geometry, moving towards a full-sized DPC demonstration, and provide a more rigorous investigation of candidate filler materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Benchmark Specification for FFTF LOFWOS Test #13

The Fast Flux Test Facility (FFTF) at the Hanford site in Washington was designed by the Westinghouse Electric Corporation for the U.S. Department of Energy. FFTF was a 400 MW thermal, oxide-fueled, liquid sodium cooled test reactor, built to assist development and testing of advanced fuels and materials for fast breeder reactors. After reaching criticality in 1980, FFTF operated until 1992, providing the U.S. Department of Energy (DOE) with the means to test fuels, materials, and other components in a fast neutron flux environment. In July 1986, a series of unprotected transients (with the plant protection system intentionally disabled) were performed in FFTF as part of the passive safety demonstration program. Among these were thirteen loss of flow without scram (LOFWOS) tests. The goals of this program included confirming the liquid metal reactor safety margins, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of specific design features. The test defined in this benchmark is LOFWOS Test #13, which was initiated at 50% power and 100% flow with the pump pony motors turned off. This benchmark specification is intended to support collaborative efforts within international partnerships on the validation of simulation tools and models in the area of Sodium-cooled Fast Reactor (SFR) safety. Validated tools and models are needed to evaluate SFR inherent safety characteristics and assess the impact of passive design features in response to accident initiators. Comparisons with experimental data and the results of safety analyses from other groups create unique opportunities to improve predictive capabilities of computational codes and methods for SFR modeling and simulation. The conditions of the LOFWOS test along with the feedback from FFTF’s limited free bow core restraint system and the novel passive reactivity control Gas-Expansion Modules (GEMs) pose a very challenging and uniquely valuable benchmark exercise.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Initial study on cross section generation requirements for a PBR equilibrium core

A Serpent model of the HTR-PM equilibrium core was developed for use in cross section prepa-ration studies in order to guide methods development for the Griffin reactor multiphysics applica-tion. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zoneswith unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluidsequilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and mod-erator temperatures for various core regions was performed with the MOOSE stochastic tools.The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluidtemperature and that the power level, heat transfer coefficient and effective conduction to neigh-boring pebbles and fluid constitute, at best, second order effects. The temperature uncertaintyrange varies from 28 K to 57 K at the core entry and exit planes, respectively, but these val-ues are probably higher. We still have to quantify the significance of these uncertainties in thepreparation of cross sections, which will be postponed for future work. In addition, we verifythat the single effective pebble approach works well for the preparation of region averaged crosssections in the infinite domain approximation. Nevertheless, there are significant discrepanciesin the cross sections when compared to the multi-pebble model. This could affect the predictionof peak values and in the depletion calculation. We conclude that is highly desirable for futurestudies with Griffin to be able to handle both the ?effective? pebble approximation and the multi-pebble approach for the various pebble burnup groups. This enables Griffin users the flexibilityto perform higher-fidelity studies. Finally, we initiate the preparation of cross sections for variouscore regions from the full core Serpent reference model. We quantify the differences in 26 groupcross sections from infinite domain models. These reference cross sections will serve to validatethe double heterogeneity, self-shielding, and spectrum-correction methods in Griffin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Supercritical CO2 Recuperators, Presented at The ASME Turbo Expo 2017, June 29, 2017, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

DOE OSTI2020

High Effectiveness, Compact, High Pressure and Low-Cost Recuperator, Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Phonon-Spin Scattering from Unpaired f-electrons in U atoms and It’s Influence on Thermal Transport in Uranium-doped Thorium Dioxide Single Crystals

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO2) is investigated. ThO2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO2 system.

36 - MATERIALS SCIENCE↗

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 OSTI2022

Assessment of thermodynamic stability of sapphire in eutectic molten chloride environment

The continued development of molten salt reactors requires highly efficient and stable instruments that can efficiently track monitor conditions of the molten salt during long-term application operation in both the fuel and coolant/heat transfer fluid loops. Sapphire (Al 2 O 3 ) fibers have shown tremendous potential by enabling radiation resistance and a broader operational range of temperature . Thermodynamic modeling using the ThermoCalc and the SGTE Molten Salts (SALT1) and Pure Substances (Pure5) database is done performed to understand the compatibility of Al 2 O 3 fibers with NaCl-MgCl 2 eutectic molten salt in the temperature range of 600 to 1000°C. The thermodynamic calculations show that sapphire fibers are not expected to be compatible with long term exposure to molten chloride salts within this temperature range. Furthermore, various pure metallic elements were evaluated as potential cladding materials for part of the Al 2 O 3 fibers. Based on the thermodynamic analysis, molybdenum (Mo) and nickel (Ni) could be effective cladding materials to enhance the stability of Al 2 O 3 in NaCl and MgCl 2 chloride salt molten bath in the desired temperature range. The presence of a metallic cladding leads to the suppression of aluminum dissolution from Al 2 O 3 , thus improving its stability in chloride salts. Additionally, it is also shown that Al 2 O 3 remains stable up to 2300 K in the presence of pre-existing Al 2 MgO 4 and Al 2 NiO 4 in the eutectic molten chloride bath environment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Idaho National Laboratory CY 2021 National Emission Standards for Hazardous Air Pollutants Analysis, Methodology and Results for Radionuclides

This report details calculations of potential dose at public receptor locations surrounding the Idaho National Laboratory (INL) Site boundary, and INL in-town facilities, from radionuclides reported to be in use and potentially emitted from INL facilities during calendar year (CY) 2021. All calculations were performed in accordance with the requirements in Code of Federal Regulations (CFR), Title 40, “Protection of the Environment,” Part 61, “National Emission Standards for Hazardous Air Pollutants (NESHAPs),” Subpart H, “National Emission Standards for Emissions of Radionuclides Other than Radon from Department of Energy Facilities” (40 CFR 61, Subpart H). Modeling methodology, model input parameters, and contribution to dose by facility, source, and radionuclide at the maximally exposed individual (MEI) location are also discussed. The information in this report supports the “National Emission Standards for Hazardous Air Pollutants - Calendar Year 2021 INL Report for Radionuclides” (INL 2022). In CY 2021, the estimated annual potential dose at the INL Site MEI location was 6.67E-02 mrem/yr, up slightly from the previous year, but far less than the regulatory standard of 10 mrem/yr (CFR 40 Part 61, Subpart H). Approximately 97% of the total dose to the INL Site MEI originated from MFC sources. Emissions from INL in-town facilities resulted in an estimated annual potential dose of 6.21E-03 mrem/yr to the MEI, down 40% from the CY 2020 estimated dose. Year-to-year variations in estimated annual dose can be attributed to adjustments in laboratory operations, changes to facility infrastructure, and variation in meteorological conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Evolution of dislocation loops and voids in post-irradiation annealed ThO 2 : A combined in-situ TEM and cluster dynamics investigation

The effect of isochronal annealing on the evolution of dislocation loop and void population in proton irradiated ThO 2 has been investigated. Post-irradiation annealing in other actinide oxides like UO 2 shows significant loop coarsening. ThO 2 samples were irradiated with 2 MeV protons up to a dose of 0.1 dpa at 600°C. Post-irradiation isochronal annealing was performed at 600, 800, 1000 and 1100°C for 1 h at each temperature using in-situ TEM. Only faulted 1/3<111> type dislocation loops were observed, and their sizes and distribution were characterized. Further, the population of self-interstitial atom (SIA) dislocation loops did not show any significant growth and coarsening. Additionally, nanometric voids were observed at annealing temperatures of 1000 and 1100°C. Using cluster dynamics (CD), we have studied the nucleation and growth of point defects and defect clusters, i.e., SIA prismatic dislocation loops and nanometric and sub-nanometric voids in proton irradiated ThO 2 . The CD model was further utilized to predict the growth and coarsening of loops and voids during isochronal annealing at the experimental and higher temperatures. The model did not predict significant SIA loop growth which closely corresponds to the TEM observations. CD predicted SIA loop coarsening is insignificant even at high annealing temperature of 1500°C because the model only considers the growth of defect clusters by absorption of like point defects, i.e., SIA loops absorb interstitials and voids absorb vacancies, and cannot account for their migration and coalescence due to elastic interaction. The CD model also predicts the evolution of nanometric voids having mean size within the error bounds of TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

High Effectiveness, Compact, High Pressure and Low Cost Recuperator for Super-Critical CO2 Power Cycles, Paper Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗