An Euler-Euler large eddy simulation (LES) of transient gas-liquid turbulence flow in a bubble column was conducted to investigate the adaptiveness of various sub-grid scale (SGS) models. The performance of three SGS models, respectively the standard Smagorinsky model, the wall-adapting local eddy (WALE) model and the Dynamic model, on the predictions of bubbly flow was studied. Comparison between SGS models for time-averaged and transient velocities, gas void fraction, power spectrum density (PSD), SGS turbulent dissipation, SGS turbulent kinetic energy, and Reynolds stress is conducted. The transient multi-scale irregular vortex structures, fluctuation characteristics of liquid motion and spiral floating motion of bubbles can be well captured by LES. The Dynamic model and WALE model with C-w = 0325 give reasonably better agreement with the time-averaged experimental data. The transient fluctuation behavior of gas-liquid turbulent flow can be correctly captured by the current Euler-Euler LES model, whatever the SGS model. The classical -5/3 law of PSD of vertical liquid velocity is hold properly in the low frequency region. The slope of PSD is around -25/3 at the high frequency region. The Smagorinsky model and WALE model are sensitive to the model constant, at least for the energy transfer of high frequency. The molecular dissipation is negligible compared to the SGS turbulent dissipation. The Dynamic model gives higher predictions for fluctuation amplitude and time-averaged values of both SGS turbulent dissipation rate profile and SGS turbulent kinetic energy profile than other SGS models. A certain calculation time is needed to reach the dynamic stable state in LES, and the specific time is related to the scale of the calculation model. (C) 2020 Elsevier B.V. All rights reserved.