A reduced model for droplet dynamics with interfacial viscosity

Cet article étend le modèle de Maffettone-Minale pour incorporer les viscosités de cisaillement et de dilatation interfaciales afin de décrire la déformation des gouttelettes en écoulement de cisaillement, et valide l'exactitude du modèle étendu qui en résulte par rapport à des simulations numériques entièrement résolues à travers divers nombres de capillarité et de Boussinesq.

L'essentiel

Imaginez une minuscule goutte d'huile flottant dans un courant d'eau. Si l'eau commence à couler plus vite, la goutte est écrasée et étirée, passant d'une sphère parfaite à une forme ovale. Les scientifiques essaient depuis longtemps de prédire exactement à quel point cette goutte va s'étirer en utilisant les mathématiques.

Pendant des décennies, ils ont utilisé une « recette » célèbre (appelée modèle de Maffettone–Minale) qui fonctionne bien pour les gouttes propres. Mais dans le monde réel, les gouttes ont souvent une « peau » composée de savon, de protéines ou d'autres molécules. Cette peau n'est pas seulement une limite ; elle possède sa propre épaisseur et sa propre adhérence, connue sous le nom de viscosité interfaciale. Imaginez que la goutte porte un pull collant et extensible.

Ce document présente une nouvelle recette améliorée (le modèle EMM ou Maffettone–Minale Étendu) qui tient compte de ce pull collant. Voici comment les auteurs ont décomposé la chose :

1. Les deux types de « collant »

Les auteurs ont réalisé que la peau de la goutte résiste au mouvement de deux manières différentes, et qu'ils devaient mesurer les deux :

• La viscosité de cisaillement (l'effet « élastique ») : Imaginez que vous essayiez de faire glisser votre main sur la surface de la goutte. Si la peau est « visqueuse au cisaillement », elle résiste à ce mouvement de glissement, comme si vous traîniez votre main dans du miel.
• La viscosité de dilatation (l'effet « respiration ») : Imaginez la goutte essayant d'étendre ou de rétrécir sa surface (comme un ballon qui se gonfle). Si la peau est « visqueuse à la dilatation », elle résiste à cet étirement ou à ce rétrécissement, comme un tissu serré et rigide qui ne veut pas s'étendre.

Le document utilise des nombres spéciaux (appelés nombres de Boussinesq) pour mesurer la force de ces deux résistances par rapport à l'épaisseur de la goutte.

2. La nouvelle recette (le modèle EMM)

Les auteurs ont pris l'ancienne recette mathématique simple et ont ajouté de nouveaux ingrédients pour gérer ces deux types de collant.

• L'objectif : Ils voulaient savoir : Jusqu'où peut-on étirer cette nouvelle recette avant qu'elle ne cesse de fonctionner ?
• La méthode : Ils n'ont pas fait que deviner. Ils ont construit une simulation informatique ultra-détaillée (comme un film en haute définition de la goutte) qui résolvait chaque infime règle de physique à partir de zéro. Cela servait de « vérité ».
• Le test : Ils ont fait fonctionner la nouvelle recette EMM aux côtés de la simulation ultra-détaillée. Ils ont comparé les résultats pour voir si la recette simple correspondait au film complexe.

3. Ce qu'ils ont trouvé

Les résultats étaient surprenants et spécifiques :

• Quand le « pull » est uniforme : Si la peau de la goutte résiste au glissement et à l'étirement de manière égale (un pull équilibré), la nouvelle recette fonctionne incroyablement bien, même lorsque la goutte est étirée de façon importante. Elle prédit avec précision la vitesse à laquelle la goutte s'étire et le temps qu'elle met pour prendre sa forme finale.
• Quand le « pull » est déséquilibré : Si la peau est très bonne pour résister au glissement mais mauvaise pour résister à l'étirement (ou vice versa), la recette simple commence à devenir un peu floue. Elle fonctionne toujours pour des écoulements doux, mais si l'écoulement devient trop fort, la recette perd en précision.
• L'effet de « ralentissement » : La découverte la plus intéressante concernait le temps. Quand la goutte possède les deux types de collant en même temps, elle met beaucoup plus de temps à changer de forme. C'est comme si la goutte restait « coincée » dans sa propre peau. Les auteurs ont trouvé que leur nouvelle recette capture parfaitement cet effet de « ralenti ».
• Le point de rupture : Si la goutte a presque aucune résistance au glissement (mais une résistance élevée à l'étirement), elle est tellement étirée qu'elle finit par se déchirer. La nouvelle recette prédit correctement que cela arrive plus tôt dans ces conditions spécifiques.

4. L'essentiel à retenir

Les auteurs ont réussi à créer un outil simple, rapide et fiable pour prédire comment les gouttes avec des « peaux collantes » se comportent dans des liquides en mouvement.

• Pourquoi c'est important : Cela évite aux scientifiques de devoir lancer des simulations informatiques massives et lentes pour chaque problème.
• Le bémol : L'outil est très précis pour les étirements de petite à moyenne taille, surtout lorsque la résistance de la peau est équilibrée. Si l'écoulement est extrêmement violent ou si la résistance de la peau est très déséquilibrée, l'outil commence à perdre sa précision, et il faut alors utiliser les simulations informatiques lourdes.

En bref, ils ont mis à jour le « calculateur de gouttes » pour gérer les peaux collantes, prouvant qu'il fonctionne très bien pour la plupart des scénarios du quotidien, tout en marquant clairement les limites où il pourrait avoir besoin d'un peu d'aide.

Résumé technique

Résumé Technique : Un modèle réduit pour la dynamique des gouttelettes avec viscosité interfaciale

Énoncé du problème
Comprendre la déformation des gouttelettes dans les écoulements hydrodynamiques est un sujet central en rhéologie. Bien que le modèle classique de Maffettone–Minale (MM) fournisse une description phénoménologique réussie pour les gouttelettes propres (celles sans viscosité interfaciale) en écoulement de cisaillement, il ne rend pas compte des complexités introduites par la viscosité interfaciale. Dans de nombreux scénarios pratiques, tels que lorsque les interfaces sont recouvertes de tensioactifs, de polymères ou de protéines, les viscosités de cisaillement ($\mu_s$) et de dilatation ($\mu_d$) de l'interface deviennent significatives. Ces viscosités introduisent une friction supplémentaire qui modifie l'équilibre des contraintes entre les fluides globaux, influençant la déformation, la relaxation et les seuils de rupture. Le défi consiste à développer un modèle d'ordre réduit qui intègre ces effets interfaciaux tout en conservant l'efficacité de calcul et la tractabilité analytique du cadre MM original.

Méthodologie
Les auteurs proposent un modèle Maffettone–Minale Étendu (EMM), qui généralise l'équation MM originale pour l'évolution d'un tenseur morphologique $S$ (décrivant la forme de la gouttelette comme un ellipsoïde).

• Extension théorique : Le modèle EMM intègre les viscosités interfaciales via les nombres de Boussinesq, $Bqs = \mu_s/\mu R$ et $Bqd = \mu_d/\mu R$. En se référant aux calculs perturbatifs pour les interfaces visqueuses (spécifiquement l'extension de Narsimhan des travaux de Rallison), les auteurs dérivent des expressions généralisées pour les coefficients du modèle $f_1$ et $f_2$. Ces coefficients, qui étaient auparavant uniquement des fonctions du rapport de viscosité $\lambda$ pour les gouttelettes propres, dépendent désormais de $\lambda$, $Bqs$ et $Bqd$.
• Validation numérique : Pour quantifier l'exactitude et la plage de validité du modèle EMM, les auteurs effectuent une comparaison systématique avec des Simulations Entièrement Résolues (FRS). Les FRS utilisent une méthode de Boltzmann sur réseau à frontière immergée pour résoudre la dynamique couplée des fluides interne et externe ainsi que l'interface déformable sans supposer de petites déformations.
• Analyse d'erreur : L'étude évalue l'erreur relative moyenne $\langle \epsilon \rangle_T$ entre la déformation temporelle $D(t)$ prédite par le modèle EMM et les résultats des FRS à travers un large espace de paramètres défini par le nombre capillaire ($Ca$) et les deux nombres de Boussinesq.

Contributions clés

1. Dérivation de coefficients généralisés : L'article fournit des expressions analytiques explicites pour les coefficients de l'EMM $f_1^{(EMM)}$ et $f_2^{(EMM)}$ qui tiennent compte à la fois des viscosités de cisaillement et de dilatation interfaciales. Ces expressions se réduisent aux coefficients MM originaux lorsque les viscosités interfaciales sont nulles.
2. Caractérisation de la dynamique transitoire : L'étude révèle que la viscosité interfaciale affecte le temps de déformation transitoire de manière non linéaire. Alors qu'une gouttelette présentant une viscosité principalement de cisaillement ou principalement de dilatation se comporte de manière similaire à une gouttelette propre concernant le temps de relaxation, une gouttelette ayant des valeurs comparables des deux viscosités subit un ralentissement dramatique de la dynamique de déformation.
3. Cartographie du domaine de validité : Par des tests numériques approfondis, les auteurs délimitent les régions spécifiques de l'espace de paramètres $(Ca, Bqs, Bqd)$ où le modèle EMM reste quantitativement précis.

Résultats

• Déformation à l'état stationnaire : Le modèle EMM prédit avec précision la déformation à l'état stationnaire $D_s$ pour une large gamme de paramètres. Le modèle capture la tendance contre-intuitive où l'augmentation de la viscosité de dilatation (avec une viscosité de cisaillement nulle) augmente légèrement la déformation, tandis que l'augmentation de la viscosité de cisaillement (ou le fait d'avoir des viscosités de cisaillement et de dilatation égales) réduit la déformation.
• Comportement transitoire : Le modèle reproduit avec succès la transition d'une relaxation monotone vers des réponses oscillatoires observées dans les FRS, particulièrement lorsque $Bqs = Bqd$. Le temps de déformation caractéristique augmente brusquement uniquement lorsque les deux viscosités interfaciales sont notables.
• Limites de précision :
  • Pour de faibles viscosités interfaciales ($Bqs, Bqd \ll 1$) et de faibles nombres capillaires ($Ca \le 0,2$), le modèle EMM reproduit les données FRS avec des erreurs typiquement inférieures à quelques pourcents.
  • La précision du modèle se dégrade plus rapidement dans les configurations de viscosité fortement asymétriques (où un nombre de Boussinesq domine) par rapport au cas équiviscous ($Bqs = Bqd$).
  • À mesure que $Ca$ augmente (par exemple, $Ca = 0,5$), les couplages non linéaires dans la contrainte interfaciale, que le modèle EMM de petite déformation approxime, entraînent des écarts plus importants. À $Ca$ élevé et faible $Bqs$, le modèle approche du seuil de rupture, où les FRS indiquent une instabilité numérique due à de grandes déformations.

Signification
L'article établit le modèle EMM comme un outil robuste et peu coûteux en calcul pour prédire la déformation des gouttelettes en présence de viscosité interfaciale. Il démontre que le modèle conserve un lien clair avec les théories classiques de petite déformation tout en offrant un contrôle indépendant sur les effets de cisaillement et de dilatation. Ce travail fournit une carte quantitative de la validité du modèle, soulignant que, bien qu'il soit très efficace pour des déformations modérées et des rapports de viscosité spécifiques, les simulations entièrement résolues restent indispensables pour les grandes déformations ou les conditions interfaciales fortement asymétriques. Les auteurs suggèrent que ce cadre pourrait servir de base à de futures extensions vers des configurations d'écoulement plus complexes et des interfaces chargées en tensioactifs, bien que de telles généralisations soient notées comme des directions futures plutôt que des accomplissements actuels.