Runge-Kutta convolution quadrature for operators arising in wave propagation
Lehel Banjai, Christian Lubich,and Jens Markus Melenk
Contact the author: Please use for correspondence this email.
Submission date: 11. Oct. 2010
published in: Numerische Mathematik, 119 (2011) 1, p. 1-20
DOI number (of the published article): 10.1007/s00211-011-0378-z
MSC-Numbers: 65R20, 65L06, 65M15
Keywords and phrases: convolution quadrature, Runge-Kutta methods, Time-domain boundary integral operators, order reduction
Download full preprint: PDF (202 kB)
An error analysis of Runge-Kutta convolution quadrature is presented for a class of non-sectorial operators whose Laplace transform satisfies, besides the standard assumptions of analyticity in a half-plane ℜs > σ0 and a polynomial bound O(sμ1) there, the stronger polynomial bound O(sμ2) in convex sectors of the form |arg s|≤ π∕2 - θ < π∕2 for θ > 0. The order of convergence of the Runge-Kutta convolution quadrature is determined by μ2 and the underlying Runge-Kutta method, but is independent of μ1.
Time domain boundary integral operators for wave propagation problems have Laplace transforms that satisfy bounds of the above type. Numerical examples from acoustic scattering show that the theory describes accurately the convergence behaviour of Runge-Kutta convolution quadrature for this class of applications. Our results show in particular that the full classical order of the Runge-Kutta method is attained away from the scattering boundary.