L'objectif de ce travail est de contribuer à la modélisation numérique des écoulements turbulents gaz-particules en conduite horizontale ou verticale non isotherme, présents dans de nombreux procédés industriels (transport pneumatique, séchage, combustion,...). Le modèle repose sur une approche eulérienne-lagrangienne permettant une description fine des mécanismes d'interactions entre les deux phases (action du fluide sur les particules (dispersion), action des particules sur le fluide (couplage « two-way ») interactions inter-particulaires (collisions)), plus ou moins influents selon les caractéristiques de l'écoulement. L'influence de la turbulence de l'écoulement gazeux sur le mouvement d'une particule est simulée par un modèle de dispersion anisotrope permettant de générer les fluctuations de vitesses et de température du fluide vu par une particule. Les développements numériques apportés au modèle en conduite ont été validés par comparaison avec les résultats expérimentaux disponibles dans la littérature. Les différents tests de sensibilité ont permis de mettre en évidence l'influence du modèle de dispersion, des collisions et de la modulation de la turbulence (actions directes et indirectes des particules sur le fluide) sur le comportement dynamique et thermique de la suspension. Le modèle est capable de prédire les échanges thermiques en présence de particules pour une large gamme d'écoulements à la fois en conduite verticale et en conduite horizontale. Cependant, des difficultés d'ordre numérique subsistent pour les simulations en couplage « two-way » en présence de très petites particules pour des taux de chargement supérieurs à 1. Ceci est lié aux problèmes de modélisation des termes de couplage entre les deux phases (en particulier les coefficients de modélisation de l'équation de transport de la dissipation, Ce2 et Ce3 ).