In This study, a computer model, based on discrete element method, is used to simulate unsteady state heat transfer at the fuser/toner and toner/coating layer interfaces during the Xerography fusing process. The model coating layers consisted of randomly arranged spherical pigment and latex particles with commercially relevant size distributions. Effects of coating characteristics, toner size, multiple toner layers, toner melting energy and toner thermal conductivity on the unsteady state heat transfer in the fusing process were investigated. Results showed that temperature variation highly depended on the toner size, toner melting energy and the fuser roll temperature. Moreover, simultaneous coupling of the stress and heat transfer indicated that the pressure exerted by fuser roll cannot significantly affect the rate of heat transfer to the toner particles.