Corrected Trapezoidal Rules for Boundary Integral Equations in Three Dimensions

by   Bowei Wu, et al.

The manuscript describes a quadrature rule that is designed for the high order discretization of boundary integral equations (BIEs) using the Nyström method. The technique is designed for surfaces that can naturally be parameterized using a uniform grid on a rectangle, such as deformed tori, or channels with periodic boundary conditions. When a BIE on such a geometry is discretized using the Nyström method based on the Trapezoidal quadrature rule, the resulting scheme tends to converge only slowly, due to the singularity in the kernel function. The key finding of the manuscript is that the convergence order can be greatly improved by modifying only a very small number of elements in the coefficient matrix. Specifically, it is demonstrated that by correcting only the diagonal entries in the coefficient matrix, O(h^3) convergence can be attained for the single and double layer potentials associated with both the Laplace and the Helmholtz kernels. A nine-point correction stencil leads to an O(h^5) scheme. The method proposed can be viewed as a generalization of the quadrature rule of Duan and Rokhlin, which was designed for the 2D Lippmann-Schwinger equation in the plane. The techniques proposed are supported by a rigorous error analysis that relies on Wigner-type limits involving the Epstein zeta function and its parametric derivatives.



There are no comments yet.


page 1

page 2

page 3

page 4


Zeta Correction: A New Approach to Constructing Corrected Trapezoidal Quadrature Rules for Singular Integral Operators

We introduce a new quadrature method for the discretization of boundary ...

High order corrected trapezoidal rules for a class of singular integrals

We present a family of high order trapezoidal rule-based quadratures for...

Nyström methods for high-order CQ solutions of the wave equation in two dimensions

We investigate high-order Convolution Quadratures methods for the soluti...

About Three Dimensional Double-Sided Dirichlet and Neumann Boundary Value Problems for the Laplacian

The orthogonality of Hilbert spaces whose elements can be represented as...

Kernel Aggregated Fast Multipole Method: Efficient summation of Laplace and Stokes kernel functions

Many different simulation methods for Stokes flow problems involve a com...

Corrected trapezoidal rules for singular implicit boundary integrals

We present new higher-order quadratures for a family of boundary integra...

High-order Corrected Trapezoidal Rules for Functions with Fractional Singularities

In this paper, we introduce and analyze a high-order quadrature rule for...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.