A refined way of solving reactor point kinetics equations for imposed reactivity insertions

We apply the concept of convergence acceleration, also known as extrapolation, to find the solution of the reactor kinetics equations (RKEs). The method features simplicity in that an approximate finite difference formulation is constructed and converged to high accuracy from knowledge of the error term. Through the Romberg extrapolation, we demonstrate its high accuracy for a variety of imposed reactivity insertions found in the literature. The unique feature of the proposed algorithm, called RKE/R(omberg), is that no special attention is given to the stiffness of the RKEs. Finally, because of its simplicity and accuracy, the RKE/R algorithm is arguably the most efficient numerical solution of the RKEs developed to date.