Exact temperature solutions and theoretical Nusselt curves, valid for Péclét numbers ranging from 1 to ∞, were obtained for thermal-entry-region heat transfer for laminar flow through concentric annuli, subject to a step jump in wall heat flux at z = 0. To allow for the effect of axial conduction, which is significant at low Péclét numbers, the inlet fluid temperature was taken to be uniform at z = −∞, and the first twenty eigenconstants were computed for the adiabatic region (−∞ < z ⩽ 0) and the heated region (0 ⩽ z < ∞), separately. By constructing two sets of orthonormal functions from the non-orthogonal eigenfunctions, the series expansion coefficients were then determined such that both the temperatures and longitudinal temperature gradients for the two regions match at z = 0. The temperature solutions corresponding to the limiting case of N pe = ∞ show excellent agreement with those reported by Lundberg et al. [4], who analyzed the entry-region problem by neglecting axial conduction.