Based on fully structured grids parallel numerical simulations of flow around a cylinder under different Reynolds number are carried out. Two-dimensional and three-dimensional models are established at the same time under specific Reynolds number, and further analyze of three-dimensional flow characteristics as well as the generated influence to overall physical quantities are presented. In order to explore efficient high Reynolds number turbulence models, a comparative research of the LES model without wall functions and the Spalart-Allmaras turbulence model is carried out. In order to improve the computational efficiency, a domain decomposition parallel computing strategy is used, and a calculation strategy that results of coarse grid was assigned to fine grid as initial field value by 3D linear interpolation is presented. Simulation results show that: Drag coefficient and Strouhal number have very good consistency with the experimental data, which verifies the correctness of the calculation method; Even if at low Reynolds number (200≤Re≤300), using a three-dimensional model is still necessary; While in the high Reynolds number stage, compared to LES model without wall functions, Spalart-Allmaras model is more applicable and more efficient.