Analytic solutions to the static and stationary boundary value field problems relative to an arbitrary configuration of parallel cylinders are obtained by using translational addition theorems for scalar Laplacian polar functions, to express the field due to one cylinder in terms of the polar coordinates of the other cylinders such that the boundary conditions can be imposed at all the cylinder surfaces. The constants of integration in the field expressions of all the cylinders are obtained from a truncated infinite matrix equation.

Translational addition theorems are available for scalar cylindrical and spherical wave functions but such theorems are not directly available for the general solution of the Laplace equation in polar coordinates. The purpose of deriving these addition theorems and applying them to field problems involving systems of cylinders is to obtain exact analytic solutions with controllable accuracies, thereby, yielding benchmark solutions to validate other approximate numerical methods.