We study the dynamics of gas-liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler-Lagrangian simulations. These simulation predictions are validated with Euler-Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard kappa-epsilon model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler-Lagrangian model is applicable at lower gas flow-rates. The Euler-Eulerian approach was found to be valid at lower as well as higher gas flow-rates. (C) 2010 Elsevier Ltd. All rights reserved.