Abstract: Numerical simulation of fluid-solid two-phase flows is reviewed with 105 references. There are two approaches (Lagrangian approach and Eulerian approach) typically used for modeling fluid-solid flows. The former treats the fluid phase as continuous phase and the solid phase as dispersed phase and predicts the particle trajectories in the fluid phase as the results of forces acting on particles.The latter treats the solid as a certain kind of continuum and solves two continuum equations for both fluid and solid phase. At first, as the basis of fluid-solid two-phase flow, several turbulence models, such as Reynolds Association Numerical Simulation (RANS), Direct Numerical Simulation (DNS), Large Eddy Simulation (LES) and Discrete Vortex Simulation (DVS), for the continuous fluid phase are summarized,and the mechanisms of RANS are discussed in detail. In the Lagrangian approach, a particle kinetic equation is presented and the importance of each force acting on a particle is discussed. For a flow with dense particles, various interaction models between particles are discussed , among which the Distinct Element Method (DEM) newly developed from rock deformation study is particularly described, and the new method shows great predicting ability. In the Eulerian approach, the basis equations for the solid and fluid phases as well as their numerical treatments are presented. To close the turbulence and phase related equations, calculation of solid phase viscosity and pressure are discussed respectively. Finally, a k-ε-k p model and model derived from kinetic theory of granular flow are especially introduced. As a complement, a discussion of one-way coupling and two-way coupling between fluid and solid phase is included. Based on the review of the advantages and disadvantages of all models, the existing problems and the future development direction of fluid-solid two phase flow are given as conclusions.