This study involves the numerical solution of the laminar heat transfer in a separating and reattaching flow by simulating the flow and heat transfer downstream of a backward‐facing step. The in‐house finite volume code has been implemented employing a hybrid differencing scheme and the SIMPLE algorithm for the pressure–velocity coupling. Three principal parameters governing heat transfer in this geometry, that is channel expansion ratio (ER), Reynolds number (Re), and Prandtl number (Pr), are systematically varied in the range ER = 1.111 to 2, Re = 1 to 200, and Pr = 0.71 to 100, and the simple correlations between these parameters have been elucidated. A series of important findings have been established by analyzing the results some of which are: (1) there is an associated shifting of the point of maximum heat transfer with respect to the flow‐reattachment point with gradually decreasing the values of ER and (2) the heat transfer enhancement increases with the increase in Pr number as a result of the compression of the thermal boundary layer and the maximum Nusselt number varies as .