By using the computational fluid dynamic software FLUENT, on the basis of the mixture model, gas-liquid two phase flow around the cylinders with equivalent pitch is simulated to obtain the drag and lift forces on the central cylinder and the flow behavior of the whole calculating domain. The results show that the root mean square value and extreme value of the lift and drag forces on the central cylinder and the frequency of the vortex shedding increase with the growth of Reynolds number. The rationality is also validated with experimental results.