L'optimisation d'une opération de mélange de poudre repose essentiellement sur un travail expérimental à l'échelle du laboratoire qui doit pouvoir être transposer aux mélangeurs de plus grandes tailles. Définir des lois d'extrapolation et améliorer notre connaissance de la dynamique du mélange est donc nécessaire. Dans ces travaux, la dynamique de mélange au sein des mélangeurs Turbula® a été étudiée en s'appuyant sur l'analyse des cinétiques de mélange et des fonctions d'autocorrélation pour plusieurs systèmes particulaires. Selon les vitesses de rotation de l'axe moteur, 3 régimes d'écoulement ont été définis et les principaux mécanismes de mélange et de ségrégation apparaissant pour chacun de ces régimes ont été identifiés en lien avec les propriétés d'écoulement des produits. Dans un deuxième temps, les qualités de mélange obtenues dans différentes tailles de mélangeurs ont été comparées sur la base du principe des similitudes afin de mettre en évidence les facteurs ayant une influence lors du passage d'une taille de mélangeur à une autre. Enfin dans le cadre d'une application industrielle, une méthodologie s'appuyant sur l'intensité de ségrégation et l'autocorrélation spatiale, a été développée pour identifier des défauts d'homogénéité au sein de plaques bipolaires composites / The optimization of a powder mixing step typically involves an experimental work at lab scale in order to be transposed to larger mixers. Defining scale-up laws and improving our knowledge of the mixing dynamics remains some of the mains industrial issues of this century. In this work, the mixing dynamics of several particulate systems has been studied in Turbula mixers thanks to the analysis of mixing kinetics and autocorrelation functions. According to the engine speed, three flow regimes have been defined. The corresponding main mixing and segregation mechanisms at play for each of these regimes have been identified in relation with the flow properties of the products. In a second phase, the qualities of the mixtures obtained in the different mixer sizes have been compared on the basis of the principle of similarities in order to shed light the factors which influencing scale-up. Finally, as part of an industrial application, a methodology has been developed using the concept of intensity of segregation and the spatial autocorrelation tools to identify heterogeneities in bipolar plates made of composite materials

Ce travail est consacré à la modélisation d'écoulements de deux fluides immiscibles et son application en microfluidique. Pour cela, nous mettons en oeuvre des méthodes Level Set actuelles permettant un suivi précis de l'interface, dont le mouvement est induit par des champs de vitesse vérifiant les équations de Stokes ou de Navier-Stokes munies d'un terme de tension de surface.<br />Dans une première partie, nous abordons la problématique du suivi d'interface et présentons en détail les composantes de la méthode Level Set. En particulier, nous détaillons les approches ENO et WENO pour discrétiser les équations de Hamilton-Jacobi ainsi que les diverses méthodes existantes de redistanciation.<br />Dans la deuxième partie, nous traitons de l'analyse et de la résolution numérique des écoulements bifluides incompressibles pilotés par la tension de surface. Après avoir décrit les modèles mathématiques ainsi que leurs discrétisations et solveurs, nous apportons une contribution nouvelle en dérivant théoriquement une condition de stabilité valable pour les nombres de Reynolds faibles à modérés, caractéristiques des configurations microfluidiques. De plus, on introduit une méthode de décomposition de l'écoulement qui permet de diminuer les temps de simulation.<br />Enfin, la troisième partie est consacrée à l'application des outils évoqués précédemment pour simuler la dynamique de gouttes dans des microcanaux. Nous présentons les résultats numériques obtenus avec d'une part, un code bidimensionnel cartésien et d'autre part, avec un code tridimensionnel axisymétrique que nous avons entièrement développés. Une bonne adéquation est obtenue relativement aux expériences microfluidiques du laboratoire LOF (Rhodia - CNRS). En particulier, avec nos simulations, nous mettons à jour différentes dynamiques de mélange au sein des gouttes.

La méthode de Boltzmann sur réseau est une formulation discrète particulière de l'équation de Boltzmann. Depuis ses débuts, il y a trente ans, cette méthode a gagné une certaine popularité, et elle est maintenant utilisée dans presque tous les problèmes habituellement rencontrés en mécanique des fluides notamment pour les écoulements multi-espèces. Dans le cadre de ce travail, une force de friction intermoléculaire est introduite pour modéliser les interactions entre les molécules de différent types causant principalement la diffusion entre les espèces. Les phénomènes de dissipation visqueuse (collision usuelle) et de diffusion moléculaire (force de friction intermoléculaire) sont séparés et peuvent être ajuster indépendamment. Le principal avantage de cette stratégie est sa compatibilité avec des optimisations de la collision usuelle et les opérateurs de collision avancés. Adapter un code mono-espèce pour aboutir à un code multi-espèces est aisé et demande beaucoup moins d’effort comparé aux précédentes tentatives. De plus, il n’ y a pas d’approximation du mélange, chaque espèce a ses propres coefficients de transport pouvant être calculés à l’aide de la théorie cinétique des gaz. En général, la diffusion et la convection sont vus comme deux mécanismes séparés : l’un agissant sur la masse d’une espèce, l’autre sur la quantité de mouvement du mélange. En utilisant une force de friction intermoléculaire, la diffusion et la convection sont couplés par l’intermédiaire la quantité de mouvement de chaque espèce. Les mécanismes de diffusion et de convection sont intimement liés dans de nombreux phénomènes physique tel que la digitation visqueuse.L’instabilité de digitation visqueuse est simulée en considérant dans un milieu poreux deux espèces dans des proportions différentes soit un mélange moins visqueux déplaçant un mélange plus visqueux. Les principaux moteurs de l’instabilité sont la diffusion et le contraste de viscosité entre les espèces. Deux stratégies sont envisagées pour simuler les effets d’un milieu poreux. Les méthodes de rebond partiel et de force de Brinkman bien que basées sur des approches fondamentalement différentes donnent dans notre cas des résultats identiques. Les taux de croissance de l’instabilité calculés à partir de la simulation coïncident avec ceux obtenus à partir d’analyses de stabilité linéaire. L’évolution de la longueur de mélange peut être divisée en deux étapes dominées d’abord par la diffusion puis par la convection. La physique de la digitation visqueuse est ainsi correctement simulée. Toutefois, les effets de diffusion multi-espèces ne sont généralement pas pris en compte lors de la digitation visqueuse de trois espèces et plus. Ces derniers ne sont pas négligeable puisque nous mettons en avant une configuration initialement stable qui se déstabilise. La diffusion inverse entraîne la digitation dont l’impact dépend de la diffusion entre les espèces. / The lattice Boltzmann method (LBM) is a specific discrete formulation of the Boltzmann equation. Since its first premises, thirty years ago, this method has gained some popularity and is now applied to almost all standard problems encountered in fluid mechanics including multi-component flows. In this work, we introduce the inter-molecular friction forces to take into account the interaction between molecules of different kinds resulting primarily in diffusion between components. Viscous dissipation (standard collision) and molecular diffusion (inter-molecular friction forces) phenomena are split, and both can be tuned distinctively. The main advantage of this strategy is optimizations of the collision and advanced collision operators are readily compatible. Adapting an existing code from single component to multiple miscible components is straightforward and required much less effort than the large modifications needed from previously available lattice Boltzmann models. Besides, there is no mixture approximation: each species has its own transport coefficients, which can be calculated from the kinetic theory of gases. In general, diffusion and convection are dealt with two separate mechanisms: one acting respectively on the species mass and the other acting on the mixture momentum. By employing an inter-molecular friction force, the diffusion and convection are coupled through the species momentum. Diffusion and convection mechanisms are closely related in several physical phenomena such as in the viscous fingering instability.A simulation of the viscous fingering instability is achieved by considering two species in different proportions in a porous medium: a less viscous mixture displacing a more viscous mixture. The core ingredients of the instability are the diffusion and the viscosity contrast between the components. Two strategies are investigated to mimic the effects of the porous medium. The gray lattice Boltzmann and Brinkman force models, although based on fundamentally different approaches, give in our case equivalent results. For early times, comparisons with linear stability analyses agree well with the growth rate calculated from the simulations. For intermediate times, the evolution of the mixing length can be divided into two stages dominated first by diffusion then by convection, as found in the literature. The whole physics of the viscous fingering is thus accurately simulated. Nevertheless, multi-component diffusion effects are usually not taken into account in the case of viscous fingering with three and more species. These effects are non-negligible as we showcase an initial stable configuration that becomes unstable. The reverse diffusion induces fingering whose impact depends on the diffusion between species.