Time discretization of an evolution equation via Laplace transforms
Artikel i vetenskaplig tidskrift, 2004

Following earlier work by Sheen, Sloan, and Thomée concerning parabolic equations we study the discretization in time of a Volterra type integro-differential equation in which the integral operator is a convolution of a weakly singular function and an elliptic differential operator in space. The time discretization is accomplished by using a modified Laplace transform in time to represent the solution as an integral along a smooth curve extending into the left half of the complex plane, which is then evaluated by quadrature. This reduces the problem to a finite set of elliptic equations with complex coefficients, which may be solved in parallel. Stability and error bounds of high order are derived for two different choices of the quadrature rule. The method is combined with finite-element discretization in the spatial variables.

memory term

quadrature error

Laplace transform

evolution equation

parallel algorithm

Författare

William McLean

University of New South Wales (UNSW)

Vidar Thomee

Chalmers, Institutionen för matematik

Göteborgs universitet

IMA Journal of Numerical Analysis

0272-4979 (ISSN) 1464-3642 (eISSN)

Vol. 24 439-463

Ämneskategorier

Beräkningsmatematik

DOI

10.1093/imanum/24.3.439