Large eddy simulations (LES) were performed to study the non-reacting flow fields of a Cambridge swirl burner. The dynamic Smagorinsky eddy viscosity model is used as the sub-grid scale turbulence model. Comparisons of experimental data show that the LES results are capable of predicting mean and root-mean-square velocity profiles. The LES results show that the annular swirling flow has a minor impact on the formation of the bluff-body recirculation zone. The vortex structures near the shear layers, visualized by the iso-surface of Q-criterion, display ring structures in non-swirling flow and helical structures in swirling flow near the burner exit. Spectral analysis was employed to predict the occurrence of flow oscillations induced by vortex shedding and precessing vortex core (PVC). In order to extract accurately the unsteady large-scale structures in swirling flow, a three-dimensional proper orthogonal decomposition (POD) method was developed to reconstruct turbulent fluctuating velocity fields. POD analysis reveals that flow fields contain co-existing helical and toroidal shaped coherent structures. The helical structure associated with the PVC is the most energetic dynamic flow structure. The latter toroidal structure associated with vortex shedding has lower energy content which indicates that it is a secondary structure. Cambridge swirl burner, large eddy simulation, proper orthogonal decomposition, vortex shedding and breakdown, precessing vortex core Citation: Zhang H D, Han C, Ye T H, et al. Large eddy simulation of unconfined turbulent swirling flow.