Search NASA⌕ Search

SEARCH · Search NASA

Results for “HTGR Modeling Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurement of Fission Product Concentration Profiles in AGR-3/4 TRISO Fuel Graphitic Matrix and Nuclear Graphites

The third Advanced Gas Reactor (AGR) irradiation experiment, AGR-3/4, was designed to investigate the migration of fission products in fuel compact graphitic matrix and reactor graphite components. Using destructive methods, radial fission product concentration profiles were measured for gamma-emitting fission products (e.g., Ag-110m, Cs-134, and Eu-154) and beta-emitting Sr-90 in irradiated graphitic and graphite components from six different AGR-3/4 irradiation capsules. These new measured concentration profiles can now be compared to non-destructive measurements and fission product transport simulations and will be used to derive new diffusivities and sorptivities to support refinement of fission product transport models and high-temperature gas-cooled reactor (HTGR) source-term analyses. Each capsule in the AGR-3/4 experiment had four fuel compacts in the middle of two concentric rings of graphitic matrix material, PCEA graphite, or IG-110 graphite. In addition to the approximately 1898 tristructural isotropic (TRISO) coated particles in each compact, there were 20 designed to fail (DTF) particles coated only in pyrocarbon so that they released fission products into the surrounding cylindrical rings of carbonaceous materials. Destructive sampling of the rings involved machining/milling material from around the circumference of the rings, collecting that material, and performing radiochemical analyses on it. Milling operations were performed in multiple steps or segments, and each segment was generally 0.508 mm (0.020 in) thick. Knowing the radial position at which each segment was milled, the volume of the milled material at each segment, and the fission product content in each segment, the radial fission product concentrations were constructed for select isotopes in each ring. Ag-110m profiles had the most variation. Some profiles were peaked at an inner or outer surface. Some were peaked at the middle of the ring wall thickness. Some increased radially outward, and some decreased radially outward. These types of variations and the fact that the measured profiles do not generally compare favorably with the transport model employed for AGR-3/4 may adversely impact the ability to extract reasonable transport parameters for this isotope. In many cases, the Cs-134 profiles decreased somewhat linearly in the outward radial direction, and in cursory comparisons, the shapes of these profiles appeared similar to those from model predictions. The step changes in concentration across the inner-outer ring gap were generally consistent with the model predictions as well. In some cases, there were local maxima in concentration at the outer surface of the rings. This suggests that fission products could have transported in the small gaps between the inner ring and the outer ring and between the outer ring and the sink ring such that some portion of a given fission product can bypass diffusion through the ring itself. The analysis of the small nubs on the outer surfaces of some of the outer rings revealed fission product concentrations in the nubs that were often higher than in the outermost segments of the rings. This further supports the hypothesis that short-circuit, gap transport occurred, causing relatively high surface concentrations on the outer surfaces of the rings. Eu-154 and Sr-90 profiles tend to have very similar shapes, suggesting that they transport via the same mechanisms. The observed profiles were indicative of a transport process where the isotopes are sorbed on the inner surface of the ring, but diffusion into the ring from that surface is quite slow. Some elevated concentrations of Sr-90 (relative to Eu-154) on the outer surface of a ring suggested that rapid, gas-gap transport of gaseous precursor Kr-90 and volatile Rb-90 could have occurred prior to their decaying to Sr-90. Overall, the Eu 154 and Sr-90 profiles were still very similar, which indicates that the transport of short-lived Sr-90 precursors is not a major effect. In some capsules, the qualitative Sr-90 behavior across the ring gaps was consistent with the model (using the available legacy Sr-90 transport parameters), but in other capsules the model was inconsistent with the measurements and seems to underestimate the amount of Sr-90 in the outer rings. The total ring Sr-90 inventories were estimated for all the rings that were subject to physical sampling. These results will be used to adjust the predicted particle and/or compact releases used in the AGR-3/4 fission product transport model. Given the different irradiation temperatures among the capsules and the rings, it was not possible to discern fundamental differences in the transport of isotopes within the different carbon materials, i.e., graphitic matrix, IG-110, or PCEA. It may be possible to do this in the course of determining transport from the concentration profiles in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Thermal-hydraulic performance of a high-temperature zigzag-channel printed circuit heat exchanger - 18563

Advanced nuclear reactors, such as high-temperature gas-cooled reactors (HTGRs) from the Generation IV nuclear systems, are designed with the capability of delivering high-pressure, high-temperature helium to power conversion units for electricity generation and industrial plants for process heat applications. The efficiency of the electricity generation and process heat applications of HTGRs is critically dependent upon an intermediate heat exchanger (IHX), which is a key component in transferring thermal energy from the primary coolant to a secondary coolant. Printed circuit heat exchangers (PCHEs) are promising to be employed in HTGR designs due to their capability for high-temperature, high-pressure applications and compactness. In the current study, a laboratory-scale zigzag-channel PCHE was fabricated and its heat transfer and pressure drop characteristics were investigated experimentally in a high-temperature helium test facility (HTHF). A commercial computational fluid dynamics (CFD) code, STARCCM+, was used to simulate the thermal-hydraulic performance of the fabricated PCHE with a simplified geometry model. Comparisons showed some differences between experimental data available for the zigzag-channel PCHE and the numerical results. Local thermal-hydraulic performance analyses indicated that fully-developed flow conditions and periodic changes were not observed, which could be attributed to the temperature variations along the flow channels, resulting in large fluid property variations. In addition, the effects of sharp-edged or rounded zigzag channel geometry were numerically evaluated. The use of rounded corners at each bend of the channels is recommended since this method can reduce the pressure drop on both the hot and cold sides of the heat exchanger while maintaining similar heat transfer capabilities. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Scoping Analysis of Pebble-Bed Reactors for the Destruction of the Transuranic Inventory of LWR Spent Nuclear Fuel

With the forecasted increase in the construction and operation of nuclear reactors, there will be a corresponding increase in the quantity of spent nuclear fuel (SNF) that requires long-term storage. In SNF, transuranic isotopes contribute the most to the long-term radiotoxicity of the fuel and pose a proliferation risk. One option that has been explored to address these issues is the removal of the transuranic isotopes from SNF and the conversion of these isotopes into transuranic fuel (TRU fuel). Here, this work sought to determine how effective a micro-modular Pebble-Bed High-Temperature Gas-Cooled Reactor (PB-HTGR); the 10-MW High Temperature Gas-cooled Test Reactor (HTR-10); and a salt-cooled small-modular pebble-bed reactor (PBR), i.e. the generic Fluoride-cooled High-temperature Reactor (gFHR), are at reducing the inventory of transuranic isotopes while still maintaining the intrinsic safety features of the PBR designs, such as negative temperature coefficients of reactivity. Optimized pebble designs utilizing TRU fuel were found for both reactors through the adjustment for the packing fraction of fuel in each pebble. The Axial Zone Equilibrium Modeling (A-ZEM) method was used in this work to help select the optimized pebble design. Once an optimized pebble design was selected and an equilibrium model was produced, the results from the deep burn (DB) HTR-10 and gFHR designs were compared to the results of two models from the literature. While both the DB gFHR and the DB HTR-10 were able to reduce the weapons-usable transuranic inventory, the performance of these reactors did not match that of the small-modular PB-HTGRs in the literature. Therefore, a need was identified for further refinement of the gFHR design using TRU fuel, as the results for this model were more promising than those of the DB HTR-10, which was strongly limited by the high leakage intrinsic to micro-modular PB-HTGRs.

Transuranic fuel↗
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 OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Thermal Hydraulics Validation Activities in the ART-GCR Campaign

The Advanced Reactor Technologies - Gas Cooled Reactor Campaign (ART-GCR) is part of the United States Department of Energy - Office of Nuclear Energy's Advanced Reactor Technologies program aimed at developing a High Temperature Gas-cooled Reactor (HTGR) which will offer enhancements in safety and efficiency. The Design, Methods, and Validation arm of the ART-GCR program provides experience and advanced tools for HTGR design and analysis. The Methods and Validation arm is developing data to validate predictions of decay heat removal through an extensive experimental campaign the Natural Convection Shutdown Heat Removal Test Facility at Argonne National Laboratory. In-core thermal hydraulics validation is being facilitated through the OECD/NEA's Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data. That benchmark contains code-to-code and code-to-data comparisons for 3 sets of HTGR thermal hydraulics phenomena. This ongoing benchmark activity has 16 participants from 9 countries developing leading insights into HTGR thermal hydraulics code validation. The ART-GCR campaign is also working on the development of state-of-the-art MOOSE-based models of the High Temperature Engineering Test Reactor (HTTR). Through the OECD/NEA's HTTR Loss of Forced Cooling (LOFC) program, the ART-GCR campaign is developing a multiphysics transient code validation benchmark based on 3 HTTR experiments that will provide an opportunity for international collaboration on multiphysics code validation.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Graphite Characterization: Baseline Variability Analysis Report

The Idaho National INL’s Baseline Material Property program provides many measurements of different material properties on a variety of nuclear graphite grades being investigated for irradiation-induced changes in material properties. The baseline material properties provide an unirradiated reference for comparison with measurements of those properties after irradiation. This study describes how uncertainty, or variance, in estimates of important material properties, like compressive strength, may depend on sample size, sample distribution and parameter estimation methods. Compressive strength data from two large collections of specimens from a billet of PCEA graphite (sample size = 230) and a billet of IG-110 graphite (sample size = 48) are used to illustrate these dependencies. The grades of graphite used represent approximate end members in terms of relative variability in strength measurements. Monte Carlo simulation and bootstrapping methods were used to examine variance in Weibull-parameter estimates as a function of sample size. Results were consistent with what previous statistical analyses have shown about the relationship between sample size and width of the confidence interval for maximum likelihood (ML) estimated parameters when samples are repeatedly drawn from a Weibull distribution. As a function of sample size, the precision of the modulus estimate increases faster than that of the characteristic strength. The modulus maximum likelihood estimates (MLEs) are biased for low sample sizes. The American Society for Testing and Materials document ASTM D7846-16 provides equations that describe the dependence of the variance of parameter MLEs on sample size. Using the same methodology, we show that the polynomials provided for variance dependence on sample size are inaccurate for sample sizes greater than ~120. We provide approximations to that relationship that illustrate relative constancy for sample size > ~10, in the logspace relationship of the modulus variance and the log transformed characteristic strength variance. Bootstrapping provided a means of comparing Weibull parameters obtained using ASME’s guidance for specimen collection to those obtained using different weighting in a case where strength is related to position in a billet and to grain orientation and a single distribution is used to characterize the billet. ASME specifies that samples be collected with equal representation in all locations and grain orientations. Parameter values obtained using ASME’s method were different from those obtained using random sampling of the complete data sets, and this was particularly true for the PCEA data. In addition, while the ASME approach yielded more conservative parameters for the PCEA data, the opposite was true for the IG-110 data. The difference observed for the PCEA data set likely reflects the larger number of against-grain specimens collected, which had generally higher strength than the with-grain specimens. This emphasizes the importance of the dependence of strength on location or grain orientation and suggests that it may be prudent to weight specimen collection toward the weaker subgroups if a conservative estimate as desired, as the equal weighting method recommended by ASME would will not always provide the most conservative Weibull parameters. Because spatial correlations commonly exist in certain grades of graphite, future work may involve developing a methodology for parameter inference that relaxes the assumption of independent specimens. ASTM and ASME codes recommend different parameter estimation methods and models of the Weibull distribution in different situations. Descriptions of the effects on parameter estimates for the example billets used in this study provide an example of how the method may affect the parameter estimate, and how—in this case—that compares to differences due to other factors, such as sample size and dependence on grain orientation or sample location.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Recent Testing of ENDF/B-VIII.0 at ORNL [Slides]

This presentation provides an overview of comparisons of ENDF/B-VII.1 and ENDF/B-VIII.0 and uses VALID suite to test >600 cases with ENDF/B-VII.1 and ENDF/B-VIII.0 in KENO V.a/KENO-VI, although a limited number of cases have been tested with TSUNAMI. It provides a results summary of KENO V.a/KENO-VI with a pooled results summary. The presentation also discusses changes for LCT experiments, MCT experiments, MST experiments, and PST experiments. A few words about graphite are also provided and the presentation states that graphite performance investigated across ENDF/B-VII.0, ENDF/B-VII.1, and ENDF/B-VIII.0 for MSR and HTGR models. A brief summary of the results is presented to make sure the Validation Committee is aware of the work. Also presented is an Eigenvalue comparison with a graphite moderated system, and comparison between an HTR-10 pebble model ENDF/B-VII.1 vs. ENDF/B-VIII.0, which provides different graphite data. In conclusion, ORNL has recently completed a comparison of ENDF/B-VIII.0 and ENDF/B-VII.1 using the VALID library. Work is also being performed for advanced reactor applications. Other efforts on-going and they are also trying to investigate more intermediate spectrum cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Multi-Scale Modeling of Thermal-Fluid Phenomena Related to Loss of Forced Circulation Transient in HTGRs

Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS), an effort had been pursued to support High Temperature Gas-cooled Reactors (HTGR) technology development and its modeling and simulation needs. There was a particular need for advanced modeling & simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially the core and lower and upper plena, during safety-related transients. In this report, two main such activities were presented: 1) system level HTGR simulations using system analysis code SAM for both normal operating and accident conditions; and 2) three-dimensional (3D) computational fluid dynamics (CFD) simulations on Deteriorated Turbulent Heat Transfer (DTHT) of interest to HTGR transients, including code benchmark and closure correlation development. In this work, Modular High Temperature Gas-cooled Reactor (MHTGR) full loop simulations were performed using System Analysis Module (SAM) code. Two approaches were used to model the reactor core; (i) unit-cell or lumped parameter model, and (ii) detailed ring model, to address the different needs under different reactor conditions. Aside from the core model changes other reactor components in the primary system, such as upper/lower plenum, heat exchanger and blower etc, were kept the same. In the unit-cell core model, the lumped parameter approach was used to model the coolant channels and the heat structures using 1-D and 2-D components, respectively. The fuel assembly heat structure was modeled based on unit-cell thickness approximation. In normal and controlled operating conditions, the calculated maximum temperatures for both coolant and heat structures were within the designed limits. The unit-cell model is, however, not able to capture the radial heat conduction through fuel assemblies, which is the dominant heat transfer mechanism during accident transients, such as pressurized conduction cooldown (PCC) event. To address this modeling challenge, ring models were developed. A simplified ring model was developed first. In this simplified model, the active core was simulated with 9 circular rings consists of 6 rings for homogenized fuel heat structure and 3 rings for gas coolant. It was quickly found that this model gives rather large simulation error as it was not able to correctly capture the solid structure representative geometries. A detailed ring model was developed to address this issue. The active core was simulated with 99 circular rings consists of 66 rings for homogenized fuel heat structure and 33 rings for gas coolant. In both ring models, six additional rings were included to represent inner reflector, outer reflector, core barrel, Reactor Pressure Vessel (RPV), RPV coolant channel, and Reactive Cavity Cooling System (RCCS). The detailed ring model was applied in two different SAM simulations; (i) normal operating conditions, and (ii) accident scenarios (more specifically, PCC event), for MHTGR design. In both simulations, the SAM calculated results for coolant and heat structure temperatures were within the reactor designed limits. The newly proposed detailed ring model approach in this work has proved to be a promising method to address the challenging modeling and simulation needs in the MHTGR design. As demonstrated in this work, this model was suitable for MHTGR simulations under both normal and accident conditions. Apart from the SAM model development for gas-cooled reactors, 3D Computational Fluid Dynamics (CFD) models were developed and validated with relaminarized Deteriorated Turbulent Heat Transfer (DTHT) MIT benchmark test. Two different CFD codes, Nek5000 and STAR-CCM+, were used for this benchmark activity. Compared with MIT data, Nek5000 shows very good agreement between its simulated and experimentally measured wall temperatures. In contrast, STAR-CCM+ under-predicted wall temperature due to over predicted turbulence in the wall region. The Nek5000 code was also used to develop wall-heat transfer correlations for laminar flow in cylindrical tube. Several simulations were performed for various Reynolds flow and wall heat fluxes. Significant deviations were found when compared the predicted heat transfer coefficient and the Sieder-Tate correlation. A new set of heat transfer correlation was proposed for laminar flow which is valid for 400<Re<2000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Automated Segmentation of Twin Boundaries in TRISO Silicon Carbide Using Deep Neural Networks

Coated particle fuels, such as the tristructural isotropic (TRISO) fuel particle, are essential for high-temperature gas reactor (HTGR) applications due to their efficiency and stability under normal and off-normal conditions. However, widespread commercialization and deployment of this technology for next-generation nuclear applications require robust quality assurance and quality control (QA/QC) methods linking fabrication, properties, and performance. Of the many important metrics for TRISO QA/QC, quantification of the silicon carbide (SiC) microstructure is critical because it correlates with fission product retention during irradiation. Previous work has shown extensive twinning of the SiC microstructure, which strongly affects microstructural metrics; however, twin grain boundaries are not expected play a significant role in fission product diffusion. This report summarizes the initial development, training, and testing of a machine learning image processing algorithm to detect twin grain boundaries in a backscattered electron image, which can be removed so that microstructural metrics can be recalculated for legacy data. Further development and deployment of this model will provide automated, scalable improvement of potential QA/QC methods for the SiC layer of TRISO particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Analysis of buffer-IPyC separation in TRISO fuel particles

During High Temperature Gas-cooled Reactor (HTGR) operation, and due to the neutron irradiation in the reactor core, damage of the nuclear fuel coating layers occurs. The mechanism of damage formation in the TRISO fuel is explored by the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, in which the debonding process between coating layers was also investigated. The purpose of this paper is to report simulation results for two models. Firstly, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted. Secondly, a debonding enabled model, where the gap between those layers is created. The simulations were performed with the Bison code. The inputs of the simulated models are based on data from the AGR-1 experiment. Further, the research included simulations on spherical and aspherical fuel types. The simulations match results obtained by the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the behavior of the gap between buffer and IPyC layers in TRISO fuel. Based on the irradiation experiments, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. In addition, the fuel sphericity does not influence the Weibull failure probability. The results obtained for spherical and aspherical fuel are nearly identical.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Europium diffusion in IG-110 nuclear graphite

Europium diffusion in graphite has been recognized to be of interest in high-temperature gas reactor safety analysis, particularly as an indicator of strontium diffusion. However, no measurements of europium diffusion coefficients have been reported in the literature. In this work, the effective diffusion coefficient of europium was measured using a time-release method. Natural europium was loaded into pre-milled unirradiated IG-110 graphite spheres using a pressurized acid digestion vessel. The time-release experiments were performed in the temperature range 1823 K – 1973 K using a SiC diffusion cell connected to an inductively-coupled plasma mass spectrometer (ICP-MS) via a He gas line. Finally, the results of this work are: D Eu, IG-110 = (1.5 x 10 -3 m 2 /s) exp $(\frac{-2.87\mathrm{x}10^{5} J/mol}{RT})$ This effective diffusion coefficient can be used to aid in predictive modelling of europium transport in HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2016

PHISICS/RELAP5-3D RESULTS FOR EXERCISES II-1 AND II-2 OF THE OECD/NEA MHTGR-350 BENCHMARK

The Idaho National Laboratory (INL) Advanced Reactor Technologies (ART) High-Temperature Gas-Cooled Reactor (HTGR) Methods group currently leads the Modular High-Temperature Gas-Cooled Reactor (MHTGR) 350 benchmark. The benchmark consists of a set of lattice-depletion, steady-state, and transient problems that can be used by HTGR simulation groups to assess the performance of their code suites. The paper summarizes the results obtained for the first two transient exercises defined for Phase II of the benchmark. The Parallel and Highly Innovative Simulation for INL Code System (PHISICS), coupled with the INL system code RELAP5-3D, was used to generate the results for the Depressurized Conduction Cooldown (DCC) (exercise II-1a) and Pressurized Conduction Cooldown (PCC) (exercise II-2) transients. These exercises require the time-dependent simulation of coupled neutronics and thermal-hydraulics phenomena, and utilize the steady-state solution previously obtained for exercise I-3 of Phase I. This paper also includes a comparison of the benchmark results obtained with a traditional system code “ring” model against a more detailed “block” model that include kinetics feedback on an individual block level and thermal feedbacks on a triangular sub-mesh. The higher spatial fidelity that can be obtained by the block model is illustrated with comparisons of the maximum fuel temperatures, especially in the case of natural convection conditions that dominate the DCC and PCC events. Differences up to 125 K (or 10%) were observed between the ring and block model predictions of the DCC transient, mostly due to the block model’s capability of tracking individual block decay powers and more detailed helium flow distributions. In general, the block model only required DCC and PCC calculation times twice as long as the ring models, and it therefore seems that the additional development and calculation time required for the block model could be worth the gain that can be obtained in the spatial resolution

Advanced Reactor Technologies↗
DOE OSTI2024

Equilibrium core modeling of a pebble bed reactor similar to the Xe-100 with SCALE

As the nuclear industry moves towards licensing and constructing advanced reactors, new attention has been focused on the advanced reactor designs that have past operational experience, such as pebble-bed high-temperature gas-cooled reactors (PB-HTGRs). Pebble-bed reactor designs have many advantages, such as their higher operating temperatures and online refueling capabilities. However, high-fidelity computational modeling of pebble-bed reactor designs, from reactor startup to operation at equilibrium, is more challenging compared to conventionally fueled reactors due to the continuous movement of the fuel pebbles through the reactor during operation. In previous work at Oak Ridge National Laboratory (ORNL), the SCALE Leap-In method for Cores at Equilibrium (SLICE) was developed around tools within the SCALE code system. This iterative method can effectively generate pebble-bed reactor zone-wise fuel inventories at equilibrium core operation within a reasonable computational time. The objective of this work was to further verify the ORNL SLICE method and to investigate the impact of considering temperature profiles during the application of the method. The SLICE method was applied to a modular high-temperature gas-cooled reactor design based upon publicly available design specifications of the Xe-100 pebble-bed reactor. Upon comparing the results from the SLICE method to published literature, the differences in the eigenvalue k effective were on the order of several hundred pcm (percent millirho). To investigate one possible cause of these differences, a study looking at the sensitivity of the full-core equilibrium k effective and discharge nuclide inventory to temperature was performed by developing equilibrium cores of two additional temperature profiles. From this temperature study, differences on the order of hundreds of pcm for the full-core equilibrium k effective , and up to 15% difference for the discharge inventories were found. In conclusion, these results indicated the strong dependence on temperature that needs to be considered for future work in equilibrium modeling of PB-HTGRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Optimization Methodology of Pebble Bed HTGR Start-Up and Running-in Strategy

In recent years interest in advanced reactor technologies has increased significantly. However, the methods used for analysis of traditional nuclear reactors are insufficient to consider all the different and varied advanced reactor designs without further development. One promising advanced reactor design is the pebble bed reactor (PBR). PBRs possess unique operational and fuel cycle features that require the development of specific analysis methodologies to adequately design and analyze the systems. There is a need in PBR research for a capability to analyze and optimize the process of transitioning from the start-up reactor core to the equilibrium core (known as the “running-in” of the reactor). The start-up of a PBR and the transition to the equilibrium core is a complex, multi-physics challenge that has not yet been well researched and has many opportunities for design, analysis, and optimization of the process. In this research, a methodology is defined to consider the potential strategies in PBR start-up and run-in to the equilibrium core. Multiple candidate software are considered, and their pros and cons are discussed for the PBR-specific application in the methodology. A preliminary software selection for the physics engine is made, and initial verification of key modules is performed. Software to support optimization of the reactor run-in strategies through reduced order modeling (ROM) and machine learning are considered. Challenges for a full implementation of the methodology are discussed. Additional code selections and verification of models relevant to the application are needed before full demonstration of the methodology can be achieved and an optimal strategy determined.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Summary Report of the FY24 DOE Contributions to the GIF VHTR CMVB

The Generation-IV Forum (GIF) Very-High-Temperature Reactor-Computational Methods Validation and Benchmark (VHTR-CMVB) initiative, involving organizations from Korea Atomic Energy Research Institute (KAERI) (South Korea), Institute of Nuclear and New Energy Technology of Tsinghua University (INET) (China), U.S. Department of Energy (DOE) (U.S.), Joint Research Centre (JRC) (Europe), and Japan Atomic Energy Agency (JAEA) (Japan), is dedicated to the verification and validation of tools for High-Temperature Gas-Cooled Reactors (HTGRs) analysis, using data shared by Computational Methods Validation and Benchmark (CMVB) signatories. For FY24, the US DOE CMVB has committed to several critical activities. Under WP1, led by the US, the integration of the High Temperature Gas Cooled Reactor - Pebble-Bed Module (HTR-PM) Phenomena Identification and Ranking Table (PIRT) into the comparative PIRT is progressing, with a new draft of the comparison tables issued earlier this year and currently being utilized by INET for their contribution. Neutronic validation efforts under WP3 include the preparation of the burnup analysis benchmark, preliminary calculations, and the development of reference models and results. In WP2, a validation exercise for hot gas mixing in the lower plenum of HTR-PM is in progress, using experimental data from INET (China) to validate modeling approaches. A model of the experimental facility has been developed using StarCCM+, with initial calculations slated for presentation at the GIF CMVB meeting this fall. Another WP2 activity focuses on validating numerical models for air-cooled Reactor Cavity Cooling System (RCCS) with experimental data from the Wisconsin Madison RCCS facility. A high-fidelity model, developed using NEK-RS, is currently being validated with available data from a low power forced convection test. These efforts are aimed at enhancing and confirming the accuracy of HTGR analysis tools, ensuring their alignment with experimental data and regulatory requirements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Methods development towards automated, physics-informed, quantitative quality control of TRISO-SiC

Tristructural-isotropic (TRISO) fuel particles have been developed as a high-performance fuel for use in high-temperature gas-cooled reactor (HTGR) systems due to their high efficiency and stability under both normal and off-normal conditions. Broader deployment of this technology in advanced nuclear applications may benefit from quantitative quality assurance and quality control (QA/QC) methods that directly link TRISO properties to downstream performance. A key TRISO property is the SiC layer microstructure, which influences fission product retention during irradiation. However, existing QA/QC for the TRISO-SiC microstructure comprises only a qualitative visual inspection; therefore, there is a clear opportunity for the development of quantitative methods for TRISO QA/QC. Here, to this end, previous work has demonstrated an image processing approach to grain boundary (GB) identification and subsequent extraction of microstructural metrics; however, extensive twinning within the SiC layer complicates such analyses because twin GBs significantly influence microstructural metrics but are not expected to contribute to fission product transport. This study presents the initial development, training, and testing of an ML-based image segmentation algorithm designed to identify and remove twin GBs from standard backscattered electron micrographs, providing an industrially applicable, quantitative, and physically meaningful QA/QC approach for the TRISO-SiC microstructure. Although pixel-wise performance metrics for the twin predictions are low, the change in grain area and the number of GB pixels after twin removal predicted by the ML workflow are within 1% of the true values calculated using crystallographic data. This suggests that the model is well capable of predicting overall twin boundary structures and grain morphology, and continued advancement of this approach could enable automated, scalable, and physics-informed QA/QC for TRISO-SiC microstructures, supporting the reliable qualification of coated particle fuels for next-generation reactor systems.

Computer vision↗