Solution of the Poisson equation by the boundary integral method

Purpose

The boundary integral method (BIM) is very attractive to practicing engineers as it reduces the dimensionality of the problem by one, thereby making the procedure computationally inexpensive compared to its peers. The principal feature of this technique is the limitation of all its computations to only the boundaries of the domain. Although the procedure is well developed for the Laplace equation, the Poisson equation offers some computational challenges. Nevertheless, the literature provides a couple of solution methods. This paper revisits an alternate approach that has not gained much traction within the community. The purpose of this paper is to address the main bottleneck of that approach in an effort to popularize it and critically evaluate the errors introduced into the solution by that method.

Design/methodology/approach

The primary intent in the paper is to work on the particular solution of the Poisson equation by representing the source term through a Fourier series. The evaluation of the Fourier coefficients requires a rectangular domain even though the original domain can be of any arbitrary shape. The boundary conditions for the homogeneous solution gets modified by the projection of the particular solution on the original boundaries. The paper also develops a new Gauss quadrature procedure to compute the integrals appearing in the Fourier coefficients in case they cannot be analytically evaluated.

Findings

The current endeavor has developed two different representations of the source terms. A comprehensive set of benchmark exercises has successfully demonstrated the effectiveness of both the methods, especially the second one. A subsequent detailed analysis has identified the errors emanating from an inadequate number of boundary nodes and Fourier modes, a high difference in sizes between the particular solution and the original domains and the used Gauss quadrature integration procedures. Adequate mitigation procedures were successful in suppressing each of the above errors and in improving the solution accuracy to any desired level. A comparative study with the finite difference method revealed that the BIM was as accurate as the FDM but was computationally more efficient for problems of real-life scale. A later exercise minutely analyzed the heat transfer physics for a fin after validating the simulation results with the analytical solution that was separately derived. The final set of simulations demonstrated the applicability of the method to complicated geometries.

Originality/value

First, the newly developed Gauss quadrature integration procedure can efficiently compute the integrals during evaluation of the Fourier coefficients; the current literature lacks such a tool, thereby deterring researchers to adopt this category of methods. Second, to the best of the author’s knowledge, such a comprehensive error analysis of the solution method within the BIM framework for the Poisson equation does not currently exist in the literature. This particular exercise should go a long way in increasing the confidence of the research community to venture into this category of methods for the solution of the Poisson equation.