Simulation of bubble populations in a gas fluidized bed