TY - JOUR
T1 - Exponential sum approximation for Mittag-Leffler function and its application to fractional Zener wave equation
AU - Lam, P. H.
AU - So, H. C.
AU - Chan, C. F.
PY - 2020/6/1
Y1 - 2020/6/1
N2 - 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 O (T2 ), where T is the interval each time step, provided that the approximation error of the MLF is sufficiently small.
AB - 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 O (T2 ), where T is the interval each time step, provided that the approximation error of the MLF is sufficiently small.
KW - Fast convolution
KW - Finite element method
KW - Fractional calculus
KW - Gaussian quadrature
KW - Mittag-Leffler function
KW - Stability analysis
UR - http://www.scopus.com/inward/record.url?scp=85081697907&partnerID=8YFLogxK
UR - https://www.scopus.com/record/pubmetrics.uri?eid=2-s2.0-85081697907&origin=recordpage
U2 - 10.1016/j.jcp.2020.109389
DO - 10.1016/j.jcp.2020.109389
M3 - 21_Publication in refereed journal
VL - 410
JO - Journal of Computational Physics
JF - Journal of Computational Physics
SN - 0021-9991
M1 - 109389
ER -