An axisymmetric boundary element method (BEM) has been developed to simulate atomization processes in a pressure-swirl atomizer. Annular ligaments are pinched from the parent sheet and presumed to breakup via the linear stability model due to Ponstein. Corrections to Ponstein's result are used to predict satellite droplet sizes formed during this process. The implementation provides a first-principles capability to simulate drop size distributions for low viscosity fluids. Results show reasonable agreement with measured droplet size distributions and the predicted SMD is 30-40% smaller than experiment. The model predicts a large number of very small droplets that cannot typically be resolved in an experimental observation of the spray. A quasi-3-D spray visualization is presented by tracking droplets in a Lagrangian fashion from their formation point within the ring-shaped ligaments. A complete simulation is provided for a case generating over 80,000 drops. (C) 2009 Elsevier Ltd. All rights reserved.