Generalization of the viscous stress tensor to the case of non-small gradients of hydrodynamic velocity: a path to numerical modeling of turbulence non-locality

Dieser Artikel verallgemeinert die Chapman-Enskog-Methode, um eine Integraldarstellung des viskosen Spannungstensors für große Geschwindigkeitsgradienten abzuleiten, was die numerische Modellierung der Nichtlokalität von Turbulenzen und von Phänomenen wie tangentialen Diskontinuitäten ermöglicht, die mit den herkömmlichen Navier-Stokes-Formulierungen nur schwer erfasst werden können.

Das Wesentliche

Das große Problem: Der „Smoothie" versus der „Sturm"

Stellen Sie sich vor, Sie versuchen vorherzusagen, wie sich eine Flüssigkeit (wie Wasser oder Luft) bewegt. Seit über einem Jahrhundert nutzen Wissenschaftler eine berühmte Regelmenge, die Navier-Stokes-Gleichungen genannt wird. Diese Regeln beruhen auf einer spezifischen Komponente, dem viskosen Spannungstensor.

Stellen Sie sich diesen Tensor als einen „Reibungsrechner" vor. In der Standardversion geht er davon aus, dass der Widerstand (die Reibung), den eine Flüssigkeit spürt, wenn Sie sie antreiben, nur davon abhängt, wie schnell sich die Flüssigkeit direkt neben dem betrachteten Punkt bewegt. Es ist, als würde man annehmen, dass der Widerstand, den Sie beim Rühren einer Tasse Kaffee spüren, ausschließlich von den Kaffeemolekülen bestimmt wird, die Ihren Löffel berühren.

Der Fehler:
Der Autor, A.B. Kukushkin, weist darauf hin, dass dieser Standard-„Reibungsrechner" versagt, wenn Dinge chaotisch werden, wie in einem turbulenten Sturm oder wenn zwei Flüssigkeitsströme mit hoher Geschwindigkeit aneinander vorbeigleiten (tangentialle Diskontinuitäten).

• Die Analogie: Stellen Sie sich eine Menschenmenge vor, die einen Flur entlanggeht. Das Standardmodell geht davon aus, dass jeder nur mit der Person direkt neben sich zusammenstößt. In einer echten Menschenmenge (Turbulenz) kann eine Person jedoch von jemandem gestoßen werden, der drei Reihen weiter hinten steht, oder eine Bewegungswelle kann sich durch den ganzen Raum fortpflanzen. Das Standardmodell ignoriert diese „Fernwirkungen".
• Das Paradoxon: Die Standardmathematik führt auch zu einem seltsamen Ergebnis: Sie legt nahe, dass die Flüssigkeit eigentlich leichter fließen sollte (niedrigere Viskosität), wenn die Teilchen häufiger kollidieren (wie in einem dichten Nebel). Dies widerspricht unserer Intuition.

Die Lösung: Den ganzen Überblick behalten

Kukushkin schlägt eine neue Methode zur Berechnung dieser Reibung vor. Anstatt nur die unmittelbare Umgebung zu betrachten, betrachtet seine neue Formel die gesamte Geschichte und den gesamten Ort der Bewegung der Flüssigkeit.

Der neue Ansatz:

1. Aufgeben der Regel der „kleinen Schritte": Die alte Mathematik (Chapman-Enskog-Methode) funktioniert nur, wenn sich die Flüssigkeit sehr langsam und glatt verändert. Kukushkin bricht diese Regel auf. Er erlaubt plötzliche, scharfe Geschwindigkeitsänderungen, wie sie in echter Turbulenz vorkommen.
2. Die „Boten"-Analogie: Anstatt nur die Flüssigkeit an einem Ort zu betrachten, stellen Sie sich vor, die Flüssigkeit ist voller winziger „Boten" (Wirbel oder Störungen).
   • Im alten Modell spricht ein Bot nur mit seinem Nachbarn.
   • Im neuen Modell von Kukushkin wird ein Bot an einem Ort geboren, fliegt durch den Raum und liefert seine Nachricht an einen weit entfernten Ort, bevor er stoppt.
3. Die Integralformel: Die neue Mathematik ist ein Integral (eine Summe über einen großen Bereich). Sie berechnet die Spannung (Reibung) an einem bestimmten Punkt, indem sie die Effekte all dieser Boten addiert, die von überall anders in der Flüssigkeit zu diesem Punkt wandern.

Warum dies wichtig ist

1. Behebung des „Paradoxons":
Indem diese Boten lange Strecken zurücklegen dürfen, korrigiert die neue Formel das seltsame Paradoxon bezüglich der Viskosität. Sie erklärt, warum Flüssigkeiten sich so verhalten, wie sie es tun, selbst wenn Teilchen häufig kollidieren. Die „langen Flüge" der Boten berücksichtigen den Widerstand auf eine Weise, die das alte Modell der „kleinen Schritte" nicht konnte.

2. Verbindung zum realen Chaos (Richardsons Gesetz):
Das Paper erwähnt eine berühmte Beobachtung, das Richardsons $t^3$-Gesetz.

• Die Analogie: Wenn Sie zwei Blätter in einen turbulenten Fluss fallen lassen, sagt das Standardmodell voraus, dass sie sich langsam voneinander entfernen (wie $t^2$). In Wirklichkeit fliegen sie jedoch viel schneller auseinander (wie $t^3$).
• Der Zusammenhang: Dieses neue „Fernwirkungs"-Modell erklärt auf natürliche Weise, warum sich Teilchen so schnell voneinander trennen. Die Boten reisen weit und schnell, tragen die Störung durch die Flüssigkeit, was der realen Beobachtung entspricht, wie sich Turbulenz ausbreitet.

3. Eine Brücke zu besseren Computersimulationen:
Derzeit müssen Computersimulationen von Turbulenzen oft „Schummeln" oder erfundene Zahlen verwenden, da die Standardmathematik an scharfen Kanten versagt (wie dort, wo ein Flügel sich von der Luft löst).

• Kukushkins neue Formel bietet eine mathematische Brücke. Sie verwandelt das „Schummeln" in eine rigorose Berechnung basierend auf ersten Prinzipien. Sie ermöglicht es Computern, Turbulenzen zu modellieren, indem sie diese Fernwirkungen summieren, anstatt nur zu raten.

Zusammenfassung auf den Punkt gebracht

Das Paper argumentiert, dass die alte Art, Flüssigkeitsreibung zu berechnen, wie der Versuch ist, ein Gespräch zu verstehen, indem man nur die Person neben sich anhört. Es verpasst das große Ganze.

Kukushkin hat ein neues Regelbuch geschrieben, das den gesamten Raum anhört. Indem es berücksichtigt, wie sich Störungen durch die gesamte Flüssigkeit ausbreiten (Nichtlokalität), leistet diese neue Mathematik Folgendes:

• Sie behebt ein logisches Paradoxon darüber, wie dick oder dünn eine Flüssigkeit sein sollte.
• Sie erklärt, warum sich Teilchen in einem Sturm so schnell voneinander entfernen.
• Sie bietet einen Weg für Computer, komplexe, chaotische Strömungen (wie Wind um ein Flugzeug oder Wasser in einem Rohr) viel genauer zu simulieren, ohne auf Raten angewiesen zu sein.

Der Autor weist darauf hin, dass dieselbe Logik schließlich auch auf Wärmeübertragung und sogar auf Plasmaphysik angewendet werden könnte, aber die Kernleistung hier ist das Umschreiben der Regeln der Flüssigkeitsreibung, um die „unordentliche" Realität der Turbulenz zu bewältigen.

Technische Zusammenfassung

Technisches Fazit: Verallgemeinerung des viskosen Spannungstensors für nicht-kleine Gradienten

Problemstellung
Der Artikel behandelt fundamentale Einschränkungen in der klassischen Navier-Stokes-Beschreibung hydrodynamischer Turbulenz, insbesondere hinsichtlich des viskosen Spannungstensors ($\sigma_{ik}$). Der konventionelle Tensor, der mittels der Chapman-Enskog-Methode aus der Boltzmann-Gleichung abgeleitet wird, beruht auf der Annahme kleiner Gradienten der hydrodynamischen Geschwindigkeit. Der Autor argumentiert, dass diese Annahme den klassischen Tensor für die Beschreibung komplexer Strömungsphänomene wie tangentialer Diskontinuitäten, abgelöster Strömungen und der in der Turbulenz inhärenten Nichtlokalität unzureichend macht.

Darüber hinaus hebt der Artikel ein theoretisches Paradoxon in der klassischen Herleitung hervor: Der dynamische Viskositätskoeffizient ($\eta$) erweist sich als umgekehrt proportional zum Teilchenstreuquerschnitt ($\sigma_{coll}$), was im Widerspruch zu intuitiven Erwartungen steht. Der Autor führt dieses Paradoxon auf die restriktive „kleine Gradienten"-Bedingung der Chapman-Enskog-Methode zurück, die zunehmend nur dann erfüllbar wird, wenn der Streuquerschnitt zunimmt, wodurch das Verhalten langlebiger lokalisierter Wirbel maskiert wird.

Bestehende Turbulenzmodelle (z. B. Spalart-Allmaras, RANS) greifen häufig auf phänomenologische Korrekturen oder lokale Differentialgleichungen zurück, die die Nichtlokalität der Turbulenz nicht erfassen können, was empirisch durch das $t^3$-Gesetz von Richardson für Paar-Korrelationen in turbulenten Medien belegt wird. Zwar wurden nicht-lokale Transportmodelle (wie Lévy-Flüge/-Pfade) vorgeschlagen, doch fehlt es diesen oft an einer rigorosen Herleitung aus ersten Prinzipien innerhalb der Hydro- und Gasdynamik.

Methodik
Der Autor schlägt eine Verallgemeinerung des viskosen Spannungstensors vor, indem er die Bedingung kleiner Geschwindigkeitsgradienten innerhalb des Chapman-Enskog-Rahmens lockert. Die Methodik verläuft wie folgt:

1. Kinematische Grundlage: Die Herleitung beginnt mit der kinetischen Boltzmann-Gleichung für ein Gas einer einzigen Sorte. Die Verteilungsfunktion $f(\mathbf{p}, \mathbf{r}, t)$ wird in einen Teil des lokalen thermodynamischen Gleichgewichts (Maxwell-Verteilung $f^{(0)}$) und einen Korrekturterm $\varphi$ zerlegt.
2. Lockerung der Gradientenbeschränkungen: Im Gegensatz zum klassischen Ansatz, bei dem räumliche Ableitungen des Korrekturterms vernachlässigt werden, behält der Autor die Ableitungen von $\varphi$ auf der linken Seite der Transportgleichung bei. Dies erkennt an, dass in Scherströmungen und lokalisierten Wirbeln die Gradienten des Korrekturterms erheblich größer sein können als die der glatten Gleichgewichtsverteilung.
3. Formulierung der Transportgleichung: Ein vereinfachtes Stoßintegral (BGK-artig) wird verwendet, was zu einer kinetischen Transfergleichung für $\varphi$ führt. Diese Gleichung weist einen Quellterm auf, der vom Gradienten der mittleren Massengeschwindigkeit abhängt, sowie einen Senkterm, der die Absorption der Träger (inverse freie Weglänge) darstellt.
4. Integral-Lösung: Die Transportgleichung wird gelöst, indem angenommen wird, dass die Bewegung der Träger signifikant schneller ist als die Bewegung des Bulk-Mediums. Die Lösung für $\varphi$ wird als Integral über Trägerbahnen von der Mediumsgrenze zum Beobachtungspunkt ausgedrückt, wobei ein exponentieller Abfallfaktor basierend auf der freien Weglänge einbezogen wird.
5. Herleitung des Spannungstensors: Der viskose Spannungstensor wird berechnet, indem das Produkt aus Korrekturterm, Gleichgewichtsverteilung und Impulsfluss über den Phasenraum integriert wird. Durch die Transformation der Bahnintegration in eine Volumenintegration über das Medium leitet der Autor eine nicht-lokale, integrale Darstellung des Spannungstensors ab.

Hauptbeiträge und Ergebnisse
Das primäre Ergebnis ist die Herleitung eines verallgemeinerten viskosen Spannungstensors $\sigma_{ik}$, der ein Integral über räumliche Koordinaten und kein lokaler Differentialoperator ist.

• Nicht-lokale Darstellung: Der neue Ausdruck für $\sigma_{ik}$ (Gleichung 2.13) hängt explizit von den Geschwindigkeitsgradienten an allen Punkten innerhalb des Mediums ab, gewichtet mit einem Kern, der den Abstand zwischen den Punkten und die freie Weglänge beinhaltet. Diese Struktur spiegelt nicht-lokale Transporttheorien wider, wie etwa das Biberman-Holstein-Modell für den Transfer resonanter Strahlung.
• Wiederherstellung der klassischen Theorie: Im Grenzfall kleiner Weglängen (kleine Gradienten) reduziert sich der Integralausdruck auf den klassischen lokalen viskosen Spannungstensor (Gleichung 1.2), wobei der Viskositätskoeffizient $\eta = P\tau$ beträgt. Dies gewährleistet die Konsistenz mit der etablierten Theorie im entsprechenden Regime.
• Auflösung des Paradoxons: Die verallgemeinerte Form löst die umgekehrte Abhängigkeit der Viskosität vom Streuquerschnitt auf. Der Autor geht davon aus, dass das Paradoxon entsteht, weil die klassische Herleitung kleine Gradienten annimmt, eine Bedingung, die mit den langlebigen lokalisierten Wirbeln unvereinbar ist, die als Träger von Störungen in turbulenten Strömungen fungieren.
• Erhaltungssätze: Die Herleitung zeigt, dass die Korrektur der Verteilungsfunktion die Gesamtzahl der Teilchen im Medium erhält und einen dynamischen Austausch zwischen dem Hauptmedium und dem Ensemble der Störungsträger darstellt.

Bedeutung und Behauptungen
Der Artikel behauptet, dass diese Arbeit einen Schritt hin zu einer nicht-lokalen Theorie der hydrodynamischen Turbulenz darstellt, die aus ersten Prinzipien abgeleitet wurde.

• Modellierung von Diskontinuitäten: Der verallgemeinerte Tensor liefert einen mathematischen Operator, der tangentialen Diskontinuitäten und Ablöseströmungen, bei denen Geschwindigkeitsableitungen unbeschränkt sind, beschreiben kann, und damit eine langjährige Lücke in der numerischen Modellierung schließt.
• Verbindung zum Gesetz von Richardson: Die nicht-lokale Struktur des abgeleiteten Tensors wird als notwendige Grundlage für die numerische Modellierung der Nichtlokalität der Turbulenz präsentiert, speziell um das empirische $t^3$-Gesetz von Richardson für Paar-Korrelationen nachzubilden.
• Validierung des Ansatzes von Kolmogorow: Der Autor schlägt vor, dass diese Verallgemeinerung den Erfolg der dimensionsanalytischen Überlegungen von Kolmogorow bei der Vorhersage von Turbulenzphänomenen erklären könnte, obwohl diese Ergebnisse historisch gesehen schwierig direkt aus der klassischen Navier-Stokes-Gleichung abzuleiten waren.
• Breitere Anwendbarkeit: Die Methodik wird als potenziell anwendbar auf Wärmeübertragung in Gasen und Flüssigkeiten sowie auf komplexere Medien wie Plasma beschrieben und bietet einen einheitlichen Ansatz für nicht-lokale Transportphänomene.

Der Artikel kommt zu dem Schluss, dass der abgeleitete Ausdruck zwar ein signifikanter theoretischer Fortschritt ist, die vollständige Demonstration seiner Fähigkeit, das Richardson-Gesetz nachzubilden, jedoch eine nachfolgende adäquate numerische Modellierung erfordert.