Abstract
Point-reactor neutron kinetics equation is an important part of neutron dynamics. It is a stiff ordinary equation, and the stiff problem would become more significant especially for fast reactor. The researches and applications of numerical methods for solving point-reactor neutron kinetics equation have obtained many achievements, but those methods are not quite ideal enough. The application of point-reactor model on real-time simulation of fast reactor core physics requires high precision and meets the needs of real-time simulation. So the numerical solution of point-reactor neutron kinetics equation still needs further study. Gauss precise time-integration method is a numerical method combining the Gauss integration method with the precise time-integration method which is used to solve exponential matrix. It is a stable explicit integration method. For solving the point-reactor neutron kinetics equation, this method can theoretically reach any high precision and own high calculation efficiency. It avoids calculating the inverse matrix by precise time-integration method, thus avoids the loss of precision, the numerical instability and the situation that the inverse matrix does not exist. From the numerical examples, the results indicate that Gauss precise time-integration method, which used to solve point-reactor neutron kinetics equation, is more efficient and of high accuracy. It can be well applied to reactor safety analysis and simulation calculation.
