English

Correction of high-order $L_k$ approximation for subdiffusion

Numerical Analysis 2023-06-27 v1 Numerical Analysis

Abstract

The subdiffusion equations with a Caputo fractional derivative of order α(0,1)\alpha \in (0,1) 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 kk (k6k\leq 6) convolution quadrature, called LkL_k 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 LkL_k approximation by the polylogarithm function or Bose-Einstein integral. To construct a τ8\tau_8 approximation of Bose-Einstein integral, the desired (k+1α)(k+1-\alpha)th-order convergence rate can be proved for the correction LkL_k scheme with nonsmooth data, which is higher than kkth-order BDFkk 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.

Keywords

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}
}
R2 v1 2026-06-24T08:32:24.649Z