Understanding the physics and chemistry involved in spray combustion, with its transient effects and the inhomogeneity of the spray is quite challenging. For efficient operation of both internal combustion and gas turbine engines, great insight into the physics of the problem can be obtained when a computational analysis is used in conjunction with either an experimental program or through published experimental data. The main area to be investigated to obtain good combustion begins with the fuel injection process and an accurate description of the mean diameter of the fuel particle, injection pressure, drag coefficient, rate shaping etc must be defined correctly. This work presents a methodology to perform the task set out in the previous paragraph and uses experimental data obtained from available literature to construct a semi-empirical numerical model for high pressure fuel injectors. A modified version of a multidimensional computer code called KIVA3V was used for the computations, with improved sub-models for mean droplet diameter, injection pressure, injection velocity, and drop distortion and drag. The results achieved show good agreement with the published in-cylinder experimental data for a Volkswagen 1.9 L turbo-charged direct injection internal combustion engine under actual operating conditions. It is crucial to model the spray distribution accurately, as the combustion process and the resulting temperature distribution and pollutant emission formation is intimately tied to the in-cylinder fuel distribution. The present scheme has achieved excellent agreement with published experimental data and will make an important contribution to the numerical simulation of the combustion process and pollutant emission formation in compression ignition direct injection engines and gas turbine engines.