Abstract
The aim of this paper is to propose a numerical strategy for computing the solution of two-dimensional time-harmonic acoustic multiple scattering problems at high-frequency. The scatterers are assumed to be circular, leading therefore to semi-analytical representation formulae of the scattered field through the solution of a large linear system of equations. Taking advantage of the special block Toeplitz structure of the matrix of the linear system, a fast iterative and preconditioned numerical method yielding large memory savings is proposed. Several numerical experiments for general configurations are presented to show the efficiency of the numerical method.