A combination of unsteady 3D CFD simulations and experimental temperature measurements was employed to determine the effect of loading patterns on the temperature distribution within steel forgings inside a gas-fired heat treatment furnace. This was aimed to obtain a more homogenous temperature distribution. Besides, a hybrid methodology using 3D numerical simulations and a high-resolution dilatometer allowed improving residence time of forgings inside the heat treatment furnace. The influence of the loading patterns and skids on temperature distribution and residence time of forgings was examined using thermal analysis four different loading patterns.Comprehensive unsteady thermal analysis of the products heating allowed quantifying the impact of skids usage and their dimensions on the extent of the uniformity of temperature distribution.The results were interpreted in terms of the inter-relationship between the skids usage, their geometry, absorbed radiation and convective heat fluxes. The analysis showed that temperature non-uniformities of up to 331 K could be produced for non-optimum loading patterns. Using the developed CFD approach it was possible to reduce the temperature non-uniformity of different sizes of blocks up to 32% via changing the loading pattern inside the furnace. Further, the slab's residence time was improved by almost 15.5% when employing the proposed hybrid approach. This approach could directly be applied to the optimization of different heat treatment cycles of forged blocks in different grades of steels.