The flow pattern in a bubble column depends upon the column diameter, height, sparger design, superficial gas velocity and the nature of gas-liquid system. In this paper, the effect of some of these parameters have been simulated using Computational Fluid Dynamics (CFD). The relationship of these parameters with the interphase force terms has been discussed. A complete energy balance has been established. Using this methodology, the flow patterns reported by Hills (1974), Menzel et al. (1990), Yao et al. (1991) and Yu and Kim (1991) have been simulated. Excellent agreement has been shown between the CFD predictions and the experimental observations. The above model has been extended to homogenization of an inert tracer. In order to confirm this model, mixing experiments were carried out in a 200 mm i bubble column. A radioactive tracer technique was used for the measurement of mixing time. Tc-99m ((99m) Tc), in the form of sodium pertechnate salt, was used as the liquid phase tracer. Good agreement has been shown between the predicted and the experimental values of mixing time. The model was further extended for the estimation of axial dispersion coefficient (D-L). Excellent agreement between the simulated and the experimental values of the axial dispersion coefficient confirms the predictive capability of the CFD simulations for the mixing process.