NT & RP Journal
SOLUTION OF THE REACTOR POINT KINETICS EQUATIONS BY MATLAB COMPUTING
......: info......: history......: editorial......: archive......: for authors......: subscription

Vol. XXX, No. 1, Pp. 01-82
March 2015
UDC 621.039+614.876:504.06
ISSN 1451-3994

Back to Contents

 


Pages: 11-17

Authors:
Sudhansu S. Singh and Dinakrushna Mohapatra

Abstract


The numerical solution of the point kinetics equations in the presence of Newtonian temperature feedback has been a challenging issue for analyzing the reactor transients. Reactor point kinetics equations are a system of stiff ordinary differential equations which need special numerical treatments. Although a plethora of numerical intricacies have been introduced to solve the point kinetics equations over the years, some of the simple and straightforward methods still work very efficiently with extraordinary accuracy. As an example, it has been shown recently that the fundamental backward Euler finite difference algorithm with its simplicity has proven to be one of the most effective legacy methods. Complementing the backward Euler finite difference scheme, the present work demonstrates the application of ordinary differential equation suite available in the MATLAB software package to solve the stiff reactor point kinetics equations with Newtonian temperature feedback effects very effectively by analyzing various classic benchmark cases. Fair accuracy of the results implies the efficient application of MATLAB ordinary differential equation suite for solving the reactor point kinetics equations as an alternate method for future applications.

Key words: point kinetics equation, backward Euler finite difference, MATLAB, ODE suite, NDF

FULL PAPER IN PDF FORMAT (1.28 MB)

Vinča Institute of Nuclear Sciences :: Designed by milas :: July 2007
Operated by acapanic :: Last updated on May, 2015