摘要

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 L-2(Omega) and H-1(Omega) 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 L-2(Omega) 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.

  • 出版日期2016-2