Time-stepping error bounds for fractional diffusion problems with non-smooth initial data

We apply the piecewise constant, discontinuous Galerkin method to discretize a fractional diffusion equation with respect to time. Using Laplace transform techniques, we show that the method is first order accurate at the \$n\$th time level \$t_n\$, but the error bound includes a factor \$t_n^{-1}\$ if we assume no smoothness of the initial data. We also show that for smoother initial data the growth in the error bound as \$t_n\$ decreases is milder, and in some cases absent altogether. Our error bounds generalize known results for the classical heat equation and are illustrated for a model problem.Comment: 22 pages, 5 figure

An hp-version discontinuous Galerkin method for integro-differential equations of parabolic type

We study the numerical solution of a class of parabolic integro-differential equations with weakly singular kernels. We use an $hp$-version discontinuous Galerkin (DG) method for the discretization in time. We derive optimal $hp$-version error estimates and show that exponential rates of convergence can be achieved for solutions with singular (temporal) behavior near $t=0$ caused by the weakly singular kernel. Moreover, we prove that by using nonuniformly refined time steps, optimal algebraic convergence rates can be achieved for the $h$-version DG method. We then combine the DG time-stepping method with a standard finite element discretization in space, and present an optimal error analysis of the resulting fully discrete scheme. Our theoretical results are numerically validated in a series of test problems

Numerical solution for a sub-diffusion equation with a smooth kernel

AbstractIn this paper we study the numerical solution of an initial value problem of a sub-diffusion type. For the time discretization we apply the discontinuous Galerkin method and we use continuous piecewise finite elements for the space discretization. Optimal order convergence rates of our numerical solution have been shown. We compare our theoretical error bounds with the results of numerical computations. We also present some numerical results showing the super-convergence rates of the proposed method

Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations

We consider an initial-boundary value problem for $\partial_tu-\partial_t^{-\alpha}\nabla^2u=f(t)$, that is, for a fractional diffusion ($-1<\alpha<0$) or wave ($0<\alpha<1$) equation. A numerical solution is found by applying a piecewise-linear, discontinuous Galerkin method in time combined with a piecewise-linear, conforming finite element method in space. The time mesh is graded appropriately near $t=0$, but the spatial mesh is quasiuniform. Previously, we proved that the error, measured in the spatial $L_2$-norm, is of order $k^{2+\alpha_-}+h^2\ell(k)$, uniformly in $t$, where $k$ is the maximum time step, $h$ is the maximum diameter of the spatial finite elements, $\alpha_-=\min(\alpha,0)\le0$ and $\ell(k)=\max(1,|\log k|)$. Here, we generalize a known result for the classical heat equation (i.e., the case $\alpha=0$) by showing that at each time level $t_n$ the solution is superconvergent with respect to $k$: the error is of order $(k^{3+2\alpha_-}+h^2)\ell(k)$. Moreover, a simple postprocessing step employing Lagrange interpolation yields a superconvergent approximation for any $t$. Numerical experiments indicate that our theoretical error bound is pessimistic if $\alpha<0$. Ignoring logarithmic factors, we observe that the error in the DG solution at $t=t_n$, and after postprocessing at all $t$, is of order $k^{3+\alpha_-}+h^2$.Comment: 24 pages, 2 figure