Correction of high-order $L_k$ approximation for subdiffusion
Abstract
The subdiffusion equations with a Caputo fractional derivative of order arise in a wide variety of practical problems, which is describing the transport processes, in the force-free limit, slower than Brownian diffusion. In this work, we derive the correction schemes of the Lagrange interpolation with degree () convolution quadrature, called approximation, for the subdiffusion, which are easy to implement on variable grids. The key step of designing correction algorithm is to calculate the explicit form of the coefficients of approximation by the polylogarithm function or Bose-Einstein integral. To construct a approximation of Bose-Einstein integral, the desired th-order convergence rate can be proved for the correction scheme with nonsmooth data, which is higher than th-order BDF method in [Jin, Li, and Zhou, SIAM J. Sci. Comput., 39 (2017), A3129--A3152; Shi and Chen, J. Sci. Comput., (2020) 85:28]. The numerical experiments with spectral method are given to illustrate theoretical results.
Cite
@article{arxiv.2112.13609,
title = {Correction of high-order $L_k$ approximation for subdiffusion},
author = {Jiankang Shi and Minghua Chen and Yubin Yan and Jianxiong Cao},
journal= {arXiv preprint arXiv:2112.13609},
year = {2023}
}