Lattice Boltmzann method (LBM) has been developed as an important technique of computational fluid mechanics. Two-dimensional turbulent flow around NACA0012 airfoil, with Re=10 5 and AOA=4 ○ (angle of attack), is simulated with "multiple-relaxation time" lattice Boltzmann method incorporated with Spalart-Allmaras (SA) turbulence model. From the viewpoint of industrial applications, the related computational aspects are discussed, such as treatments of wall boundary, far field boundary conditions, and the position of the airfoil in the domain. Meanwhile, to balance the computational loads and accuracy of the simulations, grid refinement is applied. According to the results, LBM gives pretty good predictions, and quite flexible domain size can be accepted if the far field boundary condition is well posed. Compared to 3D LES-LBM simulations, two-dimensional LBM with SA turbulent model is of lower computational load, but with good accuracy.