DefinePK

Abstract

In this work, a new numerical approximation to solve the inverse point kinetics equation is presented. The integral of the historic of the neutron population density is discretized, based on the approximation of Euler-Maclaurin formula, considering an extra term in this procedure to obtain a better approximation using two Bernoulli numbers. Reactivity is calculated for different numerical experiments, with different neutron population density forms and different time steps. In order to validate the accuracy of the proposed method, we compare the results obtained with different numerical experiments reported in the literature. The results suggest that the method shown can be used in a real time reactivity digital meter.