The adjoint double layer potential on smooth surfaces in ℝ³ and the Neumann problem

2023/09/29 by J. Thomas Beale, Beale, J. Thomas, Michael Storm +3
Computer Science · Mathematics · Physics and Astronomy · #31B10 #35J25 #65D30 #65R20 #Advanced Mathematical Modeling in Engineering #Electromagnetic Scattering and Analysis #FOS: Mathematics #Numerical Analysis (math.NA) #Numerical methods in inverse problems

paper · pdf · doi:10.48550/arxiv.2310.00188

openalex publication_date 2023/09/29 · openalex created_date 2023/10/04 · openalex updated_date 2026/07/28

Abstract

We present a simple yet accurate method to compute the adjoint double layer potential, which is used to solve the Neumann boundary value problem for Laplace's equation in three dimensions. An expansion in curvilinear coordinates leads us to modify the expression for the adjoint double layer so that the singularity is reduced when evaluating the integral on the surface. We then regularize the Green's function, with a radial parameter δ. We show that a natural regularization has error O(δ³), and a simple modification improves the error to O(δ⁵). The integral is evaluated numerically without the need of special coordinates. We use this treatment of the adjoint double layer to solve the classical integral equation for the interior Neumann problem and evaluate the solution on the boundary. Choosing δ= ch^4/5, we find about O(h⁴) convergence in our examples, where h is the spacing in a background grid.

Related

• Analytical Galerkin boundary integrals of Laplace kernel layer potentials in ℝ³
• About Three Dimensional Double-Sided Dirichlet and Neumann Boundary Value Problems for the Laplacian
• Recursive Analytical Quadrature of Laplace and Helmholtz Layer Potentials in ℝ³
• Close evaluation of layer potentials in three dimensions
• High-order close evaluation of Laplace layer potentials: A differential geometric approach