Characterisation of rough-wall drag in compressible turbulent boundary layers

Cette étude expérimentale sur les couches limites turbulentes compressibles révèle que, bien que le déficit de vitesse soit peu sensible aux transformations de vitesse, le régime pleinement rugueux présente un décalage dépendant du nombre de Mach qui est mieux corrigé par un facteur empirique lié au rapport des températures d'extrémité et de paroi, soulignant ainsi la nécessité de développer de futures transformations spécifiques aux parois rugueuses.

L'essentiel

🌪️ Le Voyage des Avions : Quand le Ciel Devient un Mur de Sable

Imaginez que vous êtes un avion. Quand vous volez lentement (en dessous de la vitesse du son), l'air autour de vous se comporte comme un fleuve tranquille. Si votre carlingue est un peu rugueuse (peut-être à cause de la saleté, de la glace ou de petits impacts), cela crée une traînée, un peu comme si vous glissiez sur du papier de verre mouillé. Les scientifiques connaissent très bien comment calculer cette résistance pour les avions lents.

Mais imaginez maintenant que vous accélérez pour voler plus vite que le son (Mach 2, Mach 3...). L'air ne se comporte plus comme un fleuve, mais comme un mur de briques solide. À ces vitesses, l'air devient "compressible" : il se comprime, chauffe et crée des ondes de choc (des murs d'air invisibles mais puissants).

Le problème ? Si votre avion a une surface rugueuse à ces vitesses supersoniques, ces petites aspérités ne créent pas seulement de la friction, elles agissent comme de minuscules obstacles qui génèrent leurs propres mini-ondes de choc. C'est comme si chaque grain de sable sur votre aile lançait sa propre petite explosion d'air. Cela ajoute une résistance supplémentaire, appelée "traînée d'onde", que les anciennes formules ne prenaient pas en compte.

🔍 L'Enquête des Scientifiques

Les chercheurs de cette étude (D. Wangsawijaya et son équipe) se sont posé une question cruciale :

"Si nous savons exactement combien une surface rugueuse freine un avion lent, pouvons-nous utiliser cette même information pour prédire combien elle freinera un avion supersonique ?"

Pour répondre, ils ont construit un tunnel à vent géant en Allemagne capable de simuler des vitesses allant du ralenti (0,3 Mach) jusqu'à des vitesses supersoniques impressionnantes (2,9 Mach, soit presque 3 fois la vitesse du son). Ils ont testé deux types de surfaces rugueuses, comme du papier de verre (un grain fin P60 et un grain grossier P24), collés sur une plaque.

🧪 Les Découvertes : Le "Choc" de la Rugosité

Voici ce qu'ils ont découvert, traduit en images simples :

1. Le vieux plan ne fonctionne plus : Les formules utilisées pour les avions lents (qui fonctionnent comme une règle à calculer parfaite) échouent complètement pour les avions rapides. Si on les utilise telles quelles, on sous-estime grandement la résistance de l'air.

2. L'effet "Mur de Choc" : À grande vitesse, la rugosité de la surface crée des ondes de choc. C'est comme si chaque grain de sable sur l'aile devenait un petit rocher dans un torrent rapide, créant des vagues qui remontent jusqu'au ciel. Plus l'avion va vite, plus ces vagues sont fortes et plus elles freinent l'avion.

3. La découverte magique (Le "Correcteur") : Les chercheurs ont essayé plusieurs façons de corriger leurs calculs.

   • Essai 1 : Changer la taille du grain de sable dans la formule. (Ça ne marche pas très bien).
   • Essai 2 : Changer la façon dont on mesure la viscosité de l'air. (Un peu mieux, mais pas parfait).
   • Essai Gagnant : Ils ont trouvé un facteur de correction basé sur la température.

   L'analogie du thermostat : Imaginez que l'air autour de l'avion est une pièce. Quand l'avion va très vite, l'air à l'avant est brûlant (à cause de la compression), mais l'air collé à la peau de l'avion (la paroi) peut être plus froid ou plus chaud selon la situation. Les chercheurs ont découvert que si l'on ajuste leur formule en fonction de ce différentiel de température entre l'air libre et la surface de l'avion, les données de tous leurs tests (lents et rapides) s'alignent parfaitement sur une seule courbe magique. C'est comme si on avait trouvé la clé universelle pour traduire la "rugosité lente" en "rugosité rapide".

🚀 Pourquoi est-ce important ?

Aujourd'hui, quand on conçoit un avion hypersonique (comme un futur avion de ligne supersonique ou un missile), on doit savoir exactement combien de carburant il faudra pour vaincre la friction de l'air. Si la surface est rugueuse (à cause de l'usure, de la glace ou de la poussière), et qu'on utilise les anciennes formules, on risque de sous-estimer la consommation de carburant ou de ne pas avoir assez de puissance.

Cette étude nous dit : "Ne regardez pas seulement la taille des grains de sable, regardez aussi la température de l'air autour d'eux."

🔮 Le Futur : Vers une Carte Précise

Bien que cette nouvelle formule (le "correcteur de température") fonctionne très bien pour les données actuelles, les chercheurs sont honnêtes : c'est encore une solution "empirique" (trouvée par l'expérience, pas encore déduite d'une théorie parfaite).

Pour l'avenir, ils veulent :

1. Mesurer directement la température et la friction sur les surfaces rugueuses (ce qui est très difficile à faire dans un tunnel à vent supersonique).
2. Créer une nouvelle "règle de transformation" spéciale pour les surfaces rugueuses, qui prendrait en compte non seulement la température, mais aussi la forme exacte des grains de sable et les ondes de choc qu'ils génèrent.

En résumé, cette recherche est un grand pas de géant pour comprendre comment les avions rapides interagissent avec des surfaces imparfaites, en nous donnant un outil pour mieux prédire leur comportement dans le ciel brûlant du supersonique.

Résumé technique

1. Problématique et Contexte

Dans les écoulements turbulents compressibles (Mach $M > 1$), la présence de rugosité de surface (due à l'ablation, aux impacts de poussière, à la glace, etc.) génère une traînée supplémentaire complexe. Contrairement aux écoulements incompressibles où la traînée de rugosité est bien caractérisée par la fonction de rugosité de Hama ($\Delta U^+$) et la rugosité équivalente en grain de sable ($k_s$) de Nikuradse, le régime compressible introduit des phénomènes supplémentaires, notamment la traînée d'onde (wave drag) générée par les ondes de choc formées devant les éléments de rugosité.

Le problème central abordé par les auteurs est le suivant : Si la rugosité équivalente $k_s$ d'une surface est connue en régime incompressible, dans quelles conditions cette valeur peut-elle être utilisée pour prédire la traînée en régime compressible ?

Les études existantes souffrent de limitations majeures :

• Elles sont souvent réalisées à un seul nombre de Mach et un seul nombre de Reynolds.
• Elles imposent souvent la rugosité à la relation de Nikuradse (incompressible) sans vérifier la validité de cette hypothèse en régime compressible.
• Il existe un manque de données expérimentales couvrant une large gamme de nombres de Reynolds ($Re_\tau$) et de Mach, avec une variation indépendante de ces paramètres.

2. Méthodologie

Les auteurs ont mené une campagne expérimentale complète dans la soufflerie trisonique intermittente (TWM) de l'Université des Forces Armées de Munich.

• Configuration expérimentale :

  • Écoulement : Couches limites turbulentes (TBL) sur des parois rugueuses.
  • Rugosité : Deux types de papier de verre (P60 et P24) utilisés comme rugosité réaliste de type "sable". La topographie a été scannée optiquement pour déterminer les paramètres physiques (hauteur $k$, distribution, etc.).
  • Plage de paramètres :
    • Nombre de Mach : $0.3 \le M \le 2.9$ (subsonique, transsonique, supersonique).
    • Nombre de Reynolds basé sur la friction : $7427 \le Re_\tau \le 30292$.
    • Variation indépendante de $Re_\tau$ à $M=2$ (en faisant varier la pression d'arrêt).
  • Mesures : Vélométrie par image de particules (PIV) 2D2C pour obtenir les profils de vitesse moyenne.
  • Limitations : La température de paroi ($T_w$) et la contrainte de cisaillement pariétal (WSS) n'ont pas été mesurées directement. Elles ont été estimées via des hypothèses de paroi adiabatique et des méthodes d'ajustement (fitting) de type Clauser modifié.

• Approche d'analyse :

  • Application de transformations de vitesse classiques (Van Driest, Howarth, Trettel & Larsson, etc.) pour "étirer" les profils de vitesse compressibles vers un équivalent incompressible.
  • Calcul du déficit de quantité de mouvement $\Delta U^+$ (fonction de Hama).
  • Comparaison avec la loi logarithmique de Nikuradse pour le régime pleinement rugueux.
  • Test de trois échelles de correction empiriques pour réconcilier les données compressibles avec le cadre incompressible.

3. Contributions Clés et Résultats

A. Indépendance vis-à-vis des transformations de vitesse

Les résultats montrent que la fonction de rugosité $\Delta U^+$ est largement insensible au choix de la transformation de vitesse utilisée (VD, TL, VO, etc.) pour corriger les effets de compressibilité. Cependant, en régime pleinement rugueux, l'ordonnée à l'origine de la loi logarithmique subit un décalage dépendant du nombre de Mach.

B. Échec de l'utilisation directe de $k_s$ incompressible

L'utilisation directe de la rugosité équivalente $k_s$ déterminée en régime incompressible ($M=0.3$) ne permet pas de faire coïncider les données supersoniques avec la loi de Nikuradse. Un décalage croissant est observé avec l'augmentation de $M$, suggérant l'apparition d'une traînée d'onde non capturée par les modèles incompressibles.

C. Évaluation de trois facteurs d'échelle (Scaling Factors)

Les auteurs ont testé trois approches pour corriger les données compressibles :

1. Rugosité équivalente incompressible ($k_{si}$) : En cherchant une étude incompressible avec une géométrie et un $\delta/k$ similaires.

   • Résultat : Succès limité. Cela ne fonctionne que pour certains cas (ex: cubes de Modesti et al.) mais échoue pour d'autres, indiquant que la simple correspondance géométrique est insuffisante.

2. Rugosité échelonnée par la viscosité ($k^* = k / \nu^*_\infty$) : Utilisation du rapport de viscosité cinématique entre l'écoulement libre et la paroi ($\nu_\infty / \nu_w$).

   • Résultat : Cela permet de faire coïncider les données expérimentales actuelles avec la loi de Nikuradse, mais l'échelle échoue à prédire correctement les données de la littérature antérieure.

3. Correction de la traînée ($\sqrt{1/F_c} \Delta U^+$) : Application d'un facteur correctif à la fonction de rugosité $\Delta U^+$ lui-même, basé sur le facteur de correction de la contrainte de cisaillement $F_c$ (dérivé de la transformation Van Driest II pour les parois lisses), qui dépend du rapport de température $T_\infty / T_w$.

   • Résultat : C'est la méthode la plus performante. Elle permet de faire coïncister de manière cohérente l'ensemble des données expérimentales actuelles (P60, P24) et une grande partie des données de la littérature (diverses rugosités, régimes supersoniques à hypersoniques) sur la loi de Nikuradse, avec une erreur inférieure à 10 %.

D. Observation physique de la traînée d'onde

Les images Schlieren instantanées et les champs de divergence de vitesse montrent clairement la formation d'une succession d'ondes de choc (ondes bow) devant les éléments de rugosité en régime supersonique. Ces ondes s'étendent jusqu'à l'écoulement libre, confirmant l'existence d'une traînée d'onde significative que les transformations de vitesse actuelles (développées pour des parois lisses) ne prennent pas en compte.

4. Signification et Perspectives

• Apport principal : L'article démontre que le cadre actuel de caractérisation de la traînée (basé sur des transformations pour parois lisses appliquées à des parois rugueuses) est insuffisant en régime compressible. Cependant, une correction empirique basée sur le rapport de température ($\sqrt{1/F_c}$) offre une solution pratique et robuste pour estimer la traînée sur une large gamme de conditions.
• Limitations : Toutes les corrections proposées sont empiriques. La méthode repose sur l'hypothèse de paroi adiabatique et sur des estimations de contrainte de cisaillement, car des mesures directes de $T_w$ et du WSS sur des surfaces rugueuses en soufflerie supersonique sont extrêmement difficiles à réaliser.
• Travaux futurs :
  • Nécessité de mesures directes de la température de paroi et de la contrainte de cisaillement sur des surfaces rugueuses.
  • Développement de transformations de vitesse spécifiques aux parois rugueuses qui intègrent explicitement les propriétés de la rugosité et les effets de la traînée d'onde, plutôt que de se contenter de corriger les données a posteriori.
  • Établissement d'une "carte" complète de la traînée en fonction des propriétés de la rugosité, du nombre de Mach et des conditions de paroi.

En conclusion, cette étude fournit une base expérimentale cruciale pour la modélisation de la traînée sur les véhicules hypersoniques et supersoniques, tout en soulignant la nécessité de dépasser les modèles de parois lisses pour comprendre pleinement la physique des couches limites rugueuses compressibles.