Numerical solution and spectrum of boundary-domain integral equation for the Neumann BVP with variable coefficient

Abstract

In this paper, a numerical implementation of a direct united boundary-domain integral equation (BDIE) related to the Neumann boundary value problem for a scalar elliptic partial differential equation with a variable coefficient is discussed. The BDIE is reduced to a uniquely solvable one by adding an appropriate perturbation operator. The mesh-based discretization of the BDIEs with quadrilateral domain elements leads to a system of linear algebraic equations (discretized BDIE). Then, the system is solved by LU decomposition and Neumann iterations. Convergence of the iterative method is discussed in relation to the distribution of eigenvalues of the corresponding discrete operators calculated numerically.