This study presents the development of a phenomenological combustion model to simulate the combustion processes in diesel engines with multi-stage fuel injection. A newly developed zero-dimensional spray propagation model and a model of spray-to-spray interaction were combined with a stochastic combustion model, which had been developed for diesel combustion calculation in the cases of single-stage injection. In this model, the combustion chamber is divided into an ambient air zone and several spray zones, where the spray formed by each injection is treated as a spray zone. Turbulent mixing, fuel evaporation, heat loss, and chemical reactions, are calculated in each spray zone respectively. A zero-dimensional spray propagation model including the spray evolution after the end of injection (EOI) and the model of interaction between the sprays from sequent injections are developed to describe the spray behavior for the case of multistage injection. Then the developed combustion model is validated against the experimental data from a single-cylinder DI diesel engine with pilot/main two-stage injection, in which the pilot injection conditions are varied with fixed main injection timing. Based on the analysis of heat release rate, entrainment rate, and the microscopic information inside the spray, such as probability density function (PDF) of equivalence ratio, the effects of wall impingement and interaction between adjacent sprays on fuel-air mixing rate and entrainment rate are formulized and employed to reproduce the measured histories of heat release rate. Reduction of fuel-air mixing rate is considered when the spray flows into the squish region after the wall impingement, which is effective to capture the measured decrease of pilot spray's heat release with advancing pilot injection timing. The effects of wall impingement of the main spray and the interaction between adjacent sprays are modeled to reproduce the heat release rate during the initial and later part of the mixing-controlled combustion. After these improvements, heat release rates of the test engine when varying the pilot injection conditions have been successfully predicted.