Two-dimensional interpolation-supplemented and Taylor-series expansion-based lattice Boltzmann method (ITLBM) is proposed to simulate flows in the non-uniform grids. The proposed method is based on the standard lattice Boltzmann method, interpolation-supplemented and Taylor-series expansion. The final formulation can be used on any mesh structure and lattice Boltzmann model. Numerical test of a two-dimensional channel flow around a square cylinder has been studied. The computational efficient and recirculation length at Re=1, 15 are obtained. Comparing the results from the ITLBM with those from the standard LBM, it has been concluded that the proposed method has good prospects in the hydrodynamic and aerodynamic engineering applications.