Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pebble bed reactor”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors

Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Conceptual Design of Passive Neutron Albedo Reactivity and Active Neutron Albedo Reactivity Instruments for Various Arrangements of Pebble Bed Reactor Pebbles

MCNP6 simulations were performed to assess the anticipated capability of the Passive Neutron Albedo Reactivity technique for nuclear safeguards measurements of irradiated pebble/pebbles. Two physical setups were examined. One with a single pebble of three different burnups and one with 27 pebbles in a cube. For all cases the sensitivity to removing all the fissile material was examined. For both of the physical setups, 4 assay cases were examined: (a) singles count rates for which the only source neutrons were the inherent neutrons in the fuel, (b) singles count rates again except this time two relatively weak AmLi sources were placed above and below the top central pebble and the combined effect of both neutron source terms were examined, (c) doubles count rates for which the only source neutrons were the inherent neutrons in the fuel, and (d) doubles count rates again except this time two relatively weak AmLi sources were placed above and below the top central pebble and the combined effect of both neutron source terms were examined. Finally, a few parameters such as Cd near the 3 He tubes were perturbed to see if the results could be improved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Utilizing Advanced Statistics to Determine Anomalistic Conditions in Pebble-Bed Reactors

Pebble-bed reactors (PBRs) utilize hundreds of thousands of fuel pebbles, which continuously circulate through the core, in lieu of traditional fuel assemblies to generate fissions and produce power. The use of unmarked fuel pebbles presents a challenge for international safeguards verification that nuclear material is not being diverted. To ensure pebble diversion is not taking place, new methods for accounting for and monitoring the pebbles should be examined to determine an appropriate methodology for performing bulk accountancy with pebbles. Here, this work examines the use of statistical methods for determining if the reactor is within a declared range of operation by examining the statistical distribution of pebble burnup as they are discharged from the core. Using this methodology, we created a model that detects diversion over 95% of the time, over multiple diversion pathways, if the reactor core maintains a constant power density during the diversion process and only falsely labels a diversion case nominal 2% of the time. For a diversion scenario where the reactor is maintained at a constant power, the statistical analysis can correctly identify if diversion is occurring over 80% of the time; however, nearly 20% of specific diversion pathways are mislabeled nominal. These results provide a basis and framework for exploring the further use of statistical methods to determine where these methods could be most useful and where additional methods, such as machine learning, could be used to capture if diversion is occurring in pebble-bed reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Deterministic High-Fidelity Neutronics Simulation of Pebble Bed Reactors Using Pebble Tracking Transport

The pebble tracking transport (PTT) algorithm offers a high-fidelity deterministic approach for neutron transport for pebble bed reactors (PBRs). This approach requires the mesh for the active-core region to consist exclusively of tetrahedral elements, where each node in the pebble-packing region represents a pebble centroid. This paper investigates the application of PTT for full-scale PBRs, considering both the isothermal and the temperature-dependent core conditions. Macroscopic cross sections are generated using Serpent 2 full-core eigenvalue simulations where pebbles are grouped into disjoint subsets using machine learning. To minimize the need for individual cross-section sets for each pebble in the core, K-means clustering is used to group pebbles by temperature and neutronic environment parameters. Here, we compare the multiplication factor and power rate distributions between PTT simulations using the Griffin reactor physics software and reference solutions from Serpent 2. Our analysis shows that a full-core, high-fidelity PTT calculation produces accurate results with minimal local (pebblewise) errors. Additionally, timing results indicate that PTT simulations converge rapidly on modern supercomputing platforms.

Griffin↗
DOE OSTI2026

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

High-fidelity Pebble Bed Reactor Depletion Based on Pebble Tracking Transport in Griffin

The pebble tracking transport (PTT) method is a high-fidelity, heterogeneous deterministic transport technique for pebble bed reactor analysis. It discretizes the broad-group neutron transport equation in space and angle with the discontinuous finite element and the discrete ordinates method, and utilizes various solving techniques, including mesh sweeping and diffusion acceleration, to provide pebble- wise reaction rates. This work presents the extension of the PTT method to enable fuel depletion capability in the Griffin code. We discuss the implementation details of the PTT-based high-fidelity depletion where isotope inventory of all individual pebbles is tracked through pre-determined pebble flow paths in the core. The implementation is verified with a generic pebble bed reactor model. Some preliminary equilibrium core results are included. Future works are also discussed.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Burnup Monitoring for Pebble Bed Reactor Systems

A pebble burnup monitoring system is a required component for domestic reactor safety and safeguards applications associated with pebble bed reactors (PBRs). One of the main requirements of a PBR burnup monitoring system is that it needs to be capable of rapid measurements to assess the burnup of each individual pebble to determine whether to recirculate it in the reactor or discard it as spent fuel. This report considers three different approaches for a burnup monitoring system for pebbles discharged from the reactor core in a pebble bed modular reactor-400: • passive gamma spectrometry measurement, • passive neutron coincidence measurement, and • active neutron counter based on the differential die-away technique. Conceptual designs have been created for each of these detectors, and preliminary analysis has been performed using Monte Carlo N-Particle and Oak Ridge Isotope Generation code simulations. The advantages and practical limitations (e.g., high radiation background) of each system were identified. Simulations suggest that each of the three measurement techniques can be successfully employed to distinguish between pebbles based on their number of passes through the core and to quantify the burnup of pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Assessment of Self-interrogation Safeguards Signatures for Pebble Bed Reactor Fuel

New reactor designs, such as pebble-bed reactors, present challenges for the safeguarding of fissile material. Due to the multi-pass fuel circulation design and the non-uniform path a pebble may traverse through a pebble-bed reactor, there will be variations in the irradiation history of spent fuel pebbles. The standard approach of estimating actinide quantities based on fission gamma spectra or neutron emission paired with depletion modeling may yield uncertainties too large for safeguards and material accountancy. This project investigated, through modeling and experimentation, the potential of neutron self-interrogation of spent fuel pebbles as an innovative method to implement materials accountability. As an example, our feasibility studies indicate that the mass of U-235, U-238, Pu-239 and Pu-241 can be predicted to 4.1, 0.86, 13 and 13 % accuracy, respectively, when measuring 100 closely packed end-of-life spent fuel pebbles over approximately 12 days using a 4π counting geometry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Equilibrium Core Model for Micro Pebble Bed Reactors Using OpenMC

Estimating the equilibrium state for pebble bed reactors (PBRs) presents complex challenges as it requires simultaneous consideration of changes in the pebbles’ movement as well as their fuel compositions. Whereas traditional approaches use multigroup diffusion codes for neutronics calculations of PBRs’ equilibrium state, the double-heterogeneity of PBRs complicates neutron cross-section generation. Continuous-energy Monte Carlo (MC) methods are better suited for detailed PBR analysis because of their natural handling of double-heterogeneity, but they demand substantially more computational resources. Here, this study introduces a novel method for efficiently estimating the equilibrium state in small and micro PBRs with reduced computational cost. The method is anticipated to accelerate the processes of core design and performing parametric studies for utilizing advanced fuel and structural materials. The HTR-10 reactor design was used for validating the method’s predictions and evaluating its computational efficiency. When compared to reference calculation values from the literature, criticality (k-effective) was predicted to be approximately within the margin of error of the MC transport calculation, average core power density (in megawatts per cubic meter) was predicted within 2.5% relative error, and maximum thermal flux (10 13 n/cm 2 .s −1 ) was predicted within 1.8% relative error. The calculated inventory of fission products and fuel composition in the equilibrium core were within 15% and 16.6%, respectively, when compared to reported values from the literature. The difference is attributed to variance in the considered values of the core temperature, which was found to significantly affect the depletion analyses.

Equilibrium core↗
DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Parameter Study of the Running-In Process for the Generic Pebble Bed Reactor (GPBR200)

The run-in period of a pebble bed reactor is complex and difficult to model given significant heterogeneity in core composition, power, and temperature. While it is understood that the initial composition of the core should eventually result in the same equilibrium core composition, the approach to equilibrium can vary significantly depending on factors such as start-up fuel enrichment and power ramp rate. To explore this, a high-fidelity model of the General Pebble Bed Reactor was used to vary power ramp schemes and start-up core compositions. It was found that both the initial core composition and power ramp rate had a significant impact on the flow rate of pebbles during early time steps, with higher ramp rates and low enrichment resulting in non-physical flow rates. Power ramp rate alone was found to dictated maximum pebble power peaking observed during the run-in process, with higher ramp rates resulting in greater peak pebble powers. Start-up fuel enrichment and power ramp were found to both impact total fuel consumption, although the impact of start-up fuel enrichment was generally secondary to ramp rate.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Pressure Drop Correlation Improvement for the Near-Wall Region of Pebble-Bed Reactors

Packed beds play an important role in several engineering fields, with their applications in nuclear energy being driven by the development of next-generation reactors utilizing pebble fuel. The random nature of a packed pebble bed creates a flow field that is complex and difficult to predict. Porous media models are an attractive option for modeling pebble-bed reactors (PBRs), as they provide intermediate fidelity results and are computationally efficient. Porous media models, however, rely on the use of correlations to estimate the effect of complicated flow features on the pressure drop and heat transfer in the system. Existing correlations were developed to predict the average behavior of the bed, but they are inaccurate in the near-wall region where the presence of the wall affects the pebble packing. This work aims to investigate the accuracy of a porous media model using the Kerntechnischer Ausschuss (KTA) correlation, the most common pressure drop correlation for PBRs compared to the high-fidelity large eddy simulation (LES). A bed of 1568 pebbles is investigated at Reynolds numbers from 625 to 10 000. The bed is divided into five concentric subdomains to compare the average velocity, friction losses, and form losses between the porous media and LES codes. The comparison between the LES simulation and the KTA correlation revealed that the KTA correlation largely underpredicts the form losses in the near-wall region, leading to an overprediction of the velocity near the wall by nearly 30%. An investigation of the form losses across the range of Reynolds numbers in the LES results provided additional insight into how the KTA correlation may be improved to better predict these spatial effects in a pebble bed. These data suggest that the form coefficient near the wall must be increased by 48% while decreasing the form coefficient of the inner bulk region of the bed by 15%. The implementation of these improvements to the KTA correlation in a porous media model produced a radial velocity profile that saw significantly improved agreement with the LES results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Generation of localized reactor point kinetics parameters using coupled neutronic and thermal fluid models for pebble-bed reactor transient analysis

The systems analysis of anticipated operating occurrences and design basis accidents for pebble-bed reactor systems requires knowledge of neutron point kinetics equations (PKE) parameters. Typically, the generation of PKE parameters is performed in a global manner using standalone neutronics calculations, without the inclusion of thermal fluid distributions. We utilize Griffin and Pronghorn for generating global and local PKE parameters which includes the use of thermal fluid distributions to account for localized effects. This work establishes a methodology for calculating PKE parameters for a pebble bed reactor with a coupled neutronics/thermal fluids analysis. PKE parameters generated on a global and local basis are compared against a diffusion solve for a typical load-following transient. Locally-generated neutron kinetic parameters are able to reduce the maximum error in the transient power level from 5% to below 1.5%; along with this, bulk temperature errors were reduced from 12 K to 4 K.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Statistical Approaches for Pebble Bed Reactor Operations and Safeguards

The design of pebble bed reactors (PBRs) and their method of operation align more closely with statistical approaches used in manufacturing and process control than traditional safeguards statistical approaches. The reason is PBRs will employ a nondestructive assay (NDA) measurement (burnup measurement system [BUMS]) that is part of the fuel handling system supporting discharge decisions in addition to the reactor code which monitor performance. The integration of these two approaches provides the opportunity to monitor reactor performance statistically for both operations and safeguards in ways not achievable using other reactor designs. For light water reactors (LWRs), knowledge about reactor code performance in predicting irradiated special nuclear material (SNM) content historically was only achieved from special measurement campaigns or from fuel reprocessing. Conversely, through statistical comparison of the BUMS with the reactor code-predicted values, PBRs can achieve this in real time. The resulting SNM distribution is an indicator of reactor performance because factors such as transit time and path of the pebble fuel through the reactor determine the plutonium production and uranium depletion. By analyzing the predicted and measured values, opportunities exist to adjust operating parameters, fuel design, and other characteristics to optimize performance and fuel utilization. From a safeguards perspective, this approach also provides the information necessary to validate declared values and evaluate whether the reactor is being operated as expected. This paper outlines statistical approaches for PBRs that can be used to support both operations and safeguards.

Ball, Cory↗
DOE OSTI2023

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗
DOE OSTI2023

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of five subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the five pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗