Abstract: In liquid-solid fluidized beds, the existence of the interstitial liquid among particles generates more energy dissipation different from that for gas-solid flows. To accurately predict the flow behavior of particles in a liquid-solid fluidization system, we simulate liquid-solid fluidization by introducing a dynamic coefficient of restitution model into the kinetic theory of granular flow. The dynamic coefficient of restitution is obtained by incorporating a semi-empirical correlation generated from a freely falling particle impacting a wet plate with a critical Stokes number, and the variable critical Stokes number is used to judge whether particles rebound after the collision. Simulated results show that the granular temperature and granular pressure are clearly affected when the dynamic coefficient of restitution model is introduced. As well, the particles velocity distribution predicted by the dynamic coefficient of restitution model is in better agreement with the experimental values compared with the constant coefficient of restitution.