Analytical solution of the fractional neutron point kinetic equations using the Mittag-Leffler function

Published: 23 August 2023| Version 1 | DOI: 10.17632/pzfbj64dvg.1


A new analytical solution of the fractional point kinetic equations is developed in the present work, which considers one group of delayed neutrons as well as a constant reactivity. One of the main novelties and theoretical contributions of the developed solution is that it does not require approximating the Inverse Laplace transform by numerical means, as other solutions do, and instead it is obtained analytically using the Green and the Mittag-Leffler functions. A set of algorithms and MATLAB codes are developed with the purpose to implement and compute the developed solution in an adequate way, which represents the major computational contribution. Numerical experiments show that the developed solution is very accurate solving short-time dynamics phenomena, being able to reproduce, in a precise way, data that is reported in literature. The development of the proposed analytic solution represents an important contribution to validate numerical methods that are currently used to solve the neutron fractional point kinetic equations.


Steps to reproduce

A brief user´s manual can be found in the following GitHub repository and in the submitted paper.


Universidad Autonoma Metropolitana Iztapalapa


Laplace Transformation, Nuclear Reactor, Fractional Calculus


Consejo Nacional de Ciencia y Tecnología