Numerical simulations based on the Eulerian-Eulerian approach have been performed in the study of interphase heat transfer in a gas solid fluidized bed. The kinetic theory of granular flow (KTGF) has been used to describe the solid phase rheology. An assessment of drag models in the prediction of heat transfer coefficients shows that no major difference is observed in the choice of the drag model used. Fluctuations of the interphase heat transfer coefficient have been found to be closely related to the bubble motion in the bed. Effects of the wall boundary condition, inlet gas velocity, initial bed height and particle size on the predicted heat transfer coefficient have also been investigated. Typical temperature profiles in the bed show that thermal saturation is attained instantaneously close to the gas distributor. Simulated results of the coefficients are in fair agreement with those reported in literature.