Direct numerical simulations of bubbly flows