Self-similar breakup of a liquid ligament with a solid particle

Cette étude démontre, par le biais de simulations numériques et de modélisation analytique, qu'une particule solide induit une dynamique de rupture universelle et auto-similaire dans les ligaments liquides en étirement, où la rupture subséquente devient indépendante de la taille de la particule et est régie par l'interaction entre l'étirement du ligament et l'instabilité de Rayleigh-Plateau.

L'essentiel

Imaginez que vous avez un long et fin filament de miel ou de sirop épais suspendu dans les airs. Si vous écartez les extrémités, le filament s'amincit de plus en plus jusqu'à ce qu'enfin, il se rompe et se brise en gouttes distinctes. C'est un spectacle courant dans la nature et la technologie, allant de la pluie tombant d'une feuille aux imprimantes à jet d'encre éjectant de minuscules gouttelettes.

Habituellement, cette rupture se produit parce que le liquide est naturellement instable ; il tend à se transformer en sphères (gouttes) pour économiser de l'énergie. Mais que se passe-t-il s'il y a un tout petit point solide — comme un grain de sable ou une particule de poussière — piégé à l'intérieur de ce filament collant ?

Cet article examine exactement ce scénario. Les chercheurs ont utilisé des simulations informatiques et des mathématiques pour voir comment une seule particule solide modifie la façon dont un filament liquide en cours d'étirement se brise.

Voici l'histoire de leurs découvertes, décomposée en concepts simples :

La Configuration : Un Fil Élastique avec un Nœud

Imaginez le filament liquide comme une longue corde élastique faite de miel. Les chercheurs ont écarté les extrémités de cette corde à une vitesse constante. À l'intérieur de la corde, ils ont placé une seule bille solide (la particule).

Au début, la corde est épaisse, et la bille est petite par rapport à la largeur de la corde. C'est comme avoir un marbre à l'intérieur d'un gros tuyau d'arrosage. Le marbre ne fait pas grand-chose ; la corde s'amincit simplement de plus en plus à mesure qu'elle s'étire, suivant un schéma prévisible.

Le Point de Bascule : Quand le « Tuyau » Rétrécit à la Taille du « Marbre »

À mesure que la corde continue de s'étirer, elle devient plus étroite. Finalement, la corde devient si fine qu'elle est presque en contact avec la surface du marbre à l'intérieur.

C'est le moment critique. L'article appelle cela le moment où le rapport entre la taille de la particule et la taille de la corde s'approche de 1. Soudain, le marbre agit comme un « nœud » ou une « bosse » dans la corde. Parce que la corde est si fine, cette bosse crée une perturbation localisée.

La Surprise : Le « Cassage Universel »

Voici la partie la plus intéressante de la découverte. Les chercheurs ont testé cela avec des marbres de différentes tailles (certains petits, d'autres grands).

• Avant le cassage : Les plus gros marbres faisaient rompre la corde plus tôt que les plus petits. Cela a du sens ; un plus grand obstacle cause des problèmes plus tôt.
• Pendant le cassage : Une fois que la corde est devenue assez fine pour toucher le marbre, quelque chose de magique s'est produit. La vitesse à laquelle la rupture finale s'est produite est devenue exactement la même, que le marbre soit petit ou grand.

Les chercheurs appellent ce comportement « auto-similaire ». C'est comme si, une fois que la corde devient assez fine pour toucher l'obstacle, la taille spécifique de l'obstacle cesse d'importer. Le liquide « oublie » la taille de la particule et suit un chemin universel et prévisible vers la rupture.

L'Analogie : L'Embouteillage

Imaginez une autoroute (le filament liquide) où les voitures s'éloignent les unes des autres, faisant s'étaler la circulation (l'étirement).

• Phase précoce : S'il y a un petit nid-de-poule (petite particule) ou un gros rocher (grosse particule) au milieu de la route, cela ne compte pas encore beaucoup car la route est large.
• Phase tardive : À mesure que la route se rétrécit jusqu'à une seule voie, tant le nid-de-poule que le rocher deviennent d'énormes obstacles.
• La Rupture : Au moment où la circulation est tellement serrée qu'elle heurte l'obstacle, la façon dont la circulation se bloque et s'arrête (la « rupture ») se produit exactement de la même manière pour le nid-de-poule et le rocher. La taille de l'obstacle ne modifie plus le timing de l'embouteillage final ; seul le fait qu'quelque chose soit présent compte.

Les Mathématiques et la Physique

Les chercheurs n'ont pas seulement observé cela se produire ; ils ont écrit une formule mathématique pour prédire exactement quand la rupture se produirait.

• Ils ont découvert que le temps de rupture dépend d'une bataille entre deux forces : l'étirement (écartant la corde) et la viscosité (à quel point le liquide est « épais » ou collant).
• Dans les liquides épais et collants (comme le miel de notre analogie), la « collantité » domine.
• Leur formule a prédit avec succès le temps de rupture, correspondant parfaitement à leurs simulations informatiques.

La Conclusion

L'article conclut que, bien qu'une particule modifie quand la rupture commence (en amincissant la corde plus rapidement près de la particule), une fois que la corde devient assez fine pour toucher la particule, l'acte final de rupture suit une règle universelle.

Dans ce régime spécifique « collant », le filament liquide se comporte comme une machine ayant un programme défini pour se rompre. Une fois que la particule est suffisamment proche pour déclencher le programme, la taille de la particule devient sans importance, et la corde se rompt de manière prévisible et auto-similaire à chaque fois.

Résumé technique

Résumé technique : Rupture auto-similaire d'un ligament liquide contenant une particule solide

Énoncé du problème
La rupture de ligaments liquides qui s'amincissent est un processus fondamental dans des applications allant de la formation de sprays et de l'impression jet d'encre à la génération de gouttelettes en microfluidique. En l'absence d'étirement externe, cette rupture est régie par l'instabilité de Rayleigh–Plateau, où les perturbations de surface croissent sous l'effet des forces capillaires. La dynamique est généralement caractérisée par le nombre d'Ohnesorge ($Oh$), qui représente le rapport entre les forces visqueuses et les forces d'inertie et de tension de surface. Bien que le comportement des ligaments propres soit bien compris, de nombreux systèmes pratiques impliquent des suspensions contenant des particules solides ou des impuretés. Des études antérieures ont examiné des jets chargés de particules avec de nombreuses particules ou des suspensions denses, démontrant des motifs de rupture altérés. Cependant, le mécanisme fondamental par lequel une seule particule solide agit comme une perturbation localisée pour modifier la dynamique de rupture et le temps de détachement d'un ligament en étirement reste mal compris. Plus précisément, il est unclear comment l'interplay entre l'étirement du ligament, la viscosité et une perturbation localisée induite par la particule dicte le régime final de rupture.

Méthodologie
Les auteurs emploient une combinaison de simulations numériques et de modélisation analytique pour investiguer ce problème.

• Configuration numérique : L'étude utilise la méthode de Boltzmann sur réseau (LBM) basée sur une approche de gradient de couleur pour des fluides biphasiques non miscibles. Le système modélise un ligament liquide cylindrique de rayon initial $R_i$ et de longueur $L$ (rapport d'aspect $L/R_i = 28.8$) entouré d'air. Une particule solide sphérique de rayon $r_p$ est placée au centre du ligament avec une condition aux limites d'adhérence (no-slip).
• Mécanisme d'étirement : Un étirement axial est imposé en appliquant des profils de vitesse linéaires aux extrémités du ligament ($u_x = \frac{2}{L}(x - L/2)u_{max}$), provoquant la sortie de masse du domaine et la diminution du rayon du ligament au cours du temps. Le taux d'étirement est caractérisé par un taux de sortie sans dimension $\tilde{\dot{\epsilon}}$.
• Paramètres : Les simulations se concentrent sur le régime visqueux où $Oh \sim 1.0$. La taille de la particule est variée par rapport au rayon initial du ligament ($r_p/R_i$), et le taux de sortie est ajusté pour contrôler l'échelle de temps d'étirement par rapport à l'échelle de temps capillaire-visqueuse ($\tau_v = \mu R_i / \sigma$).
• Approche analytique : Les auteurs dérivent une expression analytique pour le temps de détachement en modélisant le rayon du ligament comme la somme d'un rayon moyen non perturbé (gouverné par la conservation de la masse) et d'une amplitude de perturbation localisée induite par la particule. Ils supposent que la perturbation croît selon la dynamique de l'instabilité de Rayleigh–Plateau, dominée par les effets visqueux dans le régime $Oh \sim 1$.

Résultats clés

1. Dynamique en deux étapes : Le processus de rupture présente deux étapes distinctes. Initialement, lorsque le rayon du ligament $R(t)$ est significativement plus grand que le rayon de la particule ($r_p/R(t) < 1$), le ligament s'amincit de manière monotone selon la conservation de la masse, et l'influence de la particule est négligeable.
2. Début de l'auto-similarité : À mesure que l'étirement progresse et que la surface du ligament approche la particule ($\beta(t) = r_p/R(t) \sim 1$), une perturbation localisée est déclenchée. Une fois cette proximité atteinte, la croissance subséquente de la perturbation et la dynamique finale de détachement deviennent auto-similaires et indépendantes de la taille de la particule.
3. Temps de détachement universel : Bien que le temps total jusqu'à la rupture ($t_p$) dépende de la taille de la particule (les particules plus grandes déclenchent une rupture plus précoce), le temps de détachement auto-similaire ($t_{ps} = t_p - t_{exp}$), défini comme la durée depuis le début de la croissance de la perturbation jusqu'à la rupture, se réduit à une seule valeur pour différentes tailles de particules à un taux d'étirement donné.
4. Validation analytique : Les auteurs dérivent une expression analytique pour le temps de détachement ($t_{p,a}$) basée sur la compétition entre le taux d'étirement et le taux de croissance de l'instabilité capillaire visqueuse ($\lambda = \sigma / \mu R(t)$). Cette expression, $t_{p,a} = \frac{1}{\lambda_0 e^{\dot{\epsilon} t_{p,a}} + \dot{\epsilon} \ln(R_i/a_0)}$, montre un excellent accord quantitatif avec les données de simulation numérique pour diverses tailles de particules et taux de sortie.

Signification et affirmations
L'article prétend révéler un mécanisme universel par lequel les perturbations localisées contrôlent la rupture de ligaments contenant des particules solides dans le régime visqueux. La contribution principale est la démonstration que, une fois que la surface du ligament approche une particule solide, la taille spécifique de la particule devient sans importance pour la dynamique finale de rupture ; le système entre dans un régime auto-similaire universel.

Les auteurs soulignent que cette universalité découle du fait que l'échelle de temps visqueuse ($t_\mu = \mu R / \sigma$) permet aux perturbations induites par la particule de persister et de croître, contrairement au régime inertiel ($Oh \ll 1$) où la rupture se produit trop rapidement pour que de telles perturbations modifient significativement le résultat. Le cadre analytique proposé capture avec succès ce régime d'échelle universel pour $Oh \sim 1$, fournissant un outil prédictif pour le temps de détachement qui dépend du taux d'étirement et des propriétés du fluide mais est indépendant de la taille spécifique de la particule une fois la perturbation déclenchée. L'étude est limitée aux configurations de particules symétriques et stationnaires, servant de configuration minimale pour isoler ces effets, les interactions de particules plus complexes et advectées étant laissées pour des travaux futurs.