Abstract
Stiffness confinement method (SCM) can effectively alleviate the stiffness in the neutron kinetics equation and can use larger time steps to obtain the same calculation accuracy and improve the computational efficiency. However, the existing SCM is mainly used to solve two group transient diffusion equations. In this paper, the SCM is employed to solve the multigroup transient transport calculation. On the basis of the original method of characteristics (MOC) code PEACH, the transient function is added, and the PEACH-K program is developed. Based on the numerical results of the latest published OECD/NEA C5G7-TD benchmark, it shows that the PEACH-K program developed in this paper has both high computing precision and good numerical stability.