Small-amplitude harmonic oscillations of multi-cylinders are considered both experimentally and theoretically. For the theoretical model, the flow regime is separated into inner and outer regions. In the inner region, the flow is governed by the generalized Stokes boundary layer equation. In the outer region, the full Navier-Stokes equation for the steady streaming flow is solved numerically by using ADI scheme and FVM coupled with the boundary integral method. Flow visualization experiments are conducted by using the Laser Sheet Image Technique. The case of two circular cylinders and square cylinders with variable distances are chosen as a typical example. Although experimental results are based on the flow in the finite domain, both experimental and numerical results agree well qualitatively. As the separation of cylinders is increased, a numerical result shows the asymptotic convergence to a single cylinder case.