This work presents an Euler-Euler (E-E) Computational Fluid Dynamics (CFD) model developed for flow inside a rectangular air-water bubble column, in which we have also implemented population balance method (PBM) based evolution of the various bubble size classes using the Homogeneous Discrete method, Inhomogeneous Discrete method, Quadrature Method of Moments (QMOM) and Direct Quadrature Method of Moments (DQMOM) technique. The results are validated with liquid velocity field measurements obtained from the Radioactive Particle Tracking (RPT) technique. The study presents a comparative analysis of the effect of incorporating various interfacial closures like drag force, lift force and virtual mass effect in the E-E CFD coupled with PBM framework. The importance of the lift force in the transient three-dimensional simulations for predicting the time-averaged velocity profiles is highlighted. (c) 2012 Elsevier B.V. All rights reserved.