h-multigrid agglomeration based solution strategies for discontinuous Galerkin discretizations of incompressible flow problems