Abstract
We consider biperiodic weakly singular integral equations of the second kind such as they arise in boundary integral equation methods. These are solved numerically using a collocation scheme based on trigonometric polynomials. The singularity is removed by a local change to polar coordinates. The resulting operators have smooth kernels and are discretized using the tensor product composite trapezoidal rule. We prove stability and convergence of the scheme achieving algebraic convergence of any order under appropriate regularity assumptions. The method can be applied to typical boundary value problems such as potential and scattering problems both for bounded obstacles and for periodic surfaces. We present numerical results demonstrating that the expected convergence rates can be observed in practice.
Similar content being viewed by others
References
Arens, T.: Scattering by biperiodic layered media: the integral equation approach (2010). Habilitation Thesis, Karlsruhe Institute of Technology
Arens, T., Lechleiter, A., Sandfort, K., Schmitt, S.: Analysing Ewald’s method for the evaluation of Green’s functions for periodic media. IMA J. Appl. Math. 78, 405–431 (2013). doi:10.1093/imamat/hxr057
Beale, J.T.: A grid-based boundary integral method for elliptic problems in three dimensions. SIAM J. Numer. Anal. 42, 599–620 (2004)
Bruno, O., Dominguez, V., Sayas, F.: Convergence analysis of a high-order Nystrom integral-equation method for surface scattering problems. Numer. Math. 124, 603–645 (2013)
Bruno, O.P., Kunyanski, L.A.: A fast high order algorithm for the solution of surface scattering problems: basic implementation, tests and applications. J. Comput. Phys. 169, 80–110 (2001)
Bruno, O.P., Kunyanski, L.A.: Surface scattering in three dimensions: an accelerated high-order solver. Proc. R. Soc. Lond. A 457, 2921–2934 (2001)
Canino, L., Ottusch, J., Stalzer, M., Visher, J., Wandzura, S.M.: Numerical solution of the Helmholtz equation in 2D and 3D using a high-order nystroem discretization. J. Comput. Phys. 146, 627–663 (1998)
Ganesh, M., Graham, I.: A high-order algorithm for obstacle scattering in three dimensions. J. Comput. Phys. 198, 211–242 (2004)
Graham, I., Sloan, I.: Fully discrete spectral boundary integral methods for Helmholtz problems on smooth closed surfaces in \({\mathbb{R}}^3\). Numer. Math. 92(2), 289–323 (2002). doi:10.1007/s002110100343
Heinemeyer, E.: Integral equation methods for rough surface scattering problems in three dimensions. Ph.D. thesis, Georg-August-Universität Göttingen (2008)
Kress, R.: Linear Integral Equations, 2nd edn. Springer, Berlin (1999)
Kussmaul, R.: Ein numerisches Verfahren zur Lösung des Neumannschen Außenraumproblems für die Helmholtzsche Schwingungsgleichung. Computing 4, 246–273 (1969)
Martensen, E.: Über eine Methode zum räumlichen Neumannschen Problem mit einer Anwendung fr torusartige Berandungen. Acta Math. 109(1), 75–135 (1963). doi:10.1007/BF02391810
Saranen, J., Vainikko, G.: Periodic Integral and Pseudodifferential Equations with Numerical Approximation, Springer Monographs in Mathematics. Springer, Berlin (2002)
Wienert, L.: Die numerische Approximation von Randintegraloperatoren für die Helmholtzgleichung im \({\mathbb{R}}^3\). Ph.D. thesis, Univ. Göttingen (1990)
Acknowledgments
This research was supported by the Deutsche Forschungsgemeinschaft (DFG) under the grant Super-Algebraisch konvergente numerische Verfahren für biperiodische Integralgleichungen.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Arens, T., Rösch, T. A collocation method for weakly singular integral equations with super-algebraic convergence rate. Numer. Math. 134, 441–472 (2016). https://doi.org/10.1007/s00211-015-0783-9
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-015-0783-9