An $hp$-Version Discontinuous Galerkin Method for Integro-Differential Equations of Parabolic Type

Authors:
K. Mustapha;H. Brunner;H. Mustapha;D. Schötzau
Affiliations:
kassem@kfupm.edu.sa;hbrunner@math.hkbu.edu.hk;hussein.mustapha@mcgill.ca;schoetzau@math.ubc.ca
Venue:
SIAM Journal on Numerical Analysis
Year:
2011

Citing 9
Cited 1

Time discretization of an integro-differential equation of parabolic type

SIAM Journal on Numerical Analysis
Finite element methods for parabolic and hyperbolic partial integro-differential equations

Non-Linear Analysis
A posteriori error bounds and global error control for approximation of ordinary differential equations

SIAM Journal on Numerical Analysis
Discretization with variable time steps of an evolution equation with a positive-type memory term

Journal of Computational and Applied Mathematics
Numerical solution of parabolic integro-differential equations by the discontinuous Galerkin method

Mathematics of Computation
Time Discretization of Parabolic Problems by the HP-Version of the Discontinuous Galerkin Finite Element Method

SIAM Journal on Numerical Analysis
hp-Discontinuous Galerkin Time-Stepping for Volterra Integrodifferential Equations

SIAM Journal on Numerical Analysis
A second-order accurate numerical method for a fractional wave equation

Numerische Mathematik
Mixed Finite Element Approximations of Parabolic Integro-Differential Equations with Nonsmooth Initial Data

SIAM Journal on Numerical Analysis

Post-processing discontinuous Galerkin solutions to Volterra integro-differential equations: Analysis and simulations

Journal of Computational and Applied Mathematics

Quantified Score

Hi-index	0.01

Visualization

Abstract

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.