Computational Fluid Dynamics (CFD) is a powerful tool to study the dense gas-solid flow in a circulating fluidized bed. Most of the existing methods focus on the microscopic properties of individual particle. Therefore, the simulation scale is significantly limited by the huge number of individual particles. and so far the numbers of particles in most of the reported simulations are less than 10(5). The hydrodynamics behaviour of particle clustering in a dense gas-solid two-phase flow has been verified by several experimental results. The Distinct Cluster Method (DCM) was proposed in this paper by studying the macroscopic particle clustering behaviour, and comprehensive models for cluster motion, collision, break-up, and coalescence have been well developed. We model the dense two-phase flow field as gas-rich lean phase and solid-rich cluster phase. The particle cluster is directly treated as one discrete phase. The gas turbulent flow is calculated by Eulerian approach. and the particle behaviour is studied by Lagrangian approach. Using the proposed method, a three-dimension dense gas-particle two-phase flow field in a circulating fluidized bed with square-cross-section, with particle number up to 7.162 x 10(7) are able to be numerically studied, on which few results have been reported. Details on instantaneous and time-averaged distributions are obtained. Developing process of non-uniform particle distribution is visualized. These results are in agreements with experimental observations, which justified the feasibility of using the DCM method to model and simulate dense gas-solid flow in a circulating fluidized bed with large number of particle numbers. (c) 2009 Elsevier B.V. All rights reserved.