A computationally efficient spray model is presented for the simulation of transient vaporizing engine sprays. It is applied to simulate high-pressure fuel injections in a constant volume chamber and in mixture preparation experiments in a light-duty internal combustion engine. The model is based on the Lagrangian-erature. However, different values were seen for the gas-jet model constants for accurate simulations of the initial spray transient. The results show that there is a direct correlation between the predicted initial liquid-phase transient and the global gas-phase jet penetration. Model validation was also performed in engine simulations with the same set of constants. The model captured mixture preparation well in all cases, proving its suitability for simulations of transient spray injection in engines. 45 ulations, the Lagrangian-Particle/Eulerian-Fluid (LDEF, Dukowicz (1980)) approach is commonly adopted because of the many scales separating the internal injector and near-nozzle flows that affect the liquid phase development, and the consequent gas-phase turbulent flame in the combustion chamber. Because of the lack of resolution in the description of the liquid spray core, phenomeno-50 3 logical atomization models have been developed (e.g., Reitz and Bracco (1982); Reitz (1987); Reitz and Diwakar (1987); Tanner (1997); Habchi et al. (1997); Huh et al. (1998); Bianchi and Pelloni (1999); Beale and Reitz (1999); Hiroyasu (2000); Gorokhovski and Herrmann (2008b)), within the LDEF framework. This approach, although successful in a wide variety of simulations, suffers from 55 significant time-step and grid-resolution dependency, which is especially true for Reynolds-Averaged Navier-Stokes (RANS) approaches where all turbulence scales of the flow are modeled, and the grid resolution is more than one order of magnitude coarser than the characteristic injector diameter. Thus, some recent approaches have attempted to reduce the dependency of spray simula-60 tions on grid resolution and model constants by applying subcycling schemes to the Lagrangian particle step (Wang et al. (2010)), using region-of-interest instead than Computational Fluid Dynamics (CFD) grid-based collision calculations (Schmidt and Rutland (2004); Munnannur and Reitz (2009)), and using model-computed instantaneous field velocities in the near-nozzle region instead 65 of the under-resolved CFD flow fields (Ra et al. (2005); Abani et al. (2008a,b)). The last approach makes use of predictions from the theory of turbulent round jets (see, for example, Islam and Tucker (1980); Bremhorst and Hollis (1990); Abraham (1996); Iyer and Abraham (1997); Song and Abraham (2003); Singh and Musculus (2010); Musculus (2009); Liepmann and Gharib (1992)). These 70 models can predict a turbulent jet's penetration, velocity profiles, and gas concentrations within the jet, based on parameters such as an effective diameter, the densities of the gaseous environment, and the time-varying injection velocity. The models provide relevant flow properties that wou...