Abstract: In this paper,large-eddy simulation combined with multi-relaxation time lattice Boltzmann (MRT-LBM) method is proposed for solving the poor stability problem in high Reynolds number flow simulation with LBM.The implementation is discussed in the context of 9-velocity (D2Q9) MRT-LBM model in conjunction with the Smargorinsky subgrid closure model.This model is used for doing the two-dimensional numerical simulation of flow around a cylinder with a high Reynolds number.According to the experimental results,a systematic analysis of the form of shed vortex of flow around cylinder is made,and the variation of the Strouhal number and the drag coefficient with the change of Reynolds number.The experimental results show that the method can well increase numerical stability in the simulation with LBM in a high Reynolds number.