Theorema Magnum MMXIX: die Onsager-Vermutung

Die Euler-Gleichungen

beschreiben die Strömung von reibungsfreien, elastischen Flüssigkeiten und Gasen. Sie sind der Grenzfall für Viskosität 0 der Navier-Stokes-Gleichungen

.

Olga Ladyzhenskaya bewies 1959 die globale eindeutige Lösbarkeit und die Glattheit der Lösungen für die Navier-Stokes-Gleichungen auf dem R² und dem 2-dimensionalen Torus, und auch für die schwierigeren Euler-Gleichungen. Dieselbe Frage für die 3-dimensionalen Navier-Stokes-Gleichungen ist bis heute offen und ist eines der vom Clay-Institut aufgestellten Milleniumprobleme. Numerische Simulationen sind oft sehr instabil und Terence Tao bewies 2014, dass jedenfalls die Lösungen einer gemittelten Navier-Stokes-Gleichung explodieren. Das zeigte insbesondere, dass verschiedene bis dahin verfolgte Ansätze zum Beweis der Regularität der Lösungen von vornherein nicht funktionieren können, weil sie auch Regularität der Lösungen für die gemittelte Gleichung geben würden.

Ein Ansatz zur Bearbeitung unhandlicher Differentialgleichungen ist die Addition zusätzlicher stochastischer Terme und Behandlung der so entstandenen stochastischen Differentialgleichungen. Zum Beispiel konnten Hairer und Mattingly für die 2-dimensionalen Navier-Stokes-Gleichungen auf dem Torus diejenigen Störungen klassifizieren für die die Gleichungen ergodisch. Hairer entwickelte dann eine allgemeine Theorie von “Regularitätsstrukturen”, mit der man Gleichungen wie die KPZ-Gleichung oder die Navier-Stokes-Gleichungen auch bei Anfangswerten oder Störtermen niedriger Regularität als stochastische Differentialgleichungen behandeln konnte.

Für die Navier-Stokes-Gleichungen hat Kolmogorow 1941 eine statistische Theorie der Turbulenz vorgeschlagen. Für die kinetische Gesamtenergie würde man erwarten, dass für Viskosität ν>0 die Dissipation dE/dt verschwindet. Das widersprach aber physikalischen und numerischen Experimenten. Gemäß Kolmogorows Theorie sollte für ν>0 die Dissipation konstant (und positiv) sein, was durch einen stetigen Fluß von Energie von niedrigen zu hohen Frequenzen (sogenannten Kaskaden) bewirkt werden sollte.
Lars Onsager, ein bekannter Physikochemiker hatte 1949 auf eine Konsequenz aus Kolmogorows Theorie hingewiesen. Weil man im Grenzfall für ν gegen 0 die Euler-Gleichung für reibungsfreie Flüssigkeiten und Gase die vieluntersuchten Gleichungen bekommt, sollte Kolmogorows Theorie und sein 5/3-Gesetz eine gewisse Regularität und insbesondere anormale Dissipation für die schwachen Lösungen der Euler-Gleichung implizieren. Während Kolmogorows Theorie eine statistische, mathematisch schwer zu fassende war, handelte es sich bei der Onsager-Vermutung um ein mathematisches Problem für partielle Differentialgleichungen, das sich beweisen oder widerlegen ließ.

Michail Gromov hatte 1969 in seiner Dissertation ein sogenanntes h-Prinzip entwickelt: einen allgemeinen Ansatz, mit dem man für eine Reihe in der Geometrie vorkommender partieller Differential(un)gleichungen ihre Lösbarkeit entscheiden kann. Als Spezialfall seines h-Prinzips ergab sich auch ein bekannter Satz von Nash und Kuiper über die Existenz 1-mal differenzierbarer, isometrischer Einbettungen des flachen Torus in den euklidischen R3 – im Gegensatz zur mit Gauß’ Theorema Egregium leicht zu beweisenden Unmöglichkeit 2-mal differenzierbarer, isometrischer Einbettungen. Eine der von Gromov entwickelten Methoden zum Beweis von h-Prinzipien war die sogenannte “konvexe Integration”, die im Prinzip auch ein iteratives Verfahren zur Konstruktion der Lösungen liefert. Trotzdem hatte wohl 40 Jahre lang niemand versucht, eine 1-mal differenzierbare Einbettung des Torus in den R3 auch explizit zu konstruieren. Erst 2012 wurde von Borrelli-Jabrane-Lazarus-Thibert mittels konvexer Integration ein Algorithmus zu ihrer expliziten Konstruktion entwickelt. So sahen die ersten vier Iterationsschritte der Konstruktion aus:

Weil der eingebettete Torus nur einmal, nicht zweimal differenzierbar sein darf, sehen die Bilder zwangsweise etwas ‘fraktal’ aus. Die Autoren schrieben, dass in einem präziseren Sinne die Normalvektoren des Torus fraktales Verhalten zeigen.

Jean Leray wird nachgesagt, er habe manche Zeit damit verbracht, Strudel und Wirbel in der Seine an den Pfeilern des Pont Neuf zu beobachten, und sei so zur Suche nach nichtglatten Lösungen der die Strömung linear-viskoser Flüssigkeiten beschreibenden Navier-Stokes-Gleichungen motiviert worden. Eine schwache Lösung der Navier-Stokes-Gleichungen ist eine L²-Funktion, deren schwache Ableitung (im Sinne der später von Schwartz definierten Distributionen) wieder L² ist und die die schwache Version der Gleichungen erfüllt, also

für alle -Funktionen mit kompaktem Träger und verschwindender Divergenz. (Für differenzierbare f läßt sich diese schwächere Gleichung mittels partieller Integration aus der “richtigen” Gleichung herleiten.) Aus dieser Gleichung folgt die Energiegleichung und daraus wiederum, dass die Summanden auf der linken Seite (die kinetische Energie und die Dissipationsenergie) für alle Zeiten beschränkt bleiben. Die Energiegleichung wurde zum entscheidenden Argument für die Konstruktion globaler schwacher Lösungen.
Leray konstruierte schwache Lösungen der Navier-Stokes-Gleichungen in einer 1934 veröffentlichten Arbeit. Er dachte sich diese schwachen Lösungen als Fortsetzung der glatten Lösung über die Singularitäten hinaus und bezeichnete sie als “turbulente Lösungen”: ihre Singularitäten sollten die Turbulenz beschreiben.