Heat transfer problem in solid-dispersed two-phase flow is numerically studied. Temperature gradient within the finite-sized particles and inter-particle heat flux due to collisions are considered, and those effects on the flow structure and heat transfer are discussed. The interfacial flux model is extended to incorporate the heat conduction due to inter-particle contacts, based on 2-D and axisymmetric contact heat resistance solutions. The method is applied to 2-D and 3-D natural convection problems including multiple particles in a confined domain. Under high solid volume fraction conditions, the particles are observed to form densely concentrated regions, where heat flow tends to channel through the contacting points. In three-dimensional solid-dispersed flows, by decomposing the heat flux into the contributions of the convection and conduction, the change of the major heat transfer mode is studied for different solid volume fractions and conductivity ratios.