Paper Info
Title
A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane
Abstract
We consider the numerical treatment of second kind integral equations on the real line of the form $$\\phi (s) = \\psi (s) + \\int_{ - \\infty }^{ + \\infty } {K(s - t)z(t)} \\phi (t)dt, s \\in \\mathbb{R},$$ (abbreviated¿ =¿ +K z ¿) in which¿ ¿L 1(¿),z ¿L ¿ (¿), and¿ ¿BC(¿), the space of bounded continuous functions on ¿, are assumed known and¿ ¿BC(¿) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [¿A, A]) via bounds on (I ¿ K z )¿1 as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on ¿ is then analysed: in the case whenz is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases wherez is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that ifz (related to the boundary impedance in the application) takes values in an appropriate compact subsetQ of the complex plane, then the difference between¿(s) and its finite section approximation computed numerically using the iterative scheme proposed is ≤C 1[khlog(1/kh)+(1¿¿)¿1/2(kA)¿1/2] in the interval [¿¿A, ¿A] (¿<1), forkh sufficiently small, wherek is the wavenumber andh the grid spacing. Moreover this numerical approximation can be computed in ≤C 2 N logN operations, whereN = 2A/h is the number of degrees of freedom. The values of the constantsC 1 andC 2 depend only on the setQ and not on the wavenumberk or the support ofz.
Year
DOI
Venue
2002
10.1007/BF02679436
Numerische Mathematik
Keywords
Field
DocType
65r20,integral equation,collocation method,degree of freedom
Continuous function,Coefficient matrix,Mathematical analysis,Matrix (mathematics),Real line,Integral equation,Complex plane,Collocation method,Band matrix,Mathematics
Journal
Volume
Issue
ISSN
93
1
0945-3245
Citations 
PageRank 
References 
5
1.46
3
Authors
3
Name
Order
Citations
PageRank
Simon N. Chandler-Wilde111616.79
Mizanur Rahman212920.97
C. R. Ross351.46