A numerical study has been carried out to understand the effects of Unsteady Reynolds Averaged Navier-Stokes (standard k ― ε and RNG k ― ε model) and large eddy simulations (LES) on a multi-hole gasoline direct injection (GDI) system. The fuel injector considered in this study is the Spray G nozzle from the Engine Combustion Network (ECN). A blob injection model, based on empirical rate of injection (ROI) profile, is considered in this study. The latest data on spray penetrations from Engine Combustion Network is used for model validation along with experimental findings on suction velocity and local droplet diameter. The spray breakup is simulated by using the KH-RT breakup length model. The turbulence model constant Cε1, is tuned to match with the experimental data of liquid and vapor penetrations in simulations while using the standard k-ε turbulence model. On the other hand, the Kelvin-Helmholtz breakup model time constant (B1) and Rayleigh Taylor breakup length constant (Cbl) are tuned for the RNG k ― ε turbulence model. From this work it is observed that by increasing the breakup length model constants (Cbl), the radial dispersion of the spray increases, and the extent of breakup is lowered. The set of optimized model parameters used with RNG k - ε is also used for LES modeling studies with different sub-grid models. The spray penetrations with standard k ― ε turbulence (Cε1=1.44) model are reported underpredicting, and the RNG k ― ε and LES sub-grid models predicted well with the latest and recommended data from ECN. In terms of gas axial velocity comparison, the standard k-ε(Cε1=1.44) simulation setup does not perform as well as the simulation setups using RNG k-ε and LES turbulence models (with breakup parameters: Cbl = 16 and B1 = 32). However, the standard k-ε(Cε1=1.44) simulation setup perform better than the simulation setups using RNG k-ε and LES turbulence models (with breakup parameters: Cbl = 16 and B1 = 32) when it comes to predicting local droplet diameter at 15 mm downstream of the injector tip. A parametric study is also performed considering the geometry of the stepped holes in the computational domain. The rate of injection based simulation is initiated at the end of the smaller hole. The case including the stepped holes led to over-prediction compared to the case with the usual computational domain (i.e., without the stepped holes), in terms of spray penetrations, but exhibited higher levels of fluctuations in the spray morphology. Finally, parametric studies were carried out to understand the relative importance of the individual spray sub-models (breakup, evaporation and collision) and the results are conclusive that for a spray simulation the breakup models are the dominant factors.