The present study focuses on numerical modeling for droplet dispersion by turbulence and turbulence modulation by droplets. To account for the dense spray effects, modulation models, a droplet collision model, and the Reitz’s wave instability breakup model are incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. A parcel probability density function(PDF) approach is implemented to improve the efficiency in droplet dispersion calculations. The numerical results indicate that the present parcel PDF model has the capability to realistically represent turbulent dispersion in dilute and dense sprays with improved efficiency over the delta function stochastic separated flow(SSF) model. Comparative performance of the existing turbulence modulation models are discussed in detail.