Camille Carvalho

Camille Carvalho, Shilpa Khatri, Arnold D Kim

We study the evaluation of layer potentials close to the domain boundary. Accurate evaluation of layer potentials near boundaries is needed in many applications, including fluid-structure interactions and near-field scattering in nano-optics. When numerically evaluating layer potentials, it is natural to use the same quadrature rule as the one used in the Nyström method to solve the underlying boundary integral equation. However, this method is problematic for evaluation points close to boundaries. For a fixed number of quadrature points, $N$, this method incurs $O(1)$ errors in a boundary layer of thickness $O(1/N)$. Using an asymptotic expansion for the kernel of the layer potential, we remove this $O(1)$ error. We demonstrate the effectiveness of this method for interior and exterior problems for Laplace's equation in two dimensions.

Camille Carvalho, Shilpa Khatri, Arnold D. Kim

When using the boundary integral equation method to solve a boundary value problem, the evaluation of the solution near the boundary is challenging to compute because the layer potentials that represent the solution are nearly-singular integrals. To address this close evaluation problem, we apply an asymptotic analysis of these nearly singular integrals and obtain an asymptotic approximation. We derive the asymptotic approximation for the case of the double-layer potential in two and three dimensions, representing the solution of the interior Dirichlet problem for Laplace's equation. By doing so, we obtain an asymptotic approximation given by the Dirichlet data at the boundary point nearest to the interior evaluation point plus a nonlocal correction. We present numerical methods to compute this asymptotic approximation, and we demonstrate the efficiency and accuracy of the asymptotic approximation through several examples. These examples show that the asymptotic approximation is useful as it accurately approximates the close evaluation of the double-layer potential while requiring only modest computational resources.

Anne-Sophie Bonnet-Ben Dhia, Camille Carvalho, Lucas Chesnel et al.

We investigate in a $2$D setting the scattering of time-harmonic electromagnetic waves by a plasmonic device, represented as a non dissipative bounded and penetrable obstacle with a negative permittivity. Using the $\textrm{T}$-coercivity approach, we first prove that the problem is well-posed in the classical framework $H^1_{\text{loc}} $ if the negative permittivity does not lie in some critical interval whose definition depends on the shape of the device. When the latter has corners, for values inside the critical interval, unusual strong singularities for the electromagnetic field can appear. In that case, well-posedness is obtained by imposing a radiation condition at the corners to select the outgoing black-hole plasmonic wave, that is the one which carries energy towards the corners. A simple and systematic criterion is given to define what is the outgoing solution. Finally, we propose an original numerical method based on the use of Perfectly Matched Layers at the corners. We emphasize that it is necessary to design an $\textit{ad hoc}$ technique because the field is too singular to be captured with standard finite element methods.