A stochastic model that describes agglomerate growth during fluidized bed spray agglomeration is presented and numerically solved by constant volume Monte Carlo method. The methodology overcomes the difficulties of solving multivariate population balance equations and includes continuous binder addition and drying. Agglomerate formation is treated as a complex combination of consecutive and parallel micromechanisms. Due to the discrete nature of the approach, the individual role of the micro-mechanisms on the agglomeration behavior can be analyzed. The results suggest that the droplet capture mechanism governs the agglomeration speed while the maximum agglomerate diameter is ruled by the equilibrium reached between coalescence, rebound and breakage. The mechanism of deposited binder drying is found to play a negligible role on agglomerate formation because of an extremely rapid droplet consumption. The main process variables affecting each micro-mechanism have been identified showing that the liquid spraying rate affects directly the droplet capture mechanism whereas binder properties influence mainly the agglomeration and rebound interactions. The model presented in this study is able to predict qualitatively the experimentally observed response of the system as well as the general shape of the agglomerate size distribution under the variation of several process parameters, demonstrating the potential of the discrete micro-level approach. (C) 2009 Elsevier Ltd. All rights reserved.