Study on the Curved Channel Model in the Initial Core Loading of Pebble Bed High Temperature Reactors

Abstract

Continuous on-line fuel cycling is the essential feature of the pebble bed high temperature reactor (PB-HTR). The flow speed of the fuel pebbles in a PB-HTR presents a radial distribution in the reactor core, mainly due to the friction between the pebbles and the wall and the conical structure at the core bottom. In the VSOP fuel shuffling model, the simulation of unequal pebble flow speed is achieved by dividing the reactor core into some vertical flow channels with different numbers of the equal-volume regions in each channel. However, the fuel shuffling with equal-volume batches bring complexity when dealing with the change of fuel composition, such as the fuel fraction of fuel-graphite pebble mixture, during the initial core loading and early running-in phase. In this work, a curved channel model with unequal flow speed and the bottom cone is established based on the DEM simulation of pebble flow in the HTR-PM. The batch-wise fuel shuffling strategy is adapted to fit the complex situation during mixing and re-assigning the discharged fuels by employing a rounding strategy for the actual volume of fuels with similar irradiation history. The key of the adapted strategy is to divide the total number of the mixed batches with similar irradiation history by the number of flow channels, and round the quotient as the number of reloaded batches in each top region. The fuel loading process to build up the initial core, accompanied by the low-power reactor running to compensate the reactivity provided by the fresh fuels, is simulated by using the fuel shuffling model mentioned above. On the other hand, the simulation on the same process with an effective cylindrical core mesh and straight flow channels is carried out, in which dividing and rounding the batch numbers are unneeded. The results of both models are compared, indicating that the curved channel model presents less core reactivity and shorter fuel loading period than those of the cylindrical model. From the point of view of fidelity, the former is more suitable for the simulation of initial core loading process. The results in this work are important for enhancing the economy of fuel cycling of PB-HTRs.