The formation of a liquid spray emanating from a nozzle in the presence of atomizing air was studied using a computational model approach that accounted for the deformation and break up of droplets. Particular attention was given to the formation of sprays under non-swirling flow conditions. The instantaneous fluctuating fluid velocity and velocity gradient components were evaluated with the use of a probability density function (PDF)-based Langevin equation. Motions of atomized fuel droplets were analyzed, and ensemble and time averaging were used for evaluating the statistical properties of the spray. Effects of shape change of droplets, and their breakup, as well as evaporation, were included in the model. The simulation results showed that the mean-square fluctuation velocities of the droplets vary significantly with their size and shape. Furthermore, the mean-square fluctuation velocities of the evaporating droplet differed somewhat from non-evaporating droplets. Droplet turbulence diffusivities, however, were found to be close to the diffusivity of fluid point particles. The droplet velocity, concentration, and size of the simulated spray were compared with the experimental data and reasonable agreement was found. Fluids 2019, 4, 20 2 of 17 Earlier, McDonell and Samuelsen [22] provided experimental data on droplet size and velocity distributions in a simplex atomizer. Kvasnak et al. [8] examined the formation of liquid sprays in the absence of atomizing air. Their analysis included effects of primary breakup for liquid ellipsoidal droplets based on the Taylor Analogy Breakup (TAB) model developed by O'Rourke and Amsden [23], and later modified by Ibrahim [24], and Kvasnak and Ahmadi [25]. Recent reviews of the spray simulation and breakup methods were presented by Jenny et al. [5].This work was concerned with understanding the effects of droplet deformation, breakup, and evaporation on the dispersion of deforming ellipsoidal spray droplets in turbulent fuel injector flow fields with atomizing air, with a focus towards practical applications. The other goal was to develop a computational tool for the simulation of practical liquid spray fuel injectors. Kvasnak et al. [8] described a detailed numerical procedure for generating ensembles of simultaneous sample particle trajectories for estimating particle velocity and dispersion statistics. The presented study showed that, in many cases of practical importance to fuel injectors with atomizing air, the deformation time scale of the fuel droplets was comparable with their evaporation time scale. The droplet nonsphericity was also found to be an important factor for modeling the fuel injector performance. The simulation results were compared with the experimental data of McDonell and Samuelsen [22] for a simplex atomizer with atomizing air, where qualitative agreement was found.