This study is focused on the development of spray models for large eddy simulations (LESs) of high injection pressure diesel sprays. Basic governing equations are numerically solved using the Lagrangian-Eulerian approach for liquid and gas phases respectively, in the KIVA-3V code. Non-reacting simulations of evaporating spray are performed using a dynamic structure LES model in which the subgrid stress tensor is modelled with a non-viscosity tensor coefficient. This is a one equation based model, in which an extra transport equation for the subgrid kinetic energy (k) is solved. Subgrid scale energy exchange between droplets and the gas phase is identified as an important mechanism to capture the correct scaling of k, which eventually feeds back to both the spray droplets and the gas phase turbulent mixing. Hence, to account for spray-induced gas turbulence, a LES spray source model for k is developed. This model requires subfilter scale velocities, which are obtained by defiltering the filtered velocity field available from LES. Results are compared with corresponding RANS results and available experimental data for various physical spray variables such as spray penetration, gas phase penetration, droplet distribution, etc. The LES spray model was found to be important in predicting the correct liquid spray evolution and responds correctly to varying parameters such as nozzle size and ambient density.