A new spray model has been developed to improve the prediction of diesel engine combustion and emissions using the KIVA-II CFD code. The accuracy of modeling the spray breakup process has been improved by the inclusion of Rayleigh-Taylor accelerative instabilities, which are calculated simultaneously with a Kelvin-Helmholtz wave model. This model improves the prediction of the droplet sizes within a diesel spray and provides a more accurate initial condition for the evaporation, combustion, and emissions models. An improvement to the droplet drag model is also presented. This model accounts for the increased droplet drag due to the change in the droplet's shape, as well as the increase in the frontal area of the droplet. The drag model affects the breakup process locally, producing a more realistic droplet size distribution, and therefore a more accurate calculation of the vaporization process. Through the introduction of this model the prediction of the pressure, heat release, and emissions produced by single and split injections can now be accurately modeled.