Percolation of a cohesive fine particle in a static bed

Résumé technique : Percolation d'une particule fine cohésive dans un lit statique

Énoncé du problème
La percolation de particules fines à travers un lit statique de particules plus grosses est un processus fondamental dans la manipulation industrielle (par exemple, produits pharmaceutiques, traitement chimique) et les écoulements géophysiques (par exemple, transport sédimentaire, contamination des eaux souterraines). Dans les systèmes non cohésifs, le fait qu'une particule fine traverse ou reste piégée est principalement régi par des contraintes géométriques, spécifiquement le rapport des diamètres des particules grosses sur fines ($R = d_l/d_f$). Un seuil de piégeage géométrique existe à $R_t \approx 6,464$ ; lorsque $R > R_t$, il est théoriquement attendu que les fines traversent librement les gorges de pores formées par les grosses particules. Cependant, les systèmes réels impliquent souvent des forces cohésives (forces de van der Waals, électrostatiques, ponts liquides) qui modifient la dynamique de contact. Bien que des études antérieures aient abordé la percolation non cohésive et la rhéologie des écoulements cohésifs, il existe une lacune dans la compréhension de l'influence de la cohésion sur la percolation dans le régime de « criblage libre » où les contraintes géométriques seules permettraient le passage. Cette étude examine la percolation de particules fines cohésives à travers un lit statique de grosses particules où $R=7$ (au-dessus du seuil géométrique), en se concentrant sur la limite de la particule unique où les fines interagissent avec les particules du lit mais non entre elles.

Méthodologie
Les auteurs emploient des simulations par la méthode des éléments discrets (DEM) utilisant le logiciel open-source LAMMPS. Les simulations utilisent un schéma amélioré de Velocity-Verlet adapté aux grands rapports de taille.

• Configuration du système : Un lit statique de grosses particules monodisperses ($d_l = 4$ mm) est compacté avec une fraction $\phi = 0,6$. Des particules fines ($d_f \approx 0,57$ mm, $R=7$) sont libérées au repos au-dessus du lit et laissées à tomber sous l'effet de la gravité.
• Modèles de contact : Les forces de contact normales suivent le modèle élastique de Hertz. L'inélasticité est modélisée via un coefficient de restitution constant ($e_n$). Les forces tangentielles utilisent le modèle de Mindlin avec frottement de glissement ($\mu$) et résistance au roulement.
• Cohésion : Les interactions cohésives sont modélisées principalement à l'aide du modèle Derjaguin-Muller-Toporov (DMT), qui suppose une force attractive constante pendant le contact et une répulsion de Hertz. L'intensité de la cohésion est quantifiée par le nombre de Bond ($Bo = F_c / W$), défini comme le rapport de la force de décollement au poids de la particule fine. Un nombre limité de simulations utilisant le modèle Johnson-Kendall-Roberts (JKR) a également été réalisé pour comparaison.
• Paramètres : L'étude fait varier systématiquement le nombre de Bond ($Bo$), le coefficient de frottement de glissement ($\mu$), le coefficient de restitution ($e_n$) et le module de Young ($E$). La dynamique de $10^4$ particules fines indépendantes est suivie pour assurer une précision statistique.

Résultats clés

1. Piégeage non géométrique : Même lorsque $R=7$ (permettant géométriquement un passage libre), une cohésion et un frottement suffisamment élevés conduisent au piégeage des particules. À mesure que $Bo$ augmente, la distance caractéristique de percolation diminue, ce qui signifie que les particules sont piégées à des profondeurs plus faibles.
2. Distribution du piégeage : La distribution verticale des particules piégées suit une décroissance exponentielle. La distance caractéristique de piégeage ($\lambda$) sert de métrique pour le comportement de percolation. L'inverse de la distance de piégeage ($d_l/\lambda$) augmente de manière monotone avec $Bo$.
3. Rôle du frottement et de la cohésion :
   • À faible $Bo$, le frottement de glissement joue un rôle significatif dans le piégeage.
   • À fort $Bo$, la cohésion devient le mécanisme de piégeage dominant.
   • Un « coude » est observé dans la relation entre l'inverse de la distance de piégeage et le frottement. Ce coude se déplace avec $Bo$ et correspond à une transition entre les régimes où le frottement domine et où la cohésion domine.
   • La rigidité des particules ($E$) influence le piégeage : des particules plus molles (faible $E$) sont plus facilement piégées, en particulier à $Bo$ élevé.
4. Mécanisme collisionnel : Le processus de piégeage est décomposé en deux étapes séquentielles :
   • Condition d'adhésion : Lors d'une collision, une particule soit rebondit, soit adhère. L'adhésion se produit si la vitesse d'impact est inférieure à une vitesse critique ($v_{cr}$), qui dépend de $Bo$, $e_n$ et $E$.
   • Mouvement post-contact : Si la particule adhère, elle peut glisser et se détacher ou rester piégée. Une condition d'équilibre statique basée sur la balance des forces détermine si une particule ne rebondissant pas reste stationnaire. Cette condition dépend d'un coefficient de frottement redimensionné $\tilde{\mu} = \mu\sqrt{4Bo^2 - 1}$. Si $\tilde{\mu} > 1$, la particule a de fortes chances de rester stationnaire ; si $\tilde{\mu} < 1$, elle a de fortes chances de glisser et potentiellement de se détacher.
5. Modèle prédictif : Les auteurs développent un modèle collisionnel reliant la dynamique microscopique à la distance macroscopique de piégeage. L'inverse de la distance de piégeage est prédit par $1/\lambda = -k \ln(1 - P_t)$, où $k$ est la fréquence de collision et $P_t$ la probabilité de piégeage par collision (déterminée par le rapport de la vitesse d'impact à $v_{cr}$). Ce modèle prédit avec précision les distances de piégeage dans le régime dominé par les collisions (fort $Bo$).
6. Double contact : Dans les régimes à faible frottement ou cohésion, les particules glissent souvent et forment des doubles contacts (touchant simultanément deux grosses particules). Cette configuration peut conduire au piégeage même lorsque les conditions de contact unique suggéreraient un détachement, ou inversement, permettre un détachement via un rebond sur un second contact.

Signification et revendications
L'article établit que la percolation de particules fines cohésives n'est pas déterminée uniquement par l'accessibilité géométrique. Au contraire, la dynamique des interactions à l'échelle de la particule — spécifiquement l'interdépendance de la cohésion, de la restitution, du frottement et de la vitesse d'impact — peut prévaloir sur les attentes géométriques, conduisant au piégeage même lorsque les gorges de pores sont géométriquement assez larges pour le passage.

Les auteurs revendiquent fournir un cadre fondamental pour comprendre ces mécanismes dans les lits statiques, ce qui est une condition préalable pour aborder le problème plus complexe de la ségrégation granulaire dans les mélanges en écoulement de particules cohésives. En identifiant les conditions spécifiques (vitesse critique, critères d'équilibre statique) qui régissent l'adhésion et le mouvement subséquent, l'étude offre une capacité prédictive pour la distance de piégeage basée sur les propriétés matérielles fondamentales. Les résultats soulignent que dans les systèmes cohésifs, le « sort » d'une particule fine est régi par une séquence d'interactions locales (collision, adhésion, mouvement post-contact) plutôt que par la simple géométrie statique du lit.