A numerical simulation of a turbulent spray combustion flow is carried out to investigate feasibility of accurate computation for the spray combustion flows by using Large-Eddy Simulation (LES). An experimental database for a turbulent non-premixed spray flame is chosen to validate the present numerical simulation. The governing equations for the gas phases are discretized in a three-dimensional staggered grid system, and the fuel droplet motion equations are described in a Lagrangian representation. The flamelet approach is employed for a turbulent combustion model. A new thermal coupling model is proposed to estimate interaction between the gas-phase thermal energy and the dispersed-phase. The flame temperature drop caused by the thermal impact from the spray droplets is expressed properly by employing the thermal coupling model.