Long-lived high to ultra-high temperature (HT-UHT) granulitic terranes formed throughout Earth’s history. Yet, the detailed processes involved in their formation re- main unresolved and notably the sequence of appearance and duration of migmatisa- tion and granulites conditions in the orogenic cycle. These processes can be evaluated by analytical and numerical models. First, solving the steady-state heat equation allows underlining the interdependency of the parameters controlling the crustal geotherm at thermal equilibrium. Second, performing two-dimensional thermo-mechanical experiments of an orogenic cycle, from shortening to gravitational collapse, allows to consider non-steady-state geotherms and understand how deformation velocity may affect the relative timing of migmatite and granulite formation. These numeri- cal experiments with elevated radiogenic heat production and slow shortening rates allow the formation of large volumes of prograde migmatites and granulites going through the sillimanite field as observed in many HT-UHT terranes. Finally, the in- terplay between these parameters can explain the difference in predicted pressure- temperature-time paths that can be compared with the natural rock archive.