Most of today's computational fluid dynamics (CFD) calculations for gas-solid flows are carried out assuming that the solid phase is monodispersed, whereas it is well known that in many applications, it is characterized by a particle size distribution (PSD). In order to properly model the evolution of a polydisperse solid phase, the population balance equation (PBE) must be coupled to the continuity and momentum balance equations. In this work, the recently formulated direct quadrature method of moments (DQMOM) is implemented in a multi-fluid CFD code to simulate particle aggregation and breakage in a fluidized-bed (FB) reactor. DQMOM is implemented in the code by representing each node of the quadrature approximation as a distinct solid phase. Since in the multi-fluid model, each solid phase has its own momentum balance, the nodes of the DQMOM approximation are convected with their own velocities. This represents an important improvement with respect to the quadrature method of moments (QMOM) where the moments are tracked using an average solid velocity. Two different aggregation and breakage kernels are tested and the performance of the DQMOM approximation with different numbers of nodes are compared. These results show that the approach is very effective in modeling solid segregation and elutriation and in tracking the evolution of the PSD, even though it requires only a small number of scalars. (C) 2003 Elsevier B.V. All rights reserved.