In this paper, we present a method for the numerical integration of the three-dimensional Green kernel over two triangles. This method is compared with the results obtained by other algorithms. The comparison proves the efficiency of the technique in removing the numerical singularity. Furthermore, we provide a numerical algorithm for increasing the computation accuracy without refining the finite element mesh. (c) 2004 Elsevier Inc. All rights reserved.