A numerical analysis of multisized glass particles interacting with a confined bluffbody flow was performed by combining the finite-volume method for the gaseous flow with the Lagrangian approach for the particulate flow. The second-moment Reynolds-stress model was used to predict the turbulent gaseous flow in a gas-particle system, where an improved eddy-interaction model was used to predict turbulence-induced particle dispersion. The interaction between the two phases was accounted for in terms of coupling sources. Numerical predictions of two-phase mean and fluctuating velocities for different particle sizes were compared with corresponding experimental data. Reasonably good agreement was achieved for the mean properties of both the gaseous and particulate flows.