Stability analysis of the flow in a coflowing device

Cette étude démontre que dans les dispositifs de co-écoulement, l'instabilité du jet précède la déstabilisation du ménisque et que la rupture transitoire est hautement sensible aux perturbations initiales, remettant ainsi en question la validité de l'analyse de stabilité linéaire pour prédire le goutte à goutte polydisperse dans de telles configurations.

L'essentiel

Imaginez que vous essayez de créer un jet d'eau parfaitement régulier à partir d'un tuyau pour arroser votre jardin. Vous voulez que l'eau s'écoule de manière fluide sur une longue distance avant de se transformer en une pulvérisation uniforme de gouttelettes. C'est ce que les scientifiques appellent le « jet » (jetting). Cependant, il arrive parfois que l'eau commence à goutter dès la sortie du tuyau, créant une pulvérisation désordonnée et irrégulière. C'est ce qu'on appelle le « goutte-à-goutte » (dripping).

Ce document porte sur un type spécifique de dispositif de « tuyau » appelé dispositif de coflow (écoulement conjoint). Imaginez un tuyau de jardin à l'intérieur d'un tuyau plus large. Un flux extérieur de liquide en mouvement rapide pousse contre un flux intérieur plus lent, l'étirant en un cône mince et effilé (comme une forme de goutte d'eau) avant qu'il ne soit projeté sous forme de jet.

Les chercheurs ont voulu comprendre exactement quand un flux régulier se transforme en un goutte-à-goutte désordonné. Ils ont utilisé deux outils :

1. Expériences : Observer le flux réel d'un liquide dans un laboratoire.
2. Simulations informatiques : Utiliser les mathématiques pour prédire comment le liquide se comporte.

Voici la décomposition simple de ce qu'ils ont découvert et pourquoi cela est important :

1. La « boule de cristal » qui a échoué

Les scientifiques utilisent souvent une méthode appelée Analyse de Stabilité Linéaire Globale pour prédire quand un flux régulier va se transformer en goutte-à-goutte. Vous pouvez voir cette méthode comme une « boule de cristal » qui observe le flux régulier et demande : « Si je donne un petit coup à ce flux, va-t-il s'en remettre ou s'effondrer ? »

D'habitude, cette boule de cristal fonctionne bien. Elle prédit que si le flux est instable, le cône en forme de « goutte » à l'extrémité commencera à osciller et à se briser.

Mais dans ce cas précis, la boule de cristal était fausse.
Le modèle informatique (la boule de cristal) disait que le flux était stable et que le cône était parfaitement immobile. Cependant, l'expérience réelle montrait que le flux se brisait et gouttait réellement. Le modèle n'a pas vu le problème car il regardait la mauvaise chose. Il supposait que le cône en « goutte » était le point faible, alors qu'en réalité, le cône était correct ; c'était le flux mince qui en sortait qui posait problème.

2. Les « ondes fantômes » et l'explosion à court terme

Pourquoi le modèle a-t-il échoué ? Le document explique que le flux est comme un instrument de musique avec de nombreuses notes cachées (appelées modes propres ou eigenmodes).

• L'ancienne théorie : Les scientifiques pensaient que si le flux était instable, une « note sonore » spécifique (un mode propre instable) deviendrait de plus en plus forte jusqu'à ce que le flux se brise.
• La nouvelle découverte : Les chercheurs ont découvert que dans ce dispositif, toutes les « notes » essaient en fait de devenir plus silencieuses (elles sont en train de décroître). Cependant, pendant un court instant, ces notes qui s'affaiblissent peuvent interférer entre elles de manière à créer un pic d'énergie massif et temporaire.

L'analogie : Imaginez un groupe de personnes dans une pièce, essayant toutes de sortir calmement. Si elles se bousculent toutes exactement au même moment, elles peuvent créer un embouteillage chaotique et bruyant pendant une fraction de seconde avant de pouvoir enfin sortir. La boule de cristal du modèle ne regarde que le résultat à long terme (tout le monde sort calmement) et manque le chaos à court terme (l'embouteillage).

Ce chaos à court terme est ce qui provoque la rupture du flux et sa transformation en gouttelettes, même si les mathématiques disent que le flux devrait être stable.

3. La « poussée » compte

Les chercheurs ont également découvert que la façon dont on perturbe le flux est importante.

• Si vous touchez le flux juste à la pointe du cône, il ne se brisera peut-être pas.
• Si vous le touchez un peu plus bas dans le flux, il se brise beaucoup plus vite.

Cela signifie que la longueur du flux avant qu'il ne se brise n'est pas un nombre fixe écrit dans les lois de la physique pour ce montage spécifique. Elle dépend entièrement de l'endroit où la « poussée » initiale se produit. C'est comme pousser une balançoire : si vous la poussez au bon moment, elle va haut ; si vous la poussez au mauvais moment, elle bouge à peine.

4. L'observation du monde réel

Dans leurs expériences, les chercheurs ont observé ce qui se passait lorsqu'ils ralentissaient le flux du liquide intérieur :

• Flux élevé : Un long flux régulier se forme et se brise en gouttelettes uniformes loin du dispositif.
• Flux moyen : Le flux devient plus court et se brise plus près de la pointe, mais les gouttelettes sont encore majoritairement uniformes.
• Flux faible : Le flux se brise presque immédiatement, créant une pulvérisation désordonnée de gouttelettes de tailles différentes.

Le modèle informatique prédisait que la transition du « flux régulier » vers la « pulvérisation désordonnée » se produirait parce que le cône à l'extrémité commencerait à osciller. Mais l'expérience a montré que le cône restait parfaitement immobile tout le temps ! L'instabilité s'est produite dans le flux après qu'il ait quitté le cône.

L'essentiel

Ce document nous indique que pour ce type spécifique de dispositif de microfluidique, les outils mathématiques standards utilisés pour prédire la stabilité ne sont pas fiables. Ils manquent le « chaos à court terme » causé par l'interférence de différentes ondes de fluides.

Au lieu de chercher une seule « note instable » qui croît indéfiniment, nous devons comprendre comment un ensemble de « notes qui s'affaiblissent » peuvent s'entrechoquer pour provoquer une rupture soudaine. Cela change la façon dont les scientifiques doivent penser à la conception de ces dispositaux microfluidiques, car les anciennes règles ne s'appliquent pas ici.

Résumé technique

Résumé technique : Analyse de la stabilité de l'écoulement dans un dispositif en co-écoulement

Énoncé du problème
L'article étudie la stabilité hydrodynamique d'une configuration en co-écoulement, un dispositif microfluidique courant utilisé pour le « tip streaming » (formation de gouttelettes par l'extrémité), où un liquide interne est injecté à travers un capillaire entouré d'un flux coaxial externe. Bien que l'analyse de stabilité linéaire globale ait été appliquée avec succès à diverses écoulements capillaires (par exemple, jets gravitationnels, focalisation de flux) pour prédire les transitions entre les modes de jet et de gouttelettes (jetting et dripping), les tentatives précédentes n'ont pas réussi à prédire avec précision les longueurs de rupture de jet ou les conditions spécifiques de la transition vers un régime de gouttes polydisperses dans les dispositifs en co-écoulement. Le problème central abordé est la divergence entre les prédictions de l'analyse de stabilité linéaire globale et les observations expérimentales concernant la stabilité du mode de jet stationnaire et l'apparition de l'instabilité.

Méthodologie
Les auteurs emploient une approche combinée d'observation expérimentale, d'analyse de stabilité linéaire globale et de simulations numériques transitoires.

• Équations de conservation : L'écoulement est modélisé à l'aide des équations de Navier-Stokes sans dimension pour des fluides axisymétriques, incompressibles et immiscibles. Le système est caractérisé par le rapport de rayon $R$, le rapport de densité $\rho$, le rapport de viscosité $\mu$, le nombre d'Ohnesorge $Oh_i$, le nombre de Capillaire $Ca$ et le rapport de débit $Q$.
• Analyse de stabilité linéaire globale : Les auteurs supposent un écoulement de base parfaitement stationnaire et résolvent les valeurs propres de l'opérateur de Navier-Stokes linéarisé. Ils analysent le spectre pour déterminer si les perturbations croissent (instables) ou décroissent (stables) asymptotiquement.
• Simulations transitoires : Pour étudier la réponse à court terme de l'écoulement de base, les auteurs effectuent des simulations numériques directes des équations linéarisées. Ces simulations introduisent une petite perturbation initiale de faible amplitude à la surface libre et suivent l'évolution de l'interface. Crucialement, ces simulations tiennent compte de la superposition de multiples modes propres, plutôt que d'isoler un seul mode dominant.
• Validation expérimentale : Des expériences ont été menées à l'aide d'un dispositif en co-écoulement avec des propriétés géométriques et de fluides fixes ($R=6,67$, $\rho=0,962$, $\mu=36,1$, $Oh_i=0,0264$). La transition du régime de jet au régime de gouttes a été observée en faisant varier le rapport de débit $Q$.

Résultats clés

1. Échec de l'analyse de stabilité globale pour la transition jet-goutte : Les expériences révèlent qu'à mesure que le débit interne diminue, l'écoulement ne passe pas directement d'un jet stationnaire à un régime de gouttes polydisperses par la déstabilisation du ménisque. Au lieu de cela, l'écoulement passe par un état intermédiaire de micro-gouttelettes monodisperses. Dans cet état, le ménisque effilé reste stable et stationnaire, tandis que le jet émis se rompt en raison d'une instabilité convective. Par conséquent, l'analyse de stabilité globale, qui suppose un écoulement de base de jet stationnaire, ne parvient pas à prédire la transition vers les gouttes polydisperses car l'état marginalement stable réel (micro-gouttelettes monodisperses) est fondamentalement différent du jet stationnaire supposé.
2. Instabilité convective et croissance non modale : Dans le régime de jet, l'écoulement est convectivement instable. L'analyse de stabilité globale prédit que l'interface reste inchangée près du ménisque et que l'instabilité provient de modes propres parasites dépendant de la longueur du domaine. Cependant, les simulations transitoires montrent que le jet se rompt en raison de l'interférence constructive (superposition linéaire) de modes décroissants déclenchés par une perturbation initiale. Cette croissance « non modale » permet à l'énergie de la perturbation d'augmenter considérablement à court terme, entraînant la rupture du jet avant que les modes ne soient amortis.
3. Dépendance à la perturbation initiale : La longueur de rupture du jet n'est pas une propriété intrinsèque de l'écoulement de base. Les simulations démontrent que l'emplacement de la rupture est hautement sensible à la position et à la nature de la perturbation initiale. Les perturbations introduites plus près de la zone de transition ménisque-jet conduisent à une rupture plus précoce. Cela explique pourquoi les longueurs de rupture expérimentales sont souvent plus courtes que celles prédites par l'analyse de stabilité temporelle classique, qui suppose un mode dominant spécifique déclenché à l'origine du jet.
4. Effets non linéaires : Les effets non linéaires deviennent significatifs uniquement près de l'instant de la rupture de l'interface, provoquant la formation d'un étranglement et la rupture consécutive. Cependant, la croissance initiale menant à la rupture est adéquatement décrite par la superposition linéaire des modes.

Signification et conclusions
L'article conclut que l'analyse de stabilité linéaire globale classique est insuffisante pour décrire la stabilité des configurations en co-écoulement et potentiellement d'autres installations similaires. La limitation principale est que l'analyse repose sur l'hypothèse d'un écoulement de base stationnaire, qui n'existe pas au point de transition critique dans les dispositifs en co-écoulement ; l'écoulement passe en réalité par un état de micro-gouttelettes monodisperses. De plus, l'instabilité convective menant à la rupture du jet est pilotée par la nature non normale de l'opérateur linéarisé, où la superposition de modes décroissants provoque une croissance transitoire.

Les auteurs affirment que leurs conclusions remettent en question la validité de l'application de l'analyse de stabilité linéaire standard aux dispositifs en co-écoulement. Plus précisément, des critères tels que la condition de pincement de Briggs pour l'instabilité absolue peuvent être nécessaires mais restent insuffisants pour déterminer la stabilité de ces systèmes. Les résultats suggèrent que la longueur de rupture du jet ne peut être prédite uniquement par le taux de croissance d'un mode temporel dominant, car le processus est régi par l'interférence de multiples modes décroissants et dépend de la perturbation initiale. Ce travail souligne la nécessité de prendre en compte les mécanismes de croissance transitoire et non modale lors de l'analyse de la stabilité dans les configurations microfluidiques en co-écoulement.