Turbulent heat transfer is a complex phenomenon that has challenged turbulence modellers over various decades. In this regard, in the recent past, several attempts have been made for the assessment and further development/calibration of the available turbulent heat flux modelling approaches. One of the main hampering factors with respect to the further assessment of these modelling approaches is the lack of reference data. In the framework of the EU SESAME and MYRTE projects, an extensive effort has been put forward to generate a wide range of reference data, both experimental and numerical, to fill this gap. In that context, this article reports the numerical database that has been generated within these projects for various liquid metal flow configurations in different flow regimes. These high fidelity numerical data include seven different flow configurations: a wall-bounded mixed convection flow at low Prandtl number with varying Richardson number (R_i) values; a wall-bounded mixed and forced convection flow in a bare rod bundle configuration; a forced convection confined backward facing step (BFS) with conjugate heat transfer; a forced convection impinging jet for three different Prandtl fluids corresponding to two different Reynolds numbers of the fully developed planar turbulent jets; a mixed-convection cold-hot-cold triple jet configuration corresponding to R_i = 0.25; an unconfined free shear layer for three different Prandtl fluids; and a forced convection infinite wire-wrapped fuel assembly. These high-fidelity numerical databases will serve the further development of turbulent heat transfer models by providing unique, new and detailed data for the thermal-hydraulic behaviour of liquid metals in various flow configurations.