In this work, two powerful methods are combined for bubble column simulations namely, the volume of fluid and the discrete bubble model. While the former method takes care of the free surface, the discrete bubble model tracks and handles the dynamics of the dispersed bubbles. The hybrid model presented in this work is verified and validated with existing established experimental results. A model parameter study for bubble break-up and coalescence is performed to find the optimum values of the model parameters, i.e. the critical Weber number and the coalescence calibration factor.