Wave Particle Turbulent Simulation of Spatially Developing Round Jets Using a Non Equilibrium Transport Model with a Mixing Length Characteristic Time Closure

In dieser Arbeit wird die neu entwickelte Wave-Particle Turbulent Simulation (WPTS) unter Verwendung eines auf der Prandtlschen Mischungshypothese basierenden Zeitabschlusssystems auf räumlich sich entwickelnde Rundstrahlen angewendet und erfolgreich validiert.

Das Wesentliche

Das Rätsel der wirbelnden Luft: Wie man Turbulenzen „sieht“, ohne alles zu zählen

Stellen Sie sich vor, Sie stehen an einer Autobahn an einem windigen Tag. Die Luft bewegt sich nicht einfach nur gleichmäßig von A nach B. Stattdessen gibt es riesige Luftmassen, die wirbeln, sich gegenseitig hochschleudern und in winzigen, chaotischen Wirbeln tanzen. In der Wissenschaft nennen wir das Turbulenz.

Das Problem für Ingenieure (die zum Beispiel Flugzeuge oder Triebwerke bauen) ist: Diese Wirbel sind so winzig und so schnell, dass selbst die stärksten Supercomputer der Welt kaum in der Lage sind, jeden einzelnen kleinen Wirbel exakt zu berechnen. Es ist, als wollten Sie ein riesiges, wildes Feuer beobachten, indem Sie jedes einzelne glühende Staubkorn zählen – das dauert viel zu lange und kostet zu viel Energie.

Die alte Methode: Die „Honig-Taktik“

Früher hat man versucht, das Chaos zu vereinfachen, indem man so tat, als wäre die Luft nicht wie Gas, sondern wie zäher Honig. Man hat einfach eine „künstliche Zähigkeit“ (die sogenannte Viskosität) hinzugefügt, um die kleinen Wirbel mathematisch „glattzubügeln“. Das funktioniert oft, aber es ist ungenau, weil es die echte, wilde Bewegung der Luft ignoriert.

Die neue Idee: Die „Welle-Teilchen-Party“ (WPTS)

Die Forscher in dieser Arbeit haben einen neuen Weg gefunden, den sie WPTS nennen. Stellen Sie sich das wie eine Party in einem Club vor:

1. Die Welle (Das Grundrauschen): Das ist die Musik und die allgemeine Stimmung im Raum. Sie ist überall präsent, gleichmäßig und lässt sich leicht beschreiben. Das entspricht der „glatten“ Strömung der Luft.
2. Die Teilchen (Die wilden Tänzer): Das sind die Leute, die mitten auf der Tanzfläche völlig ausrasten. Sie bewegen sich unvorhersehbar, springen hin und her und beeinflussen die Leute um sie herum. In der Simulation sind das „stochastische Teilchen“.

Anstatt zu versuchen, jeden einzelnen Tänzer in der ganzen Welt zu berechnen, konzentriert sich das Modell nur dort auf die „Teilchen“, wo es wirklich wild zugeht (also dort, wo die Turbulenz hoch ist). Sobald die Party in einer Ecke des Raumes ruhig wird, „verschwinden“ die wilden Tänzer einfach und es bleibt nur noch die ruhige Musik (die Welle) übrig. Das spart unglaublich viel Rechenleistung!

Der Clou: Die „Prandtl-Regel“ (Der Misch-Abstand)

Die Forscher haben diesem Modell nun eine neue Regel gegeben, basierend auf einer alten Idee von Prandtl. Man kann sich das so vorstellen: Ein wilder Tänzer (ein Luftwirbel) rennt nicht einfach nur ziellos herum. Er rennt eine bestimmte Strecke, bevor er wieder mit der Masse zusammenstößt und seine Energie abgibt.

Die Forscher haben eine Formel entwickelt, die berechnet, wie weit dieser „Tänzer“ rennen kann, je nachdem, wie weit er vom Zentrum des Luftstrahls (dem Jet) entfernt ist. Je weiter man nach außen kommt, desto eher wird die Party ruhig.

Was haben sie bewiesen?

Sie haben dieses Modell an einem „Luftstrahl“ (einem Jet) getestet – so ähnlich wie das Abgas aus einer Düse.

• Das Ergebnis: Selbst wenn sie die Auflösung der Simulation niedrig hielten (also nicht jeden winzigen Punkt berechneten), konnte das Modell die Form und die Kraft des Luftstrahls fast perfekt vorhersagen.
• Der Clou bei hohen Geschwindigkeiten: Selbst als sie die Geschwindigkeit massiv erhöhten (von Reynolds-Zahl 5.000 auf 20.000), blieb das Modell stabil und präzise.

Zusammenfassend

Anstatt zu versuchen, das Chaos mit purer Gewalt (Rechenpower) zu bezwingen, nutzt diese neue Methode eine kluge Mischung: Sie nutzt eine ruhige „Welle“ für die Ordnung und schickt gezielt „Teilchen“ los, um das Chaos zu simulieren. Das ist schneller, effizienter und verhält sich viel mehr wie die echte Welt.

Technische Zusammenfassung

Technische Zusammenfassung: Wave–Particle Turbulent Simulation (WPTS) von räumlich entwickelnden Rundstrahlströmungen

1. Problemstellung

Die Simulation turbulenter Strömungen stellt eine der größten Herausforderungen der numerischen Strömungsmechanik dar. Während die direkte numerische Simulation (DNS) hochpräzise Ergebnisse liefert, ist sie aufgrund des enormen Rechenaufwands bei hohen Reynolds-Zahlen (Re) praktisch nicht anwendbar. Herkömmliche RANS-Modelle (Reynolds-Averaged Navier–Stokes) sind zwar effizient, basieren jedoch meist auf Gleichgewichtsannahmen (Eddy-Viscosity-Modelle), die nicht-lokale Transporteffekte und Nichtgleichgewichtszustände nur unzureichend erfassen können. Das Ziel dieser Arbeit ist es, eine effiziente Methode zu entwickeln, die auf groben Gittern arbeitet, aber dennoch die wesentlichen physikalischen Merkmale der Turbulenz – insbesondere den Nichtgleichgewichtstransport – präzise abbildet.

2. Methodik

Die Autoren verwenden die Wave–Particle Turbulent Simulation (WPTS), einen multiskaligen Ansatz, der auf der kinetischen Gastheorie basiert.

• Wave-Particle-Konzept: Die Strömung wird in zwei Komponenten zerlegt:
  • Wave-Komponente (Eulerisch): Repräsentiert die großskaligen, durch das Gitter aufgelösten Strukturen mittels der Navier-Stokes-Gleichungen.
  • Particle-Komponente (Lagrange): Repräsentiert die unaufgelösten kleinskaligen Turbulenzen (z. B. turbulente kinetische Energie, TKE) durch stochastische Partikel.
• Neuartige Zeit-Abschlussmodellierung ($\tau_t$): Der entscheidende Beitrag ist die Kopplung von WPTS mit einem neuen Modell für die turbulente Charakteristikzeit $\tau_t$, das direkt aus Prandls Mischungshypothese abgeleitet wurde. Anstatt einer einfachen Viskositätsanpassung wird die Mischungslänge $l$ als Funktion des axialen Abstands $x$ modelliert ($l = \sqrt{D x + b_{ml} l}$). Die turbulente Viskosität $\nu_t$ wird über die Deformationsrate (Strain Rate Tensor $S$) definiert.
• Nichtgleichgewichtstransport: Im Gegensatz zu klassischen Modellen bewegen sich die Partikel über eine endliche Distanz, bevor sie mit dem Hintergrundfeld interagieren. Dies ermöglicht die natürliche Abbildung von nicht-lokalen Effekten.

3. Wesentliche Beiträge

• Entwicklung eines Mischungslängen-basierten $\tau_t$-Modells: Ein maßgeschneidertes Modell für Jet-Strömungen, das die räumliche Entwicklung der Turbulenz berücksichtigt.
• Einheitliches Framework: Die Methode kann nahtlos zwischen laminaren und turbulenten Regimen wechseln, da die Partikelkonzentration in laminaren Zonen automatisch gegen Null geht.
• Skalierbarkeit: Die Einführung einer Reynolds-Zahl-abhängigen Skalierungsstrategie für die Modellkoeffizienten und die Gitterauflösung.

4. Ergebnisse

Die Validierung erfolgte anhand von räumlich entwickelnden Rundstrahlströmungen (Round Jets) bei zwei Reynolds-Zahlen:

• Re = 5.000: Das WPTS-Modell reproduzierte die Selbstähnlichkeit (Self-Similarity) der Strömung exzellent. Der berechnete Abfall der Zentrumgeschwindigkeit ($B_u = 5,55$) lag sehr nah an den Referenzwerten aus DNS- und Experimentaldaten (ca. 5,50 bis 6,06). Die Profile der Reynolds-Spannungen zeigten eine hervorragende Übereinstimmung und Kollaps der Profile an verschiedenen axialen Positionen.
• Re = 20.000: Trotz der deutlich höheren Reynolds-Zahl und der Verwendung eines vergleichsweise groben Gitters (ca. $10^6$ Zellen) lieferte das Modell präzise Ergebnisse. Der Abfallwert $B_u$ von 6,18 wich nur um etwa 2 % vom experimentellen Referenzwert ab.
• Skalierungsregeln: Die Autoren stellten fest, dass die Modellierung für höhere Re-Zahlen durch eine Skalierung der Gittergröße ($\propto Re^{-0,20}$) und der Modellkoeffizienten ($\propto Re^{-0,40}$) stabilisiert werden kann.

5. Bedeutung der Arbeit

Die Arbeit demonstriert, dass die WPTS-Methode ein vielversprechendes Werkzeug für die industrielle Anwendung darstellt. Sie schließt die Lücke zwischen der hohen Rechenlast von LES/DNS und der geringen Genauigkeit klassischer RANS-Modelle. Durch die Kombination kinetischer Theorie mit der klassischen Mischungshypothese wird eine Methode geschaffen, die sowohl physikalisch fundiert (Nichtgleichgewichtstransport) als auch numerisch effizient (grobe Gitter) ist und eine robuste Vorhersage über verschiedene Reynolds-Zahlen hinweg ermöglicht.