Scale-resolving simulations and data-driven modal analysis of turbulent transonic buffet cells on infinite swept wings

Cette étude démontre que le buffet transsonique sur des ailes à flèche infinie résulte de la superposition d'un mouvement bidimensionnel du choc et d'instabilités tridimensionnelles de séparation, ces dernières n'émergeant de manière dominante que lorsque la séparation moyenne à l'endroit du choc est suffisamment prononcée.

L'essentiel

🌤️ Le Problème : Le "Tremblement" de l'Avion

Imaginez que vous êtes dans un avion qui vole très vite, presque à la vitesse du son. À un certain moment, l'air qui passe sur l'aile devient turbulent et crée une onde de choc (comme un petit mur d'air invisible). Parfois, ce mur d'air commence à osciller d'avant en arrière, comme une corde de guitare qu'on pince.

Ce phénomène s'appelle le "buffet". C'est comme si l'avion avait le hoquet ou tremblait de manière incontrôlable. Cela peut être dangereux pour la structure de l'avion et très inconfortable pour les passagers.

Les scientifiques savaient déjà deux choses :

1. Ce tremblement a une composante 2D (comme une vague qui va et vient d'avant en arrière sur toute l'aile en même temps).
2. Mais il y a aussi une composante 3D (comme des bulles ou des cellules qui se forment et voyagent le long de l'aile, de gauche à droite).

Le mystère, c'est : Comment ces deux mouvements interagissent-ils ? Et surtout, comment la forme de l'aile (son angle de balayage, comme une aile d'oiseau en vol) change-t-elle la donne ?

🔬 L'Expérience : Un Laboratoire Numérique Géant

Pour comprendre cela, les chercheurs (David Lusher et Andrea Sansica) ont construit un "laboratoire virtuel" immense.

• Le Défi : Simuler l'air autour d'une aile infinie (sans bout) est très difficile. Les ordinateurs habituels sont trop petits pour voir les détails fins de la turbulence. C'est comme essayer de voir les grains de sable sur une plage avec une paire de jumelles de mauvaise qualité.
• La Solution : Ils ont utilisé des superordinateurs ultra-puissants (comme le Fugaku au Japon) pour faire des simulations d'une précision incroyable. Ils ont créé une "aile infinie" virtuelle avec une largeur équivalente à 3 fois la longueur de l'aile. C'est comme si on regardait une section d'aile d'avion de taille réelle, mais dans un monde où l'aile ne finit jamais.

Ils ont testé deux scénarios principaux :

1. L'aile "calme" (Angle d'attaque faible) : L'air colle bien à l'aile, il y a peu de turbulence.
2. L'aile "agitée" (Angle d'attaque fort) : L'air commence à se décoller de l'aile, créant de gros tourbillons.

🎢 Les Découvertes : La Danse des Bulles

Voici ce qu'ils ont découvert, avec des analogies simples :

1. Quand l'air est calme (Angle faible)

Même si l'aile est large, le tremblement reste très 2D. C'est comme une vague qui monte et descend uniformément sur toute la surface.

• L'analogie : Imaginez une corde de guitare que vous secouez. Tout le mouvement se fait d'avant en arrière. Il n'y a pas de vagues qui voyagent le long de la corde.
• Résultat : Même si on incline l'aile (on ajoute du "balayage"), le tremblement reste simple et prévisible.

2. Quand l'air est agité (Angle fort)

C'est là que la magie opère. Dès que l'air commence à se décoller de l'aile (comme quand un avion décolle trop vite), des cellules de buffet apparaissent.

• L'analogie : Imaginez maintenant que vous secouez la corde de guitare, mais qu'en plus, des vagues commencent à voyager le long de la corde, de gauche à droite. Ces vagues sont les "cellules de buffet".
• Le rôle de l'inclinaison (Sweep) : Quand l'aile est inclinée (comme sur un avion de ligne moderne), ces vagues ne font pas que voyager ; elles accélèrent et changent de rythme. Plus l'aile est inclinée, plus les vagues voyagent vite et changent de fréquence.

3. Le Secret : La Séparation de l'Air

La découverte clé est que ces vagues 3D ne peuvent pas exister si l'air reste collé à l'aile.

• L'analogie : C'est comme si vous aviez besoin d'une "piscine de boue" (la séparation de l'air) pour que des bulles puissent se former et voyager. Sans cette boue, vous n'avez que de l'eau calme.
• Conclusion : Pour que le tremblement 3D (le plus dangereux et complexe) apparaisse, il faut que l'air se décolle de l'aile au niveau du mur de choc.

🚀 Pourquoi c'est important ?

Avant cette étude, les scientifiques pensaient que le tremblement 2D (la corde de guitare) et le tremblement 3D (les vagues voyageuses) étaient des choses totalement séparées ou qu'il fallait des ailes très complexes pour voir les vagues.

Cette recherche montre que :

1. C'est une seule et même danse : Le tremblement de l'avion est un mélange de ces deux mouvements.
2. L'inclinaison est la clé : Sur les vrais avions (qui ont des ailes inclinées), le mouvement 3D devient dominant. C'est pour cela qu'on ne voit presque jamais le simple mouvement 2D sur les vrais avions : il est "noyé" par le mouvement 3D plus rapide et plus énergique.
3. Prédire le danger : En comprenant comment ces vagues voyagent et accélèrent avec l'inclinaison, les ingénieurs peuvent mieux concevoir les ailes pour éviter ces tremblements dangereux.

En résumé

Imaginez que l'air autour d'une aile est comme une foule de gens.

• Parfois, la foule bouge tous ensemble d'avant en arrière (le mouvement 2D).
• Mais si la foule devient trop dense et commence à se bousculer (la séparation), des groupes de personnes commencent à courir le long de la rangée (le mouvement 3D).
• Si la rangée est inclinée, ces groupes de coureurs accélèrent et deviennent le mouvement principal, effaçant le mouvement lent de la foule entière.

Cette étude nous donne la carte exacte de comment ces "groupes de coureurs" se comportent, ce qui aidera à construire des avions plus sûrs et plus confortables.

Résumé technique

1. Problématique et Contexte

Le buffet transsonique turbulent est un phénomène d'interaction onde de choc/couche limite (SBLI) auto-entretenue, critique pour l'intégrité structurelle et les enveloppes de vol des aéronefs. Ce phénomène se caractérise par deux instabilités principales :

• Une instabilité bidimensionnelle (2D) : Oscillation de l'onde de choc principale dans la direction de la corde (nombre de Strouhal $St \sim 0.05 - 0.1$).
• Une instabilité tridimensionnelle (3D) : Formation de "cellules de buffet" (bulles de séparation périodiques le long de l'envergure) se propageant dans la direction de l'envergure ($St \sim 0.2 - 0.4$).

Jusqu'à présent, les simulations haute fidélité (LES/DNS) de ces phénomènes étaient limitées à des géométries à faible rapport d'aspect (AR < 0.25) en raison du coût computationnel. Ces dimensions insuffisaient pour capturer les structures 3D complètes (longueur d'onde $\lambda_z \approx 1-1.5c$). De plus, la relation entre les modes 2D et 3D, ainsi que le rôle de la séparation de l'écoulement moyen dans l'émergence des cellules 3D, restaient mal compris, en particulier sur des ailes balayées.

2. Méthodologie

Les auteurs ont mené une étude numérique de haute fidélité utilisant des simulations implicites de grandes tourbillons (ILES) sur des ailes infinies (périodiques dans la direction de l'envergure).

• Configuration : Ailes infinies basées sur le profil NASA-CRM (65% de la demi-envergure), avec un rapport d'aspect élevé AR = 3 (largeur d'envergure $L_z = 3c$).
• Paramètres de vol : Nombre de Mach transsonique $M_\infty = 0.72$ et nombre de Reynolds $Re_c = 0.5 \times 10^6$.
• Variables étudiées :
  • Angles d'attaque ($\alpha$) : Deux cas comparés : $\alpha = 5^\circ$ (séparation moyenne minimale) et $\alpha = 6^\circ$ (séparation moyenne importante au niveau du choc).
  • Angles de balayage ($\lambda$) : Une gamme de $0^\circ$ à $35^\circ$ (angles typiques des avions commerciaux).
• Ressources computationnelles : Utilisation de supercalculateurs (GPU NVIDIA V100 et CPU Fujitsu A64FX) pour des simulations massives impliquant environ 8 milliards de points de maillage par cas.
• Analyse des données : Application de la Décomposition Orthogonale Propre Spectrale (SPOD) sur les champs de pression et de vitesse pour extraire les modes cohérents dominants et leurs énergies spectrales.

3. Contributions Clés

1. Échelle sans précédent : Première application d'ILES haute fidélité sur des ailes infinies avec un rapport d'aspect suffisant (AR=3) pour résoudre pleinement les cellules de buffet 3D, au-delà des limites des études précédentes (souvent limitées à AR < 0.25).
2. Lien entre modes stationnaires et voyageurs : Établissement d'un lien direct entre le mode de séparation 3D quasi-stationnaire à basse fréquence observé sur les ailes non balayées et le mode de cellule de buffet voyageur à fréquence intermédiaire observé sur les ailes balayées.
3. Rôle de la séparation moyenne : Identification de la séparation moyenne de l'écoulement au niveau du choc comme condition nécessaire à l'émergence de la dynamique 3D dominante.
4. Validation de modèles d'échelle : Vérification de la loi d'échelle de la vitesse de convection des cellules de buffet en fonction de l'angle de balayage.

4. Résultats Principaux

A. Cas à faible séparation ($\alpha = 5^\circ$)

• Comportement quasi-2D : Même avec un rapport d'aspect élevé et un balayage jusqu'à $25^\circ$, l'écoulement reste essentiellement bidimensionnel au niveau du choc.
• Absence de cellules 3D majeures : Aucune cellule de buffet 3D macroscopique n'est observée le long du choc. De faibles perturbations 3D intermittentes apparaissent uniquement près du bord de fuite lors des phases de forte séparation instantanée, mais elles ne sont pas couplées à l'oscillation du choc.
• Dominance du mode 2D : Le mode d'oscillation du choc ($St \approx 0.08$) domine totalement, indépendamment de l'angle de balayage.

B. Cas à forte séparation ($\alpha = 6^\circ$)

• Émergence des cellules 3D : Une séparation moyenne importante au niveau du choc déclenche l'apparition de cellules de buffet 3D prononcées avec une longueur d'onde caractéristique $\lambda_z \approx 1 - 1.5c$.
• Effet du balayage sur la fréquence :
  • À $\lambda = 0^\circ$, un mode 3D quasi-stationnaire à très basse fréquence ($St \approx 0.02$) est observé.
  • Dès qu'un balayage est imposé, ce mode se transforme en un mode voyageur dans la direction de l'envergure.
  • La fréquence de ce mode 3D augmente de manière monotone avec l'angle de balayage, passant de $St \approx 0.06$ à $St \approx 0.35$ pour $\lambda = 35^\circ$.
• Indépendance de la longueur d'onde : La longueur d'onde des cellules 3D reste constante ($\approx 1.5c$) quelle que soit l'angle de balayage, étant dictée par le rapport d'aspect du domaine computationnel.
• Compétition des modes : À des angles de balayage élevés ($\lambda \ge 25^\circ$), l'énergie du mode 3D devient comparable, voire supérieure, à celle du mode 2D. Le mode 2D reste insensible au balayage, tandis que le mode 3D gagne en énergie et en fréquence.

C. Vitesse de convection

La vitesse de convection des perturbations le long du choc suit une loi d'échelle en fonction de l'angle de balayage : $V_c \propto \tan(\lambda)$. Les résultats de cette étude ($V_c \approx 0.78 \tan(\lambda)$) sont en excellent accord avec les modèles théoriques et les simulations URANS précédentes, malgré les différences de géométrie et de nombre de Reynolds.

5. Signification et Implications

• Compréhension fondamentale : L'étude démontre que le buffet transsonique sur des ailes réelles (balayées) résulte de la superposition d'un mouvement de choc 2D et d'instabilités 3D pilotées par la séparation.
• Explication de l'absence de mode 2D sur les ailes finies : Sur les configurations réalistes (ailes finies, balayage), le mode 3D devient dominant à haute fréquence et énergie, masquant souvent le mode 2D pur observé sur les profils 2D ou les ailes non balayées étroites. Cela explique pourquoi les études expérimentales sur des avions complets ne voient souvent que le mode 3D.
• Condition critique : La séparation moyenne de l'écoulement au niveau du choc est identifiée comme la condition sine qua non pour l'apparition de la dynamique 3D dominante.
• Validation des méthodes : La capacité de la SPOD appliquée à des données ILES à distinguer les modes 2D et 3D, même lorsqu'ils se chevauchent spectralement, valide l'approche pour l'analyse de flux complexes.

En conclusion, ce travail fournit une référence haute fidélité pour la compréhension du buffet 3D, reliant les observations expérimentales et les simulations basse fidélité à la physique fondamentale résolue, et ouvre la voie à des stratégies de contrôle plus ciblées sur les mécanismes de séparation.