While slow granular flows have been an area of active research in recent years, heat transfer in flowing particulate systems has received relatively little attention. We employ a computational technique that couples the discrete element method (DEM), computational fluid dynamics (CFD), and heat transfer calculations to simulate realistic heat transfer in a rotary kiln. To maintain simplicity, while simulating the cylindrical kiln, we use a non-uniform grid in our code. Different materials, particle sizes, and rotation speeds are used to track the transition from convection-dominated heat transfer to conduction-dominated heat transfer. At low particle conductivities, the heat transfer is dominated by gas-solid conduction: however, at higher particle conductivities solid-solid conduction plays a more important role. Moreover, our results suggest that the rate of change of the average bed temperature can display a transition as the conductivity of the interstitial medium is increased. At low interstitial transport rates, such as in vacuum, high conductivity, high heat capacity particles get heated most rapidly, but with increased interstitial transport coefficients, lower heat capacity material may get heated faster despite lower values of conductivity. (c) 2008 Elsevier Ltd. All rights reserved.