🔬 condensed matter

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

Cet article présente un modèle de simulation à base de particules qui prédit avec succès la viscosité et la microstructure de suspensions denses de sphérocylindres granulaires sous écoulement de cisaillement, révélant comment ces propriétés évoluent en fonction de la fraction volumique et du rapport d'aspect tout en corroborant des données expérimentales limitées.

Auteurs originaux : Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

Publié 2026-07-31

📖 5 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez un monde où ce que vous versez, pompez ou mélangez n'est pas seulement de l'eau ou du sable, mais une soupe chaotique de minuscules bâtonnets rigides flottant dans un liquide. C'est le domaine de la rhéologie, la science de la façon dont les choses s'écoulent. Vous avez probablement déjà vu cela en action : pensez à une rivière obstruée par des arbres tombés formant un énorme « embâcle », ou à la bouillie épaisse et collante utilisée pour fabriquer du papier. Dans ces situations, les particules sont si encombrées qu'elles s'entrechoquent constamment, et leur forme importe tout autant que leur nombre. Si vous essayez de remuer un seau de billes, elles roulent les unes à côté des autres facilement. Mais si vous essayez de remuer un seau de longs bâtons fins, ils s'emmêlent, se verrouillent, et soudain, tout le mélange se comporte comme un solide. Les scientifiques comprennent depuis longtemps comment les boules rondes se comportent dans ces foules, mais les règles pour les particules longues en forme de tiges ont été un peu mystérieuses, surtout lorsqu'elles sont très serrées. Comprendre cela est crucial car ces mélanges de « bâtonnets collants » apparaissent partout, de la fabrication industrielle à l'écoulement de la roche fondue à l'intérieur des volcans.

Dans cet article, une équipe de chercheurs a construit une simulation informatique super intelligente pour servir de laboratoire virtuel pour ces suspensions de bâtonnets encombrées. Au lieu de mélanger de vrais produits chimiques dans un laboratoire, ils ont créé un monde numérique rempli de milliers de « sphérocylindres » — des particules qui ressemblent à des capsules ou des pilules (un cylindre avec une hémisphère à chaque extrémité). Ils ont fait flotter ces capsules dans un fluide virtuel, puis ont commencé à les « cisaillement », ce qui est une façon élégante de dire qu'ils ont entraîné la couche supérieure du fluide latéralement tout en maintenant la couche inférieure immobile, tout comme si l'on étale du beurre sur une tartine. Le but était de voir comment l'épaisseur (viscosité) du mélange changeait à mesure qu'ils emballaient davantage de bâtonnets dans l'espace et à mesure qu'ils rendaient les bâtonnets plus longs.

La simulation a révélé des comportements fascinants. D'abord, lorsque le cisaillement a commencé, le mélange est devenu incroyablement épais pendant une fraction de seconde — un « pic de viscosité » — avant de se stabiliser. Cela s'est produit parce que les bâtonnets étaient initialement éparpillés dans des directions aléatoires, se cognant les uns aux autres comme dans un mosh pit. Une fois le flux lancé, les bâtonnets ont commencé à s'aligner, comme des soldats marchant en formation, et le mélange est devenu plus facile à remuer. Les chercheurs ont découvert que plus ils emballaient de bâtonnets (la fraction volumique, $\phi$ ), plus le mélange devenait épais, atteignant finalement un point où il se bloquait (jamming) et cessait de couler complètement. Curieusement, plus les bâtonnets étaient longs (le rapport d'aspect, $A$ , allant de 1 pour une sphère jusqu'à 20 pour de très longs bâtons), plus le point de blocage arrivait tôt. Par exemple, avec des bâtonnets de longueur 5, le point de blocage se situait autour d'une fraction volumique de 0,57, mais pour des bâtonnets de longueur 20, le mélange se bloquait à une densité beaucoup plus faible d'environ 0,41.

L'équipe a également découvert que les bâtonnets ne se contentaient pas de s'aligner ; ils s'alignaient spécifiquement avec le flux. Ils ont mesuré un « paramètre d'ordre » ( $S$ ) pour voir à quel point l'alignement était efficace. À faible densité, les bâtonnets oscillaient et tournaient en boucles (connues sous le nom d'orbites de Jeffery), mais à mesure que la foule devenait plus dense, ils se verrouillaient en une ligne droite le long de la direction du flux. Cependant, si la foule devenait trop dense, les bâtonnets commençaient en fait à être frustrés et ne pouvaient plus rester parfaitement alignés, ce qui faisait chuter légèrement l'ordre. Les chercheurs ont également décomposé les forces en jeu, montrant que lorsque le mélange était fluide, la propre résistance du liquide importait le plus, mais qu'à mesure qu'il devenait encombré, les chocs physiques et la friction entre les bâtonnets devenaient la force dominante, prenant le relais pour rendre le mélange épais.

En créant un modèle qui équilibre la physique du contact, de la friction et de la traînée fluide, les auteurs ont fourni un nouvel outil pour prédire comment ces mélanges complexes se comportent. Ils n'ont pas seulement deviné ; ils ont exécuté des milliers de simulations avec différents nombres de particules, différentes longueurs de bâtonnets et différentes densités de compactage. Leurs résultats suggèrent que bien que l'idée de base selon laquelle « plus de matière égale un fluide plus épais » soit vraie, la forme des particules change radicalement le point auquel le mélange se verrouille. Ce travail ne prétend pas résoudre tous les mystères de la mécanique des fluides, mais il offre une base simulée solide pour comprendre comment gérer ces fluides complexes remplis de bâtonnets dans le monde réel, de la prévention des embâcles dans les rivières à la conception de meilleures bouillies industrielles.

Résumé technique : Rhéologie de suspensions denses de sphérocylindres granulaires par simulation à base de particules

Énoncé du problème
Les suspensions denses de particules granulaires en forme de bâtonnets sont omniprésentes dans les contextes naturels et industriels, allant des embâcles dans les rivières aux suspensions cristallines dans le magma et aux réseaux de fibres dans la fabrication du papier. Bien que la rhéologie des suspensions sphériques soit bien caractérisée, particulièrement à travers l'unification des suspensions non-browniennes et des écoulements granulaires secs par Boyer et al., le comportement des particules allongées reste incomplètement compris. Le degré de liberté supplémentaire associé à l'orientation des particules introduit des phénomènes microstructuraux complexes, tels que les orbites de Jeffery et l'alignement induit par le cisaillement, qui modifient considérablement la viscosité effective et la réponse de contrainte. Les données expérimentales existantes pour les régimes denses (hautes fractions volumiques de solides, $\phi$ ) sont limitées en raison des difficultés de préparation et de mesure. De plus, les approches de simulation précédentes font face à des compromis : la dynamique de Stokes (SD) est coûteuse en calcul pour les systèmes denses en raison des inversions de matrices requises pour l'hydrodynamique à n-corps, tandis que les modèles de la méthode des éléments discrets (DEM) simplifient souvent la géométrie des particules (par exemple, des chaînes de sphères) ou négligent les forces hydrodynamiques nécessaires. Il existe un besoin pour un modèle à base de particules calculablement traitable qui capture fidèlement la physique microscopique dominante des suspensions de sphérocylindres non-browniennes denses sous cisaillement.

Méthodologie
Les auteurs ont développé un modèle de simulation de particules pour des suspensions denses de sphérocylindres non-pliables et non-browniens (cylindres avec des extrémités hémisphériques) sous un écoulement de cisaillement simple. Le modèle opère dans un régime où l'inertie des particules est présente mais négligeable.

Équations du mouvement : Les trajectoires des particules sont mises à jour à l'aide de l'algorithme Velocity-Verlet, intégrant les équations de mouvement de translation et de rotation. Ces équations tiennent compte des forces de contact par paire, des forces de lubrification par paire et des forces hydrodynamiques à corps unique (traînée, portance et couple).
Hydrodynamique : Les forces à corps unique sont modélées par des solutions analytiques pour les sphéroïdes prolatés en tant qu'approximation pour les sphérocylindres. Cela inclut la traînée dépendante de l'orientation et la portance induite par le cisaillement. Le couple hydrodynamique reproduit les orbites de Jeffery.
Contact et Lubrification :
- Contact : Modélisé comme étant de type Hooké avec friction de Coulomb, incluant des composantes ressort-amortisseur normales et tangentielles.
- Lubrification : Un nouveau schéma a été développé pour calculer les forces de lubrification entre les sphérocylindres. Pour traiter la discontinuité de courbure à la jonction extrémité-arbre, le modèle calcule simultanément trois composantes d'interaction (arbre-arbre, extrémité-arbre et extrémité-extrémité). Celles-ci sont pondérées en fonction de la position relative des particules pour assurer des transitions fluides entre les types d'interactions, résolvant ainsi les discontinuités présentes dans les modèles précédents (par exemple, Yamane et al., Butler et Shaqfeh).
Pas de temps dynamique : Une innovation méthodologique clé est le calcul d'un pas de temps dynamique. Le pas de temps est déterminé par la plus petite échelle de temps caractéristique parmi les interactions de contact, les forces de lubrification et les couples de lubrification. Cela garantit la stabilité numérique tout en évitant des pas de temps inutilement petits, améliorant ainsi l'efficacité de calcul.
Configuration de la simulation : Le système utilise des conditions aux limites périodiques dans les directions d'écoulement et de gradient, et des conditions de Lees-Edwards pour le cisaillement. Un mélange bidisperse de rayons (1:1,4) est utilisé pour empêcher la cristallisation. Les simulations couvrent des rapports d'aspect ( $A$ ) allant de 1 à 20 et des fractions volumiques ( $\phi$ ) allant de dilué à proche du blocage (jamming).

Résultats clés

Comportement de la viscosité : Le modèle prédit une augmentation systématique de la viscosité à l'état stationnaire ( $\eta_s$ $η_{s}$ ) avec la fraction volumique ( $\phi$ $ϕ$ ) et le rapport d'aspect ( $A$ $A$ ).
- Réponse transitoire : Lors du démarrage du cisaillement, un pic de viscosité est observé, attribué à la distribution initiale aléatoire des orientations des particules. Ce pic diminue à mesure que les particules s'alignent, menant à un état stationnaire. L'amplitude du pic par rapport à l'état stationnaire augmente avec $\phi$ .
- État stationnaire : La viscosité à l'état stationnaire suit une divergence de loi de puissance $\eta_s \sim (1 - \phi/\phi_m)^{-\varsigma}$ . La fraction volumique de blocage $\phi_m$ diminue à mesure que $A$ augmente (passant d'environ $0,57$ pour $A=5$ à des valeurs plus basses pour des $A$ plus élevés). Le modèle confirme l'indépendance du taux de la viscosité à l'état stationnaire, ce qui est cohérent avec le cadre des suspensions granulaires.
- Contributions de force : À faible $\phi$ , la viscosité est dominée par la contrainte fluide et les forces de lubrification. À mesure que $\phi$ augmente, les forces de contact (normales et tangentielles) deviennent dominantes, dépassant finalement les contributions de lubrification autour de $\phi \approx 0,5$ .
Microstructure et Alignement :
- Paramètre d'ordre : Le paramètre d'ordre orientationnel ( $S$ ) augmente avec $\phi$ jusqu'à un point critique, après quoi il diminue légèrement en raison de la frustration géométrique empêchant un alignement parfait.
- Direction d'alignement : Les particules s'alignent préférentiellement selon la direction de l'écoulement (axe $x$ ), avec un petit angle d'inclinaison ( $\phi_v$ ) par rapport au plan écoulement-gradient. Cet angle d'inclinaison augmente avec $\phi$ jusqu'à un seuil, puis diminue.
- Effet du rapport d'aspect : Des rapports d'aspect plus élevés entraînent un alignement plus fort à de plus faibles fractions volumiques, mais entraînent également un point de blocage plus bas. Pour des $A$ très élevés, l'angle d'inclinaison devient négligeable.
- Transitoire vs État stationnaire : Le pic de viscosité initial est corrélé au temps nécessaire pour que le système passe d'un état aléatoire à un état aligné par le cisaillement.

Signification et Revendications
L'article affirme fournir un outil fondamental pour explorer la mécanique des fluides, la microstructure et la rhéologie des suspensions denses de bâtonnets granulaires. En modélisant avec succès la transition des régimes dilués à denses et en capturant l'interaction entre les forces hydrodynamiques, les contacts frictionnels et la géométrie des particules, le modèle corrobore les données de rhéologie expérimentale limitées disponibles pour les bâtonnets granulaires.

Les auteurs déclarent que leur travail étend la compréhension de la rhéologie des suspensions granulaires (précédemment établie pour les sphères) à la large classe des particules en forme de bâtonnets. Le modèle offre une plateforme pour des études futures visant à cartographier la dépendance du point de blocage vis-à-vis des détails microstructuraux, étudier les contraintes normales et explorer le comportement dans des conditions inhomogènes. L'article présente sa contribution de manière modeste comme une étape vers une description complète de la mécanique des fluides de ces matériaux, laissant la quantification complète de la relation entre le point de blocage, l'alignement et la friction à des travaux futurs. Il ne prétend pas résoudre directement des problèmes de traitement industriel spécifiques, mais fournit plutôt l'infrastructure de simulation pour étudier la physique sous-jacente de ces systèmes.

Articles similaires