Comprendre la dynamique des bulles est fondamental pour une vaste gamme de processus industriels.De réacteurs chimiques et des cuves de fermentation aux usines de traitement des eaux usées et aux opérations de récupération du pétrole, la formation, la croissance, le mouvement et la coalescence des bulles de gaz influencent directement l'efficacité, la sécurité et la qualité du produit.La dynamique des fluides computationnels (CFD) est apparue comme un outil indispensable pour modéliser ces phénomènes multiphasés complexes, permettant aux ingénieurs de simuler le comportement des bulles dans des conditions d'exploitation réalistes sans les frais et les limites des expériences physiques.

Les fondamentaux de la dynamique des bulles dans les fluides industriels

La dynamique des bulles englobe l'ensemble du cycle de vie d'une bulle de gaz dans un milieu liquide. Le processus commence par la nucléation, où le gaz dissous sort de solution à un site de nucléation – souvent une imperfection de surface ou une microbulle préexistante. Une fois formée, les bulles se développent à mesure que le gaz se diffuse en elles ou que la pression diminue. La bulle augmente alors en raison de la flottabilité, se déforme sous l'influence de la tension de surface et des contraintes visqueuses, et peut interagir avec d'autres bulles par coalescence ou rupture.

Les paramètres physiques clés qui régissent la dynamique des bulles comprennent la densité et la viscosité des liquides, le coefficient de tension de surface, la solubilité des gaz et les champs de pression et de température locaux. Les nombres sans dimension Eötvös (Eo), Reynolds (Re) et Morton (Mo) sont couramment utilisés pour caractériser les régimes de forme des bulles, de la cape sphérique à ellipsoïdale à la cape sphérique.

Importance industrielle : Dans un bioréacteur aérobie, le transfert d'oxygène des bulles de gaz à la culture liquide est l'étape limitante de la croissance cellulaire. Dans un réacteur à colonne de bulles pour la synthèse Fischer-Tropsch, la distribution de la taille des bulles affecte directement le transfert de masse et la sélectivité de la réaction.

Pourquoi la dynamique des fluides calculateurs pour la modélisation des bulles?

La mesure expérimentale de la dynamique des bulles dans les réacteurs industriels opaques, à haute pression ou à haute température est extrêmement difficile. L'instrumentation, comme les caméras à grande vitesse, les sondes optiques et les capteurs de conductivité, fournit des données ponctuelles limitées et peut perturber le débit. CFD offre une approche complémentaire ou alternative qui fournit des informations sur la vitesse, la distribution de phase et les caractéristiques des bulles en plein champ, résolues dans le temps.

Le CFD résout les équations qui régissent le mouvement du fluide — les équations Navier-Stokes — pour la phase liquide, avec des modèles supplémentaires pour tenir compte de la présence de la phase gazeuse. Selon l'approche de modélisation, le gaz peut être traité comme une interface distincte (bulbe résolue) ou comme une phase dispersée avec des propriétés moyennes. Le choix du modèle dépend de l'échelle d'intérêt : de la formation de bulles individuelles à un orifice (échelle de millimètre) à un réacteur entier (échelle de mètre).

Les logiciels CFD modernes tels que ANSYS Fluent, STAR-CCM+ et OpenFOAM comprennent des modules multiphases dédiés qui peuvent gérer la dynamique de bulles avec des degrés de fidélité variables. La disponibilité croissante de ressources de calcul haute performance (HPC) a permis de simuler des flux industriels à grande échelle avec des millions de cellules de calcul et des interfaces de type bulle.

Approches numériques pour la dynamique des bulles dans le CFD

Plusieurs méthodes de CFD existent pour modéliser les flux gaz-liquides avec des bulles. Chacune a des forces et des limites en ce qui concerne la précision, le coût de calcul et l'applicabilité à différents régimes de débit.

Modèle eulérien-eulérien (deux rainures)

Dans cette approche, le gaz et le liquide sont traités comme des continuums interpénétrateurs, chacun avec sa propre fraction de volume. Le modèle résout deux ensembles d'équations de conservation (masse et moment) couplées par des termes d'échange interphase. Un modèle de bilan démographique (PBM) est souvent ajouté pour suivre la distribution de la taille des bulles, en tenant compte de la nucléation, de la croissance, de la coalescence et de la rupture.

Eulerian-Lagrangian (modèle de bulles de béton, DBM)

Ici, la phase liquide est résolue dans un cadre eulérien, tandis que les bulles individuelles sont suivies dans un cadre lagrangien. Chaque bulle est assignée une position, vitesse, taille et forme (souvent supposé sphérique pour la simplicité). Des modèles pour la traînée, le levage, la masse virtuelle et la dispersion turbulente sont appliqués à chaque bulle. DBM peut fournir des statistiques détaillées sur les trajectoires et les interactions de bulles, mais devient calculativement intense pour les comptes de bulles élevés (généralement limité à <106 bulles).

Méthode du volume de fluide (VOF)

Cette méthode résout explicitement la forme de bulle, la déformation et les événements de coalescence/découpe. VOF est idéal pour étudier la formation de bulles à un orifice, la montée de bulles dans des colonnes stagnantes ou l'interaction avec les obstacles. Cependant, elle nécessite une résolution très fine du maillage près de l'interface (souvent avec un raffinement adaptatif du maillage) et est généralement limitée à un nombre modéré de bulles (de dizaines à centaines) en raison du coût élevé du calcul.

Méthode de définition de niveau

Comme pour VOF, la méthode de réglage de niveau utilise une fonction de distance signée pour suivre l'interface. Elle offre une reconstruction plus fluide de l'interface, ce qui est bénéfique pour calculer avec précision les forces de tension de surface. Associée à une étape de réinitialisation, elle peut gérer des changements topologiques comme la coalescence.

Modèles hybrides

Les développements récents comprennent des approches hybrides qui combinent des méthodes d'interface résolues pour les grandes bulles avec un modèle de phase dispersé pour les petites bulles. Par exemple, le modèle de densité interfaciale algébrique (AIAD) ou le concept de flux biphasé généralisé (GENTOP) peut changer dynamiquement entre les formulations VOF et eulériennes en fonction des conditions de débit locales.

Équations et considérations de modélisation

Quelle que soit l'approche, le débit en phase liquide est décrit par les équations de Navier-Stokes incompressibles ou compressibles, souvent avec modélisation de turbulence. Pour les écoulements en phase liquide, l'équation de moment comprend les termes d'échange de moment interphases – drag, lifting, lubrification de mur, et dispersion turbulente. Des modèles de fermeture précis pour ces termes sont critiques.

La modélisation de la turbulence dans les flux de bulles est particulièrement difficile. La présence de bulles peut soit supprimer ou améliorer la turbulence en phase liquide selon la taille de la bulle et la fraction vide. Les approches courantes comprennent le modèle standard k-ε avec des termes de source supplémentaires pour la turbulence induite par la bulle, ou des modèles plus avancés comme les modèles de stress de Reynolds (RSM) ou la simulation de grande taille Eddy (LES). LES fournit une meilleure résolution de turbulence à grande échelle mais à un coût de calcul beaucoup plus élevé.

La tension de surface joue un rôle dominant dans les petites bulles. Le saut de pression de Laplace à travers l'interface est pris en compte par le modèle de Force de surface continu (CSF) ou le modèle de Force de surface Sharp (SSF).

Les schémas numériques doivent être soigneusement choisis pour maintenir la netteté de l'interface et éviter la diffusion numérique. Les schémas de différenciation compressive (par exemple, HRIC, CICSAM) sont souvent utilisés dans le VOF, tandis que l'intégration de temps de plus grande ordre est nécessaire pour la dynamique transitoire des bulles.

Mesh Generation pour simulations de bulles

Pour les méthodes d'interface résolues (VOF, jeu de niveaux), un minimum de 10-20 cellules par diamètre de bulle est typique, et le raffinement local autour des interfaces de déformation est essentiel. Les techniques de raffinement adaptatif du maillage (AMR) sont largement utilisées pour maintenir la résolution tout en maintenant les coûts computationnels gérables. Pour les modèles eulérien-eulérien et lagrangien, le maillage doit résoudre les échelles moyennes de flux et de turbulence, mais les bulles ne sont pas explicitement capturées; au contraire, les modèles sous-réseaux tiennent compte de la physique à l'échelle bulle.

Dans les géométries industrielles, qui peuvent comprendre des chicanes, des espargers, des turbines et des échangeurs de chaleur, les mailles hexaédriques sont préférées pour leur précision et leur efficacité. Les mailles polyédriques offrent une flexibilité pour les géométries complexes.

Études de cas : Applications industrielles de la DFC de bulle

Réacteurs de la colonne de bulles

Les simulations CFD peuvent prédire le maintien des gaz, la distribution de la taille des bulles, les schémas de circulation des liquides et les coefficients de transfert de masse. Dans une étude utilisant l'approche eulérienne-eulérienne avec un modèle de bilan de population, les chercheurs ont pu optimiser la conception de la sparger pour obtenir une distribution uniforme des bulles, augmentant l'efficacité du transfert d'oxygène de 30 % par rapport à une conception conventionnelle. La simulation a également révélé des zones mortes où la recirculation des bulles a conduit à une réduction du mélange, qui a été atténuée par la modification du rapport d'aspect de la colonne.

Réservoirs d'aération pour traitement des eaux usées

Dans les procédés de boues activées, les diffuseurs de bulles fines fournissent de l'oxygène pour la dégradation microbienne de la matière organique. La modélisation CFD aide à déterminer le placement optimal du diffuseur, le débit d'air et la taille des bulles pour maximiser le transfert d'oxygène tout en minimisant la consommation d'énergie.

Récupération du pétrole et débit multiphasique dans les pipelines

Dans le domaine de l'ingénierie pétrolière, les bulles de gaz peuvent se former en raison de la chute de pression pendant l'extraction (écoulement en deux phases du gaz). La méthode CFD est utilisée pour prédire les régimes de débit (boucle, limace, annulaire) et pour concevoir des séparateurs et des pipelines qui évitent l'accumulation et le lissage de gaz. La méthode VOF a été appliquée pour simuler la formation de bulles à un orifice dans un tuyau horizontal sous haute pression, révélant que la tension de surface et la viscosité ont des effets plus forts à des débits plus faibles.

Mélange et réacteurs chimiques avec des réservoirs stirred

Les simulations CFD utilisant la technique de maille coulissante ou de cadre de référence multiple (MRF) pour la turbine à turbine, combinées à un modèle multiphase eulérien-eulérien, permettent de prédire le nombre de puissance, le maintien en place du gaz et la taille des bulles en fonction de la vitesse de la turbine à turbine à turbine à turbine à Rushton.

Les défis de la modélisation industrielle des bulles

Malgré des progrès importants, plusieurs défis demeurent :

  • Physique à plusieurs échelles: La formation de bulles se produit à l'échelle micrométrique (millimètres), tandis que la performance du réacteur dépend des schémas d'écoulement macroscopiques.
  • Modèle de coalescence et de rupture:[ Les modèles pour ces phénomènes sont souvent empiriques et validés uniquement pour des conditions spécifiques. Prévoir l'évolution de la distribution de la taille des bulles dans les flux turbulents complexes nécessite des noyaux plus robustes.
  • Coûts d'utilisation: Les simulations de haute fidélité (interface résolue, ERP, équilibre de population) sont encore trop coûteuses pour le design industriel courant.
  • Les données de validation:[ Les ensembles de données expérimentales complets dans des conditions industrielles pertinentes (haute pression, fluides non néotoniens, géométries importantes) sont rares et limitent l'étalonnage du modèle.
  • Contamination et surfactants:[ Les fluides industriels contiennent souvent des agents de surface qui suppriment la coalescence et modifient l'hydrodynamique des bulles. La modélisation de ces effets nécessite un couplage avec les équations de transport pour la cinétique de concentration et d'adsorption de l'interface des surfactants.

Orientations futures et innovations

Le champ de la dynamique des bulles CFD évolue rapidement, en raison des progrès des méthodes numériques, de la puissance de calcul et des techniques expérimentales.

Apprentissage automatique et modèles data-driven

Des approches fondées sur les données sont en cours de développement pour remplacer ou augmenter les modèles de fermeture empiriques. Les réseaux neuronaux peuvent prédire les grains de traînée, de levage et de coalescence à partir de données de simulation de haute fidélité.

Calcul à haute performance et accélération du GPU

L'avènement de résolveurs accélérés GPU (p. ex. dans OpenFOAM, ANSYS Fluent et les paquets CFD commerciaux) a permis de réaliser des simulations de milliers de bulles complètement résolues. L'informatique à l'échelle exascale permettra des simulations globales à l'échelle du réacteur qui résolvent les bulles individuelles, éliminant ainsi le besoin de modèles sous-réseaux dans de nombreux cas.

Couplage multiphysique

Les modèles CFD couplés qui résolvent simultanément le transport des espèces, la cinétique des réactions et la population de bulles deviennent plus fréquents. Par exemple, dans un réacteur à colonne à bulles gaz-liquide pour la capture du carbone avec des solvants d'amine, le modèle CFD doit tenir compte de l'augmentation de bulles, de la dissolution, de la chaleur de réaction et de la dégradation des solvants.

Méthodes de délimitation et de coupe immergées

La manipulation de géométries mobiles complexes (par exemple, bulles déformantes, chicanes flexibles) est un défi pour les mailles ajustées par le corps. Les méthodes de bordures immergées permettent des simulations sur les grilles cartésiennes, simplifient le maillage et permettent une manipulation efficace des changements topologiques pendant la coalescence et la rupture des bulles.

Directives pratiques à l'intention des professionnels de l'industrie

Lorsqu'on entreprend une étude sur la dynamique des bulles dans le cadre d'un processus industriel, il faut tenir compte des étapes suivantes :

  1. Définir l'objectif : Êtes-vous intéressé par la distribution de la taille des bulles, le transfert de masse, le temps de mélange ou le régime de débit ?
  2. Sélectionnez le modèle CFD approprié en fonction de la taille de bulle prévue, de la fraction de vide et du budget de calcul. Commencez par un modèle plus simple (p. ex., eulérien-eulérien avec PBM) et raffinez au besoin.
  3. Identifier avec précision les propriétés physiques : densité, viscosité, tension de surface (y compris la température et la concentration), et solubilité dans les gaz.
  4. Valider en fonction des données expérimentales disponibles provenant de la littérature ou des essais contrôlés en laboratoire selon une géométrie représentative.
  5. Effectuer une étude de sensibilité à la grille et évaluer la dépendance en temps et en étape, en particulier pour les simulations transitoires avec suivi VOF ou Lagrangian.
  6. Utiliser des modèles de turbulence appropriés : pour les fractions à vide élevé, envisager des modifications de turbulence induites par bulles.
  7. Pour les modèles d'équilibre de population, commencer par un faible nombre de bacs de taille (p. ex., 10-15) et augmenter progressivement; assurer la stabilité du schéma de discrétisation (p. ex., méthode des classes, méthode des moments quadrature).
  8. Tirer parti du calcul parallèle et envisager le raffinement adaptatif des mailles pour les méthodes d'interface résolues afin de réduire les temps d'exécution.

Conclusion

En prenant en compte la physique complexe de la formation, de l'augmentation et de l'interaction des bulles, le CFD permet aux ingénieurs de concevoir des équipements plus efficaces, de réduire les coûts opérationnels et d'améliorer la performance environnementale. Bien que les défis liés à la modélisation à plusieurs échelles, aux dépenses de calcul et à la validation persistent, les développements continus dans les méthodes numériques, l'apprentissage des machines et l'informatique à haute performance dépassent constamment ces obstacles.