The paper discusses the latest improvements obtained in performing heat transfer calculations by RANS turbulence models when dealing with fluids at supercritical pressure.An algebraic heat flux model (AHFM) is adopted as an advanced tool for calculating the turbulent Prandtl number distribution to be used in the energy equation. Though maintaining a simple gradient approach, the proposed model manages to obtain interesting results when dealing with temperatures spanning from low to supercritical values.This is due to the introduction of a correlation for defining one of the relevant AHFM parameters. As stated in previous works, in fact, single fixed constant values could not be sometimes sufficient for dealing with very different operating conditions such as the ones occurring with supercritical fluids. In order to make the relation suitable for different fluids, a dependence on a non-dimensional quantity which proved to be relevant by parallel work is assumed.Some sensitivity analyses are also performed, showing some interesting capabilities in reproducing a sort of threshold behaviour when working in transition regions.Buoyancy induced phenomena are much better captured than in past attempts, though incomplete accuracy is observed for some boundary conditions.