Population balance modelling of bubble columns under the heterogeneous regime