Deriving closures for bubbly flows using direct numerical simulations by