A high-fidelity numerical simulation of jet breakup and spray formation from a complex diesel fuel injector has been performed. A full understanding of the primary atomization process of diesel fuel injection has not been achieved for several reasons, including the difficulties in accessing the optically dense region. Due to recent advances in numerical methods and computing resources, high-fidelity simulations of realistic atomizing flows are currently feasible, providing a new mechanism to study the jet breakdown process. In the present study, a novel volume-of-fluid (VOF) method coupled to a stochastic Lagrangian spray (LSP) model is employed to simulate the atomization process. A common rail fuel injector is modeled by a nozzle geometry provided by the engine combustion network (ECN). The working conditions correspond to a single 90 µm orifice JP-8 fueled injector operating at 90 bar and 373 K and releasing into a 100% nitrogen, 29 bar, 300 K ambient with a Rel= 16, 071 and Wel= 75, 334, placing the liquid jet in the atomization breakup regime. The experimental dataset from Army Research Lab (ARL) is used for validation and the Kelvin-Helmholtz/Rayleigh-Taylor (KH-RT) breakup model (Reitz & Bracco 1979) is used for verification, both in terms of spray angle. Droplet distributions of the simulated spray are provided for future experimental comparisons and secondary atomization simulations using LSP modeling.