In this work the quadrature method of moments (QMOM) is implemented in a commercial computational fluid dynamics (CFD) code (FLUENT) for modeling simultaneous aggregation and breakage. Turbulent and Brownian aggregation kernels are considered in combination with different breakage kernels (power law and exponential) and various daughter distribution functions (symmetric, erosion, uniform). CFD predictions are compared with experimental data taken from other work in the literature and conclusions about CPU time required for the simulations and the advantages of this approach are drawn. (C) 2003 Elsevier Ltd. All rights reserved.