Abstract Thermal-hydraulics is recognized as a key safety challenge in the development of liquid metal cooled reactors. At nominal operating conditions, the Prandtl number of liquid metals which are used as primary coolants, such as lead and sodium, is very low: typically of the order of 0.025–0.001. Obtaining an accurate prediction of the turbulent heat transfer at such a low Prandtl number is not an easy task for the standard turbulence models and has challenged the modellers over several decades. In the framework of the EU SESAME project, an effort has been put forward to assess and/or further develop/calibrate different turbulent heat flux closures. In this regard, the present article reports an assessment of four different turbulent heat flux closures for applications involving low-Prandtl fluids. These closures include: (i) the Reynolds analogy based on a constant turbulent Prandtl number (ii) a four-equation explicit algebraic heat flux model (AHFM) (ii) a three-equation implicit AHFM called AHFM-NRG and (iv) a non-linear second-order heat flux model called Turbulence Model for Buoyant Flows (TMBF). The performance of these turbulence models has been assessed in three different test cases against high-fidelity numerical reference data been generated within the SESAME project. The three test cases are: a natural Rayleigh-Benard convection flow, a mixed convection planar channel flow and a forced convection impinging jet flow. The shortcomings of the classical Reynolds analogy approach for low-Prandtl fluids in all flow regimes are highlighted; hence, more advanced and well-calibrated closures are recommended.