Understanding unsteady effects of the Gasoline Direct injection (GDi) process, due to short injection duration, plays a major role in the analysis of the mixture formation and so combustion efficiency in the spark ignition engines. Focusing on the phase when the needle is at its maximum lift, there still are some uncertainties. For instance, the time evolution of the upstream pressure, together with the detailed geometry of the needle and the seat, shape of the rate of injection could affect the uniformity of the flow between orifices. Experimentally addressing these issues nowadays remains challenging, if not impossible, so Computational Fluid Dynamics (CFD) are used. Homogeneous Relaxation Model (HRM) is employed to consider the mass exchange between liquid and vapor phases of the fuel inside the nozzle. Different Large Eddy Simulations (LES) sub-grid models are employed in order to account for the unsteady effects of turbulence and analyze the predictive capabilities in resolving the spatial-temporal scales. Results are validated against experimental data. Vortices are generated inside the counterbore and enhance the liquid dispersion near the injector. Spray parameters, after a proper time-averaging, also accurately match the results reported in the literature. Root mean square velocity (Urms) fluctuation has higher levels of spatially located fluctuations in areas with significant levels of small-scale turbulence and noticeable liquid-vapor mixture.