We study super-Planckian near-field heat exchanges for multilayer hyperbolic metamaterials using exact S-matrix calculations. We investigate heat exchanges between two multilayer hyperbolic metamaterial structures. We show that the superPlanckian emission of such metamaterials can either come from the presence of surface phonon-polaritons modes or from a continuum of hyperbolic modes depending on the choice of composite materials as well as the structural configuration.In the last few years several fascinating experiments have demonstrated that for small separation distances compared with the thermal wavelength the thermal radiation exchanged between two hot bodies out of thermal equilibrium increases dramatically compared with what we observe at large distances and can even exceed the well-known Stefan-Boltzmann law by orders of magnitude 1-6 . Accordingly, thermal emission is in that case also called super-Planckian emission emphasizing the possibility to go beyond the classical black-body theory. There are several promising applications of super-Planckian emitters ranging from thermal imaging 7-9 and thermal rectification/management 10-12 to near-field thermophotovoltaics 13-17 . This has triggered many studies on the possibilities of tailoring and controlling the super-Planckian radiation spectrum by means of designing the material properties 18-23 , using phase-change materials 24 or 2D systems such as graphene 25,26 , for instance.Recently, it was shown that hyperbolic metamaterials can lead to broad-band photonic thermal conductance inside the material itself 27 and between two hyperbolic materials only separated by a vacuum gap 28 . Further Nefedov et al. considered nanorod-like structures