Exponential sum approximation for Mittag-Leffler function and its application to fractional Zener wave equation

Research output: Journal Publications and ReviewsRGC 21 - Publication in refereed journalpeer-review

8 Scopus Citations
View graph of relations

Related Research Unit(s)


Original languageEnglish
Article number109389
Journal / PublicationJournal of Computational Physics
Online published9 Mar 2020
Publication statusPublished - 1 Jun 2020


The Mittag-Leffler function (MLF) is fundamental to solving fractional calculus problems. An exponential sum approximation for the single-parameter MLF with negative input is proposed. Analysis shows that the approximation, which is based on the Gauss-Legendre quadrature, converges uniformly for all non-positive input. The application to modelling of wave propagation in viscoelastic material with the finite element method is also presented. The propagation speed of the solution to the approximated wave equation is proved to have the same upper bound as the original. The discretized scheme, based on the generalised alpha method, involves only a single matrix inverse whose size is the degree of freedom of the geometric model per time step. It is proved that the scheme is unconditionally stable. Furthermore, the solution converges to the true solution at (T), where T is the interval each time step, provided that the approximation error of the MLF is sufficiently small.

Research Area(s)

  • Fast convolution, Finite element method, Fractional calculus, Gaussian quadrature, Mittag-Leffler function, Stability analysis