An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces

Johan Helsing,Anders Karlsson

DOI: https://doi.org/10.1016/j.jcp.2014.04.053

2016-10-05

Abstract:A high-order accurate, explicit kernel-split, panel-based, Fourier-Nyström discretization scheme is developed for integral equations associated with the Helmholtz equation in axially symmetric domains. Extensive incorporation of analytic information about singular integral kernels and on-the-fly computation of nearly singular quadrature rules allow for very high achievable accuracy, also in the evaluation of fields close to the boundary of the computational domain.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: to develop a high - order accurate, explicit kernel - splitting panel - based Nyström discretization scheme for solving integral equations related to the Helmholtz equation in axisymmetric domains. Specifically, the goals of the paper include:

1. **Improve numerical accuracy**: By introducing analytical information and real - time calculation of almost - singular integration rules, extremely high achievable accuracy can be achieved, especially in field evaluation near the boundary of the computational domain.

2. **Handle nearly - singular and weakly - singular oscillating kernels**: Compared with previous work, this method uses 16 - order analytical product integration and the transformed kernel - splitting technique instead of the pre - calculated Kolm - Rokhlin integral formula. This method enables more accurate handling of near - singular situations.

3. **Solve specific physical problems**: This method is particularly suitable for the interior Neumann problem and is of great significance for solving Neumann Laplace eigenfunctions with homogeneous boundary conditions (such as acoustic resonance problems). In addition, it also involves problems such as electromagnetic scattering.

4. **Provide a flexible and efficient computational framework**: By calculating weights in real - time in the post - processor, it provides the ability to perform extremely accurate field evaluations throughout the computational domain, and can maintain high accuracy even near surfaces where ordinary integral equation techniques encounter difficulties.

5. **Combine multiple computational techniques**: The paper comprehensively uses multiple computational techniques, such as Fourier transform, convolution, semi - analytical methods, etc., to ensure an efficient and accurate discretization scheme.

In summary, this paper aims to provide an efficient and high - precision numerical solution method for integral equations of the Helmholtz equation in axisymmetric domains, especially suitable for physical problems requiring high - precision field evaluation.