Predicting turbulent heat transfer in high-speed flows dominated by shock waves poses significant challenge to computational fluid dynamics (CFD) methods [1, 2]. The majority of the turbulence models are based on the gradient diffusion hypothesis, where the turbulent heat flux is modeled in terms of turbulent conductivity and turbulent Prandtl number (PrT ). Conventional turbulence models in Reynolds-averaged Navier-Stokes simulations use a constant PrT value of 0.89; they often overpredict the wall heat transfer in shock-boundary layer interaction (SBLI). A few variable PrT models proposed in literature improve the heat flux prediction, but are computationally intensive [3, 4].