A comprehensive study on confined swirling flows in an operational gas-turbine injector was performed by means of large-eddy simulations. The formulation was based on the Favre-filtered conservation equations and a modified Smagorinsky treatment of subgrid-scale motions. The model was then numerically solved by means of a preconditioned density-based finite-volume approach. Calculated mean velocities and turbulence properties show good agreement with experimental data obtained from the laser-Doppler velocimetry measurements. Various aspects of the swirling flow development (such as the central recirculating flow, precessing vortex core and Kelvin–Helmholtz instability) were explored in detail. Both co- and counter-rotating configurations were considered, and the effects of swirl direction on flow characteristics were examined. The flow evolution inside the injector is dictated mainly by the air delivered through the primary swirler. The impact of the secondary swirler appears to be limited.