Numerical stability of Grunwald-Letnikov method for time fractional delay differential equations

This paper is concerned with the numerical stability of time fractional delay differential equations (F-DDEs) based on Grunwald-Letnikov (GL) approximation for the Caputo fractional derivative. In particular, we focus on the numerical stability region and the Mittag-Leffler stability. Using the boundary locus technique, we first derive the exact expression of the numerical stability region in the parameter plane, and show that the fractional backward Euler scheme is not tau(0)-stable, which is different from the backward Euler scheme for integer DDE models. Secondly, we prove the numerical Mittag-Leffler stability for the numerical solutions provided that the parameters fall into the numerical stability region, by employing the singularity analysis of generating function. Our results show that the numerical solutions of F-DDEs are completely different from the classical integer order DDEs, both in terms of tau(0)-stability and the long-time decay rate. Numerical examples are given to confirm the theoretical results.

Automated summary

This paper investigates the numerical stability of F-DDEs based on the Grunwald-Letnikov approximation for the Caputo fractional derivative, focusing on the numerical stability region and Mittag-Leffler stability. Using the boundary locus technique, the exact expression of the numerical stability region in the parameter plane is derived. The study shows that the numerical solutions of F-DDEs differ significantly from classical integer order DDEs in terms of stability and long-time decay rate.