Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model

Abstract

In the last few years, reaction–diffusion models associated with discrete fractional calculus have risen in prominence in scientific fields, not just due to the requirement for numerical simulation but also due to the described biological phenomena. This work investigates a discrete equivalent of the fractional reaction–diffusion glycolysis model. The discrete fractional calculus tool is introduced to the discrete modeling of diffusion problems in the Caputo-like delta sense, and a fractional discretization diffusion model is described. The local stability of the equilibrium points in the proposed discrete system is examined. We additionally investigate the global stability of the equilibrium point by developing a Lyapunov function. Furthermore, this study indicates that the L1 finite difference scheme and the second-order central difference scheme can successfully preserve the characteristics of the associated continuous system. Finally, an equivalent summation representing the model’s numerical formula is shown. The diffusion concentration is further investigated for different fractional orders, and examples with simulations are presented to corroborate the theoretical findings.