Using statistical learning to close two-fluid multiphase flow equations for bubbly flows in vertical channels