A number of droplet breakup models have been developed to predict the diesel spray. The capabilities of droplet deformation and breakup models such as TAB, ETAB, DDB and APTAB models are evaluated in modeling the non-evaporating diesel sprays injected into atmosphere. New methods are also suggested that take into account the non- spherical shape of droplets and the reduced drag force by the presence of neighbouring droplets. The KIVA calculations with standard ETAB, DDB, and APTAB models predict well the spray tip penetrations of the experiment, but overestimate the Sauter mean Diameter(SMD) of droplets. The calculation with non spherical droplets injected from the nozzle shows very similar results to the calculation with spherical droplets. The drag coefficient which is linearly increased with the time after start of injection during the breakup time gives the smaller SMD that agrees well with the experimental result.