Abstract: An LB-based gas-solid two-phase model with two-way coupling is developed considering feedback forcing of particles in evolution equation of fluid particles. Smagorinsky subgrid model is also introduced in simulation of flow field with high Reynolds numbers. Classic particle-laden flow over a backward facing step is simulated and velocity profiles of gas phase and particles (considering one-way coupling and two-way coupling respectively) are compared with experimental results. The results of two-way coupling LB model are obviously better than these of one-way coupling LB model. Furthermore, preferential concentration of particles with different Stokes numbers (St) is investigated. It is found that small particles (St~0(0.1)) show better following behaviors with gas phase and are uniformly distributed in the flow field. Particles with moderate Stokes numbers (St~0(1)) are hard to be entrained into the vortex and show strong preferential concentration. On the other hand, large particles (St~0(10)) can enter into the vortex because of great inertial and are distributed more uniformly in flow field.