The three dimensional Neumann Green's function for general surfaces: singular asymptotics and boundary integral methods

Abstract

We present an asymptotic analysis and high-order boundary integral method for the three-dimensional Neumann Green's function in general geometries. The Neumann Green's function is a fundamental quantity which arises in numerous fields of science and engineering. In the application of singular perturbation methods to strongly localized reactions and diffusive transport, the Green's function plays the key role in mediating global dynamics. However, this essential quantity can only be determined in closed form for a limited set of geometries. The Green's function for the Laplacian is an elliptic problem with a Dirac forcing term. Accurate resolution of the solution requires a careful decomposition into a singular and a regular part. The bulk scenario is where the source is placed off surface and the singularity is given by the free-space function. In the surface case, where the source is placed at a curved point on the boundary, we use asymptotic analysis to determine a three-term singularity structure. With explicit knowledge of these singularities, we develop a high-order boundary integral method for the determination of the remaining regular part. To resolve the singular boundary data, our integral method uses a custom discretization with Duffy patches near the source. We validate our method using several test cases in which closed form solutions can be developed, including spheres, prolate spheroids and constructed domains. We demonstrate the applicability of our method to address some open problems in narrow capture theory.