To accurately simulate industrial sized bubble columns, the accurate prediction of the turbulent structures is important. Currently used RANS and LES models cannot capture the dynamics of the bubble columns or require a high grid resolution, respectively. In this work the Detached Eddy Simulation (DES) was used to combine the advantages of RANS and LES approaches. The DES method was based on Spalart–Allmaras, k−ϵ and k−ω SST turbulence models and was used to simulate gas–liquid flow for a bubble column with a square cross-section. The results are compared with experimental data and Large Eddy Simulation (LES) results, based on the Vreman and Smagorinsky subgrid-scale models. Profiles of turbulence kinetic energy and average axial liquid and...