Numerical solution analysis of fractional point kinetics and heat exchange in nuclear reactor

The paper presents the neutron point kinetics and heat exchange models for the nuclear reactor. The models consist of a nonlinear system of fractional ordinary differential and algebraic equations. Two numerical algorithms are used to solve them. The first algorithm is application of discrete Grünwald–Letnikov definition of the fractional derivative in the model. The second involves building an analog scheme in the FOMCON Toolbox in MATLAB environment where Oustaloup recursive filter is defined. The simulation results for input step reactivity are discussed and compared. The impact of the model parameters is examined.