Thermal properties of complex suspension flows are investigated using numerical computations. The objective is to develop an efficient and accurate computational method to investigate heat transport in suspension flows. The method presented here is based on solving the lattice Boltzmann equation for the fluid phase, as it is coupled to the Newtonian dynamics equations to model the movement of particles and the energy equation to find the thermal properties. This is a direct numerical simulation that models the free movement of the solid particles suspended in the flow and its effect on the temperature distribution. Parallel implementations are done using MPI (message passing interface) method. Convective heat transfer in internal suspension flow (low solid volume fraction, φ<10%), heat transfer in hot pressing of fiber suspensions and thermal performance of particle filled thermal interface materials (high solid volume fraction, φ>40%) are investigated. The effects of flow disturbance due to movement of suspended particles, thermo-physical properties of suspensions and the particle micro structures are discussed.