Internal flow and concentration in neighbouring evaporating binary droplets and rivulets

Cette étude numérique et théorique examine comment la proximité de gouttes ou de rivières évaporantes binaires induit un effet de blindage et des asymétries dans les écoulements et les profils de concentration, dont l'évolution dépend de paramètres clés tels que l'angle de contact, la distance inter-éléments et, dans le cas des gouttes, les forces de tension superficielle (nombre de Marangoni).

L'essentiel

🌧️ La Danse des Gouttes : Quand l'évaporation crée des courants secrets

Imaginez que vous avez laissé tomber deux petites gouttes d'un mélange (comme de l'eau et un peu de sirop) l'une à côté de l'autre sur une table. Vous savez que l'eau va s'évaporer, mais saviez-vous que ces gouttes ne sont pas de simples réservoirs passifs ? À l'intérieur, il se passe une véritable danse liquide, un ballet de courants invisibles qui changent tout au fur et à mesure que l'eau disparaît.

Les scientifiques de l'Université de Twente (Pim Dekker, Duarte Rocha et leurs collègues) ont décidé de regarder sous le capot de ce phénomène pour comprendre comment deux gouttes voisines s'influencent mutuellement.

1. Le problème : Les gouttes qui se "protègent" mutuellement

Quand une seule goutte s'évapore, elle perd de l'eau uniformément (ou presque). Mais quand deux gouttes sont proches, elles créent un effet de bouclier.

• L'analogie : Imaginez deux personnes qui parlent fort dans une pièce calme. Si elles sont loin l'une de l'autre, tout le monde les entend bien. Si elles sont très proches, leurs voix se mélangent et créent une zone de bruit intense entre elles. De la même façon, les gouttes saturent l'air entre elles en vapeur d'eau.
• La conséquence : L'air entre les gouttes devient "saturé", ce qui ralentit l'évaporation dans cette zone. C'est comme si les gouttes se protégeaient mutuellement de l'air sec. Ce ralentissement n'est pas uniforme : il est plus fort au centre (entre les gouttes) et plus faible sur les côtés.

2. Le moteur invisible : La tension de surface

C'est ici que ça devient fascinant. Ces gouttes sont des mélanges (eau + un ingrédient qui ne s'évapore pas, comme du 1,2-hexanediol).

• Le mécanisme : L'eau s'évapore plus vite sur les bords que dans le centre (à cause du bouclier). Cela crée une différence de concentration : le centre devient plus "sirupeux" (riche en ingrédient non volatil) que les bords.
• L'effet Marangoni : Pensez à la tension de surface comme à une peau élastique qui recouvre la goutte. Cette peau est plus tendue là où il y a plus de sirop. La peau va donc tirer le liquide des zones "lâches" (les bords) vers les zones "tendues" (le centre).
• Le résultat : Cela crée un courant interne puissant qui fait tourner le liquide comme un tourbillon, emportant les particules avec lui.

3. La découverte clé : La goutte vs Le ruisseau

Les chercheurs ont comparé deux formes :

1. Les gouttes (3D) : Comme des dômes de savon.
2. Les rivières (2D) : Comme de longs rubans de liquide (ce qu'ils appellent des "rivulets").

Ils ont cherché à trouver le point mort (ou point de stagnation) : l'endroit précis à l'intérieur de la goutte où le courant s'arrête et où les deux tourbillons se rencontrent. C'est comme le centre d'une rivière où l'eau ne coule ni vers l'avant ni vers l'arrière.

La grande différence trouvée :

• Pour les rivières (2D) : Le point mort est très prévisible. Il dépend de la distance entre les rivières et de l'angle qu'elles forment avec la table. Curieusement, la force du courant (Marangoni) n'a presque aucune influence sur sa position. C'est comme si le courant était bloqué dans un couloir étroit et ne pouvait pas faire de détours.
• Pour les gouttes (3D) : C'est beaucoup plus compliqué ! Le point mort bouge selon la force du courant. Pourquoi ? Parce que les gouttes sont rondes, le liquide peut tourner dans toutes les directions (y compris autour de l'axe, comme une toupie). Cette liberté de mouvement crée des tourbillons complexes qui déplacent le point mort de manière imprévisible.

4. Pourquoi est-ce important ?

Cette recherche n'est pas juste de la théorie. Elle aide à comprendre des technologies du quotidien :

• L'impression jet d'encre : Pour éviter que les gouttes d'encre ne laissent des traces en forme d'anneau (l'effet "tache de café").
• Le refroidissement par pulvérisation : Pour refroidir des moteurs ou des serveurs informatiques de manière efficace.
• Les pesticides : Pour que les produits chimiques restent bien répartis sur les feuilles.

En résumé

Les chercheurs ont créé des modèles mathématiques et des simulations pour prédire comment ces gouttes se comportent. Ils ont découvert que :

1. Les gouttes voisines se gênent mutuellement, créant des courants internes asymétriques.
2. Pour les formes allongées (rivulets), la physique est simple et prévisible.
3. Pour les gouttes rondes, la physique devient un jeu de billard complexe où la force du courant change tout.

C'est une victoire pour la science : en comprenant ces petits tourbillons invisibles, nous pouvons mieux contrôler les technologies qui utilisent des gouttes, de l'imprimante de votre bureau au système de refroidissement de votre voiture.

Résumé technique

1. Problématique et Contexte

L'évaporation de gouttes ou de ruisselets (rivulets) binaires voisins est un phénomène critique dans de nombreuses applications industrielles telles que l'impression jet d'encre, le refroidissement par pulvérisation et la délivrance de pesticides. Contrairement aux gouttes isolées, la proximité de ces corps liquides induit un effet d'ombrage (shielding effect) : la concentration en vapeur dans la phase gazeuse augmente localement, réduisant le taux d'évaporation individuel et modifiant la dynamique des écoulements internes.

Dans le cas de mélanges binaires (par exemple, eau et 1,2-hexanediol), l'évaporation sélective d'un composant crée des gradients de tension de surface, générant des écoulements de Marangoni. Ces écoulements peuvent être beaucoup plus intenses que l'effet « tache de café » (flux vers le bord de contact) et peuvent soit l'atténuer, soit l'amplifier. L'objectif principal de cette étude est de comprendre comment la symétrie des champs d'écoulement et de concentration est brisée par la proximité, et d'identifier les paramètres clés gouvernant cette asymétrie : l'angle de contact ($\theta$), la distance inter-objets ($b$) et le nombre de Marangoni ($Ma$).

2. Méthodologie

Les auteurs ont adopté une approche multi-échelle combinant simulations numériques directes et modélisation théorique simplifiée :

• Simulations Transitoires 2D (Ruisselets) : Pour contourner la complexité computationnelle d'un problème 3D complet, les auteurs ont d'abord résolu les équations complètes (Navier-Stokes, advection-diffusion, diffusion de vapeur) pour deux ruisselets voisins en utilisant une méthode éléments finis Lagrangien-Eulerien (via le code pyoomph). Ces simulations ont permis de valider les hypothèses de base.
• Modèle Quasi-Stationnaire : Compte tenu de l'évaporation lente des mélanges eau/1,2-hexanediol, les auteurs ont introduit un modèle quasi-stationnaire. Ce modèle suppose que l'écoulement interne atteint un état stable à chaque instant, permettant de négliger les termes transitoires d'inertie et de simplifier les équations de transport de masse.
• Approximation Lubrification : Pour réduire le problème à deux dimensions et permettre une analyse analytique, une approximation de lubrification a été appliquée (valable pour de petits angles de contact). Ce modèle intègre la dispersion de Taylor-Aris pour tenir compte du mélange vertical de la concentration.
• Flux Évaporatif : Des expressions analytiques pour le flux évaporatif local ont été utilisées, tenant compte de l'effet d'ombrage entre les objets voisins, bien que des approximations soient nécessaires pour les angles de contact non nuls.

3. Résultats Clés

A. Validation des Modèles

Les simulations transitoires complètes montrent que le modèle quasi-stationnaire reproduit avec une grande précision les caractéristiques de l'écoulement et de la concentration, à condition que le nombre de Marangoni soit suffisamment élevé pour que le mélange induit par Marangoni domine le transport par l'effet « tache de café ».

B. Comportement des Ruisselets (2D)

Pour les ruisselets, l'asymétrie de l'écoulement est quantifiée par la position du point de stagnation ($x_0$) à l'interface, là où les deux tourbillons internes se rejoignent.

• Indépendance vis-à-vis de $Ma$ : De manière surprenante, pour les ruisselets, la position du point de stagnation est indépendante du nombre de Marangoni. Elle ne dépend que de l'angle de contact ($\theta$) et de la distance inter-ruisselets ($b$).
• Analytique : Une solution analytique a été dérivée montrant que le déplacement du point de stagnation est déterminé par l'équilibre entre le flux évaporatif local et la géométrie du ruisseau, sans dépendance explicite de la force de Marangoni dans ce cas 2D.
• Évolution temporelle : Au cours de l'évaporation, alors que l'angle de contact diminue, le point de stagnation migre vers le centre du ruisseau.

C. Comportement des Gouttes (3D)

L'application du modèle de lubrification aux gouttes révèle des comportements différents dus à la géométrie axisymétrique :

• Dépendance à $Ma$ : Contrairement aux ruisselets, la position du point de stagnation dans les gouttes dépend fortement du nombre de Marangoni.
• Mécanisme : Cette différence s'explique par la présence d'un écoulement azimutal supplémentaire dans les gouttes. L'augmentation de $Ma$ modifie la distribution de la concentration de manière non linéaire, créant des gradients verticaux de vitesse significatifs et déplaçant le point de stagnation différemment selon qu'il est mesuré à l'interface ou au substrat.
• Influence de l'angle de contact : Pour les gouttes, l'asymétrie dépend également de l'angle de contact, même avec des approximations de flux pour des gouttes plates.

4. Contributions Principales

1. Identification des paramètres de contrôle : Démonstration que l'asymétrie des écoulements internes est principalement contrôlée par $\theta$, $b$ et $Ma$, mais avec des dépendances radicalement différentes entre géométries 2D (ruisselets) et 3D (gouttes).
2. Modélisation simplifiée : Développement et validation d'un modèle quasi-stationnaire et d'une approximation de lubrification capables de prédire avec précision les champs de concentration et d'écoulement sans le coût computationnel des simulations transitoires 3D complètes.
3. Résultat analytique pour les ruisselets : Dérivation d'une expression analytique pour le déplacement du point de stagnation dans les ruisselets, prouvant son indépendance vis-à-vis de la force de Marangoni dans ce régime 2D.
4. Distinction 2D/3D : Mise en évidence du rôle crucial de la dimensionnalité (écoulement azimutal) dans la réponse des écoulements aux forces de Marangoni.

5. Signification et Perspectives

Ce travail fournit un cadre théorique et numérique essentiel pour la prédiction du comportement de mélanges binaires évaporants dans des configurations denses, fréquentes en impression jet d'encre et en dépôt de revêtements.

• Implications pratiques : La compréhension de l'asymétrie permet de mieux contrôler le dépôt de solutés (évitant ou favorisant l'effet tache de café) en ajustant la distance entre les gouttes ou la composition du fluide.
• Limites et travaux futurs : Le modèle actuel suppose un écoulement quasi-stationnaire, ce qui exclut les phases finales où l'effet « tache de café » redevient dominant (singularité temporelle). De plus, l'absence d'expressions analytiques précises pour le flux évaporatif de gouttes/ruisselets avec des angles de contact non nuls reste un défi. Les auteurs suggèrent que la généralisation de ces flux pour des angles de contact finis est une priorité pour les recherches futures.

En résumé, l'article établit que la géométrie (2D vs 3D) dicte la sensibilité des écoulements internes aux forces de Marangoni, offrant des outils de modélisation robustes pour optimiser les procédés industriels impliquant l'évaporation de gouttes multiples.