In this study, a three-dimensional model based on RANS, slender-body theory and Newton-Euler dynamics is established to study the number concentration, one of the most important fluidization characteristics of cylindrical particles. Also, the effects of interaction between cylindrical particles are taken into account by introducing the rigid collision dynamics. To validate the model, the fluidization experiments of cylindrical particles in a cold-state fluidized bed are carried out. The number concentration characteristics of cylindrical particles are obtained from computational fluid dynamics (CFD) simulation. It is found that cylindrical particles arriving at the exit of the riser the earliest come from the near-wall regions, the horizontal transfer of so many cylindrical particles from the radial centre regions to the near-wall regions is evident. Meanwhile, there is no distinct relationship between the number concentration and inlet wind velocity.