A hybrid reduced order method for modelling turbulent heat transfer problems