Flow behavior of solid phases is simulated by means of Eulerian-Eulerian in a liquid-solid fluidized bed with modified drag model based on energy-minimization multi-scale (EMMS) method. The modified EMMS drag coefficient is characterized by the treatment of the particle-rich dense phase and the liquid-rich dilute phase as the two interpenetrating continua. It was shown that the modified EMMS drag coefficient can predict reasonably the solid concentration profiles in a liquid-solid fluidized bed. The distributions of solid velocity, granular temperature and granular pressure are predicted. The phenomenon of back-mixing near the wall is found in the liquid-solids fluidized beds. (C) 2014 The Society of Powder Technology Japan. Published by Elsevier B.V. and The Society of Powder Technology Japan. All rights reserved.