Fully Turbulent Wakes at Low Reynolds Numbers: the Case of the Thin Flat Plate

Diese Arbeit zeigt durch direkte numerische Simulation und experimentellen Vergleich auf, dass die Nachlaufströmung hinter einer dünnen zweidimensionalen Flachplatte bei einer relativ niedrigen Reynolds-Zahl von 400 voll turbulent wird und statistische sowie spektrale Merkmale aufweist, die von Turbulenznachläufen bei höheren Reynolds-Zahlen ununterscheidbar sind, ein Übergangspfad, der sich signifikant von dem kanonischer kreisförmiger oder quadratischer Zylinder unterscheidet.

Das Wesentliche

Stellen Sie sich vor, Sie halten ein dünnes, flaches Stück Pappe (wie eine Spielkarte) in einem starken Wind. Während der Wind auf die Karte trifft, erzeugt er eine unordentliche, wirbelnde Spur aus Luft hinter ihr, die man „Nachlauf“ (Wake) nennt. Lange Zeit glaubten Wissenschaftler, dass der Wind sehr schnell blasen müsste, damit dieser Nachlauf wirklich chaotisch und „turbulent“ wird, oder dass das Objekt eine bestimmte Form wie ein runder Rohrkörper oder ein quadratischer Block haben müsste.

Dieses Paper erzählt eine andere Geschichte. Die Forscher entdeckten, dass die Luft hinter einer dünnen, flachen Platte bei einer viel geringeren Windgeschwindigkeit vollkommen chaotisch und turbulent wird, als man es erwartet hätte. Tatsächlich geschieht dies bei einer Geschwindigkeit, bei der die Luft hinter anderen Formen noch relativ ruhig und geordnet ist.

Hier ist eine Aufschlüsselung ihrer Ergebnisse unter Verwendung einfacher Analogien:

1. Die Überraschung beim „Turbulenz-Schwellenwert“

Stellen Sie sich Turbulenz wie eine belebte Tanzfläche vor.

• Der alte Glaube (Runde Zylinder): Wenn ein runder Pfosten im Wind steht, beginnt die Luft dahinter als ein ruhiger, rhythmischer Tanz (ein Hin- und Herwiegen). Es bedarf viel Energie (hohe Geschwindigkeit), bevor die Tänzer anfangen, gegeneinander zu stoßen, wild herumzuspringen und ein chaotisches Durcheinander zu verursachen (Turbulenz). Dieser Übergang geschieht langsam über einen weiten Bereich von Geschwindigkeiten.
• Die Neuentdeckung (Dünne Flachplatte): Die Forscher fanden heraus, dass die Tanzfläche bei einer dünnen, flachen Platte fast augenblicklich vom ruhigen Tanz zum wilden Moshpit wird. Schon bei einer relativ geringen Windgeschwindigkeit (Reynolds-Zahl 400) ist die Luft hinter der Platte bereits vollkommen chaotisch. Sie durchläuft nicht die langsamen, rhythmischen Phasen, die runde Pfosten durchlaufen. Sie springt direkt zur Party.

2. Wie sie es bewiesen haben

Um sicherzugehen, dass sie sich das nicht nur einbildeten, handelten die Forscher wie Detektive, die Tatorte vergleichen.

• Die Simulation (Das virtuelle Labor): Sie nutzten Supercomputer, um den Wind zu simulieren, der auf die Platte bei niedrigen Geschwindigkeiten trifft (Re 150 und Re 400).
• Der Realwelt-Test (Der Windkanal):ung Sie betrachteten auch reale Experimente, bei denen der Wind viel schneller wehte (Re 12.500 und Re 19.700).
• Die Übereinstimmung: Als sie die langsame Computersimulation (Re 400) mit den Hochgeschwindigkeits-Experimenten aus der realen Welt verglichen, passten die Muster perfekt zusammen. Die „Fingerabdrücke“ der Turbulenz – wie sich die Luft bewegte, wie viel Energie sie hatte und wie sie wirbelte – waren identisch.
• Die Kontrollgruppe: Als sie die Simulation bei einer noch niedrigeren Geschwindigkeit betrachteten (Re 150), waren die Muster völlig anders. Die Strömung befand sich noch in der „ruhigen“ Phase, noch nicht im Chaos. Dies bewies, dass der Übergang zum Chaos irgendwo zwischen 150 und 400 stattfindet.

3. Der „Fingerabdruck“ der Turbulenz

Woran erkennt man, ob eine Strömung wirklich turbulent ist? Das Paper sucht nach spezifischen „Lebenszeichen“ in den Daten:

• Das Energiespektrum (Der Klang des Windes): In einer ruhigen Strömung ist die Energie auf wenige spezifische Töne konzentriert (wie eine Flöte, die einen einzelnen Ton spielt). In einer turbulenten Strömung klingt es wie weißes Rauschen oder statisches Rauschen, wobei die Energie über ein riesiges Spektrum an Frequenzen verteilt ist. Die Forscher fanden heraus, dass das „Geräusch“ des Windes hinter der Platte bei Re 400 bereits voller dieses chaotischen Rauschens war, genau wie in den Hochgeschwindigkeits-Experimenten.
• Die „Intermittenz“ (Der gelegentliche Schrei): In einer wirklich turbulenten Strömung wirbelt die Luft nicht nur sanft; sie weist plötzliche, intensive Ausbrüche von Geschwindigkeit und Rotation auf. Die Forscher fanden diese „Schreie“ in den Daten bei Re 400, aber sie waren bei Re 150 nicht vorhanden.

4. Warum ist das anders?

Das Paper legt nahe, dass der Grund für diesen plötzlichen Sprung die Form des Objekts ist.

• Runde/Quadratische Objekte: Wenn Wind auf ein rundes oder quadratisches Objekt trifft, wirkt der hintere Teil des Objekts wie ein Schild, der die Luftströmung dahinter stabilisiert. Es braucht viel Energie, um diese Stabilität zu brechen.
• Die dünne Platte: Weil die Platte so dünn ist, gibt es kein „Hinterteil“, das die Luft abschirmt. Die Druckschwankungen (das Drücken und Ziehen der Luft) sind von Beginn an direkt mit den wirbelnden Wirbeln verbunden. Es ist, als würde man versuchen, einen Bleistift auf seiner Spitze zu balancieren im Vergleich zu einem Bowlingball; der Bleistift (die dünne Platte) ist von Natur aus instabil und kippt viel schneller in das Chaos.

Das Fazit

Dieses Paper verändert unser Verständnis davon, wie Luft um flache Objekte strömt. Es beweist, dass dünne Flachplatten bereits bei überraschend niedrigen Geschwindigkeiten voll turbulente Nachläufe erzeugen, viel niedriger als runde oder quadratische Objekte. Der Übergang zum Chaos ist für diese Formen kein langsamer, gradueller Prozess; es ist ein plötzlicher, fundamentaler Umschlag, der sehr früh im Geschwindigkeitsbereich stattfindet.

Die Forscher haben nicht diskutiert, wie dies auf den Bau von Brücken, das Design von Autos oder medizinische Geräte Anwendung findet. Sie konzentrierten sich strikt darauf, zu beweisen, dass dieses Phänomen existiert und wie die Physik der Luftströmung von dem abweicht, was wir bisher angenommen hatten.

Technische Zusammenfassung

Technische Zusammenfassung: Vollständig turbulente Nachlaufströmungen bei niedrigen Reynolds-Zahlen: Der Fall der dünnen Flachplatte

Problemstellung
Die numerische Simulation der nominell zweidimensionalen (2D) Strömung an einer dünnen Flachplatte, die senkrecht zu einer gleichförmigen Anströmung steht, hat sich historisch als schwierig erwiesen, da CFD-Studien widersprüchliche Ergebnisse hinsichtlich der Nachlaufdynamik bei niedrigen Reynolds-Zahlen ($Re < 1000$) geliefert haben. Während experimentelle Daten darauf hindeuten, dass integrale Parameter (wie etwa Drag-Koeffizienten und Strouhal-Zahlen) im 3D-Regime schwach von $Re$ abhängen, haben Simulationen sehr unterschiedliche Nachlaufverhaltensweisen berichtet. Eine zentrale Frage bleibt: Geht der Nachlauf einer dünnen Flachplatte unmittelbar nach dem Einsetzen der 3D-Instabilitäten in einen vollständig turbulenten Zustand über, ähnlich wie bei kanonischen Zylinder-Nachläufen, oder folgt er einem distinkten, mehrstufigen Übergangspfad? Die Autoren untersuchen, ob der Nachlauf einer dünnen Flachplatte bereits bei einer relativ niedrigen Reynolds-Zahl von $Re = 400$ vollständig turbulent ist – ein Regime, in dem kanonische kreisförmige und quadratische Zylinder-Nachläufe typischerweise eine starke $Re$-Abhängigkeit und nicht-turbulente Fluktuationscharakteristika aufweisen.

Methodik
Die Studie verwendet einen facettenreichen Ansatz, der Direct Numerical Simulations (DNS) mit experimenteller Validierung kombiniert:

1. Numerische Simulationen: Es wurden zwei unabhängige DNS-Codes eingesetzt, um die Robustheit zu gewährleisten:

   • Spektrale/hp-Elementmethode (SEM): Implementiert in Nektar++, verwendet für die primäre hochgenaue DNS bei $Re = 400$ (DNS400) und $Re = 150$ (DNS150).
   • Finite-Volumen-Methode (FVM): Implementiert in OpenFOAM, verwendet für Gitterkonvergenzstudien und Kreuzvalidierung (DNS400V).
   • Das Rechengebiet wurde gekürzt und mit No-Slip-Bedingungen an der Platte, Freistrom-Einlassbedingungen und spezifischen Auslassbedingungen gemäß Dong et al. versehen. In der Spannweitenrichtung wurde Periodizität erzwungen.
   • Die Gitterkonvergenz wurde streng getestet, um sicherzustellen, dass die Kolmogorov-Längenskala ($\eta$) mit $\Delta x/\eta$-Werten aufgelöst wird, die vergleichbar mit oder besser als die bestehende Literatur für ähnliche Geometrien sind.

2. Experimentelle Daten: Zwei experimentelle Datensätze wurden mittels zeitaufgelöster stereoskopischer Particle Image Velocimetry (PIV) in einem Windkanal an der University of Calgary gewonnen:

   • EX12K: $Re = 12.500$
   • EX20K: $Re = 19.700$
     Diese Datensätze repräsentieren vollständig turbulente Regime und dienen als Benchmark für die Niedrig-$Re$-Simulationen.

3. Analyse: Die Studie vergleicht globale integrale Größen, mittlere Geschwindigkeitsfelder, Reynolds-Spannungen, turbulente kinetische Energie-Budgets (TKE), Energiespektren und statistische Eigenschaften von Geschwindigkeitsgradienten (speziell Nicht-Gaussianität) über die verschiedenen $Re$-Fälle hinweg.

Wesentliche Beiträge und Ergebnisse

• Frühes Einsetzen der Turbulenz: Die primäre Erkenntnis ist, dass der Nachlauf einer dünnen Flachplatte bei $Re = 400$ vollständig turbulent wird. Dies liegt signifikant niedriger als der Übergang, der bei kanonischen kreisförmigen ($Re \approx 10^5$) oder quadratischen ($Re \approx 6000$) Zylinder-Nachläufen beobachtet wird.
• Quantitative Übereinstimmung mit High-$Re$-Experimenten: Die DNS-Ergebnisse bei $Re = 400$ (DNS400) zeigen eine bemerkenswerte Übereinstimmung mit den experimentellen Daten bei $Re = 12.500$ und $Re = 19.700$ (EX12K, EX20K) hinsichtlich:
  • Globaler Größen: Der mittlere Drag-Koeffizient ($C_D$) und die Strouhal-Zahl ($St$) sind langsam variierende Funktionen von $Re$für$Re \ge 400$.
  • Strömungsstruktur: Mittlere Stromlinien, Geschwindigkeitsprofile und Reynolds-Spannungsverteilungen ($u'^2, v'^2, w'^2$) sind innerhalb der Messunsicherheit statistisch ununterscheidbar zwischen $Re = 400$ und den High-$Re$-Experimenten.
  • Energiebudgets: Die räumliche Verteilung der TKE-Produktion und der Advektionsterme in der $k$-Transportgleichung ist über die Fälle $Re = 400$ und die High-$Re$-Fälle hinweg konsistent.
• Turbulente Signaturen bei niedrigem $Re$: Der Fall $Re = 400$ weist definitive Merkmale der Turbulenz auf, die im Fall $Re = 150$ fehlen:
  • Energiespektren: Die fluktuierenden Geschwindigkeitskomponenten zeigen Breitbandspektren mit einem Inertialbereich, der näherungsweise der $f^{-5/3}$-Skalierung (Kolmogorovs Gesetz) folgt, was auf eine Energiekaskade hindeutet. Im Gegensatz dazu zeigen die $Re = 150$ Spektren Energietransfers zwischen nur wenigen Moden.
  • Nicht-Gaussianität: Die Wahrscheinlichkeitsdichtefunktionen (PDFs) der fluktuierenden Dehnungsrate ($\Sigma$) und Rotationsrate ($\Omega$) weisen schwere Enden (heavy tails) auf, was Intermittenz und seltene Ereignisse hoher Intensität signalisiert, die charakteristisch für Turbulenz sind. Der $Re = 150$ Fall zeigt Gaußsche Verteilungen.
• Distinkter Übergangspfad: Der Übergang bei $Re = 150$ beinhaltet 3D-Instabilitäten, besitzt jedoch nicht die kohärenten, voll turbulenten Strukturen, die bei $Re = 400$ zu sehen sind. Der $Re = 150$ Nachlauf weist einen schwächeren Rückstrom, eine langsamere Nachlauferholung und eine nicht-monotone Entwicklung der Nachlaufbreite auf, was sich grundlegend vom Verhalten kanonischer Zylinder unterscheidet, bei denen 3D-Instabilitäten (Moden A, B, C) über einen weiten $Re$-Bereich vor der vollen Turbulenz evolvieren.

Bedeutung und Behauptungen
Das Paper argumenttiert, dass der Pfad zum Übergang zur Turbulenz für eine dünne Flachplatte fundamental anders ist als für kanonische Zylinder. Während Zylinder-Nachläufe eine langwierige Sequenz von Instabilitäten (Modus A, B, C) über einen weiten Reynolds-Zahlenbereich durchlaufen, bevor sie einen vollständig turbulenten Zustand erreichen, geht der Nachlauf der dünnen Flachplatte fast unmittelbar nach dem Einsetzen der 3D-Instabilitäten (zwischen $Re = 150$ und $Re = 400$) in einen vollständig turbulenten Zustand über.

Die Autoren legen nahe, dass dieser Unterschied auf das Fehlen eines stromabwärts gelegenen Hinterkörpers (Afterbody) der Geometrie zurückzuführen sein könnte. Bei Zylindern stabilisiert der Hinterkörper die Strömung und entkoppelt Druckfluktuationen vom Ablösepunkt. Im Gegensatz dazu sind bei dünnen Flachplatten Druckfluktuationen direkt mit der Wirbelbildungsregion gekoppelt, was potenziell einen früheren und schnelleren Übergang zur Turbulenz vorantreibt. Diese Ergebnisse helfen, die Diskrepanzen in der bisherigen Literatur bezüglich der Sensitivität von Flachplatten-Nachläufen gegenüber Aspektverhältnissen und Reynolds-Zahlen aufzuklären, indem sie feststellen, dass das "vollständig turbulente" Regime für diese Geometrie bei einem viel niedrigeren $Re$ beginnt, als bisher für kanonische Blöcke angenommen wurde.