An Eulerian granular multiphase model with a drag coefficient based on the energy minimization multiscale (EMMS) model was used to perform a three-dimensional (3D), full-loop, time-dependent simulation of hydrodynamics of a 150 MW, circulating fluidized bed (CFB) boiler. Simulation results were presented in terms of the pressure profile around the whole loop of solids circulation, profiles of solids volume fraction and solids vertical velocity, as well as the non-uniform distribution of solid fluxes into two parallel cyclones. (C) 2010 Elsevier B.V. All rights reserved.