A comprehensive computational model based on the Eulerian-Eulerian approach was developed to simulate gas-liquid flows in a stirred vessel. A separate submodel was developed to quantitatively understand the influence of turbulence and presence of neighboring bubbles on drag acting on bubbles. This submodel was used to identify an appropriate correlation for estimating the interphase drag force. The standard k-epsilon turbulence model was used to simulate turbulent gas-liquid flows in a stirred vessel. A computational snapshot approach was used to simulate motion of the standard Rushton turbine in a fully baffled vessel. The computational model was mapped onto FLUENT4.5, a commercial CFD solver. The model predictions were compared with the previously published experimental data of Bombac and co-workers. The model was used to simulate three distinct flow regimes in gas-liquid stirred vessels: vortex clinging (VC), alternating small cavities (S33), and alternating large cavities (L33). The predicted results show reasonably good agreement with the experimental data for all three regimes. The computational model and results discussed in this work would be useful for understanding and simulating gas holdup distribution and flow regimes in stirred vessels. (c) 2006 American Institute of Chemical Engineers.