Global convergence and local superconvergence of first-kind Volterra integral equation approximations

We present a comprehensive convergence analysis for discontinuous piecewise polynomial approximations of a first-kind Volterra integral equation with smooth convolution kernel, examining the attainable order of (super-) convergence in collocation, quadrature discontinuous Galerkin (QDG) and full dis...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:IMA journal of numerical analysis 2012-07, Vol.32 (3), p.1117-1146
Hauptverfasser: Brunner, H., Davies, P. J., Duncan, D. B.
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We present a comprehensive convergence analysis for discontinuous piecewise polynomial approximations of a first-kind Volterra integral equation with smooth convolution kernel, examining the attainable order of (super-) convergence in collocation, quadrature discontinuous Galerkin (QDG) and full discontinuous Galerkin (DG) methods. We introduce new polynomial basis functions with properties that greatly simplify the convergence analysis for collocation methods. This also enables us to determine explicit formulae for the location of superconvergence points (i.e., discrete points at which the convergence order is one higher than the global bound) for all convergent collocation schemes. We show that a QDG method, which is based on piecewise polynomials of degree m and uses exactly m + 1 quadrature points and nonzero quadrature weights, is equivalent to a collocation scheme, and so its convergence properties are fully determined by the previous collocation analysis and they depend only on the quadrature point location (in particular, they are completely independent of the accuracy of the quadrature rule). We also give a complete analysis for QDG with more than m + 1 quadrature points when the degree of precision (d.o.p.) is at least 2m + 1. The behaviour (but not the approximation) is the same as that for a DG scheme when the d.o.p. is at least 2m + 2. Numerical test results confirm that the theoretical convergence rates are optimal.
ISSN:0272-4979
1464-3642
DOI:10.1093/imanum/drr029