Neutronics Analysis on High-Temperature Gas-Cooled Pebble Bed Reactors by Coupling Monte Carlo Method and Discrete Element Method

Abstract

The High-Temperature Gas-Cooled Pebble Bed Reactor (HTG-PBR) is notable in the advanced reactor realm for its online refueling capabilities and inherent safety features. However, the multiphysics coupling nature of HTG-PBR, involving neutronic analysis, pebble flow movement, and thermo-fluid dynamics, creates significant challenges for its development, optimization, and safety analysis. This study focuses on the high-fidelity neutronic modelling and analysis of HTG-PBR with an emphasis on achieving an equilibrium state of the reactor for long-term operations. Computational approaches are developed to perform high-fidelity neutronics analysis by coupling the superior modelling capacities of the Monte Carlo Method (MCM) and Discrete Element Method (DEM). The MCM-based code OpenMC and the DEM-based code LIGGGHTS are employed to simulate the neutron transport and pebble movement phenomena in the reactor, respectively. To improve the computational efficiency to expedite the equilibrium core search process, the reactor core is discretized by grouping pebbles in axial and radial directions with the incorporation of the pebble position information from DEM simulations. The OpenMC model is modified to integrate fuel circulation and fresh fuel loading. All of these measures ultimately contribute to a successful generation of an equilibrium core for HTG-PBR. For demonstration, X-energy’s Xe-100 reactor—a 165 MW thermal power HTG-PBR—is used as the model reactor in this study. Starting with a reactor core loaded with all fresh pebbles, the equilibrium core search process indicates the continuous loading of fresh fuel is required to sustain the reactor operation after 1000 days of fuel depletion with depleted fuel circulation. Additionally, the model predicts 213 fresh pebbles are needed to add to the top layer of the reactor to ensure the keff does not reduce below the assumed reactivity limit of 1.01.

Keywords

• high-temperature gas-cooled pebble bed reactor
• Monte Carlo Method
• discrete element method
• equilibrium state
• neutronics analysis