The characteristics of flow and heat transfer in a bubbling fluidized bed are investigated by means of computational fluid dynamics (CFD). To simulate two-phase flow for the gas and solid flows, Eulerian-Eulerian approach is applied. Attention is paid for a heat transfer from the wall to fluidized bed by bubbling motion of the flow. From the result, it is confirmed that heat transfer is promoted by chaotic bubbling motion of the flow by enhancement of mixing among solid particles. In particular, the vortical flow motion around gas bubble plays an important role for the mixing and consequent heat transfer. Discussion is made for the time and space averaged Nusselt number which shows peculiar characteristics corresponding to different flow regimes.