Abstract: The cross-flow induced vibration of a row of cylinders at sub-critical Reynolds number Re = 2.67104 is simulated based on the surface vorticity method and a fluid-structure interaction model incorporating the effects of cylinder motions and displacements. The computational results of the vorticity map, fluid force, vibration response, and shedding frequency are presented. The non-uniform flow behind a rigid cylinder row with a small pitch ratio T/D = 1.5 characterized by narrow and wide near wake, multiple frequencies, and the vortex shedding behind a rigid cylinder row with T/D = 2.0 are replicated properly by simulation results. Furthermore, the FIV and fluid elastic instability of a fully flexible cylinder row at T/D = 1.5 is studied. The critical reduced velocity of the cylinder row at SG = 1.29 is also calculated.