Numerical analysis is made on the turbulent heat transfer with suspension of solid particles in circular tube with constant heat flux. The mean motion of suspending particles in mixture is treated as the secondary gas flow with virtual density and viscosity. Our modeling of turbulent transport phenomena of suspension flow is based on this assumption and conventional mixing length theory. This paper gives the evidence that the mixing length models can be extended to close the governing equations for two phase turbulent flow with solid boundary at a first order level. Results on Nusselt numbers obtained by analytical treatments are compared with available experimental data and discussed. They suggest that the most important parameters of two phase turbulent heat transfer phenomena are relative particle diameter to pipe diameter, gas-solid loading ratio, and specific heat of suspending material.