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

Diese Studie zeigt mittels impliziter Large-Eddy-Simulationen und modaler Analysen an unendlich gepfeilten Flügeln, dass transsonisches Buffeting aus einer Überlagerung von zweidimensionaler Stoßbewegung und dreidimensionalen, durch mittlere Strömungsablösung getriebenen Instabilitäten besteht, wobei die Entstehung dominanter 3D-Buffet-Zellen eine signifikante Ablösung an der Stoßstelle voraussetzt.

Das Wesentliche

Das große Problem: Der "Flügel-Zittern" (Transonic Buffet)

Stellen Sie sich vor, ein Flugzeug fliegt sehr schnell, fast so schnell wie der Schall. In diesem Bereich passiert etwas Seltsames: Die Luft strömt über den Flügel, wird an einer Stelle plötzlich langsamer (ein "Schock"), und dann wieder schneller.

Das Problem ist, dass dieser Schock nicht ruhig bleibt. Er wackelt hin und her, wie ein nervöses Herz. Das nennt man Buffet.

• Die Folge: Das ganze Flugzeug fängt an zu vibrieren. Das ist unangenehm für die Passagiere und kann dem Flugzeug auf Dauer schaden.

Früher dachten Forscher, dieses Wackeln sei nur eine einfache, gerade Linie von vorne nach hinten (wie ein Seil, das auf und ab schwingt). Aber neuere Studien zeigten: Es ist viel chaotischer! Es gibt auch 3D-Muster, die sich quer über den Flügel ausbreiten, wie Wellen auf einem See oder wie eine Welle, die durch ein Stadion läuft. Diese nennt man "Buffet-Zellen".

Das Experiment: Ein unendlicher Flügel im Supercomputer

Um dieses Chaos zu verstehen, haben die Forscher (David Lusher und Andrea Sansica von JAXA) einen riesigen digitalen Versuchsaufbau gebaut.

1. Der "Unendliche" Flügel: In der echten Welt sind Flügel endlich (sie haben eine Spitze). Das macht es schwer zu messen, was genau passiert. Die Forscher haben also einen unendlichen Flügel simuliert. Stellen Sie sich einen Flügel vor, der links und rechts ins Unendliche reicht, ohne Enden. So können sie sich nur auf das Wesentliche konzentrieren.
2. Der Supercomputer: Um die winzigen Wirbel in der Luft zu sehen, brauchten sie einen extrem leistungsfähigen Computer. Sie haben so viele Datenpunkte berechnet, dass es fast unmöglich ist, das zu verstehen (etwa 8 Milliarden Gitterpunkte!). Es war, als würde man versuchen, jeden einzelnen Wassertropfen in einem Ozean zu zählen, um zu verstehen, wie eine Welle entsteht.
3. Die Drehung (Sweep): Sie haben den Flügel in verschiedenen Winkeln gedreht (von gerade bis stark geneigt, wie bei modernen Passagierflugzeugen).

Was haben sie herausgefunden? (Die Entdeckungen)

Die Forscher haben zwei verschiedene Szenarien getestet, ähnlich wie bei einem Auto, das entweder sanft oder hart bremst:

Szenario 1: Der sanfte Bremsvorgang (Wenig Ablösung)

• Die Situation: Der Flügel ist nur leicht geneigt. Die Luft strömt relativ ruhig über die Oberfläche.
• Das Ergebnis: Hier ist das Wackeln sehr vorhersehbar. Es ist fast wie ein Metronom. Der Schock wackelt nur auf und ab (vor und zurück), aber er bleibt gerade. Es gibt keine wilden 3D-Muster, die quer über den Flügel laufen. Selbst wenn der Flügel geneigt ist, bleibt das Wackeln "zweidimensional".
• Analogie: Stellen Sie sich einen geraden Stock vor, der auf und ab wippt. Egal, wie Sie ihn halten, er wippt nur auf und ab.

Szenario 2: Der harte Bremsvorgang (Starke Ablösung)

• Die Situation: Der Flügel wird steiler geneigt. Die Luft kann der Oberfläche nicht mehr folgen und löst sich ab (wie Wasser, das von einer rutschigen Treppe fällt).
• Das Ergebnis: Hier wird es wild! Es entstehen die berühmten Buffet-Zellen.
  • Das Muster: Stellen Sie sich vor, Sie laufen über einen Feldweg. Plötzlich bilden sich Löcher im Boden, die sich quer über den Weg ausbreiten. Diese Löcher wandern von links nach rechts über den Flügel.
  • Der Einfluss der Drehung: Je mehr der Flügel geneigt ist (Sweep), desto schneller wandern diese "Löcher" quer über den Flügel.
  • Die Frequenz: Bei wenig Drehung wandern sie langsam. Bei viel Drehung (wie bei echten Flugzeugen) wandern sie sehr schnell und werden so dominant, dass sie das einfache "Auf-und-Ab"-Wackeln fast komplett übertönen.

Die große Erkenntnis: Warum wir das vorher nicht sahen

Früher haben Forscher oft nur schmale Flügel simuliert (wie ein schmales Brett). Das war wie ein Fenster, durch das man nur einen kleinen Ausschnitt sieht. In diesen schmalen Fenstern konnten die großen 3D-Muster (die Zellen) nicht entstehen, weil ihnen der Platz fehlte.

Die Forscher haben nun einen riesigen Flügel (drei Mal so breit wie lang) simuliert. Erst auf dieser großen Fläche konnten die "Buffet-Zellen" ihre volle Kraft entfalten.

Die wichtigste Botschaft:
Das chaotische 3D-Wackeln entsteht nur, wenn die Luft an der Stelle, wo der Schock ist, auch wirklich abreißt (sich löst).

• Wenn die Luft klebt (wenig Ablösung) -> Nur einfaches Wackeln (2D).
• Wenn die Luft abreißt (viel Ablösung) -> Wildes 3D-Chaos mit wandernden Zellen.

Zusammenfassung für den Alltag

Stellen Sie sich vor, Sie halten ein Taschentuch im Wind.

• Wenn der Wind sanft weht und das Tuch straff ist, wackelt es nur leicht auf und ab (2D-Modus).
• Wenn der Wind stark ist und das Tuch schlaff wird, fängt es an, sich in wilden, wandernden Wellen zu bewegen, die sich über die ganze Breite des Tuches ausbreiten (3D-Modus).

Die Studie zeigt uns, dass bei modernen Flugzeugen (die stark geneigte Flügel haben) diese "schlaffen" 3D-Wellen das Hauptproblem sind. Sie sind so stark, dass sie das einfache Wackeln fast unsichtbar machen. Das hilft Ingenieuren, Flugzeuge zu bauen, die sicherer sind und weniger vibrieren, indem sie genau wissen, wann und warum diese "Wellen" entstehen.

Kurz gesagt: Um das Wackeln von Flugzeugen zu verstehen, muss man nicht nur auf das Auf-und-Ab schauen, sondern auch auf die wilden Querwellen, die entstehen, wenn die Luft nicht mehr sauber über den Flügel strömt. Und je mehr der Flügel geneigt ist, desto schneller und stärker werden diese Querwellen.

Technische Zusammenfassung

Titel:

Auflösungssimulierungen und datengesteuerte Modalanalyse von turbulenten transonischen Buffeting-Zellen auf unendlich gepfeilten Flügeln

1. Problemstellung und Hintergrund

Transonisches Buffeting ist eine Form der selbstunterhaltenen Wechselwirkung zwischen Stoßwelle und Grenzschicht (SBLI), die bei hohen Unterschall-Mach-Zahlen und moderaten bis hohen Anstellwinkeln auftritt. Es manifestiert sich durch zwei Hauptinstabilitäten:

1. Zweidimensionale (2D) Moden: Niederfrequente, chordale Oszillationen der Hauptstoßwelle (Strouhal-Zahl $St \approx 0.05 - 0.1$).
2. Dreidimensionale (3D) Moden: Spannungsmodulierte Ablösungs-/Wiederanlegezellen, sogenannte "Buffeting-Zellen" (Buffet cells), die typischerweise bei höheren Frequenzen ($St \approx 0.2 - 0.4$) auftreten und sich spannungsweise ausbreiten.

Bisherige hochauflösende Simulationen (LES/DNS) waren aufgrund des hohen Rechenaufwands auf Flügel mit sehr geringer Spannweite (Seitenverhältnis $AR = L_z/c \approx 0.05 - 0.25$) beschränkt. Diese Dimensionen reichten nicht aus, um die beobachteten 3D-Buffeting-Zellen mit Wellenlängen von $\lambda_z \approx 1 - 1.5$ Sehnenlängen ($c$) vollständig aufzulösen. Zudem ist unklar, unter welchen Bedingungen diese 2D- und 3D-Instabilitäten koexistieren und wie die Flügelpfeilung ($\lambda$) diese Dynamik beeinflusst.

2. Methodik

Die Autoren führten eine Reihe von Implicit Large-Eddy Simulations (ILES) durch, die auf dem Open-Source-Code OpenSBLI basieren.

• Konfiguration: Unendlich lange, gepfeilte Flügel (spannungsperiodische Randbedingungen) basierend auf dem NASA-CRM-Profil (65% Spannweitenstation).
• Geometrie: Seitenverhältnis bis zu $AR = 3$($L_z = 3c$) und Pfeilwinkel $\lambda$ im Bereich von $0^\circ$ bis $35^\circ$.
• Strömungsbedingungen: Transonische Mach-Zahl $M_\infty = 0.72$ (basierend auf der Sehne) und Reynolds-Zahl $Re_c = 0.5 \times 10^6$.
• Untersuchte Fälle: Zwei Anstellwinkel wurden verglichen, um den Einfluss der mittleren Strömungsablösung zu untersuchen:
  • $\alpha = 5^\circ$: Minimale mittlere Ablösung an der Stoßwelle.
  • $\alpha = 6^\circ$: Große mittlere Ablösung an der Stoßwelle.
• Rechenleistung: Die Simulationen nutzten ca. $8 \times 10^9$ Gitterpunkte pro Fall und erforderten ca. 648.000 GPU-Stunden (NVIDIA V100) bzw. $3.8 \times 10^8$ CPU-Kernstunden.
• Datenanalyse: Die gewonnenen hochfideligen Zeitreihen wurden mittels Spektraler Eigenervektor-Zerlegung (SPOD) analysiert, um dominante Frequenzen und kohärente Strukturmuster zu extrahieren.

3. Wichtige Beiträge und Ergebnisse

A. Einfluss der mittleren Strömungsablösung ($\alpha = 5^\circ$ vs. $6^\circ$)

• Fall $\alpha = 5^\circ$ (Minimale Ablösung): Die Strömung verhält sich im Wesentlichen zweidimensional. Die Stoßwelle oszilliert chordal ($St \approx 0.08$) ohne signifikante spannungsweise Modulation. Zwar treten intermittierende Ablösungsblasen nahe dem Hinterkantenbereich auf, diese interagieren jedoch schwach mit der Stoßwelle und bilden keine makroskopischen Buffeting-Zellen. Die 3D-Energie ist hier um drei Größenordnungen geringer als die der 2D-Stoßmoden.
• Fall $\alpha = 6^\circ$ (Große Ablösung): Hier entstehen ausgeprägte 3D-Buffeting-Zellen mit einer charakteristischen Wellenlänge von $\lambda_z \approx 1 - 1.5c$. Die mittlere Strömungsablösung an der Stoßwelle wird als notwendige Bedingung identifiziert, damit dominante 3D-Buffeting-Dynamiken entstehen können.

B. Einfluss der Pfeilung ($\lambda$)

• 2D-Stoßmode: Die Frequenz der chordalen Stoßoszillation ($St \approx 0.08$) ist weitgehend unempfindlich gegenüber der Pfeilung.
• 3D-Buffeting-Mode:
  • Bei unswepten Flügeln ($\lambda = 0^\circ$) wurde ein quasi-stationärer, niederfrequenter 3D-Ablösungsmodus ($St \approx 0.02$) beobachtet.
  • Mit zunehmender Pfeilung wandelt sich dieser Modus in eine spannungsweise wandernde Welle um.
  • Die Frequenz dieser 3D-Mode steigt monoton mit dem Pfeilwinkel an (von $St \approx 0.06$ bei $\lambda=5^\circ$ bis $St \approx 0.35$ bei $\lambda=35^\circ$).
  • Die Wellenlänge $\lambda_z$ bleibt dabei unabhängig vom Pfeilwinkel konstant und wird durch das Seitenverhältnis des Rechengebiets bestimmt.
• Dominanz: Bei hohen Pfeilwinkeln (typisch für kommerzielle Flugzeuge, $\lambda \approx 35^\circ$) übersteigt die Energie der 3D-Buffeting-Mode die der 2D-Stoßmode. Dies erklärt, warum auf realen, endlichen Flügeln oft nur die 3D-Phänomene beobachtet werden, während die 2D-Stoßoszillationen verschwinden.

C. Konvektionsgeschwindigkeit

Die Konvektionsgeschwindigkeit der Buffeting-Zellen entlang der Stoßfront folgt einer Skalierungsgesetz-Beziehung in Abhängigkeit vom Pfeilwinkel:
$$V_c \approx 0.78 \cdot \tan(\lambda) \cdot U_{\infty}$$
Dies steht in guter Übereinstimmung mit früheren Studien (z.B. Plante et al.) und experimentellen Daten, trotz signifikanter Unterschiede in Reynolds-Zahl und Geometrie.

4. Bedeutung und Fazit

Diese Studie liefert den ersten hochfideligen Nachweis (mittels ILES und SPOD), der die Verbindung zwischen dem auf unswepten Flügeln beobachteten quasi-stationären 3D-Modus und dem auf gepfeilten Flügeln dominierenden wandernden Buffeting-Modus herstellt.

Kernaussagen:

1. Transonisches Buffeting auf unendlichen Flügeln resultiert aus der Superposition einer 2D-Stoßbewegung und einer ablösungsgetriebenen 3D-Instabilität.
2. Eine signifikante mittlere Strömungsablösung an der Stoßwelle ist die Voraussetzung für das Auftreten dominanter 3D-Buffeting-Zellen.
3. Die Pfeilung wandelt den 3D-Modus von einer stehenden in eine wandernde Welle um und verschiebt dessen Frequenz in einen höheren Bereich, wodurch die 3D-Dynamik bei typischen Flugzeugpfeilwinkeln dominiert.
4. Die Ergebnisse widerlegen die Annahme, dass 2D-Stoßoszillationen auf endlichen Flügeln vollständig verschwinden; sie werden vielmehr durch die stärkere 3D-Instabilität überlagert und maskiert.

Die Arbeit unterstreicht die Notwendigkeit von hochauflösenden Simulationen mit ausreichender Spannweite, um die komplexe Wechselwirkung zwischen Stoßwellen und 3D-Ablösungsstrukturen korrekt zu erfassen, und liefert wichtige Erkenntnisse für die Vorhersage und Kontrolle von Buffeting in der Luftfahrt.