A methodology for simulating the effect of wide range of compositional varieties of international market diesel fuels on its spray combustion is proposed and experimentally validated. Spray combustion properties of market diesel fuels from various countries were investigated in a constant volume vessel while the compositions of those tested fuels were analyzed through GCxGC and H-NMR. It was then found that ignition delay times, flame region and luminous intensity were significantly affected within the compositional variations. For its numerical simulation, a numarical evaporation surrogate consisting of 22 compounds, covering the basic molecular types and a wide range of carbon numbers of complex real fuels, is developed. Each compound in the evaporation surrogate is then imaginarily converted to a reaction surrogate consisting of n-cetane, methylcyclohexane and 1,2,4-trimethyl benzene so that the reaction surrogate could take over the combustion propensity of each compound in evaporation surrogates. Then, in CFD, light components evaporate at lower temperature than heavy components, and thus an upstream fuel/air mixture contains higher proportion of light components while a downstream one includes higher proportion of heavy components, which is consistent with the distillation curve of original market fuel. Thus, the ignition event is calculated based on local composition and temperature with a kinetic model. The simulations were performed using commercial software for various market fuels, and were validated with the above experimental results. It is concluded that, for experimentally observed magnitude relationships in ignition delay times, flame region and luminous intensity between different fuels, the proposed method reproduces 90% of all cases for ignition delay times, 60% of flame region and 70% for luminous intensity, respectively.