Warum das Navier-Stokes-Problem schwierig ist

Viele der Gleichungen, denen man in der Physik zuerst begegnet, sind linear: Verdoppelt man die Eingabe, verdoppelt sich die Antwort. Navier-Stokes ist nicht so.

Navier-Stokes? Nichtlinear. Die Geschwindigkeit des Fluids beeinflusst ihre eigene Änderungsrate, was bedeutet, dass das Fluid sich selbst antreibt. Stellen Sie sich vor, vorhersagen zu wollen, wohin eine Menschenmenge gehen wird, wenn die Bewegung jeder einzelnen Person davon abhängt, was alle um sie herum tun, und das, was diese Personen tun, wiederum von allen um sie herum abhängt, sich also endlos nach außen fortsetzt. Genau auf diese Situation blickt man hier.

Der Schuldige ist der Selbstwechselwirkungsterm $(u \cdot \nabla)u$. Er erzeugt Rückkopplungsschleifen, in denen kleine Störungen zu großen anwachsen, und er ist der Grund, warum Turbulenz in Fluiden so außerordentlich komplex ist (siehe Teilprobleme für mehr dazu).

Die konvektive Nichtlinearität $(u \cdot \nabla)u$ ist das fundamentale Hindernis. In der Wirbelformulierung $\omega = \nabla \times u$ wird die Gleichung zu

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$

Der Wirbelstreckungsterm $(\omega \cdot \nabla)u$ hat kein Vorzeichen. Er kann die Vortizität unbegrenzt verstärken. In 2D verschwindet dieser Term (da $\omega$ ein Skalar senkrecht zur Strömung ist), weshalb die globale Regularität in 2D bekannt ist (Ladyzhenskaya, 1969). In 3D ist die Wirbelstreckung der wichtigste Kandidatenmechanismus für Blow-up in endlicher Zeit.

Hier liegt der Kern: Die Nichtlinearität ist quadratisch in $u$. Die $H^1$-Energieabschätzung liefert $\|\nabla u\|_{L^2}$, aber die Kontrolle von $(u \cdot \nabla)u$ in $L^2$ erfordert typischerweise stärkere Informationen wie $u \in L^\infty$ zusammen mit $\nabla u \in L^2$, oder eine äquivalente Kontrolle auf kritischer Skala; die Energieklasse allein liefert dies nicht.

Die Navier-Stokes-Gleichungen besitzen eine Skalierungssymmetrie. Vergrößert man eine Lösung im Nahbereich, macht alles im richtigen Maß kleiner und schneller, erhält man wieder eine vollkommen gültige Lösung. Diese Symmetrie ist mathematisch natürlich, analytisch jedoch gefährlich.

Warum? Die einzige Größe, die wir zuverlässig kontrollieren können, ist die gesamte Energie des Fluids, und sie liegt auf völlig der falschen Skala: Sie sagt uns etwas über das Gesamtbild, aber absolut nichts darüber, was auf den mikroskopischen Skalen geschieht, auf denen sich ein Blow-up tatsächlich bilden würde.

Das ist, als würde man den gesamten Stromverbrauch einer Stadt überwachen, um einen einzelnen Funken zu entdecken. Nützlich? Sicher. Fein genug aufgelöst? Nicht annähernd. Diese Lücke ist das ganze Problem.

Unter der natürlichen Skalierung $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$ ist der kritische Sobolev-Raum $\dot{H}^{1/2}(\mathbb{R}^3)$ (äquivalent $L^3$). Eine Norm ist skalenkritisch, wenn sie unter dieser Transformation invariant ist. Wird eine kontrollierte Norm beim Hineinzoomen kleiner, ist die Gleichung relativ zu dieser Kontrolle superkritisch: Die Schranke liegt unterhalb der Skala, die nötig wäre, um Konzentration zu sehen.

Die Energieungleichung liefert Kontrolle über $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. Auf dem Energieniveau ist das 3D-Problem superkritisch:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

Für $\lambda \to \infty$ (Hineinzoomen) werden diese Schranken schwächer relativ zu den invarianten Größen $\|u\|_{L^3}$ und $\|u\|_{\dot{H}^{1/2}}$. Der nichtlineare und der viskose Term skalieren gemeinsam; das Problem ist, dass die bekannte A-priori-Abschätzung unterhalb der kritischen Skala liegt. Von der Energieklasse zu einer kritischen Norm zu gelangen, ist die zentrale Schwierigkeit.

Betrachten Sie einen Fluss. Die Fluidbewegung wird chaotisch. Turbulent. Große Wirbel zerfallen in kleinere, diese in noch kleinere, und so kaskadiert es bis hinab zu mikroskopischen Skalen, auf denen die Viskosität schließlich alles glättet.

Kolmogorov beschrieb diese Energiekaskade im Jahr 1941, und die Navier-Stokes-Gleichungen erfassen sie auf elegante Weise. Doch Folgendes lässt Forschende nachts wachliegen: Was, wenn sich Energie schneller in immer kleineren Bereichen konzentriert, als die Viskosität sie dissipieren kann? Das ist ein Blow-up.

Kann das tatsächlich passieren? Oder gewinnt die Viskosität immer? Das ist die offene Frage, Punkt. Zur physikalischen Brücke von der Reynolds-Zahl zu diesem Bild kleiner Skalen siehe Reynolds-Zahl, Turbulenz und warum kleine Skalen wichtig sind.

Kolmogorovs K41-Theorie sagt ein Energiespektrum $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ im Inertialbereich $k_f \ll k \ll k_\eta$ voraus, wobei $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ die Kolmogorov-Dissipationswellenzahl ist. Konstanter Energiefluss über die Skalen hinweg. Das ist das klare Bild.

Die Regularitätsfrage stellt etwas Düstereres: Kann diese Kaskade entarten? Kann die Dissipationsskala $k_\eta^{-1}$ in endlicher Zeit gegen null schrumpfen, mit $\|\nabla u\|_{L^2} \to \infty$ während die Gesamtenergie endlich bleibt?

Onsagers Dissipationsanomalie-Vermutung (1949) besagt in gewissem Sinne ja: Im Grenzfall verschwindender Viskosität $\nu \to 0$ bleibt Energiedissipation bestehen, und schwache Lösungen der Euler-Gleichungen können Energie dissipieren. Bestätigt für Hölder-Exponenten unter $1/3$ (Isett, 2018; Buckmaster et al., 2018). Doch was das für die Regularität der Navier-Stokes-Gleichungen bedeutet, bleibt wirklich unklar. Für eine Intuition auf Ebene der Regime siehe Reynolds-Zahl, Turbulenz und warum kleine Skalen wichtig sind.

Druck in den Navier-Stokes-Gleichungen ist eigenartig. Er ist überhaupt keine unabhängige Variable; die Geschwindigkeit bestimmt ihn vollständig über eine einzige Nebenbedingung: Die Flüssigkeit ist inkompressibel, kann also nicht zusammengedrückt werden.

Das macht den Druck nichtlokal. Im inkompressiblen Modell beeinflusst eine Änderung der Geschwindigkeit in einem Bereich das Druckfeld global, weil der Druck zu diesem Zeitpunkt durch das gesamte Geschwindigkeitsfeld bestimmt wird.

Für die Analysis ist das verheerend. Man kann nicht untersuchen, was an einem einzelnen Punkt geschieht, ohne die gesamte Flüssigkeit zugleich zu berücksichtigen. Lokales Argumentieren? Vergessen Sie es.

Die Inkompressibilitätsbedingung $\nabla \cdot u = 0$ bestimmt den Druck durch die Poisson-Gleichung

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j)$$

also $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$, unter Beteiligung von Riesz-Transformationen (singulären Integraloperatoren). Der Druck ist eine nichtlokale Funktion der Geschwindigkeit, und diese Nichtlokalität ist das zentrale Hindernis für punktweise oder räumlich lokale Abschätzungen.

Standardargumente mit dem Maximumprinzip versagen hier: Auch wenn der viskose Term $\nu \Delta u$ dissipativ ist, kann der Druckgradient $-\nabla p$ Energie aus entfernten Bereichen konzentrieren. Die Caffarelli-Kohn-Nirenberg-Theorie behandelt dies über lokale Energieungleichungen auf parabolischen Zylindern, doch aus ihnen punktweise Regularität zu gewinnen bleibt der schwierige Schritt.

Für die 2D-inkompressiblen Navier-Stokes-Gleichungen in den Standardkonfigurationen sind globale glatte Lösungen bekannt; dies wurde in klassischen Arbeiten, unter anderem von Ladyzhenskaya (1969), gezeigt. Siehe warum 2D einfacher ist.

Drei Dimensionen? Alles bricht auseinander, und der Grund läuft auf einen Mechanismus hinaus: Wirbelstreckung. In 2D können Wirbel rotieren und verschmelzen, aber sie können sich nicht strecken. In 3D kann die Flüssigkeit Wirbelröhren erfassen und sie immer dünner und dünner ziehen, sodass sich potenziell der letzte Rest Energie in einem unendlich dünnen Filament konzentriert.

Kann diese Streckung in endlicher Zeit ins Unendliche ausbrechen, oder greift die Viskosität immer ein? Das ist die Millionen-Dollar-Frage. Buchstäblich.

Die Dichotomie ist scharf.

2D: Die Vortizität $\omega$ ist ein Skalar, der $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$ erfüllt. Das Maximumprinzip liefert $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$, BKM impliziert globale Regularität, und der Wirbelstreckungsterm $(\omega \cdot \nabla)u$ ist identisch null.

3D: Die Vortizität $\omega \in \mathbb{R}^3$ erfüllt $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$, wobei der Streckungsterm $(\omega \cdot \nabla)u$ $|\omega|$ superlinear verstärken kann (formal $\sim |\omega|^2$ über Biot-Savart). Kein Maximumprinzip verfügbar.

Die Enstrophie $\|\omega\|_{L^2}^2$ erfüllt

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

Die Kontrolle von $\|\nabla u\|_{L^\infty}$ erfordert $\omega \in L^\infty$, was die Kontrolle von $\|\nabla u\|_{L^\infty}$ erfordert. Zirkulär. Keine bestehende Technik hat diese Schleife durchbrochen.

Dieser Artikel ist Teil von Das Problem.

Diese Hindernisse haben Mathematiker dazu veranlasst, das Problem in Teilprobleme zu zerlegen und spezialisierte Ansätze für jedes davon zu entwickeln.

Zum Kontext, warum der viskose Term hilft, aber nicht ausreicht, siehe Euler vs. Navier-Stokes.

Dieser Artikel ist Teil von Das Problem.

Das Regularitätsproblem zerfällt in handhabbare Teilstücke; siehe Teilprobleme. Die analytischen Werkzeuge, die entwickelt wurden, um diese Hindernisse anzugehen, werden behandelt in Ansätze. Warum hilft der viskose Term, schließt die Lücke aber dennoch nicht? Euler vs. Navier-Stokes.