Modelling of buoyancy-driven circulating bubbly flow