Large eddy simulation of spray combustion was performed with Dynamic SGS model for compressible turbulence. SGS model of combustion reaction was also introduced to take account for contribution of SGS components to combustion reaction rate term. The SGS models were evaluated by comparing results of LES with DNS, which governing equations were not filtered and which had high resolution. A partition method for phase coupling terms and two types methods for partition of gaseous phase variables were applied to compare the results of simulations which have different grid resolutions. The results of LES were greatly improved with the increase of spatial resolution and the contribution of SGS stress decreased with increase of the resolution. However, LES with coarse resolution also could show good agreement with the DNS, in spite of the low resolution. On the other hand, the SGS diffusivity increased with the increase of spatial resolution. From the above results, we could conclude that the proposed partition method for phase coupling terms and partition type a were appropriate for LES of spray combustion.