Les écoulements de mélanges de deux fluides incompressibles, miscibles et ayant des densités très différentes sont gouvernés par les équations de Navier-Stokes non-homogène, couplées à une équation de convection-diffusion décrivant l'évolution de la composition du mélange. Nous proposons un algorithme associant la méthode des caractéristiques, pour la discrétisation des termes de transport, à la méthode des éléments finis avec une adaptation automatique de maillage, et démontrons que la solution de cet algorithme converge vers la solution exacte lorsque les pas de discrétisation tendent vers zéro. La robustesse de cet algorithme permet d'obtenir les premiers résultats de simulations numériques directes d'écoulements d'échange à très forte différence de densité, et de les valider par comparaison avec des expériences. Des écoulements de nuages denses sur des pentes sont simulés, et permettent d'analyser l'influence de la différence de densité sur les écoulements d'avalanches.

Cette thèse est consacrée aux écoulements de grains secs en régime dense. Nous nous restreigons aux géométries de cisaillement simple, sur plan incliné ou en conduite verticale, autorisant l'étude d'écoulement stationnaires et uniformes. Nous mesurons les grandeurs permettant de discuter la loi de comportement : profils de vitesse, de compacité, et distribution des forces de contact.<br />En conduite verticale, le profil de vitesse est mesuré à la paroi par imagerie rapide, et à l'intérieur de l'écoulement par imagerie par résonnance magnétique nucléaire. On met en évidence la localisation du cisaillement près des parois, et l'invariance de la forme du profil de vitesse avec le débit, signe d'un régime d'écoulement quasi-statique. On a aussi étudié l'influence de la rugosité sur le profil des vitesses, et montré le caractère intermittent de ce type d'écoulement.<br />Sur le plan incliné, la simulation numérique discrète (méthode de dynamique des contacts) permet d'étudier des écoulements d'assemblées de disques polydisperses frottants et complètement inélastiques. Nous avons mis en évidence la plage des écoulements stationnaires et uniformes, limité par une épaisseur d'arrêt dépendant de l'inclinaison. Nous avons observé un profil de compacité constante, avec une décroissance de la compacité moyenne lorsque l'inclinaison augmente. Le profil de vitesse est convexe, et montre une zone d'influence du socle sur les dix premières couches et une zone centrale en accord avec la prédiction du modèle de Bagnold. La vitesse moyenne suit une loi d'échelle en accord avec l'observation de Pouliquen. L'étude du réseau de contact met en évidence des anisotropies importantes ainsi qu'une statistique étalée des forces de contacts. Enfin nous avons montré que le tenseur des contraintes se restreint essentiellement à la contribution des forces de contacts, et que la pression est isotrope.<br />Ces observations conduisent à une discussion du rôle relatif des collisions, du frottement, de la compacité, du socle et du réseau de contacts, et une comparaison qualitative des prédictions de plusieurs modèles rhéologiques est proposée. L'analyse micromécanique des problèmes couplés de frottement et rotation permet de comprendre la rotation moyenne des grains, et la réduction notable du frottement effectif entre grains.