CFD and FE tools are intensively adopted by engine manufacturers in order to prevent thermo-mechanical failures reducing time- and cost-to market. The capability to predict correctly the physical factors leading to damages is hence essential for their application in the industrial practice. This is even more important for last generation SI engines, where the more and more stringent need to lower fuel consumption and pollutant emissions is pushing designers to reduce engine displacement in favor of higher specific power, usually obtained by means of turbocharging. This brings to a new generation of SI engines characterized by higher and higher adiabatic efficiency and thermo-mechanical loads.A recent research highlighted the different behavior of the thermal boundary layer of such engines operated at high revving speeds and high loads if compared to the same engines operated at low loads and revving speeds or even engines with a lower specific power. This means that CFD heat transfer models proposed and validated in the past decades on these last may not be predictive when applied to high specific power engines. This is why an alternative formulation was proposed in a previous work, for the estimation of the heat transfer in in-cylinder CFD simulations of high performance turbocharged SI engines.Nevertheless, for both the proposed alternative thermal wall function and the other ones available in literature, there are essential limitations due to the dimensionless distance y+. In fact, even if heat transfer models provide a further formulation for low y+ (viscous sub-layer), industrial practice seldom makes use of turbulence models enabling the integration up to the wall and low Reynolds approaches are even less used because of computational costs. Therefore in the present paper the authors aim to analyze critically the use of thermal wall functions along with high Reynolds turbulence models for the prediction of heat transfer in CFD in-cylinder simulations for different values of y+.