Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

Cette étude utilise des simulations numériques d'un réservoir agité contenant des liquides miscibles pour démontrer que, bien que le temps de mélange soit corrélé positivement au nombre de Richardson, une loi d'échelle exponentielle dérivée et fondée sur les nombres de puissance, de Froude et de Richardson permet de rassembler toutes les données sur une seule courbe maîtresse pour des différences de densité intermédiaires à grandes.

L'essentiel

Imaginez que vous essayez de mélanger une cuve géante contenant deux liquides différents : un sirop lourd et épais au fond, et un jus plus léger et plus fluide au-dessus. Vous y plongez une énorme pale tournante (une turbine) pour les agiter ensemble.

Dans le monde réel, c'est une tâche courante dans les usines produisant tout, des médicaments au traitement des eaux usées. Mais voici le hic : comme les liquides ont des poids différents (densités), le lourd veut rester au fond et le léger veut flotter au-dessus. Cela crée une « lutte » entre la pale tournante qui tente de les mélanger et la gravité qui tente de les maintenir séparés.

Ce papier est comme une histoire policière où des scientifiques ont utilisé de puissantes simulations informatiques pour déterminer exactement combien de temps il faut pour mélanger parfaitement ces deux liquides, et comment prédire ce temps sans avoir à construire un réservoir physique et à effectuer des tests coûteux à chaque fois.

Le Déroulement : Une Cuisine d'Essai Numérique

Les chercheurs ont construit une version virtuelle d'un réservoir de mélange industriel standard.

• Le Réservoir : C'est un grand cylindre avec des parois et quatre ailettes verticales (chicanes) pour empêcher le liquide de simplement tourner en rond comme une rivière paresseuse.
• La Pale : Une lame tournante au centre.
• Les Liquides : Ils ont simulé un mélange 50/50 d'un liquide lourd et d'un liquide plus léger. Ils n'ont pas utilisé de vrais produits chimiques ; ils les ont simplement traités comme des fluides « lourds » et « légers » ayant la même épaisseur (viscosité).
• La Méthode : Au lieu d'utiliser des équations mathématiques standard, ils ont employé une astuce ingénieuse appelée la Méthode de Boltzmann sur Réseau. Imaginez cela comme simuler le liquide non pas comme une masse continue, mais comme des milliards de minuscules billes de billard invisibles qui rebondissent et entrent en collision. Cela leur a permis de voir exactement comment la turbulence (le tourbillon chaotique) se comportait.

La Grande Question : À quelle vitesse pouvons-nous mélanger ?

L'objectif principal était de trouver une « formule magique » pour prédire le temps de mélange.

• Les Variables : Ils ont modifié deux choses principales :
  1. La vitesse de rotation de la pale (nombre de Reynolds) : Une rotation plus rapide signifie généralement plus de turbulence et un mélange plus rapide.
  2. La différence de poids (nombre de Richardson) : Si les liquides ont presque le même poids, ils se mélangent facilement. Si l'un est beaucoup plus lourd, la gravité lutte contre le mélange, créant des couches difficiles à briser.

La Découverte : La Bataille « Gravité contre Rotation »

Les chercheurs ont trouvé des motifs intéressants :

1. Quand la Gravité n'a pas d'Importance (Même Poids) :
   Si les deux liquides ont exactement le même poids, le temps de mélange est étonnamment constant. Peu importe la vitesse à laquelle vous faites tourner la pale (dans une certaine plage), le « temps de mélange sans dimension » (une manière élégante de dire « combien de tours de pale sont nécessaires ») reste constant à environ 20 tours. C'est comme une loi de la nature : une fois que l'eau est suffisamment brassée, la faire tourner plus vite ne la fait pas mélanger plus vite en termes de tours de pale.

2. Quand la Gravité Résiste (Poids Différents) :
   Dès que les liquides ont des poids différents, le lourd veut rester au fond. Plus la différence de poids est grande, plus il est difficile de mélanger.

   • La Tendance : Plus les poids sont différents, plus le temps de mélange est long.
   • La Surprise : Si vous maintenez la « différence de poids » constante et que vous faites simplement tourner la pale plus vite, le temps de mélange ne diminue pas toujours. Parfois, tourner plus vite fait en réalité qu'il faut plus de temps pour atteindre un point de mélange spécifique.
   • Pourquoi ? Imaginez que le liquide lourd est comme une couverture épaisse. Si vous faites tourner la pale trop vite, vous créez beaucoup d'énergie, mais le liquide lourd forme un « couvercle » stable que le liquide léger ne peut pas pénétrer. L'énergie est gaspillée à faire tourbillonner la couche supérieure tandis que la couche inférieure reste scellée. C'est comme essayer de remuer une marmite de soupe où les légumes lourds se sont déposés en un bloc solide au fond ; faire tourner la cuillère plus vite éclabousse simplement le bouillon au-dessus sans briser le bloc de légumes.

La Solution : Une Nouvelle « Courbe Maîtresse »

La plus grande réussite de l'équipe a été de créer une formule unique et simple qui combine tous ces facteurs. Ils ont découvert que si l'on considère le temps de mélange à travers le prisme de trois nombres spécifiques (Puissance, Froude et Richardson), tous leurs points de données désordonnés se regroupent sur une seule courbe exponentielle lisse.

Imaginez-le ainsi : Auparavant, les ingénieurs devaient deviner ou effectuer des centaines de tests pour voir comment un nouveau liquide se mélangerait. Maintenant, ils ont une « recette ». Si vous leur indiquez la différence de poids et la vitesse de rotation prévue, cette formule prédit le temps de mélange avec une grande précision.

La Conclusion

Le papier conclut que pour ces réservoirs industriels spécifiques :

• La turbulence est la clé : Une fois que le liquide est complètement brassé, le comportement de mélange est prévisible.
• La gravité est le patron : Si les liquides ont des densités différentes, la gravité crée une « stratification » (mise en couches) qui résiste au mélange.
• Plus rapide n'est pas toujours mieux : Dans les systèmes avec de fortes différences de densité, augmenter simplement la vitesse du moteur ne garantit pas un mélange plus rapide ; cela crée parfois une séparation plus stable.

Les auteurs fournissent cette nouvelle formule pour aider les ingénieurs à concevoir de meilleurs processus de mélange sans avoir besoin de construire des prototypes coûteux au préalable. Ils prévoient de tester cette formule sur différentes formes de réservoirs et types de pales à l'avenir, mais pour l'instant, elle fonctionne parfaitement pour les réservoirs standards qu'ils ont simulés.

Résumé technique

Résumé technique : Mise à l'échelle sans dimension pour des différences de densité intermédiaires à élevées dans un réservoir agité

Énoncé du problème
Le mélange de liquides miscibles constitue une opération unitaire critique dans les procédés industriels, allant de la production pharmaceutique au traitement des eaux usées. Bien que le mélange turbulent dans les systèmes monophasiques et homogènes soit bien caractérisé, la présence de différences de densité entre les fluides introduit une stratification pilotée par la flottabilité. Ce phénomène crée des régions stratifiées verticalement où les forces gravitationnelles s'opposent au transport convectif, supprimant potentiellement le mélange vertical et augmentant le temps nécessaire pour atteindre l'homogénéité. Bien que des corrélations empiriques existent pour le mélange de traceurs en phase unique et pour des cas stratifiés spécifiques, un cadre généralement valable qui prenne systématiquement en compte les effets gravitationnels, inertiels, visqueux et de flottabilité à travers différents régimes d'écoulement et géométries reste insaisissable. Cette étude aborde le défi de la prédiction des temps de mélange dans des réservoirs agités munis de chicanes contenant deux liquides miscibles présentant des différences de densité intermédiaires à élevées.

Méthodologie
Les auteurs ont utilisé la Dynamique des Fluides Numérique (CFD) avec la méthode de Boltzmann sur réseau (LBM) pour simuler le mélange turbulent dans un réservoir agité muni de chicanes. Le système consistait en un rapport volumique 50/50 de deux liquides miscibles ayant des viscosités dynamiques identiques mais des densités variables ($\rho_l$ et $\rho_h$), où le fluide plus léger était initialement stratifié au-dessus du fluide plus lourd.

• Approche numérique : Les simulations ont résolu les équations de Navier-Stokes incompressibles en utilisant un schéma de réseau D3Q19. Pour gérer les nombres de Reynolds élevés ($Re > 10^5$), un modèle de simulation des grandes échelles (S-L LES) de type Smagorinsky-Lilly a été mis en œuvre pour modéliser la turbulence sous-maille.
• Transport scalaire : Un modèle de transport scalaire basé sur le flux de masse a été adopté pour suivre la concentration du liquide plus léger. Ce modèle a couplé le champ scalaire à l'équation de quantité de mouvement via une force volumique gravitationnelle dépendante de la concentration, permettant la simulation des effets de flottabilité résultant des variations locales de densité.
• Espace de simulation : Une étude paramétrique systématique a été menée en faisant varier la vitesse de rotation de l'agitateur ($N$) et la densité de la phase légère ($\rho_l$). Cela a abouti à une large gamme de nombres sans dimension : des nombres de Reynolds ($Re$) allant de 50 000 à 125 000 et des nombres de Richardson ($Ri$) allant de 0 à 14.
• Validation : La précision numérique a été vérifiée par une analyse de l'indice de convergence de maillage (GCI), qui a estimé l'incertitude numérique à environ 2,39 % sur le maillage de production. De plus, les résultats de simulation ont été validés par rapport à des données expérimentales de traceurs issues de la littérature, montrant une bonne concordance dans la réponse des sondes et les définitions du temps de mélange.

Résultats clés
L'étude a quantifié le temps de mélange ($t_m$), défini comme le temps nécessaire pour atteindre 90 % d'homogénéité, et l'a exprimé sous forme sans dimension ($t^*_m = N \cdot t_m$).

1. Référence homogène ($Ri = 0$) : En l'absence de différences de densité, le temps de mélange sans dimension s'est effondré vers une valeur constante de $t^*_m \approx 20$ sur toute la plage de nombres de Reynolds investiguée. Cela confirme que l'écoulement est pleinement turbulent et que le mélange inertiel domine, ce qui est cohérent avec la littérature établie pour les réservoirs munis de chicanes.
2. Effet de la stratification ($Ri > 0$) : À mesure que le nombre de Richardson augmentait, indiquant de plus grandes différences de densité, le temps de mélange sans dimension augmentait systématiquement. Les données suivaient une tendance exponentielle claire sur tous les nombres de Reynolds.
3. Loi d'échelle proposée : Les auteurs ont dérivé une courbe maîtresse qui projette toutes les données de simulation sur une seule relation d'échelle exponentielle :
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Où $N_p$ est le nombre de puissance, $Fr$ est le nombre de Froude et $Ri$ est le nombre de Richardson. Le préfacteur de 20 représente la référence inertielle, tandis que le terme exponentiel capture l'effet stabilisateur de la flottabilité.
4. Comportement non monotone : En analysant les temps de mélange dimensionnels ($t_m$) à $Ri$ constant, l'étude a observé une dépendance non monotone vis-à-vis de $Re$ pour les nombres de Richardson plus élevés ($Ri \geq 8$). Bien que l'augmentation de $Re$ (et donc de la vitesse de l'agitateur) réduise initialement le temps de mélange en renforçant la turbulence, des augmentations ultérieures conduisent à une stratification de densité plus forte (due au couplage des paramètres dans la configuration de la simulation), ce qui supprime le transport vertical. Cette interaction se traduit par un temps de mélange minimum à des nombres de Reynolds intermédiaires.
5. Structure de l'écoulement : L'analyse du champ de vitesse a révélé qu'à haut $Ri$, la phase légère peut devenir « isolée » de l'apport énergétique de l'agitateur, des structures verticales apparaissant principalement dans le liquide plus dense. L'analyse de la densité spectrale de puissance (DSP) a indiqué que, bien que la flottabilité affecte l'organisation spatiale de l'écoulement, la dynamique pilotée par l'agitateur domine toujours la cascade d'énergie.

Importance et revendications
L'article revendique fournir une corrélation prédictive compacte pour les temps de mélange turbulent dans des systèmes initialement stratifiés par densité, dans le régime pleinement turbulent. La loi d'échelle dérivée, basée sur l'influence combinée des nombres de puissance, de Froude et de Richardson, offre une expression analytique simple pour les ingénieurs et les experts en applications.

Les auteurs notent que, à leur connaissance, aucune corrélation comparable pour des systèmes stratifiés dans le régime pleinement turbulent n'a été rapportée précédemment. L'importance de ce travail réside dans sa capacité à rationaliser la conception des procédés de mélange en réduisant la dépendance à des essais réels coûteux. Le modèle proposé permet d'estimer les temps de mélange dans des systèmes où les effets de flottabilité sont significatifs, capturant la compétition complexe entre la turbulence accrue et la stratification renforcée. Les auteurs suggèrent modestement que les travaux futurs consisteront à tester cette corrélation sur différentes géométries de réservoirs et types d'agitateurs, y compris les agitateurs inclinés, afin d'évaluer la robustesse de la mise à l'échelle proposée.