Communications on Applied Mathematics and Computation ›› 2025, Vol. 7 ›› Issue (3): 865-884.doi: 10.1007/s42967-024-00401-1

Previous Articles     Next Articles

Sharp Error Analysis for Averaging Crank-Nicolson Schemes with Corrections for Subdiffusion with Nonsmooth Solutions

Baoli Yin1, Yang Liu1, Hong Li1   

  1. School of Mathematical Sciences, Inner Mongolia University, Hohhot 010021, Inner Mongolia, China
  • Received:2023-09-30 Revised:2024-02-16 Accepted:2024-03-15 Online:2025-09-20 Published:2025-05-23
  • Supported by:
    This work is supported by the Autonomous Region Level High-Level Talent Introduction Research Support Program in 2022 (No. 12000-15042224 to B.Y.), the National Natural Science Foundation of China (No. 12201322 to B.Y., 12061053 to Y.L., and 12161063 to H.L.), and the Natural Science Foundation of Inner Mongolia, China (No. 2020MS01003 to Y.L. and 2021MS01018 to H.L.).

Abstract: Thanks to the singularity of the solution of linear subdiffusion problems, most time-stepping methods on uniform meshes can result in O(τ) accuracy where τ denotes the time step. The present work aims to discover the reason why some type of Crank-Nicolson schemes (the averaging Crank-Nicolson(ACN) scheme)for the subdiffusion can only yield O(τ α) (α < 1) accuracy, which is much lower than the desired. The existing well-developed error analysis for the subdiffusion, which has been successfully applied to many time-stepping methods such as the fractional BDF-p (1 ≤ p ≤ 6), requires singular points to be out of the path of contour integrals involved. The ACN scheme in this work is quite natural but fails to meet this requirement. By resorting to the residue theorem, some novel sharp error analysis is developed in this study, upon which correction methods are further designed to obtain the optimal O(τ 2O(τ 2) accuracy. All results are verified by numerical tests.

Key words: Subdiffusion, Uniform mesh, Crank-Nicolson scheme, Convolution quadrature, Singularity

CLC Number: