Includes abstract.Includes bibliographical references.The intimate contact achieved between the gas and liquid phases in bubble columns, coupled with the inherent efficient mixing these reactors offer, yield excellent heat and mass transfer characteristics. These attributes have been exploited commercially for decades, however, due to the complexity of the underlying hydrodynamics, the prediction of bubble columns based on empirical models can be unreliable outside of the operating ranges used to fit these models. Computational Fluid Dynamics (CFD) has emerged as an attractive tool for simulating these reactors and is based on numerically approximating the fundamentally based Navier-Stokes equations on a discretized domain. The application of CFD has become more practical as the cost of computational resources has declined and has lead to the establishment of three distinct modelling approaches which have been evaluated for the purpose of bubble column simulation in a number of research papers over the past two decades. Here the Euler-Euler approach has been recommended for the simulation of large scale columns, however, this approach is based on the most assumptions and yields the least amount of flow field information. The Euler-Lagrange approach treats bubbles as discrete particles which allows for the incorporation of a deterministic bubble size distribution and the direct consideration of heat and mass transfer effects. The most fundamental approach, Direct Numerical Simulation (DNS), predicts flow properties at the bubble scale, however, is extremely computationally expensive and is therefore only practically applicable to the investigation of a very small number of bubbles. The objective of this study is to contribute to the simulation of gasliquid flow interaction occurring in bubble columns by proposing a novel technique for simulating bubble scale flow information at a significantly reduced computational expense. For this purpose, it is proposed to predict the micro-flow fields around individual bubbles, within an Euler-Lagrange framework, with an algebraic model termed the Bubble Cell Model (BCM). The high gradient regions around individual bubbles are thereby accounted for with an algebraic flow model that can be rapidly evaluated as opposed to the two-phase partial differential Navier-Stokes equations, thereby reducing the numerical complexity of the problem. Since no such flow models currently exist and accuracy and fast evaluation are imperative, a statistical approach to the construction of the BCM is justified