Coupled DEM-LBM method for the free-surface simulation of heterogeneous suspensions