Neutron Point Kinetics Model With A Distributed-Order Fractional Derivative

In this paper, the solutions of an extended form of the Fractional-order Neutron Point Kinetics (FNPK) equation in terms of Caputo-time derivatives of the same order are investigated. Instead of using a Caputo derivative, a distributed-order fractional derivative in the Caputo sense was employed in the term of the FNPK equation which is multiplied by the reactivity. This term plays an important role in the description of neutron kinetics during the start-up, shutdown, and steady-state processes in nuclear reactors. The extended (DFNPK) model was solved using the beta, normal, bimodal and Dirac delta distributions to investigate their effect on the transient state solutions of the neutron density. Regardless of the distribution used, the most significant finding is that a destabilizing effect on the neutron density is induced when the mode (or the instant of application of the Dirac delta) of the distribution tends to one while maintaining the orders of the Caputo-time derivatives constant. What defines the destabilizing effect are large magnitude oscillations, a rapid decay, and an oscillation-free steady state with a monotonic increase that is parallel to but somewhat above the trend determined by the FNPK equation. The extended model is anticipated to be effective for modeling neutron density dispersion in a highly heterogeneous medium that may be described using distributed derivatives.