English

A discontinuous Galerkin method for time fractional diffusion equations with variable coefficients

Numerical Analysis 2015-11-03 v1

Abstract

We propose a piecewise-linear, time-stepping discontinuous Galerkin method to solve numerically a time fractional diffusion equation involving Caputo derivative of order μ(0,1)\mu\in (0,1) with variable coefficients. For the spatial discretization, we apply the standard piecewise linear continuous Galerkin method. Well-posedness of the fully discrete scheme and error analysis will be shown. For a time interval~(0,T)(0,T) and a spatial domain~Ω\Omega, our analysis suggest that the error in L2((0,T),L2(Ω))L^2\bigr((0,T),L^2(\Omega)\bigr)-norm is of order O(k2μ2+h2)O(k^{2-\frac{\mu}{2}}+h^2) (that is, short by order μ2\frac{\mu}{2} from being optimal in time) where kk denotes the maximum time step, and hh is the maximum diameter of the elements of the (quasi-uniform) spatial mesh. However, our numerical experiments indicate optimal O(k2+h2)O(k^{2}+h^2) error bound in the stronger L((0,T),L2(Ω))L^\infty\bigr((0,T),L^2(\Omega)\bigr)-norm. Variable time steps are used to compensate the singularity of the continuous solution near t=0t=0.

Keywords

Cite

@article{arxiv.1511.00163,
  title  = {A discontinuous Galerkin method for time fractional diffusion equations with variable coefficients},
  author = {K. Mustapha and B. Abdallah and K. M. Furati and M. Nour},
  journal= {arXiv preprint arXiv:1511.00163},
  year   = {2015}
}
R2 v1 2026-06-22T11:33:53.076Z