An Eulerian/Eulerian multiphase flow model coupled with a population balance model is used as the basis for numerical simulation of atomization in swirling flows. The objective of this exercise is to develop a methodology capable of predicting the local point-wise drop size distribution in a spray, such as would be measured by the Phase Doppler Particle Analyzer (PDA). Model predictions are compared to experimental measurements of particle size distributions in an airblast atomizer spray to demonstrate good qualitative and quantitative agreement. It is observed that the dependence of velocity on drop size inherent in a multiphase description of the drop cloud appears necessary to capture some features of the experimental data. Using this model, we demonstrate the relative contributions of secondary atomization and transport to the variation observed in the downstream spray drop size distribution.