The steady, two-dimensional and incompressible flow of power-law fluids across an unconfined isothermal heated circular cylinder is investigated numerically to ascertain the effect of temperature-dependent viscosity on the flow and forced convection heat transfer phenomena. Extensive numerical results elucidating the variation of the heat transfer characteristics and drag coefficient on the severity of temperature dependence of viscosity (0 <= b <= 0.5), power law index (0.6 <= n <= 1.6), Prandtl number (1 <= Pr <= 100) and Reynolds number (1 <= Re <= 30) are presented. The coupled momentum and energy equations are expressed in the stream function/vorticity formulation and solved using a second-order accurate finite difference method to determine the local and surface-averaged Nusselt numbers, the drag coefficient, and to map the flow domain in terms of the temperature and flow fields near the cylinder. The variation of viscosity with temperature is shown to have a substantial effect on both the local and surface-averaged values of the Nusselt number. As expected, the results also suggest that the rate of heat transfer shows positive dependence on the Reynolds number and Prandtl number. Furthermore, stronger the dependence of viscosity on the temperature, the greater is the enhancement in the rate of heat transfer. Finally, all else being equal, shear-thinning fluid behaviour facilitates heat transfer while the shear-thickening behaviour has deleterious effect on heat transfer. (C) 2010 Elsevier Ltd. All rights reserved.