Error estimates for approximations of distributed order time fractional diffusion with nonsmooth data

Abstract

In this work, we consider the numerical solution of a distributed order subdiffusion model, arising in the modeling of ultra-slow diffusion processes. We develop a space semidiscrete scheme based on the Galerkin finite element method, and establish error estimates optimal with respect to data regularity in L2(Ω) and H1(Ω) norms for both smooth and nonsmooth initial data. Further, we propose two fully discrete schemes, based on the Laplace transform and convolution quadrature generated by the backward Euler method, respectively, and provide optimal L2(Ω) error estimates, which exhibits exponential convergence and first-order convergence in time, respectively. Extensive numerical experiments are provided to verify the error estimates for both smooth and nonsmooth initial data.