Nuclear Science and Engineering / Volume 25 / Number 2 / June 1966 / Pages 183-188
Technical Paper / dx.doi.org/10.13182/NSE66-A17735
Articles are hosted by Taylor and Francis Online.
A numerical solution of the reactor kinetics equations is given in terms of a difference equation involving an exponential matrix. It is shown that this equation yields a local discretinzation error of the order of the time step-size cubed. The exponential of the original equation is approximated by rational matrix functions and a bound for the error resulting from such approximations is established. Three particular rational matrix functions are proposed, and the two problems involving ramp changes in reactivity are solved using them.