Convolution spline approximations for time domain boundary integral equations

Penny J. Davies, Dugald Black Duncan

Research output: Contribution to journalArticlepeer-review

10 Citations (Scopus)
128 Downloads (Pure)


We derive a new \convolution spline" approximation method for convolution Volterra integral equations. This shares some properties of convolution quadrature, but instead of being based on an underlying ODE solver is explicitly constructed in terms of basis functions which have compact support. At time step tn = nh > 0, the solution is approximated in a \backward time" manner in terms of basis functions j by u(tn􀀀t) Pnj=0 un􀀀jj(t=h) for t 2 [0; tn]. We carry out a detailed analysis for B-spline basis functions, but note that the framework is more general than this. For B-splines of degree m 1 we show that the schemes converge at the rate O(h2) when the kernel is
suciently smooth. We also establish a methodology for their stability analysis and obtain new stability results for several non-smooth kernels, including the case of a highly oscillatory Bessel function kernel (in which the oscillation frequency can be O(1=h)). This is related to convergence analysis for approximation of time domain boundary integral equations (TDBIEs), and provides evidence that the new convolution spline approach could provide a useful time-stepping mechanism for TDBIE problems. In particular, using
compactly supported basis functions would give sparse system matrices.
Original languageEnglish
Pages (from-to)369-410
Number of pages42
JournalJournal of Integral Equations and Applications
Issue number3
Publication statusPublished - 31 Oct 2014


  • Convolution quadrature
  • Volterra integral equations
  • time dependent boundary integral equations


Dive into the research topics of 'Convolution spline approximations for time domain boundary integral equations'. Together they form a unique fingerprint.

Cite this