This work evaluates the performance of the Sigma-Y Eulerian atomization model at reproducing the internal structure of a diesel spray in the near-field. In the study, three different computational domains have been used in order to perform 3D and 2D coupled simulations, where the internal nozzle flow and external spray are modeled in one continuous domain, and 2D decoupled simulations, where only the external spray is modeled. While the 3D simulation did the best job of capturing the dense zone of the spray, the 2D simulations also performed well, with the coupled 2D simulation slightly outperforming the decoupled simulation. The similarity in results between the coupled and the decoupled simulation show that internal and external flow calculations can be performed independently. In addition, the use of spatially averaged nozzle outlet conditions, in the case of an axisymmetric (single-hole) convergent nozzle, leads to a slightly worse near-field spray predictions but to an accurate far-field ones. Finally, a novel constraint on turbulent driven mixing multiphase flows is introduced which prevents the slip velocity from exceeding the magnitude of the turbulent fluctuations through a realizable Schmidt number. This constraint increased model stability, allowing for a 4x increase in Courant number. (C) 2015 Elsevier Ltd. All rights reserved.