This study proposes a CFD model to simulate large-scale bubble columns operating in different flow regimes. Transient 3-D simulations were performed employing a commercial code (ANSYS Fluent), and the numerical results were compared with available experimental data. The superficial gas velocity ranges between 0.0037 m/s and 0.2 m/s, covering both the mono-dispersed and pure-heterogenous flow regimes, where bubbles coalescence and breakup were modelled. The results have been critically analysed, and the discrepancies between the numerical and experimental results have been deeply commented on, setting the stage for future improvements.