A new method for solution of thermal fields.around arrays of large pipes is described. It is shown that increased interaction between members of the array can be matched by appropriate splitting and variation of intensities of sources, sinks and dipoles with a view to maintain boundary conditions. Rapidly converging series of elementary analytic functions provide high accuracy for a wide range of array packing. In this context a specific example of a system comprising a wall heated on one side and cooled by array of pipes is considered