This paper describes a Large Eddy Simulation (LES) investigation into flow fields in a model gas turbine combustor equipped with a swirl burner. A probability density function was used to describe the interaction physics of chemical reaction and turbulent flow as liquid fuel was directly injected into the combustion chamber and rapidly mixed with the swirling air. Simulation results showed that heat release during combustion accelerated the axial velocity motion and made the recirculation zone more compact. As the combustion was taking place under lean burn conditions, NO emissions was less than 10 ppm. Finally, the effects of outlet contraction on swirling flows and combustion instability were investigated. Results suggest that contracted outlet can enhance the generation of a Central Vortex Core (CVC) flow structure. As peak RMS of velocity fluctuation profiles at center-line suggested the turbulent instability can be enhanced by CVC motion, the Power Spectrum Density (PSD) amplitude also explained that the oscillation at CVC position was greater than other places. Both evidences demonstrated that outlet contraction can increase the instability of the central field.