A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Ce papier présente un nouveau solveur couplé par blocs Volume de Fluide et Frontière Immergée (BC-VOF-IB) implémenté dans OpenFOAM qui surmonte les instabilités numériques causées par les forts contrastes de viscosité pour simuler avec précision le mélange de polymères à surface libre dans des extrudeurs industriels partiellement remplis.

L'essentiel

Imaginez que vous essayez de mélanger un immense bac de miel épais et collant avec de l'air à l'intérieur d'une machine tournante. C'est essentiellement ce qui se produit dans le mélange industriel des polymères, où des entreprises comme Pirelli doivent mélanger du plastique fondu avec des additifs pour fabriquer des pneus, des dispositifs médicaux ou des pièces automobiles. L'objectif est d'obtenir un mélange parfaitement homogène afin que le produit final soit résistant et uniforme.

Cependant, simuler ce processus sur un ordinateur est un cauchemar pour les mathématiciens et les ingénieurs. Voici pourquoi, et comment cet article le résout, en utilisant des analogies simples :

Le Problème : La Lutte « Miel Épais contre Air Léger »

Dans ces machines, vous avez deux fluides très différents :

1. La Fonte de Polymère : Extrêmement épaisse, collante et lente (comme du miel froid).
2. L'Air : Très léger et rapide.

Lorsque vous essayez de simuler comment ces deux éléments interagissent à l'intérieur d'une machine dotée de vis tournantes, les programmes informatiques standards se perdent. C'est comme essayer de calculer le mouvement d'un escargot et d'une voiture de course sur la même piste en utilisant le même ensemble de règles. L'ordinateur tente de faire des pas minuscules, très minuscules, pour empêcher l'« escargot » (le plastique épais) de se déplacer trop vite, ce qui rend la simulation incroyablement lente — prenant parfois des jours pour simuler quelques secondes de mélange en temps réel.

De plus, les machines possèdent des pièces complexes et tournantes (des vis) qui se déplacent à l'intérieur d'un conteneur fixe. Traditionnellement, pour simuler cela, vous devez construire un maillage numérique (une grille de petits cubes) qui enveloppe parfaitement les vis tournantes. Alors que les vis tournent, cette grille doit constamment se remodeler elle-même, ce qui revient à essayer de tricoter un pull pendant que la personne qui le porte court un marathon. C'est désordonné, difficile et sujet aux erreurs.

La Solution : Une Nouvelle « Grille Intelligente » et une « Approche d'Équipe »

Les auteurs de cet article ont développé une nouvelle méthode pour exécuter ces simulations en utilisant un logiciel appelé OpenFOAM. Ils ont combiné deux techniques puissantes :

1. La Méthode des Frontières Immergées (L'Astuce du « Mur Fantôme »)
Au lieu de remodeler la grille pour qu'elle s'adapte aux vis tournantes, ils ont maintenu la grille fixe et rigide (comme un bloc de glace solide). Ils ont ensuite dit à l'ordinateur : « Hé, il y a une vis qui tourne à l'intérieur de ce bloc de glace. »

• L'Analogie : Imaginez une piscine avec une grille fixe de carreaux au fond. Au lieu de déplacer les carreaux pour s'adapter à un nageur, vous dites simplement à l'eau : « Ne passe pas à travers le nageur. » L'ordinateur utilise les mathématiques pour créer un « mur fantôme » autour de la vis, forçant le fluide à s'écouler autour sans jamais avoir besoin de reconstruire la grille. Cela rend la gestion de formes complexes et mobiles beaucoup plus facile.

2. La Méthode du Volume de Fluide (VOF) (L'Astuce de la « Peinture de Suivi »)
Pour voir où se termine le plastique épais et où commence l'air, ils utilisent une « peinture » qui remplit les cellules.

• L'Analogie : Imaginez que la grille informatique est un damier en 3D. Certaines cases sont à 100 % de plastique, d'autres à 100 % d'air, et certaines sont un mélange. L'ordinateur suit la quantité de « peinture plastique » dans chaque case pour visualiser la surface du liquide.

3. Le Schéma Couplé par Blocs (Le « Huddle d'Équipe »)
C'est la percée la plus importante. Dans les simulations standard, l'ordinateur résout la vitesse du fluide dans les directions X, Y et Z l'une après l'autre, comme trois personnes qui parlent à tour de rôle. Lorsque le fluide est super épais (comme le polymère), cette approche de « tour de rôle » provoque un plantage de la simulation ou un ralentissement extrême, car le fluide épais couple toutes les directions ensemble de manière très serrée.

Les auteurs ont changé cela pour une approche Couplée par Blocs.

• L'Analogie : Au lieu que trois personnes parlent à tour de rôle, elles se regroupent toutes et résolvent le problème ensemble exactement au même moment. En traitant le mouvement dans toutes les directions comme une seule équipe géante et interconnectée, l'ordinateur peut gérer la différence massive entre le plastique épais et l'air léger sans se bloquer.

Les Résultats : De Heures à Minutes

L'équipe a testé sa nouvelle méthode sur deux scénarios :

1. Un Canal en Forme de Dog-Bone : Un cas de test où du plastique épais est injecté dans un canal étroit et sinueux.

   • Ancienne Méthode : Le programme informatique standard plantait ou prenait 7 heures pour simuler quelques secondes car il était forcé de faire des pas minuscules.
   • Nouvelle Méthode : Leur nouvelle méthode de « Huddle d'Équipe » a terminé le même travail en seulement 16 minutes et n'a pas planté, même lorsque le plastique devenait extrêmement épais.

2. Machines Industrielles Réelles : Ils ont simulé des extrudeuses à vis unique et à vis doubles réelles (les machines utilisées pour fabriquer des granulés de plastique).

   • Ils ont montré avec succès comment le plastique remplit la machine, comment la pression s'accumule et comment l'air est expulsé.
   • Ils ont prouvé que leur méthode de « Mur Fantôme » fonctionne aussi bien que l'ancienne méthode difficile de remodelage de la grille, mais beaucoup plus rapidement et plus facilement à mettre en place.

Et Après ?

L'article conclut que c'est une étape majeure pour l'industrie. Il comble le fossé entre les mathématiques académiques et les besoins réels des usines. Cependant, les auteurs notent que leur modèle actuel suppose que la température reste constante (isotherme). En réalité, le mélange du plastique génère de la chaleur, ce qui modifie l'épaisseur du plastique. L'ajout des effets de température et de comportements de plastique « extensible » plus complexes constituent les prochaines étapes pour les recherches futures.

En résumé : Ils ont créé un moyen plus rapide et plus stable d'observer sur un ordinateur comment le plastique épais se mélange à l'air dans des machines tournantes, transformant un processus qui prenait autrefois des heures en un processus qui prend des minutes, sans avoir besoin de reconstruire le monde numérique à chaque fois qu'une vis tourne.

Résumé technique

Résumé Technique : Une Méthode de Frontière Immergée Volume de Fluide pour le Mélange Industriel de Polymères

Énoncé du Problème
La simulation numérique des procédés de mélange de polymères présente des défis significatifs dus à l'interaction complexe entre la rhéologie non newtonienne, la dynamique de surfaces libres multiphasiques et les géométries complexes en mouvement. Les mélangeurs industriels, tels que les extrudeuses à vis simple (SSE) et à vis double (TSE), fonctionnent souvent dans des conditions de remplissage partiel où des polymères fondus très visqueux coexistent avec de l'air. Simuler ces scénarios nécessite de capturer l'interface entre les phases tout en gérant des frontières solides en rotation.

Les approches standards de la Dynamique des Fluides Numérique (CFD) rencontrent deux obstacles principaux dans ce contexte :

1. Complexité Géométrique : La génération de maillages adaptés (conformes) pour des vis en rotation à géométrie complexe est coûteuse en calcul et nécessite souvent des techniques de maillage dynamique qui peuvent entraîner une distorsion du maillage.
2. Instabilité Numérique : Le contraste extrême de viscosité entre les polymères fondus et l'air (souvent dépassant $10^5$ à $10^8$) provoque de graves instabilités numériques dans les solveurs Volume de Fluide (VOF) ségrégés standards. Plus précisément, le traitement explicite du gradient de vitesse transposé dans le terme de diffusion visqueuse impose des contraintes sévères sur le pas de temps ($\Delta t \propto h^2$), rendant les simulations d'écoulements à haute viscosité prohibitives en termes de calcul ou instables.

Méthodologie
Ce travail propose un cadre numérique novateur intégrant une approche de capture d'interface Volume de Fluide (VOF) avec une méthode de Frontière Immergée (IB) non conforme, implémentée au sein de la bibliothèque de volumes finis OpenFOAM.

• Approche Frontière Immergée (IB) : La méthode utilise une grille de fond fixe où les frontières solides (par exemple, les vis en rotation) sont représentées par des surfaces STL triangulées. Un opérateur d'interpolation par Moindres Carrés Pondérés (WLS) est employé pour imposer les conditions aux limites (Dirichlet pour la vitesse, Neumann pour la fraction volumique) sur les cellules intersectées par la frontière immergée. Cela élimine le besoin de génération de maillage conforme et permet la gestion de lignes de contact mobiles le long de géométries complexes en rotation.
• Formulation VOF : L'interface est suivie en résolvant une équation de transport pour la fraction volumique de phase ($\alpha$). La méthode emploie l'algorithme MULES (Limiteur Universel Multidimensionnel avec Solution Explicite) pour maintenir la netteté et la bornitude de l'interface, incluant un terme de vitesse de compression pour contrer la diffusion numérique.
• Schéma Couplé par Blocs (BC) : Pour résoudre les problèmes de stabilité découlant des forts contrastes de viscosité, les auteurs introduisent un schéma couplé par blocs pour l'équation de quantité de mouvement. Contrairement aux solveurs ségrégés standards qui traitent les composantes de vitesse séquentiellement, l'approche BC traite le terme complet de diffusion visqueuse $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ de manière implicite. Cela implique de coupler les trois composantes de vitesse au sein d'un seul système linéaire. La discrétisation du terme de gradient tangentiel utilise une méthode de Surface Finie (FA) combinée à une interpolation WLS sur des stencils étendus (impliquant des voisins par point plutôt que seulement des voisins par face), adaptée des applications en élasticité [11].
• Algorithme de Résolution : Le solveur emploie un algorithme PIMPLE modifié (combinant SIMPLE et PISO) pour gérer le couplage pression-vitesse. Les contraintes IB sont intégrées dans le système algébrique, modifiant les coefficients de matrice et les termes sources pour satisfaire les conditions aux limites sur les surfaces immergées tout en maintenant la conservation de la masse.

Contributions Clés

1. Développement de Solveurs BC-VOF-IB : L'article introduit le solveur UCInterIbFoam, qui combine la méthode IB non conforme avec l'approche VOF couplée par blocs. Il s'agit de la première implémentation d'une telle méthode couplée au sein du cadre OpenFOAM-10 pour des écoulements biphasiques.
2. Amélioration de la Stabilité par Couplage par Blocs : Le travail démontre que le traitement entièrement implicite de la diffusion visqueuse relâche considérablement les contraintes de stabilité du pas de temps. Cela permet la simulation d'écoulements de polymères à haute viscosité avec des ratios de viscosité allant jusqu'à $10^8$ sans que la simulation ne plante en raison de violations du nombre de Courant.
3. Gestion des Lignes de Contact Mobiles : Le cadre résout avec succès les difficultés géométriques et algorithmiques des lignes de contact mobiles où l'interface biphasique intersecte les frontières immergées, une exigence critique pour les dispositifs de mélange partiellement remplis.

Résultats Numériques
Les méthodes proposées ont été validées et testées sur deux catégories de problèmes :

• Cas de Référence (Moulage par Injection en forme d'os de chien) :

  • Une comparaison entre les solveurs ségrégés (interFoam), couplés par blocs conformes (UCInterFoam) et couplés par blocs non conformes (UCInterIbFoam) a été effectuée sur un cas d'injection de fluide à loi de puissance à haute viscosité.
  • Le solveur ségrégé standard n'a pas réussi à produire des résultats stables dans des délais raisonnables en raison de restrictions sévères sur le pas de temps ($\Delta t \approx 10^{-6}$ s).
  • Les solveurs couplés par blocs ont atteint des solutions stables avec des pas de temps allant jusqu'à $10^{-3}$ s, réduisant le temps de calcul de 7 heures à 16 minutes (conforme) et 58 minutes (non conforme).
  • La solution IB non conforme a montré un accord qualitatif avec la solution conforme, avec de légères divergences dans les pics de viscosité près des restrictions géométriques, attribuées à la nature non conforme du maillage.

• Applications Industrielles (SSE et TSE) :

  • Des simulations ont été menées sur des extrudeuses à vis simple (SSE) simplifiées et des extrudeuses à vis double (TSE) réalistes dans des conditions de remplissage complet et partiel (alimentation en starvation).
  • Les solveurs ont capturé avec succès l'évolution de la surface libre, les taux de remplissage et les champs de pression.
  • Les résultats étaient conformes aux attentes physiques : les scénarios de remplissage complet ($Q_{in} > Q_0$) présentaient des gradients de pression négatifs de l'entrée à la sortie, tandis que les scénarios de remplissage partiel ($Q_{in} < Q_0$) montraient des augmentations de pression près de la sortie.
  • La méthode IB s'est avérée efficace pour gérer la géométrie complexe d'engrènement de la TSE sans nécessiter de régénération dynamique du maillage.

Signification et Revendications
L'article revendique que le cadre BC-VOF-IB proposé représente une étape significative vers le rapprochement entre la recherche académique en CFD et les exigences du traitement industriel des polymères. En surmontant les instabilités numériques associées aux forts contrastes de viscosité et en éliminant le besoin d'une génération complexe de maillages adaptés pour les géométries en rotation, la méthode permet des prédictions physiquement cohérentes des champs de vitesse et de pression dans des dispositifs de mélange partiellement remplis pertinents pour l'industrie.

Les auteurs reconnaissent que le travail actuel est limité aux procédés isothermes et à la rhéologie newtonienne généralisée. Ils déclarent que les développements futurs doivent inclure les effets thermiques (équation de l'énergie) et des modèles constitutifs viscoélastiques pour capturer pleinement la physique du mélange industriel des polymères. Cependant, le cadre actuel établit une base rigoureuse pour ces extensions et fournit un outil robuste pour simuler des écoulements à surface libre dans des géométries complexes en mouvement.