Small amplitude in-phase oscillations of multi-cylinders are considered both numerically and experimentally. Flow field is separated into inner and outer regions. The basic unsteady solution is obtained analytically and the secondary flow, termed as steady streaming flow, can be obtained numerically by using Finite Volume Code with Panel Method. The Particle Induced Velocimetry, one of whole field measurements, is introduced for comparison with numerical flow visualization quantitatively. Among the algorithms for PIV, the Three Step Vector Searching Technique is applied to reduce CPU time. Small but non-zero damping coefficient, that is important in lightly damped system can be obtained with varying number of bodies and distances.