Numerical solutions of 2-D laminar flow over a backward-facing step using the lattice Boltzmann equation method (LBEM) are presented in this article. Unlike conventional numerical schemes based on macroscopic continuum equation (mass conservation and Navier-Stokes) discretisation, the LBEM is based on microscopic models and mesoscopic kinetic equations. The simulations were validated for a wide range of Reynolds numbers (100 ≤ Re ≤ 1,000), comparing them to previous studies. Several flow features, such as primary and secondary vortex location at the bottom and top of the wall, respectively, were investigated regarding Reynolds number. Two typical classes of boundary condition were implemented in the LBEM model: the Drichlet condition at the inlet flow (parabolic speed profile) and the Newman condition at the outlet flow (zero gradient speed). The results showed that the LBEM gave accurate results over a wide range of Reynolds number; these were compared with other numerical methods and experimental data.