The newly proposed one-equation Wray-Agarwal (WA) turbulence model coupled with Euler-Euler (E-E) approach was used to simulate the gas-liquid flow in a square bubble column, which further validated the accuracy and expanded the application of the WA model comparing to a preliminary study (Ouyang et al., 2019). All simulation results of the WA model were compared with two widely used turbulence models of standard k-ε and SST k-ω and literature data. The numerical results of the time-averaged axial liquid and gas velocity were validated with experimental data, and it was found that the WA model showed overall the best agreement. Profiles of turbulent kinetic energy, turbulent dissipation rate, turbulent viscosity and gas holdup were investigated at three different heights. Results for the first time validate the WA model based on the turbulent characteristic values of gas-liquid flow in the bubble column, which shows its promising prospective in multiphase flow simulation.