Numerical simulations and kinetic theory for simple shear motion of bubbly suspensions