Spray formation and break-up are crucial phenomena for mixture formation inside diesel engines, both for combustion control and pollutant formation. Since the emission restrictions have become more and more severe in the last years, many studies have been conducted in order to improve diesel injection. Numerical simulations have proven to be reliable in producing results in a faster and cheaper way than experimental measures. The recent great progresses in computer science, then, have allowed to reach great accuracy in the simulations.In this work, a novel methodology based on Boltzmanns Kinetic Theory is applied to diesel injection. Lattice Boltzmann BGK (LBGK) provides and alternative method for solving fluid-dynamic problems and allows even superior accuracy as compared to conventional CFD.The multiphase approach used in this paper to study spray formation and primary is based on the works by Shan and Chen and their successive modifications.The results are compared to those predicted by means of the most accredited models in literature.