Characterisation of rough-wall drag in compressible turbulent boundary layers

Die Studie untersucht, wie sich die Rauigkeit in kompressiblen turbulenten Grenzschichten über einen weiten Mach-Zahlen-Bereich verhält, und zeigt, dass zwar die Geschwindigkeitstransformation wenig Einfluss hat, jedoch eine Mach-zahl-abhängige Korrektur für den logarithmischen Bereich notwendig ist, wobei ein temperaturbasierter Korrekturfaktor die konsistentesten Ergebnisse liefert.

Das Wesentliche

Titel: Wenn rauer Boden im Wind fliegt – Wie Wissenschaftler den Luftwiderstand bei hohen Geschwindigkeiten verstehen

Stellen Sie sich vor, Sie fahren mit dem Fahrrad. Wenn Sie eine glatte Straße nehmen, gleitet der Wind sanft an Ihnen vorbei. Aber wenn Sie über einen Weg mit vielen kleinen Steinen oder Sand fahren, wird es schwieriger. Der Wind prallt auf die Steine, wirbelt herum und zieht Sie zurück. Das ist der Luftwiderstand (oder in der Physik: der Drag).

In der Luftfahrt ist das ähnlich, nur viel extremer. Wenn ein Flugzeug oder eine Rakete sehr schnell fliegt (sogar schneller als der Schall), wird die Luft nicht nur „gedrückt", sondern sie verhält sich wie ein elastischer Stoff, der sich zusammenzieht und ausdehnt. Das nennt man kompressible Strömung.

Dieser wissenschaftliche Artikel untersucht genau dieses Problem: Wie berechnet man den Widerstand, wenn eine Oberfläche rau ist und das Flugzeug dabei extrem schnell fliegt?

Hier ist die Geschichte der Forschung, einfach erklärt:

1. Das alte Regelbuch (für langsame Autos)

In der Welt der langsamen, „inkompressiblen" Strömung (wie bei einem normalen Auto) haben Wissenschaftler eine perfekte Formel gefunden. Sie nennen sie die „Sandkorn-Regel".

• Die Idee: Egal wie die rauhe Oberfläche aussieht (ob Sandpapier, Schotter oder Rost), man kann sie so tun, als wäre sie aus gleich großen Sandkörnern gemacht.
• Die Messung: Man misst, wie viel Energie die Luft verliert, wenn sie über diese Sandkörner strömt. Das Ergebnis ist eine Zahl, die man „Sandkorn-Rauhigkeit" ($k_s$) nennt.
• Das Problem: Diese Regel funktioniert super, wenn man langsam fährt. Aber was passiert, wenn man mit Überschallgeschwindigkeit fliegt?

2. Das neue Problem: Der „Schockwellen-Effekt"

Wenn ein Flugzeug mit Überschall fliegt, passiert etwas Magisches (und Ärgerliches) an den rauen Stellen:
Stellen Sie sich vor, Sie laufen durch einen dichten Wald. Wenn Sie langsam laufen, weichen die Bäume aus. Wenn Sie aber sprinten, stoßen Sie gegen die Bäume und erzeugen eine Art „Druckwelle" vor sich.
Genau das passiert bei rauen Oberflächen im Überschall: Die kleinen Unebenheiten (wie Sandkörner) erzeugen winzige Schockwellen. Diese Wellen kosten zusätzliche Energie. Das alte Regelbuch kennt diese Wellen nicht! Es sagt nur: „Rauheit ist Rauheit." Aber im Überschall ist Rauheit + Schockwellen = viel mehr Widerstand.

3. Der Experimentier-Labor-Versuch

Die Autoren dieses Papers (Wissenschaftler aus Großbritannien und Deutschland) haben sich einen riesigen Windkanal geschnappt. Sie haben zwei Arten von Sandpapier (grob und feiner) auf eine Platte geklebt und Luft mit verschiedenen Geschwindigkeiten (von ganz langsam bis fast dreifacher Schallgeschwindigkeit) darüber geblasen.

Ihre große Frage war:
„Wenn wir die Sandkorn-Regel aus der langsamen Welt kennen, können wir sie einfach auf die schnelle Welt übertragen? Wenn ja, müssen wir dann etwas an der Formel ändern?"

4. Die drei Versuche, das Problem zu lösen

Die Wissenschaftler haben drei verschiedene „Rezepte" (Korrekturfaktoren) ausprobiert, um die alte Formel für die schnelle Welt anzupassen:

• Rezept A: Der „Sandkorn-Tausch"

  • Idee: Wir nehmen einfach die Sandkorn-Größe aus der langsamen Welt und hoffen, dass sie auch im Überschall passt.
  • Ergebnis: Fehlanzeige. Das funktioniert nur bei sehr niedrigen Geschwindigkeiten. Bei hohen Geschwindigkeiten stimmt die Rechnung nicht mehr.

• Rezept B: Der „Luft-Dicke"-Trick

  • Idee: In der schnellen Welt wird die Luft anders „dick" (viskös). Vielleicht müssen wir die Rauheit mit der Dicke der Luft an der Spitze der Sandkörner vergleichen?
  • Ergebnis: Besser, aber immer noch nicht perfekt. Es funktionierte bei einigen Daten, aber bei anderen nicht.

• Rezept C: Der „Temperatur-Hebel" (Der Gewinner)

  • Idee: Wenn Luft sehr schnell über eine Wand strömt, wird sie heiß. Die Temperatur an der Wand ist anders als in der freien Luft. Die Wissenschaftler haben einen Faktor eingeführt, der genau dieses Temperatur-Verhältnis berücksichtigt.
  • Analogie: Stellen Sie sich vor, Sie messen den Widerstand eines Schuhs. Im Sommer (heiß) ist der Schuh weich und formbar, im Winter (kalt) ist er steif. Wenn Sie die Temperatur ignorieren, ist Ihre Messung falsch. Dieser Faktor „korrigiert" die Messung basierend auf der Hitze, die durch die Geschwindigkeit entsteht.
  • Ergebnis: Das war der Treffer! Wenn man diesen Temperatur-Faktor anwendet, passen fast alle Daten – sowohl die neuen Experimente als auch alte Daten aus anderen Studien – wieder in die alte, bekannte Formel.

5. Was bedeutet das für die Zukunft?

Die Wissenschaftler haben herausgefunden, dass man den Widerstand von rauen Oberflächen im Überschall berechnen kann, ABER man muss einen „Temperatur-Korrektur-Filter" über die Formel legen.

Warum ist das wichtig?

• Flugzeuge: Wenn ein Flugzeug durch Sand oder Eis fliegt, wird die Oberfläche rau. Das erhöht den Widerstand und kostet Treibstoff. Mit dieser neuen Formel können Ingenieure besser vorhersagen, wie viel Kraft sie brauchen.
• Raumfahrt: Bei Raketen, die in die Atmosphäre eintreten, ist die Oberfläche oft uneben. Hier ist die Berechnung von Hitze und Widerstand überlebenswichtig.

Fazit in einem Satz

Die alten Regeln für rauen Boden funktionieren auch im Überschall, wenn man sie mit einem „Temperatur-Thermometer" kombiniert, das die Hitze der schnellen Luft berücksichtigt. Ohne diesen Trick würde die Rechnung scheitern, weil die Luft bei hohen Geschwindigkeiten nicht nur strömt, sondern auch „schreit" (Schockwellen) und sich erhitzt.

Die Forscher hoffen nun, dass zukünftige Experimente direkt messen können, wie heiß die Wand wird und wie stark der Widerstand ist, um noch genauere „Karten" für Ingenieure zu zeichnen.

Technische Zusammenfassung

Titel: Charakterisierung des Rauigkeitswiderstands in kompressiblen turbulenten Grenzschichten

Autoren: D. D. Wangsawijaya et al.
Zielzeitschrift: Journal of Fluid Mechanics (in Prüfung)

---

1. Problemstellung und Hintergrund

In der Strömungsmechanik ist der Widerstand (Drag) durch Oberflächenrauigkeit in inkompressiblen turbulenten Grenzschichten (TBL) gut verstanden. Er wird typischerweise durch die Hama-Funktion $\Delta U^+$ (Impulsdefizit) und die äquivalente Sandkornrauigkeit $k_s$ charakterisiert, die in der Nikuradse-Beziehung (1933) mit dem logarithmischen Geschwindigkeitsprofil verknüpft sind.

In kompressiblen Strömungen (z. B. bei Überschallflugzeugen) ist die Situation komplexer:

• Neben dem Druckwiderstand der Rauigkeitselemente entsteht durch Schockwellen, die an den Rauigkeitsspitzen entstehen, ein zusätzlicher Wellenwiderstand.
• Der aktuelle Ansatz zur Charakterisierung besteht darin, die kompressiblen Geschwindigkeitsprofile zunächst mittels Transformationen (z. B. Van Driest) in ein äquivalentes inkompressibles Profil zu überführen. Anschließend wird $\Delta U^+$ berechnet und $k_s$ abgeleitet.
• Das Kernproblem: Es ist unklar, unter welchen Bedingungen ein $k_s$-Wert, der in inkompressibler Strömung gemessen wurde, auf kompressible Strömungen übertragen werden kann. Bisherige Studien basierten oft auf Einzelpunkten (einzelne Mach-Zahl und Reynolds-Zahl), was die Rauigkeit erzwungenermaßen in die inkompressible Nikuradse-Beziehung zwang, ohne die Mach-Zahl-Abhängigkeit zu berücksichtigen.

2. Methodik

Die Studie basiert auf einer umfangreichen experimentellen Untersuchung im Windkanal der Universität der Bundeswehr München.

• Versuchsaufbau: Ein intermittierender Blow-down-Windkanal (Trisonic Wind Tunnel) ermöglichte Mach-Zahlen von $0.3 \le M \le 2.9$ und Reynolds-Zahlen bis $Re_\tau \approx 30.000$.
• Oberflächen: Es wurden zwei Arten von Sandpapier-Rauigkeiten untersucht (P60 und P24 Körnung), die realistische, abriebbedingte Oberflächen nachahmen. Eine glatte Wand diente als Referenz.
• Messverfahren:
  • PIV (Particle Image Velocimetry): Zweidimensionale, zeitaufgelöste Geschwindigkeitsmessungen im $x-y$-Bereich zur Erfassung der Grenzschichtprofile.
  • Schlieren-Visualisierung: Zur Sichtbarmachung von Schockwellen vor den Rauigkeitselementen.
  • Temperaturabschätzung: Da direkte Wandtemperaturmessungen nicht möglich waren, wurden adiabatische Wandbedingungen angenommen und Temperaturprofile über die Crocco-Beziehung und Sutherland-Formel berechnet.
• Datenanalyse:
  • Anwendung von fünf verschiedenen Kompressibilitätstransformationen (HO, VD, BR, TL, VO), um die kompressiblen Profile in inkompressible äquivalente ($U_I, y_I$) zu überführen.
  • Bestimmung des Impulsdefizits $\Delta U^+_I$ durch Anpassung an das logarithmische Gesetz.
  • Schätzung der Wandreibungsspannung ($C_f$) mittels einer modifizierten Clauser-Anpassung.

3. Wichtige Beiträge und Ergebnisse

A. Unabhängigkeit von der Transformation

Die Ergebnisse zeigen, dass das Impulsdefizit $\Delta U^+_I$ weitgehend unempfindlich gegenüber der Wahl der spezifischen Geschwindigkeitstransformation ist. Das bedeutet, dass die Wahl der Transformationsfunktion (z. B. Van Driest vs. Trettel & Larsson) für die Bestimmung des Rauigkeitswiderstands weniger kritisch ist als die Berücksichtigung der Mach-Zahl-Effekte selbst.

B. Mach-Zahl-Abhängigkeit im vollrauhigen Regime

Im vollrauhigen Regime ($\Delta U^+ > 7$) wurde ein signifikanter Effekt beobachtet:

• Bei subsonischen und transsonischen Strömungen ($M \le 0.8$) kollabieren die Daten gut auf die klassische Nikuradse-Beziehung.
• Im Überschallbereich ($M = 2$ und $2.9$) verschiebt sich der Schnittpunkt des logarithmischen Gesetzes nach oben. Dies deutet auf einen zusätzlichen Widerstandsanteil hin, der mit der Mach-Zahl skaliert und vermutlich durch die Wellenwiderstände (Schockwellen vor den Rauigkeitsspitzen) verursacht wird.

C. Evaluation von Korrekturfaktoren

Die Autoren testeten drei empirische Skalierungsansätze, um die Daten aus kompressiblen Strömungen auf die inkompressible Nikuradse-Kurve zurückzuführen:

1. Äquivalente inkompressible $k_s$: Versuche, $k_s$ aus inkompressiblen Referenzstudien zu übernehmen, zeigten nur begrenzten Erfolg. Die Übereinstimmung hing stark von der genauen Geometrie und den Strömungsbedingungen ab und war nicht universell übertragbar.
2. Viskositäts-Skalierung ($k^* = k / \nu^+_\infty$): Eine Skalierung der physikalischen Rauigkeitshöhe $k$ mit dem Verhältnis der kinematischen Viskosität im Freistrom zur Viskosität an der Wand ($\nu_\infty / \nu_w$) führte zu einer guten Kollapsierung der eigenen experimentellen Daten. Bei Daten früherer Studien aus der Literatur funktionierte dieser Ansatz jedoch nicht konsistent.
3. Korrektur des Impulsdefizits ($\sqrt{1/F_c} \cdot \Delta U^+_I$): Dies war der erfolgreichste Ansatz.
   • Der Faktor $F_c$ ist eine Korrektur für den Hautreibungskoeffizienten glatter Wände, die von der Temperaturverhältnis $T_\infty / T_w$ abhängt (abgeleitet aus der Van Driest-II-Transformation).
   • Durch Multiplikation des gemessenen $\Delta U^+_I$ mit $\sqrt{1/F_c}$ kollabierten sowohl die eigenen experimentellen Daten als auch ein großer Teil der Daten aus früheren Studien (verschiedene Rauigkeitstypen, Mach-Zahlen bis hypersonisch) konsistent auf die Nikuradse-Beziehung.

4. Signifikanz und Schlussfolgerungen

• Herausforderung der aktuellen Modelle: Der bestehende Rahmen zur Widerstandscharakterisierung (Transformation glatter Wände + Hama-Funktion) berücksichtigt nicht die durch Rauigkeit erzeugten Schockwellen. Die Schlieren-Bilder der Studie zeigen deutlich, dass diese Schockwellen bis in den Freistrom reichen und einen signifikanten zusätzlichen Widerstand verursachen.
• Empirische Lösung: Der vorgeschlagene Korrekturfaktor $\sqrt{1/F_c}$ bietet derzeit die konsistenteste Möglichkeit, den Widerstand kompressibler, rauer Grenzschichten unter Verwendung etablierter inkompressibler Modelle zu schätzen.
• Limitationen: Alle getesteten Korrekturfaktoren sind rein empirisch. Da keine direkten Messungen der Wandreibung oder Wandtemperatur an den rauen Oberflächen durchgeführt werden konnten, bleiben Unsicherheiten bestehen.
• Zukunftsausblick: Die Studie plädiert für die Entwicklung einer maßgeschneiderten Transformation für raue Wände, die nicht nur die Kompressibilität, sondern auch die spezifischen Rauigkeitseigenschaften und die daraus resultierenden Wellenwiderstände explizit modelliert. Dies erfordert zukünftige Experimente mit direkten Messungen von Wandtemperatur und Wandreibungsspannung.

Fazit: Die Arbeit liefert einen wichtigen Schritt zur Quantifizierung des Rauigkeitswiderstands in kompressiblen Strömungen und zeigt, dass eine einfache Übertragung inkompressibler $k_s$-Werte ohne Mach-Zahl-Korrektur unzureichend ist. Der vorgeschlagene Korrekturfaktor basierend auf dem Temperaturverhältnis ist ein vielversprechender empirischer Ansatz für die Ingenieurspraxis, bis physikalisch fundiertere Modelle entwickelt sind.