A series of numerical simulations is undertaken to study a highly swirling turbulent flow generated by rotor-stator interaction in a swirl generator. The purpose is to assess the applicability of different turbulence models in swirling flow with a high level of unsteadiness and a significant production and dissipation of turbulence in the flow away from the wall. Nine turbulence models are compared: four high-Reynolds URANS, two low-Reynolds URANS and three hybrid URANS-LES. These are the standard k − , SST k − ω, realizable k − , RNG k − , Launder-Sharma k − , Lien-Cubic k − , delayed DES Spalart-Allmaras, DDES SST k − ω and improved DDES-SA. The URANS models are capable of capturing the main unsteady feature of this flow, the so-called helical vortex rope, which is formed by the strong centrifugal force and an on-axis recirculation region. However, the size of the on-axis recirculation region is overestimated by the URANS models. Although the low-Reynolds URANS formulations resolve the boundary layers in the runner and the draft tube more accurately, they still encounter difficulties in predicting the main flow features in the adverse pressure gradient in the draft tube. It is shown that a more detailed resolution, which is provided by the hybrid URANS-LES methods, is necessary to capture the turbulence and the coherent structures. The flow contains a strong disintegration of the vortex rope which is predicted well by the hybrid RANS-LES models. The hybrid methods also capture the blade wakes better than the other models, elucidating the wake interaction with the vortex rope. The frequency of the vortex rope is predicted well and the total turbulence (resolved and modeled), suggested by DDES-SA, corresponds reasonably well to the experimental results.