Modelling of bubble columns operating at various flow conditions: effect of different options of simulation and different codes