Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Diese Studie präsentiert die erste direkte numerische Simulation verteilter Rauheit an einem stumpfen Zylinder bei Mach 6 und zeigt auf, dass die Anordnung der Rauheit den Übergangsmechanismus bestimmt, indem sie entweder sinuose Streifenmoden oder zweidimensionale T-S-Wellen fördert, wobei letztere durch eine interne akustische Rückkopplungsschleife, die eine externe Erregung überflüssig macht, einzigartig aufrechterhalten werden.

Das Wesentliche

Stellen Sie sich ein Raumschiff vor, das mit dem Sechsfachen der Schallgeschwindigkeit durch die Atmosphäre rast. Um die intensive Hitze zu überleben, ist seine Oberfläche mit einem speziellen Material bedeckt, das langsam verbrennt (abläuft), um das Schiff zu schützen. Doch während dieses Material verbrennt, hinterlässt es keine perfekt glatte Oberfläche, sondern eine unebene, raue Textur, die fast wie Sandpapier wirkt.

Dieses Paper ist eine hochgeschwindigkeitsbasierte Computersimulation, die eine einfache, aber entscheidende Frage stellt: Wie verwandeln diese winzigen Unebenheiten auf der Oberfläche eine glatte, geordnete Luftströmung in eine chaotische, turbulente Luftströmung?

Hier ist die Aufschlüsselung ihrer Ergebnisse unter Verwendung alltäglicher Analogien:

1. Der Aufbau: Der „Sandpapier“-Zylinder

Die Forscher bauten ein digitales Modell eines stumpfen Zylinders (wie die Nase einer Rakete), der mit Mach 6 fliegt. Anstatt einer glatten Oberfläche überzogen sie ihn mit winzigen, künstlichen „Höckern“ (Rauheit), um die sandähnliche Textur nachzuahmen, die durch das Verbrennen des Materials hinterlassen wird.

Sie testeten drei verschiedene Arten, diese Höcker anzuordnen:

• Ausgerichtet: Wie Soldaten, die in perfekten Reihen und Spalten stehen.
• Versetzt: Wie eine Ziegelwand, bei der die Höcker in einer Reihe gegenüber der Reihe dahinter versetzt sind.
• Zufällig: Wie Kieselsteine, die ohne Muster auf einem Gehweg verstreut liegen.

2. Die alte Theorie vs. die neue Entdeckung

Lange Zeit glaubten Wissenschaftler, dass der Übergang zur Turbulenz durch einen „langsamen Aufbau“ von Energie verursacht wurde, ähnlich wie eine Schaukel an Höhe gewinnt, wenn man sie über die Zeit genau richtig anschubst. Dies wird als „transientes Wachstum“ bezeichnet.

Das Ergebnis des Papers:
Die Simulation zeigte, dass diese Theorie des „langsenden Energieaufbaus“ hier nicht wirklich erklärt, was passiert. Die Höcker auf der Oberfläche haben die Energie nicht nur langsam verstärkt, sondern wirkten wie Destabilisatoren. Sie nahmen die Luftströmung und machten sie sofort instabil, indem sie sie in eine spezifische Art von Welle verwandelten, die sehr schnell wächst.

Man kann es sich so vorstellen: Die alte Theorie dachte, die Höcker würden einen Domino erst sanft anstoßen, damit er umfällt. Die neue Entdeckung zeigt, dass die Höcker den Domino tatsächlich treten, was dazu führt, dass er sofort gegen den nächsten prallt.

3. Die zwei Arten von „Wellen“

Je nachdem, wie die Höcker angeordnet waren, reagierte die Luftströmung auf zwei verschiedene Arten:

• Die „Schlange“ (Sinuous Mode): Wenn die Höcker ausgerichtet waren (perfekte Reihen), begann die Luftströmung, sich seitlich hin und her zu winden wie eine Schlange. Dies ist ein sehr spezifisches, organisiertes Wackeln.
• Die „Flache Welle“ (Tollmien-Schlichting oder T-S Wellen): Wenn die Höcker versetzt oder zufällig angeordnet waren, begann die Luftströmung, in einem flachen, 2D-Wellenmuster auf und ab zu kräuseln. Dies ist eine klassische Wellenart, die normalerweise in viel langsamerer, niedriggeschwindigkeits-Luft vorkommt, was in dieser Hochgeschwindigkeitsumgebung überraschend war.

Die zentrale Erkenntnis: Die Anordnung der Höcker bestimmte, welchen „Tanz“ die Luft aufführte. Der „Schlangentanz“ geschah bei ausgerichteten Höckern, während der „flache Wellentanz“ bei den anderen auftrat.

4. Das „Haarnadel“-Finale

Sobald diese Wellen stark genug waren, lösten sie die finale Phase des Absturzes aus. Die stetigen „Streifen“ (Streaks) der Luft, die durch die Höcker erzeugt wurden (welche wie lange, unsichtbare Bänder langsamer Luft sind), verdrehten und rissen plötzlich in Haarnadel-Wirbel (hairpin vortices) auf.

Stellen Sie sich ein Gummiband vor, das straff gezogen ist. Plötzlich verdreht es sich und bildet eine Schlaufe, die wie eine Haarnadel aussieht. Diese Schlaufen sind die Geburtsstunde der Turbulenz. Sobald diese Haarnadeln entstehen, bricht die glatte Luft vollständig in Chaos zusammen, und die Hitze auf der Oberfläche des Raumschiffs steigt dramatisch an.

5. Die Überraschung der „Echo-Kammer“

Eine der faszinierendsten Entdeckungen war, wie die Turbulenz in den versetzten und zufälligen Fällen überhaupt erst entstand.

Normalerweise denken Wissenschaftler, dass man einen externen „Anstoß“ benötigt (wie eine Windböe oder eine Vibration), um diese Wellen zu starten. Aber die Simulation zeigte etwas Selbsterhaltendes:

1. Die Turbulenz beginnt an einem Punkt am Zylinder.
2. Da die Luft hinter der Stoßwelle langsamer als der Schall ist (subsonisch), wandert der Lärm dieser Turbulenz wie ein Echo rückwärts stromaufwärts.
3. Dieser „Echo“ trifft auf den glatten Teil der Oberfläche vor der Turbulenz und regt die Luft dort an, wodurch neue Turbulenz entsteht.
4. Dies erzeugt eine Feedback-Schleife: Turbulenz erzeugt Lärm, der Lärm wandert zurück, und der Lärm erzeugt neue Turbulenz.

Es ist, als ob ein Mikrofon sein eigenes Lautsprecherausgangssignal aufnimmt und eine kreischende Rückkopplung erzeugt – aber in diesem Fall ist das „Kreischen“ die Luft, die turbulent wird.

Zusammenfassung

Dieses Paper nutte einen Supercomputer, um den Luftstrom über einem unebenen, Hochgeschwindigkeits-Zylinder zu beobachten. Sie fanden heraus, dass:

• Das Muster der Höcker entscheidet, wie genau die Luft turbulent wird.
• Die alte Idee des „langsamen Energieaufbaus“ nicht der Hauptschuldige ist; stattdessen destabilisieren die Höcker direkt bestimmte Wellen.
• Diese Wellen wachsen, bis sie sich in „Haarnadel“-Formen verdrehen und die Luft chaotisch werden lassen.
• In einigen Fällen erzeugt die Turbulenz ihr eigenes „Echo“, das rückwärts wandert, um den Prozess von vorn zu beginnen, ohne dass Hilfe von außen nötig ist.

Dies hilft Ingenieuren zu verstehen, dass die winzigen, zufälligen Unebenheiten, die durch die brennenden Hitzeschilde hinterlassen werden, nicht bloß kleine Unvollkommenheiten sind, sondern die primären Architekten dafür, wie und wann die Oberfläche eines Raumschiffs gefährlich heiß wird.

Technische Zusammenfassung

Technische Zusammenfassung: Verteilte rauheitsinduzierte Transition an einem stumpfen Körper bei Mach 6

Problemstellung
Die Vorhersage des laminaren-turbulenten Übergangs (LTT) über Hochgeschwindigkeits-Geometrien mit stumpfen Formen bleibt eine kritische Herausforderung für das Ingenieursdesign, insbesondere für atmosphärische Eintrittsfahrzeuge, bei denen die Oberflächenablation eine verteilte Rauheit erzeugt. Während empirische Korrelationen existieren (z. B. Hollis 2017, 2019, 2025), mangelt es ihnen an einer klaren physikalischen Grundlage. Frühere theoretische Bemühungen, wie die Reshotko & Tumin (2004) Korrelation, stützen sich auf Mechanismen des transienten Wachstums, stoßen jedoch an Grenzen: Sie setzen eine lineare Skalierung der Störungen mit der Rauheitshöhe voraus (gültig nur für kleine $k/\delta$), benötigen Anpassungskonstanten für experimentelle Daten und haben Schwierigkeiten, den Übergang in Fällen zu erklären, in denen die Rauheitshöhen die Grenzschichtdicke überschreiten ($k/\delta > 1$). Darüber hinaus untersuchten frühere numerische Untersuchungen oft Rauheitspatches statt vollständig verteilter Rauheit oder beobachteten keine Verstärkung der Modeninstabilität, wodurch die spezifischen Mechanismen, die den Übergang in voll verteilten, hyperschall-basierten Strömungen an stumpfen Körpern antreiben, unklar blieben.

Methodik
Diese Studie präsentiert die erste direkte numerische Simulation (DNS) eines stumpfen Zylinders in Mach-6-Querströmung mit Rauheitselementen, die über die gesamte Oberfläche verteilt sind. Die Simulationen nutzen den CHAMPS-Solver, um die kompressiblen Navier-Stokes-Gleichungen (CNSE) für ein ideales Gas zu integrieren. Die Geometrie ist ein Zylinder mit einem Durchmesser von 0,1524 m, modelliert basierend auf den experimentellen Bedingungen von Hollis (2017) (Run 38, Test 6975), mit einer Freistrom-Machzahl von 5,98.

Die Rauheit wird als sinusförmige Verteilung mit einer festen Amplitude ($A \approx 35\,\%$ der Grenzschichtdicke) und einer Wellenlänge entsprechend der nominalen Sandkornrauheit modelliert. Drei verschiedene Phasenkonfigurationen in der Spanweite wurden untersucht:

1. Gestaffelt (Staggered): Reihen sind um $180^\circ$ phasenverschoben.
2. Ausgerichtet (Aligned): Reihen sind in Phase ($0^\circ$ Verschiebung).
3. Zufällig (Random): Reihen weisen zufällig verteilte Phasenverschiebungen auf.

Zur Analyse des Transitionsprozesses setzten die Autoren Folgendes ein:

• DNS: Um die vollständige nichtlineare Entwicklung der Strömung zu erfassen.
• Lineare Störungsgleichungen (LDE): Gelöst auf zeitgemittelten Mittelfeldern aus der DNS, um Instabilitätsmechanismen mittels Impulserregung zu isolieren.

• Lineare Stabilitätstheorie (LST): Angewendet auf spanwise-gemittelte Profile zum Vergleich.

• Spektralanalyse: Schnelle Fourier-Transformationen (FFT) von Geschwindigkeits- und Druckfeldern zur Identifizierung dominanter Frequenzen und Wellenzahlen.

Kernergebnisse

• Mittelfeld und Streifen (Streaks): Alle Konfigurationen erzeugten stromlinienparallele Geschwindigkeitsstreifen via eines Lift-up-Mechanismus, der durch gegenläufige Wirbelpaare (CVPs) getrieben wird. Die ausgerichtete Konfiguration erzeugte die stärksten stationären Streifen aufgrund der konsistenten Phasenlage, die Fluid mit geringem Massestrom höher in die Grenzschicht hebt. Die gestaffelten und zufälligen Konfigurationen erzeugten schwächere, komplexere Streifenmuster. Entgegen den Vorhersagen der optimalen transienten Wachstums-Theorie zeigten die Streifenamplituden ein geringes relatives Wachstum (weniger als ein Faktor 2) und blieben im gesamten Gebiet nichtlinear, was die Nützlichkeit des linearen transienten Wachstums als Vorhersagewerkzeug für dieses Regime infrage stellt.

• Instabilitätsmechanismen: Der Transitionsprozess wurde durch konvektive Instabilitäten angetrieben, wobei der spezifische Modentyp von der Phasenlage der Rauheit abhing:

  • Ausgerichteter Fall: Dominiert von einem fundamentalen sinuosen Streifenmodus. Die Tollmien-Schlichting (T-S) Wellen wurden nur schwach destabilisiert.
  • Gestaffelte und Zufällige Fälle: Dominiert von 2D T-S Wellen (varicose Mode). Diese Konfigurationen zeigten eine signifikante Destabilisierung von T-S Wellen, was zu einem schnellen exponentiellen Wachstum führte.
  • Breakdown: In allen Fällen lösten die dominanten Instabilitätswellen, sobald sie eine ausreichende Amplitude erreichten, die Bildung von Haarpfropfwirbeln (Hairpin Vortices) in den stationären Streifen aus, was zum Breakdown und zur Turbulenz führte.

• Transitionsorte: Die DNS sagte Transitionsorte ($\Theta_{tr}$) deutlich weiter stromabwärts voraus (ca. $11^\circ$–$17^\circ$) als modernste empirische Korrelationen (die $6^\circ$–$13^\circ$ vorhersagten). Der ausgerichtete Fall trat unter den Rauheitsfällen am frühesten auf, da die zur Initiierung des LTT erforderliche Erregung gegeben war, während die gestaffelten und zufälligen Fälle später transitierten.

• Feedback-Mechanismus: Eine neuartige Erkenntnis in den gestaffelten und zufälligen Fällen war die Identifizierung einer Rückkopplungsschleife. Die in der Subsonic-Region des Zylinders erzeugte Turbulenz streut akustische Wellen, die stromaufwärts propagieren. Diese langsamen akustischen Wellen regten T-S Wellen auf der gegenüberliegenden Seite des Zylinders (über die Stagnationslinie hinweg) an ihrem neutralen Punkt an. Dieser Mechanismus, der keine externe Erregung benötigt, hält den Transitionsprozess aufrecht.

Bedeutung und Ansprüche
Die Autoren behaupten, dass diese Arbeit ein neues physikalisches Verständnis der rauheitsinduzierten Transition an hyperschallbasierten stumpfen Körpern liefert und die bisherige Abhängigkeit von der transienten Wachstums-Theorie infrage stellt. Die wesentlichen Beiträge sind:

1. Ablehnung des transienten Wachstums als primärer Treiber: Die Studie zeigt, dass für voll verteilte Rauheit mit $k/\delta \approx 0,35$ die nichtlineare Entwicklung von Streifen und die beobachteten geringen Verstärkungsverhältnisse die optimale transiente Wachstums-Theorie zu einem unzureichenden Prädiktor für den tatsächlichen Transitionsprozess machen.
2. Destabilisierung von Modeninstabilitäten: Die Arbeit etabliert, dass verteilte Rauheit Modeninstabilitäten (speziell T-S Wellen) effektiv destabilisieren kann, selbst an stumpfen Körpern, bei denen die glatte Wandströmung typischerweise eigenmodenstabil ist. Die Rauheitskonfiguration bestimmt, ob der Übergang durch Streifenmoden oder T-S Wellen getrieben wird.
3. Entdeckung eines akustischen Feedback-Mechanismus: Die Identifizierung eines selbstaufrechterhaltenden Feedback-Mechanismus, bei dem durch Turbulenz erzeugte Akustik T-S Wellen auf der gegenüberliegenden Seite des Körpers anregt, bietet eine potenzielle Erklärung für den Übergang in Subsonic-Regionen von Hyperschallströmungen – ein Phänomen, das zuvor durch lineare Mechanismen nicht erklärt werden konnte.
4. Sensitivität gegenüber der Phasenlage: Die Ergebnisse unterstreichen die extreme Sensitivität des LTT gegenüber der relativen Phasenlage von Rauheitselementen, was darauf hindeutet, dass einfache Topologien das Verhalten komplexerer, zufälliger Rauheitsmuster eingrenzen können.

Die Autoren kommen zu dem Schluss, dass auch wenn der vorgeschlagene akustische Feedback-Mechanismus einer weiteren strengen Beweisführung bedarf, die Offenlegung dieser Instabilitätsmechanismen durch hochauflösende DNS auf einen grundlegenden Wandel im Verständnis der den LTT auf hyperschallbasierten stumpfen Körpern verantwortlichen Störungen hindeutet.