Abstract: Numerical simulation of spray combustion in liquid propellant rocket engine with PISO algorithm was implemented.A second-order upwind scheme developed by the authors was combined with thePISO algorithm to ensure the numerical accuracy.The Navier-Stokes equation in cylindrical coordinates was transferred into general curvilinear coordinate system and discretized by finite-volume method with non-staggered grid.The body-fitted grid generated by TTM method was used in the calculation.A stochastic spray model was employed to determine the droplet initial conditions,and the chemical reaction rate was given by Arrhenius model.The results show that the numerical method and physical models are proper for simulation of spray combustion processes.