Numerical Study on the Effects of Wall Wettability on Heat Transfer and Flow Characteristics during Flow Boiling in Microchannels

Résumé technique : Étude numérique des effets de la mouillabilité des parois sur le transfert de chaleur et les caractéristiques d'écoulement lors de l'ébullition en microcanaux

Énoncé du problème
La miniaturisation des puces à haute performance a entraîné une augmentation significative de la densité de flux thermique local, dépassant souvent 1000 W/cm², ce que les méthodes de refroidissement traditionnelles peinent à dissiper. L'ébullition en écoulement dans les microcanaux est une solution prometteuse en raison de ses coefficients de transfert de chaleur élevés et de sa capacité d'absorption de la chaleur latente. Cependant, les performances de ces systèmes sont fortement influencées par la dynamique des bulles et l'évolution des régimes d'écoulement. Un facteur critique, bien que complexe, dans ce processus est la mouillabilité des parois. Alors que la littérature existante a exploré les surfaces hydrophiles, hydrophobes et à mouillabilité mixte, les mécanismes couplés régissant la manière dont les conditions de mouillabilité extrêmes affectent l'interaction entre le comportement des bulles, la stabilité de l'écoulement, l'efficacité du transfert de chaleur et l'irréversibilité thermodynamique (génération d'entropie) restent insuffisamment compris. Plus précisément, il est nécessaire de quantifier les compromis entre l'amélioration du transfert de chaleur et la stabilité de l'écoulement sous différents régimes de mouillabilité.

Méthodologie
Cette étude emploie un modèle numérique tridimensionnel basé sur la méthode du Volume de Fluide (VOF) pour simuler l'ébullition en écoulement dans un microcanal rainuré. La simulation utilise le modèle de Force de Surface Continue (CSF) pour capturer les effets de tension superficielle et un modèle d'évaporation-condensation (basé sur Lee et al.) pour gérer le changement de phase.

• Modèle physique : La géométrie consiste en un microcanal avec des rainures triangulaires sur les parois latérales, conçu pour étudier la nucléation et la croissance des bulles. Le fluide de travail est de l'eau désionisée et le domaine solide est du silicium.
• Conditions de mouillabilité : Deux conditions de surface extrêmes sont comparées : une paroi hydrophile (angle de contact $\theta = 30^\circ$) et une paroi hydrophobe ($\theta = 150^\circ$).
• Approche numérique : Les équations de conservation (masse, quantité de mouvement et énergie) sont résolues à l'aide de la méthode des volumes finis avec l'algorithme PISO pour le couplage pression-vitesse. Une étude d'indépendance de maillage a confirmé qu'un maillage d'environ 275 000 éléments offre une précision suffisante. Le modèle a été validé par rapport à des données expérimentales, montrant une divergence maximale de moins de 8 % dans l'élévation de la température de la base de la puce.
• Métriques d'analyse : L'étude évalue la dynamique des bulles, les régimes d'écoulement, les fluctuations de température de sortie, le nombre de Nusselt (Nu), la chute de pression et la génération d'entropie (composantes thermique et frictionnelle).

Principales contributions et résultats

1. Dynamique des bulles et évolution du régime d'écoulement :

   • Surfaces hydrophiles ($\theta = 30^\circ$) : Les interactions liquide-solide fortes provoquent l'adhérence des bulles à la paroi et leur coalescence, formant des structures de film liquide. Cela favorise une transition vers un écoulement annulaire mais entrave le détachement des bulles.
   • Surfaces hydrophobes ($\theta = 150^\circ$) : Les bulles se détachent rapidement et sont transportées dans le flux principal. Cela accélère l'accumulation de vapeur et conduit à une transition plus précoce de l'écoulement de type bulle vers l'écoulement de type bouchon (slug flow) et l'écoulement annulaire.

2. Performance thermique :

   • Efficacité du transfert de chaleur : La surface hydrophobe démontre une performance de transfert de chaleur supérieure. Le détachement rapide des bulles facilite l'évacuation de la chaleur latente, maintenant des températures de paroi plus basses et des nombres de Nusselt plus élevés. À l'inverse, la surface hydrophile souffre d'une accumulation de chaleur due aux bulles attachées, entraînant des températures de paroi plus élevées et un déclin continu du nombre de Nusselt.
   • Stabilité de la température : Bien que les surfaces hydrophobes améliorent le transfert de chaleur moyen, elles présentent une instabilité de la température de sortie plus faible par rapport aux surfaces hydrophiles. La surface hydrophile montre des fluctuations de température plus importantes en raison de l'entrave à la dissipation de chaleur par les films liquides persistants.

3. Caractéristiques d'écoulement et instabilité :

   • Chute de pression : Le canal hydrophobe subit des chutes de pression nettement plus élevées et des oscillations de pression plus prononcées. Ceci est attribué à la fraction volumique de vapeur plus élevée et à la décharge périodique des bulles modifiant la structure locale de l'écoulement.
   • Champs de vitesse : Le canal hydrophobe présente des perturbations de flux locales plus fortes et des gradients de vitesse plus élevés en raison de la dynamique active des bulles, tandis que le canal hydrophile présente des variations de vitesse plus fluides et plus faibles.

4. Analyse de la génération d'entropie :

   • Thermique vs Frictionnelle : L'étude distingue la génération d'entropie thermique (induite par les gradients de température) de la génération d'entropie frictionnelle (induite par la dissipation visqueuse).
   • Mécanisme dominant : La génération d'entropie thermique est le contributeur dominant à la génération totale d'entropie. La surface hydrophile génère une entropie thermique nettement plus élevée en raison des grands gradients de température causés par la détérioration du transfert de chaleur.
   • Compromis : Bien que la surface hydrophobe génère plus d'entropie frictionnelle en raison des gradients de vitesse intenses, sa génération totale d'entropie reste inférieure à celle de la surface hydrophile grâce à la réduction substantielle de l'irréversibilité thermique.

Signification et affirmations
Les auteurs affirment que cette étude fournit des orientations théoriques précieuses pour la conception de la mouillabilité et l'optimisation des systèmes thermohydrauliques dans les dissipateurs thermiques à microcanaux sous des flux thermiques élevés. Les résultats mettent en évidence un compromis fondamental : alors que les surfaces hydrophobes améliorent considérablement le transfert de chaleur par ébullition et réduisent l'irréversibilité thermique, elles introduisent une plus grande résistance à l'écoulement et des fluctuations de pression. En revanche, les surfaces hydrophiles, bien qu'offrant des chutes de pression plus stables, souffrent d'une grave détérioration du transfert de chaleur et de pertes thermodynamiques plus élevées.

L'étude conclut que l'optimisation de la mouillabilité nécessite de trouver un équilibre entre l'amélioration du transfert de chaleur par changement de phase et la pénalité de l'augmentation de l'instabilité de l'écoulement. Les résultats suggèrent que pour les applications privilégiant la performance thermique et la minimisation de l'entropie, les surfaces hydrophobes sont avantageuses, à condition que le système puisse accommoder la chute de pression et les oscillations de flux associées. L'étude ne propose pas d'applications futures spécifiques ou de validations expérimentales au-delà de la portée de l'investigation numérique actuelle, mais souligne la nécessité de comprendre ces mécanismes couplés pour une conception efficace de la gestion thermique.