In many biological processes, increasing the rate of transport of a limiting nutrient can enhance the rate of product formation. In aerobic fermentation systems, the rate of oxygen transfer to the cells is usually the limiting factor. A key factor that influences oxygen transfer is bubble size distribution. The bubble sizes dictate the available interfacial area for gas-liquid mass transfer. Scale-up and design of bioreactors must meet oxygen transfer requirements while maintaining low shear rates and a controlled flow pattern. This is the motivation for the current work that captures multiphase hydrodynamics and simultaneously predicts the bubble size distribution. Bubbles break up and coalesce due to interactions with turbulent eddies, giving rise to a distribution of bubble sizes. These effects are included in the modeling approach by solving a population balance model with bubble breakage and coalescence. The population balance model was coupled to multiphase flow equations and solved using a commercial computational fluid mechanics code FLUENT 6. Gas holdup and volumetric mass transfer coefficients were predicted for different superficial velocities and compared to the experimental results of Kawase and Hashimoto (1996). The modeling results showed good agreement with experiment. (C) 2004 Elsevier Ltd. All rights reserved.