Développement d'un modèle Euler-Lagrange robuste pour la simulation des écoulements solide-liquide dans les opérations de mélange

Les opérations de mélange solide-liquide en cuves agitées jouent un rôle clef dans de nombreux procédés, que ce soit dans la fabrication de produits alimentaires, pharmaceutiques, cosmétiques, ou pour promouvoir l'homogénéité de suspensions, ce qui est particulièrement vital au bon fonctionnement des réacteurs chimiques employant un catalyseur solide. Malgré leur importance indéniable pour l'industrie chimique et les efforts considérables qui ont été déployés afin de mieux les comprendre, la conception et l'optimisation de ces opérations demeurent un grand défi. En effet, la quasi-totalité de la littérature se concentre sur le régime d'opération pleinement turbulent, malgré le fait que de nombreux procédés industriels soient opérés en régimes laminaire ou transitoire. Malheureusement, la littérature sur ces derniers régimes d'opération est quasi-inexistante. De plus, bien que le régime d'opération turbulent ait fait l'objet d'un grand nombre d'études, celles-ci se sont principalement concentrées sur la prédiction, à l'aide de corrélations empiriques ou semi-empiriques, de la vitesse nécessaire pour la suspension complète des particules (Njs), c'est-à dire la vitesse d'agitation nécessaire pour suspendre toutes les particules hors du fond de la cuve. Cependant, de nombreux procédés pourraient être opérés dans des conditions différentes telles que la suspension partielle ou complètement homogène. Le premier cas permettrait d'économiser de l'énergie et d'éviter des contraintes trop fortes sur l'agitateur, tandis que le second assurerait une cinétique de réaction et une qualité de produit nettement mieux contrôlée. Dans ces deux situations, la connaissance de Njs n'est que d'une aide limitée. Pour mieux opérer et concevoir ces unités, il est nécessaire de pouvoir prédire la distribution et la dispersion des particules solides ainsi que les patrons d'écoulement au sein de la cuve. L'étude des systèmes de mélange solide-liquide représente une difficulté importante, car l'opacité des suspensions solide-liquide limite fortement la mesure de variables locales telles que les profils de concentration. Ainsi, la majeure partie des travaux expérimentaux n'ont mesuré que des paramètres globaux, tels que le couple sur l'agitateur, la fraction de particules suspendues ou Njs. La mécanique des fluides numérique (CFD), complémentaire à l'expérience, permet quant à elle d'investiguer à la fois les paramètres globaux, mais aussi ce qui se passe localement en tout point de la cuve. Cependant, la modélisation d'écoulements multiphasiques renferme de nombreux défis compte tenu de l'interaction multiéchelle (de temps et d'espace) entre les phases. Les nombreux modèles capables de modéliser ces types d'écoulements sont décrits dans cette thèse, dans l'objectif de faire ressortir leurs forces ainsi que leurs limites. Parmi ceux-ci, il est montré que seuls les modèles à deux fluides, où la phase solide est modélisée comme un second fluide, ont été utilisés exhaustivement pour aborder le mélange solide-liquide. Cependant, ces modèles souffrent d'une incapacité à bien décrire les régimes rapides et denses d'écoulement granulaire (comportement de Burnett et de super Burnett) ainsi que de plusieurs difficultés à saturer la concentration de solide lorsque les particules sont à leur fraction maximale d'empilement. La CFD-DEM, une famille de modèles relativement récents qui combinent la CFD pour la phase fluide et la méthode des éléments discrets (DEM) pour les particules solides permet quant à elle de décrire la dynamique de la phase granulaire avec un grand degré de précision et a largement fait ses preuves dans l'étude de milieux solide-gaz. Cependant, cette méthode n'a jamais été employée rigoureusement pour l'étude d'écoulement solide-liquide dans des géométries complexes telles que des mélangeurs. Pour que ceci soit possible, de nombreux développements mathématiques sont nécessaires afin de s'assurer que le schéma soit stable, qu'il converge (en temps et en espace) et qu'il soit capable de simuler des géométries complexes en mouvement tels que les agitateurs. La conception d'un tel modèle et son application à l'étude de la dynamique du mélange solide-liquide et de la mise en suspension de particules solides est l'objectif principal de cette thèse. En premier lieu, un schéma volume fini de type Pressure Implicit with Splitting of Operator (PISO) pour résoudre les équations de Navier-Stokes moyennées volumiquement (VANS), nommé PISO-VANS, est établi. La résolution de ces équations est une partie essentielle d'un modèle CFD-DEM applicable à des écoulements concentrés. Afin de vérifier la cohérence du schéma PISO-VANS pour la résolution des équations VANS, une méthodologie basée sur la méthode des solutions manufacturées est développée afin d'établir des cas tests analytiques permettant d'effectuer des tests de convergence numérique. Ces tests, les premiers de ce genre, démontrent que le schéma proposé converge à la précision désirée, c'est-à-dire qu'il est bien de second ordre en temps et en espace. Ensuite, cette méthodologie est employée pour vérifier un nouveau schéma permettant de résoudre les équations VANS avec la méthode de Boltzmann sur réseau (LBM). Ce schéma est basé sur un nouvel opérateur de collision. Il est démontré, à l'aide d'une analyse de Chapmann-Enskogg, que la formulation proposée permet de retrouver les équations VANS. Cet opérateur est le premier permettant de résoudre les équations VANS avec la LBM lorsque la fraction volumique n'est pas constante dans l'espace, une capacité essentielle pour l'étude d'écoulements polyphasiques ou dans des milieux poreux. Dans la troisième partie de ce travail, une nouvelle méthode de condition immergée semiimplicite permettant de modéliser des corps rigides en rotation est développée. Cette méthode est conçue pour bien s'harmoniser avec le schéma PISO et pour être fonctionnelle sur un maillage non structuré polyédrique tout en demeurant parallèle. Cette méthode est tout d'abord vérifiée sur des cas tests académiques tels que l'allée de von Karman derrière un cylindre, ainsi qu'un écoulement de Taylor-Couette entre deux cylindres. Il est montré que le schéma peut bien reproduire les vortex de von Karman, mais qu'il dégrade l'ordre du schéma volume fini de 2 à 1.33 dans le cas de l'écoulement de Taylor-Couette. Le schéma est finalement validé expérimentalement et comparé à d'autres méthodes numériques permettant de simuler des géométries en rotation. Un accord quasi parfait est obtenu. La quatrième partie de ce travail résulte directement de la combinaison du schéma volume fini PISO-VANS avec la méthode de conditions immergées afin de simuler le mélange solideliquide, du démarrage au régime permanent, à l'aide de la CFD-DEM. Différentes stratégies de couplage entre les phases sont testées et il est montré que, contrairement au cas gaz-solide, un couplage explicite est préférable, car il atténue les erreurs de moyenne volumique. Il est aussi démontré qu'un modèle de rhéologie est nécessaire afin de considérer la dissipation à l'échelle inférieure à une taille de maille. Le modèle complet est ensuite validé qualitativement, en comparant des profiles particuliers d'écoulements obtenus au début de la suspension des particules, et quantitativement, à travers la comparaison avec l'expérience de la fraction de particules suspendues mesurée par la technique de pression de jauge (PGT). À nouveau, un excellent accord est obtenu entre les données expérimentales et les résultats du modèle. Dans la cinquième partie de ce travail, le modèle CFD-DEM est modifié afin de permettre l'étude des écoulements turbulents à l'aide de la simulation aux grandes échelles turbulentes (LES). Suite à une validation avec l'expérience, deux nouvelles techniques de mesure de la fraction de particules suspendues, l'analyse lagrangienne de fraction suspendue (LSFA) et l'analyse de fraction de décorrélation (DFA), sont introduites. Les résultats issus de ces méthodes sont ensuite comparés à ceux obtenus numériquement et expérimentalement par la méthode de pression de jauge. Il est montré que ces deux méthodes sont pratiquement aussi précises que la PGT, mais qu'elles sont aussi plus versatiles, car elles peuvent être appliquées à toutes les géométries et ne nécessitent pas de simulations sur des larges plages de vitesse d'agitation. Dans la sixième partie de ce travail, le modèle CFD-DEM est utilisé pour étudier en détail le mélange solide-liquide en régime laminaire et transitoire. Notamment, l'impact du dégagement de l'agitateur et de la présence de chicanes sur la fraction de particules suspendues et la dynamique du mélange solide-liquide est établi. Il est montré que de réduire le dégagement au fond de l'agitateur permet de prévenir l'apparition de zone mortes à haute vitesse. De surcroît, une étude de sensibilité sur les paramètres de la DEM est effectuée et montre que seule la friction entre les particules joue un rôle significatif sur la dynamique de la phase solide. Finalement, une brève discussion permet de résumer les résultats obtenus et de donner de nombreuses pistes de travaux pouvant faire suite à ce qui fut développé dans le cadre de cette thèse.