Development and Verification of Neutron Kinetics Parameter Calculation Program Applicable to Subcritical Systems of Long-Life Reactors

Abstract

The fundamental prompt neutron delayed constant α, the effective delayed neutron fraction β eff , and the effective neutron generation time Λ eff are important for accurate measurements of keff in subcritical systems. In this paper, based on the open-source program OpenMC, a module with the function of calculating the α-eigenvalue is developed based on the k-α iterative method. Based on the nodal expansion method (NEM), an adjoint neutron flux calculation program was developed, which in turn used the idea of coupling the k-α iterative method with the NEM to obtain the neutron kinetics parameters (α, β eff , and Λ eff ) of the subcritical system and allowed for a more accurate keff under the subcritical system as compared to the traditional k-mode program. In this paper, the results of the α-eigenvalues calculated by the program are validated using the Godiva benchmark and the MUSE-4 device as computational objects. The accuracy of the results for β eff and Λ eff was verified with the KUCA (Kyoto University Critical Assembly) benchmark. The results of the keff calculations using k-mode and α-mode are compared with the ADS benchmark. The results show that the calculation program (α-mode) can obtain keff under the subcritical system more accurately, which makes up for the shortcomings of the traditional program, which cannot be calculated under the deep subcritical system. The relative error between the calculated and reference values was less than 750 pcm at deep subcriticality (keff = 0.899). The error was less than 400 pcm at the deeper subcritical system (keff = 0.952), which can satisfy the demand for neutron dynamics parameters and keff calculation in the subcritical system. The program can also continue to be upgraded in the later stage in the calculation method and other aspects.