Population balance modelling of bubble columns under the heterogeneous flow regime