Development of Digital Reactor High-fidelity Neutronics Code SHARK
ZHANG Hongbo,ZHAO Chen,PENG Xingjie,ZHAO Wenbo,CHEN Zhang,LI Qing,YU Yingrui,GONG Zhaohu,ZENG Wei,LIU Kun,RAO Junjie,WANG Bo
DOI: https://doi.org/10.7538/yzk.2021.youxian.0902
2022-01-01
Abstract:The concept of digital nuclear reactor, which is featured with high-resolution modeling and highfidelity multiphysics simulation, has been proposed and developed in recent years with the dramatic development of highperformance computers. It has special and important meanings for the design improvement and mechanistic understanding of nuclear energy system. In the digital reactor system, the highfidelity core physics neutronics simulator is a key component. In this paper, a newly developed highfidelity neutronics code SHARK (Simulationbased Highfidelity Advanced Reactor physics Kit) was introduced. This is a onestep heterogeneous transport code developed by Nuclear Power Institute of China (NPIC). It adopts constructive solid geometry (CSG) method to model the problem geometry, which enhances its geometry flexibility including square and hexagonal lattices with different kinds of fuel pellets. In methodlogies, the code is currently characterized by a 45group neutron library, subgroup selfshielding method, and both 2D/1D method of characteristics (MOC) and quasi-3D MOC as transport solvers. For the multigroup library, the selected energy group structure was proved to be accurate enough and effective for light water reactor (LWR) applications. The isotopes in the library covered fuels, moderator and coolant, poison absorbers, structural materials and their depletion chain daughters. Reaction types including fission, capture, (n, 2n)/(n, 3n) and decays were considered. For the resonance selfshielding, an improved subgroup method based on equivalence crosssection tables was utilized. This method kept the advantages of conventional methods in the sense of geometry flexibility and distribution effects, and gained better efficiencies on computational costs. Resonance scattering and interferences were also properly treated. For the heterogeneous transport solvers, 2D/1D MOC was applied with a discrete ordinate (SN) finite difference axial solver. To improve the stability of neutron transport calculations, quasi-3D MOC was also developed, which did not contain a transverse leakage term and gained stability. For acceleration and parallelization, optimally diffusive coarse mesh finite difference method (odCMFD) was used to improve the convergence stability, and an MPI plus OpenMP mixed parallel scheme was implemented. The overall calculation frame was established by the languages of Python and C++. The code used objectoriented design concept and had good maintainability and user friendliness. The steadystate verification results were described in the paper including some VERA progression benchmarks and macro BEAVRS problem. The preliminary benchmarking results demonstrate the methodology and code developments in terms of reactivities, radial and axial power distributions and control worth. In the future, performance and functionality will be improved continuously.