Search

Reaction–diffusion systems have a broad variety of applications, particularly in biology, and it is well known that fractional calculus has been successfully used with this type of system. However, analyzing these systems using discrete fractional calculus is novel and requires significant research in a diversity of disciplines. Thus, in this paper, we investigate the discrete-time fractional-order Lengyel–Epstein system as a model of the chlorite iodide malonic acid (CIMA) chemical reaction. With the help of the second order difference operator, we describe the fractional discrete model. Furthermore, using the linearization approach, we established acceptable requirements for the local asymptotic stability of the system's unique equilibrium. Moreover, we employ a Lyapunov functional to show that when the iodide feeding rate is moderate, the constant equilibrium solution is globally asymptotically stable. Finally, numerical models are presented to validate the theoretical conclusions and demonstrate the impact of discretization and fractional-order on system dynamics. The continuous version of the fractional-order Lengyel–Epstein reaction–diffusion system is compared to the discrete-time system under consideration. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]