Numerical simulation of gas-solid flow behaviors in a rectangular fluidized bed is carried out three dimensionally by the discrete element method (DEM).Euler method and Lagrange method are employed to deal with the gas phase and solid phase respectively.The collided force among particles,striking force between particle and wall,drag force,gravity,Magnus lift force and Saffman lift force are considered when establishing the mathematic models.Soft-sphere model is used to describe the collision of particles.In addition,the Euler method is also used for modeling the solid phase to compare with the results of DEM.The flow patterns,particle mean velocities,particles’ diffusion and pressure drop of the bed under typical operating conditions are obtained.The results show that the DEM method can describe the detailed information among particles,while the Euler-Euler method cannot capture the micro-scale character.No matter which method is used,the diffusion of particles increases with the increase of gas velocity.But the gathering and crushing of particles cannot be simulated,so the energy loss of particles’ collision cannot be calculated and the diffusion by using the Euler-Euler method is larger.In addition,it is shown by DEM method,with strengthening of the carrying capacity,more and more particles can be schlepped upward and the dense suspension upflow pattern can be formed.However,the results given by the Euler-Euler method are not consistent with the real situation. Numerical simulation of gas-solid flow behaviors in a rectangular fluidized bed is carried out three dimensionally by the discrete element method (DEM). Euler method and Lagrange method are employed to deal with the gas phase and solid phase respectively. Collided force among particles , striking force between particle and wall, drag force, gravity, Magnus lift force and Saffman lift force are considered when establishing the mathematic models. soft-sphere model is used to describe the collision of particles. In addition, the Euler method is also used for modeling the solid phase to compare with the results of DEM. flow patterns, particle mean velocities, particles’ diffusion and pressure drop of the bed under typical operating conditions are. The results show that the DEM method can describe describe the detailed information among particles, while the Euler-Euler method can not capture the micro-scale character. No matter which method is used, the diffusion of particles increases with the incrementation ase of gas velocity.But the gathering and crushing of particles can not be simulated, so the energy loss of particles’ collision can not be calculated and the diffusion by using the Euler-Euler method is larger.In addition, it is shown by DEM method, with strengthening of the carrying capacity, more and more particles can be schlepped upward and the dense suspension upflow pattern can be formed. However, the results given by the Euler-Euler method are not consistent with the real situation.