The paper is devoted to the solution of Laplace equation by the boundary element method. The coupling between a finite element solution inside a bounded domain and a boundary integral formulation for an exterior infinite domain can be performed by producing a "stiffness" or "impedance matrix" which is equivalent to the interaction coming from BEM, when coupled with FEM stiffness or impedance matrix. It is shown in a first step that the use of classical Green's functions for plane domains can lead to impedance matrices which have eigenvalues of different signs, which is physically and numerically unsatisfying and also to singular impedance matrices, corresponding to the classical degenerate scale problem. Avoiding the degenerate scale problem is classically overcome by adding to the Green's function a constant which is large compared to the size of the 2D domain. However, it is shown that this constant affects the solution of exterior problems in the case of non-null resultant of the normal gradient at the finite boundary. It becomes therefore important to precise the value of this constant related to a characteristic length introduced into the Green's function. Using a "slender body theory" allows to show that for long cylindrical domains with a given cross section, the characteristic length is asymptotically equal to the length of the cylindrical domain.Comparing numerical or analytical 3D and 2D solutions on circular cylindrical domains confirms this result for circular cylinders.