The heat transfer between a vertical pipe wall and turbulent gas-particle flow is numerically investigated according to the Eulerian-Lagrangian approach and the k-epsilon turbulence model. The particles are introduced homogeneously into the simulation volume by a unique technique referred to as an artificial feeding volume. The numerical code using additional computer programs is validated with available experimental results for the constant heat flux boundary condition. An average deviation of about 4% and a maximum deviation of about 7% were attained from the numerical predictions for various particle and pipe diameters. The effect of the geometrical parameters and the flow parameters on the gas/particle temperature, the convection heat transfer coefficient between the wall and the gas-particle mixture, and the thermal entry length were studied. An increase in particle diameter ( loading ratio approximate to 0.5) extended the thermal entry length and decreased the bulk mixed temperature, particle temperature, and convection heat transfer coefficient. Increasing the pipe diameter led to a significant reduction in bulk mixed temperature and thermal entry length, in addition to a decrease in particle temperature and Nusselt number. Increasing the loading ratio up to 2.36 led to a reduction in wall temperature and bulk mixed temperature, in addition to an increase in the convective heat transfer coefficient and thermal entry length.