The ambitious performance goals set by the aerodynamic as well as the turbomachinery industries call for more accurate, scale-resolving simulation tools. A viable path towards industrial-level scale-resolving simulations of flows at high Reynolds numbers is through wall-modeled large-eddy simulation (WMLES). This talk will discuss the best practice in WMLES of problems in which heat transfer plays a role, covering topics including the gird resolution, the LES/wall-model matching location, the low Mach number limit, and the turbulent Prandtl number. Special attention is given to problems at high Mach numbers. Our results will answer the outstanding question: why our wall model works well in the context of WMLES yet the temperature transformations that give rise to our wall model fail at collapsing data.