Numerical solution of a generalized boundary value problem for the modified Helmholtz equation in two dimensions

We propose numerical schemes for solving the boundary value problem for the modified Helmholtz equation and generalized impedance boundary condition. The approaches are based on the reduction of the problem to the boundary integral equation with a hyper-singular kernel. In the first scheme the hyper-singular integral operator is treated by splitting off the singularity technique whereas in the second scheme the idea of numerical differentiation is employed. The solvability of the boundary integral equation and convergence of the first method are established. Exponential convergence for analytic data is exhibited by numerical examples. (C) 2021 International Association for Mathematics and Computers in Simulation (IMACS). Published by Elsevier B.Y. All rights reserved.