A numerical investigation was conducted on the transient behavior of a hydrodynamically, fully developed, laminar flow of power-law fluids in the thermally developing entrance region of circular ducts taking into account the effect of viscous dissipation but neglecting the effect of axial conduction. In this regard, the unsteady state thermal energy equation was solved by using a finite difference method, whereas the steady state thermal energy equation without wall heat flux was solved analytically as the initial condition of the former. The effects of the power-law index and wall heat flux on the local Nusselt number and thermal entrance length were investigated. Moreover, the local Nusselt number of steady state conditions was correlated in terms of the power-law index and wall heat flux and compared with literature data, which were obtained by an analytic solution for Newtonian fluids. Furthermore, a relationship was proposed for the thermal entrance length. (C) 2008 Elsevier Ltd. All rights reserved.