The evaluation of the mass transfer rates and the fluid-dynamics aspects of bubble columns are strongly affected by the intrinsic poly-dispersity of the gas phase, namely the different dispersed bubbles are usually distributed over a certain range of size and chemical composition values. In our previous work, gas-liquid systems were investigated by coupling Computational Fluid Dynamics with mono-variate population balance models (PBM) solved by using the quadrature method of moments (QMOM). Since mass transfer rates depend not only on bubble size, but also on bubble composition, the problem was subsequently extended to the solution of multi-variate PBM (Buffo et al. 2013). In this work, the conditional quadrature method of moments (CQMOM) is implemented in the open-source code OpenFOAM for describing bubble coalescence, breakage and mass transfer of a realistic partially aerated rectangular bubble column, experimentally investigated by Diaz et al.(2008). Eventually, the obtained results are here compared with the experimental data availabl