Numerical solution of two-dimensional mixed problems with variable coefficients by the boundary-domain integral and integro-differential equation methods

This paper presents a formulation of the boundary-domain integral equation (BDIE) and the boundary-domain integro-differential equation (BDIDE) methods for the numerical solution of two-dimensional mixed boundary-value problems (BVP) for a second-order linear elliptic partial differential equation (PDE) with variable coefficients. The methods use a specially constructed parametrix (Levi function) to reduce the BVP to a BDIE or BDIDE. The numerical formulation of the BDIDE employs an approximation for the boundary fluxes in terms of the potential function within the domain cells; therefore, the solution is fully described in terms of the variation of the potential function along the boundary and domain. Linear basis functions localised on triangular elements and standard quadrature rules are used for the calculation of boundary integrals. For the domain integrals, we have implemented Gaussian quadrature rules for two dimensions with Duffy transformation, by mapping the triangles into squares and eliminating the weak singularity in the process. Numerical examples are presented for several simple problems with square and circular domains, for which exact solutions are available. It is shown that the present method produces accurate results even with coarse meshes. The numerical results also show that high rates of convergence are obtained with mesh refinement.