ARCHER - a new Three-Dimensional method of characteristics neutron transport code for Pebble-bed HTR with coarse mesh finite difference acceleration

Abstract

A new deterministic code ARCHER using Three-Dimensional (3-D) Method of Characteristics (MOC) has been developed by Institute of Nuclear and new Energy Technology (INET), Tsinghua University, aiming to obtain high-fidelity transport solution of pebble-bed High Temperature gas-cooled Reactor (HTR) with explicit pebble-bed geometry. However, 3-D MOC usually suffers from the slow convergence rate without efficient acceleration methods. In this work, the Coarse Mesh Finite Difference (CMFD) based on regular mesh is implemented to accelerate 3-D MOC calculation, which is available for the problems with cylinder and cuboid geometry. The complex geometry of reactor core with randomly stacked spherical fuel elements and irregular helium flow area, as well as the topological relationships between fine Flat Source Region (FSR) mesh for MOC and coarse CMFD mesh, are the key but challenging tasks for high-fidelity transport solution. Constructive Solid Geometry (CSG) and Boolean operation are constructed to describe problem geometry, which has the ability to deal with complex topological relationships between meshes. In addition, tree grid structure is utilized to reduce the time complexity of ray tracing and obtain the grid mapping relationship between coarse mesh and fine mesh. Spatial domain decomposition parallel based on MPI is utilized to alleviate the memory pressure and pursue high computational performance. Ray parallel is further achieved through OpenMP with shared memory. Several problems have been used to verify the geometric modeling capability and computational performance of this new 3-D MOC code for pebble-bed HTR.