The heat transfer and hydrodynamics of particle flows in stirred tanks are investigated numerically in this paper by using a coupled CFD−DEM method combined with a standard k-e turbulence model. Particle−fluid and particle−particle interactions, and heat transfer processes are considered in this model. The numerical method is validated by comparing the calculated results of our model to experimental results of the thermal convection of gas-particle flows in a fluidized bed published in the literature. This coupling model of computational fluid dynamics and discrete element (CFD−DEM) method, which could calculate the particle behaviors and individual particle temperature clearly, has been applied for the first time to the study of liquid-solid flows in stirred tanks with convective heat transfers. This paper reports the effect of particles on the temperature field in stirred tanks. The effects on the multiphase flow convective heat transfer of stirred tanks without and with baffles as well as various heights from the bottom are investigated. Temperature range of the multiphase flow is from 340 K to 350 K. The height of the blade is varied from about one-sixth to one-third of the overall height of the stirred tank. The numerical results show that decreasing the blade height and equipping baffles could enhance the heat transfer of the stirred tank. The calculated temperature field that takes into account the effects of particles are more instructive for the actual processes involving solid phases. This paper provides an effective method and is helpful for readers who have interests in the multiphase flows involving heat transfers in complex systems.