Circulating fluidized beds (CFBs) have applications in many industrial processes like fluid catalytic cracking, coal gasification, coal combustion, biomass gasification and chemical looping. Consequently, the process have been studied intensively for many researchers seeking to understand the complex gas-solid flow encountered in circulating fluidized beds. In this sense, the Computational Fluid Dynamics (CFD) tools are very useful to study gas-solid flows, but this approach directly depends on the correct choice of the mathematical models. There are many numerical studies on CFBs at low solids fluxes and only a few studies for high solids flux, which is very common in many applications, therefore, this regime was chosen to conduct the simulations in this work. The results were compared with experimental data on a riser with an internal diameter of 76 mm and a height of 10 m. Radial and axial profiles were computed using a two-phase 3-D computational fluid dynamics model, with four different correlations for the drag between the phases, Gidaspow, Syamlal-O’Brien, and other two models that account for particle clusters, the energy-minimization multi-scale (EMMS) and the Four-Zone model. The results indicates that all the correlations predicts the solids concentration, found experimentally, in the dilute and developed flow regions but, in the region of highest solids concentration, especially at the level close to the inlet, the Gidaspow and Syamlal-O’Brien correlations not properly represented the solids concentration. The EMMS and Four-Zone models improved the results on these regions showing the influence of these models especially at the bottom of the riser, where the solids concentration is higher.