Numerical solution of an evolution equation with a positive-type memory term

Abstract

We study the numerical solution of an initial-boundary value problem for a Volterra type integro-differential equation, in which the integral operator is a convolution product of a positive-definite kernel and an elliptic partial-differential operator. The equation is discretised in space by the Galerkin finite-element method and in time by finite differences in combination with various quadrature rules which preserve the positive character of the memory term. Special attention is paid to the case of a weakly singular kernel. Error estimates are derived and numerical experiments reported.