Modeling the fuel injection process in modern gasoline direct injection engines plays a principal role in characterizing the in−cylinder mixture formation and subsequent combustion process. Flash boiling, which usually occurs when the fuel is injected into an ambient pressure below the saturation pressure of the liquid, is characterized by fast breakup and evaporation rates but could lead to undesired behaviors such as spray collapse, which significantly effects the mixture preparation. Four mono−component fuels have been used in this study with the aim of achieving various flashing behaviors utilizing the Spray G injector from the Engine Combustion Network (ECN). The numerical framework was based on a Lagrangian approach and was first validated for the baseline G1 condition. The model was compared with experimental vapor and liquid penetrations, axial gas velocity, droplet sizes and spray morphology and was then extended to the flash boiling condition for iso−octane, n−heptane, n−hexane, and n−pentane. A good agreement was achieved for most of the fuels in terms of spray development and shape, although the computed spray morphology of pentane was not able to capture the spray collapse. Overall, the adopted methodology is promising and can be used for engine combustion modeling with conventional and alternative fuels.