Application of stiffness confinement method within variational nodal method for solving time-dependent neutron transport equation

Abstract

The stiffness confinement method (SCM) is implemented in the variational nodal method (VNM) framework to solve time-dependent multi-group neutron transport equations. With frequency-transform, the time-dependent neutron equation is decoupled from the precursor equation and transformed into a dynamic eigenvalue problem which can be solved with the in-house VNM code VITAS. An improved time-source iteration method within the VNM-SCM framework is proposed to accelerate constructing response matrices and reduces memory usage. Based on the VNM-SCM scheme, the kinetics module of VITAS is developed for the Cartesian and hexagonal-z geometry and verified by several benchmark problems.Numerical results demonstrate that VNM-SCM can perform accurate and efficient transient calculations with large time-step sizes. For the two-dimensional TWIGL benchmark, the VNM-SCM produces a maximum core power error of -0.57% with the time-step size of 20 ms. For the three-dimensional LMW benchmark, the VNM-SCM results in a maximum core power error of 0.95% with the time-step size of 2 s. The VNM-SCM also yields an accurate neutron flux prediction for the hexagonal-z geometry problem. In addition, the effect of the angular expansion is investigated systematically. Examinations with the three-dimensional LMW benchmark display an error reduction in the core power error from 0.131% to 0.011% from P1 to P3 interfacial angular expansion, with a substantial computing time increase of 6.14 hours. The study indicates that the high-order angular expansion significantly increases the computing time but contributes little to the accuracy of the core power. Besides, it was found that the initial steady-state calculation dominates the shape error of the flux in the transient calculation.