Inkompressible vs. kompressible Navier-Stokes-Gleichungen

Die Navier-Stokes-Gleichungen sind eine Familie von Systemen. Der Unterschied zwischen inkompressibler und kompressibler Strömung ist nicht kosmetisch. Er verändert die Unbekannten, die Mathematik und die offenen Probleme.

Die physikalische Trennung: Dichte, die sich ändert, vs. Dichte, die sich nicht ändert

„Inkompressibel vs. kompressibel“ läuft auf die Dichte hinaus. Bleibt sie konstant, oder ändert sie sich?

Probieren Sie es aus. Füllen Sie eine Spritze mit Wasser und drücken Sie den Kolben. Das Wasser bewegt sich, aber unter Alltagsbedingungen lässt es sich nicht merklich komprimieren. Wasser widersetzt sich der Kompression so stark, dass die Behandlung als inkompressibel eine ausgezeichnete Näherung ist. Inkompressibel. Füllen Sie nun dieselbe Spritze mit Luft und verschließen Sie das Ende. Drücken Sie den Kolben hinein, und Sie spüren, wie die Luft nachgibt: dieselbe Luftmasse ist nun in ein kleineres Volumen gepackt, während sie unter Ihrem Daumen komprimiert wird. Das ist kompressible Strömung.

Bei inkompressibler Strömung ist die Dichte $\rho$ im gesamten Fluid konstant, und jedes winzige Fluidelement behält sein Volumen, während es sich durch den Raum bewegt. Kompressible Strömung ist anders. Die Dichte wird zu einer Variablen, die sich von Ort zu Ort und von Moment zu Moment ändern kann. Luft um ein Strahltriebwerk, Gas bei einer Explosion, die Atmosphäre auf großen Skalen: alles kompressibel, alles von Dichtevariationen getrieben.

Warum ist das wichtig? Weil diese Unterscheidung neu formt, wie die Navier-Stokes-Gleichungen aussehen, was sie vorhersagen und wie schwierig sie zu analysieren und zu lösen sind.

Die Annahme der Inkompressibilität setzt das Dichtefeld $\rho$ im gesamten Strömungsgebiet auf eine positive Konstante. Physikalisch bedeutet dies, dass jedes materielle Volumenelement unter der Strömungsabbildung sein Volumen erhält. Im kompressiblen Fall wird $\rho(x,t)$ zu einer vollwertigen Unbekannten, die von ihrer eigenen Entwicklungsgleichung bestimmt wird.

Das sind nicht zwei Notationen für dasselbe System. Sie unterscheiden sich in der Anzahl der Unbekannten, in der Struktur der Nebenbedingungsgleichungen und im Charakter des Drucks. Das inkompressible System hat $d+1$ skalare Unbekannte ($u$ und $p$ in $\mathbb{R}^d$); das kompressible System hat typischerweise $d+2$ primäre Felder (Geschwindigkeit, Dichte und innere Energie oder Temperatur), wobei der Druck über eine Zustandsgleichung bestimmt wird.

Die Unterscheidung verändert auch den PDE-Typ des Drucks: elliptisch im inkompressiblen Modell, wo der Druck durch eine momentane räumliche Lösung bestimmt wird, gegenüber gemischt hyperbolisch-parabolisch im kompressiblen Fall, wo sich Druckstörungen mit endlicher Geschwindigkeit ausbreiten. Das ist keine technische Kleinigkeit. Es bestimmt, wie sich Information durch das Fluid bewegt.

Die inkompressiblen Navier-Stokes-Gleichungen

Die inkompressiblen Navier-Stokes-Gleichungen beschreiben Fluide, deren Dichte konstant ist. Sie sind die Version, die im Clay-Millennium-Problem erscheint, und die Version, auf die sich diese Website konzentriert.

Das System hat zwei Teile. Die Impulsgleichung:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

und die Inkompressibilitätsbedingung:

$$\nabla \cdot u = 0$$

Die Bedingung $\nabla \cdot u = 0$ besagt, dass das Geschwindigkeitsfeld divergenzfrei ist: Fluid häuft sich nirgends an und verdünnt sich nirgends. Was in ein winziges Gebiet hineinströmt, muss mit derselben Rate wieder herausströmen. Diese eine Bedingung ersetzt die gesamte Dichtegleichung. Die Dichte ändert sich nicht, also braucht man keine Gleichung, um sie zu verfolgen.

Der Druck spielt hier eine besondere Rolle. Er wird nicht durch ein thermodynamisches Gesetz (wie das ideale Gasgesetz) bestimmt. Stattdessen passt er sich überall instantan an, um die Strömung divergenzfrei zu halten. Mathematisch löst $p$ eine Poisson-Gleichung, die aus der Bedingung abgeleitet wird. Druckänderungen breiten sich unendlich schnell aus. In inkompressibler Strömung gibt es keine „Schallgeschwindigkeit“.

Das inkompressible Navier-Stokes-System hat zwei unbekannte Felder: Geschwindigkeit $u$ und Druck $p$. Diese Einfachheit täuscht. Der nichtlineare Term $(u \cdot \nabla)u$ macht das System in drei Dimensionen dennoch extrem schwierig.

Das inkompressible Navier-Stokes-System auf $\mathbb{R}^3$ mit kinematischer Viskosität $\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

Die divergenzfreie Bedingung $\nabla \cdot u = 0$ ist eine punktweise Nebenbedingung, keine Entwicklungsgleichung. Sie kodiert lokale Volumenerhaltung: Die Strömungsabbildung $\Phi_t$ erfüllt $\det(D\Phi_t) = 1$ für alle $t$.

Wendet man den Divergenzoperator auf die Impulsgleichung an und verwendet die Inkompressibilität, erhält man die Druck-Poisson-Gleichung:

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

Dies ist eine elliptische Gleichung für $p$ zu jedem festen Zeitpunkt. Der Druck ist keine unabhängige thermodynamische Variable; in der Standardinterpretation der PDE wirkt er als Lagrange-Multiplikator, der die divergenzfreie Bedingung erzwingt, global und instantan durch das Geschwindigkeitsfeld bestimmt. Information breitet sich über den Druck mit unendlicher Geschwindigkeit aus, ein struktureller Unterschied zum kompressiblen System, der sich nicht übertünchen lässt.

Die Unbekannten sind $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ und $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$. In der Clay-Formulierung (Fefferman, 2000) ist $u_0 \in C^\infty(\mathbb{R}^3)$ divergenzfrei, und die Frage ist, ob $u$ in $C^\infty(\mathbb{R}^3 \times [0,\infty))$ mit beschränkter Energie bleibt.

Die kompressiblen Navier-Stokes-Gleichungen

Die kompressiblen Navier-Stokes-Gleichungen beschreiben Strömungen, bei denen die Dichte variiert. Größeres System. Mehr Unbekannte. Mehr Gleichungen.

Man hat weiterhin eine Impulsgleichung, aber nun tritt die Dichte $\rho$ explizit auf:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

Die Nebenbedingung $\nabla \cdot u = 0$ entfällt. An ihre Stelle tritt eine Kontinuitätsgleichung die verfolgt, wie sich die Dichte entwickelt:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Das besagt, dass Masse erhalten bleibt: Die Dichte ändert sich, weil die Strömung Fluidelemente komprimiert oder expandiert.

Das System benötigt außerdem eine Energiegleichung und eine Zustandsgleichung, eine thermodynamische Beziehung wie $p = \rho R T$ (das ideale Gasgesetz), die den Druck mit Dichte und Temperatur verknüpft. Der Druck ist nicht länger ein passiver Durchsetzer einer Nebenbedingung. Er hat seine eigene Physik, seine eigene Dynamik, und er breitet sich mit endlicher Geschwindigkeit aus: der Schallgeschwindigkeit.

Das kompressible System ist unverzichtbar für Aerodynamik bei hohen Geschwindigkeiten, astrophysikalische Gasdynamik, Verbrennung und jede Strömung, bei der Dichteänderungen eine Rolle spielen. Es ist jedoch ein wirklich anderes mathematisches Objekt als die inkompressiblen Gleichungen. Mehr Unbekannte, mehr Gleichungen, eine völlig andere PDE-Struktur.

Das kompressible Navier-Stokes-System koppelt die Geschwindigkeit $u(x,t)$, Dichte $\rho(x,t)$, Druck $p(x,t)$ und spezifische innere Energie $e(x,t)$ (oder Temperatur $\theta$). In Erhaltungsform:

Kontinuität: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Impuls: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

Energie: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

wobei $E = e + \tfrac{1}{2}|u|^2$ die gesamte spezifische Energie ist, $\tau$ der viskose Spannungstensor ist (für ein Newtonsches Fluid $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$ mit Volumenviskosität $\lambda$), und $\kappa$ die Wärmeleitfähigkeit ist.

Der Abschluss erfordert eine Zustandsgleichung, z. B. $p = (\gamma - 1)\rho e$ für ein ideales Gas mit adiabatischem Exponenten $\gamma$.

Im kompressiblen System breiten sich akustische Störungen mit endlicher Geschwindigkeit aus, mit der charakteristischen Schallgeschwindigkeit $c = \sqrt{\partial p / \partial \rho |_s}$. Dies steht im Gegensatz zur elliptischen, instantanen Druckbestimmung im inkompressiblen Modell. Die kompressiblen Gleichungen erlauben Stoßwellen, Verdünnungswellen und Kontaktunstetigkeiten, die in inkompressiblen Strömungen kein Analogon haben.

Die Mach-Zahl: Wann spielt Kompressibilität eine Rolle?

Wann spielt Kompressibilität eine Rolle? Eine Zahl entscheidet: die Mach-Zahl.

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ ist die Strömungsgeschwindigkeit. $c$ ist die Schallgeschwindigkeit. Ihr Verhältnis sagt aus, wie schnell sich die Strömung im Vergleich zu der Geschwindigkeit bewegt, mit der sich Druckstörungen durch das Medium ausbreiten können; dieser Vergleich bestimmt, ob man Dichteänderungen gefahrlos vernachlässigen kann oder ob sie die Physik dominieren.

Wenn $\text{Ma} < 0.3$, density changes by less than about 5%. Incompressible equations work. Air in a room, water in a pipe, wind around a building: all low-Mach flows where pressure disturbances travel so much faster than the flow itself that density barely budges.

Oberhalb von $\text{Ma} \approx 0.3$ beginnt die Kompressibilität spürbar zu werden, und bei etwa $\text{Ma} \approx 1$ erreicht man den transsonischen Bereich, in dem lokale Überschallzonen auftreten und sich Stoßwellen bilden. Kampfjets. Raketendüsen. Wiedereintretende Raumfahrzeuge.

Kein binärer Schalter. Die meisten alltäglichen Fluidströmungen und das Clay-Millennium-Problem liegen klar im Niedrig-Mach-Bereich, in dem die inkompressiblen Gleichungen gelten.

Die Mach-Zahl $\text{Ma} = |u|/c$ parametrisiert die Bedeutung der Kompressibilität, wobei $c = \sqrt{\partial p / \partial \rho |_s}$ die isentrope Schallgeschwindigkeit ist. Formal entstehen die inkompressiblen Gleichungen als Niedrig-Mach-Grenzwert des kompressiblen Systems.

Die asymptotische Entwicklung in Potenzen von $\text{Ma}^2$ (siehe Klainerman & Majda, 1981, 1982; Schochet, 1986) zeigt, dass für $\text{Ma} \to 0$ bei geeigneten Anfangsdaten die kompressiblen Lösungen gegen die inkompressible Lösung konvergieren. In üblichen Niedrig-Mach-Nondimensionalisierungen mit gut vorbereiteten Daten schreibt man den Druck oft als nahezu räumlich uniformen thermodynamischen Hintergrund plus eine kleinere dynamische Korrektur, die im Grenzwert die Inkompressibilitätsbedingung durchsetzt.

Dies ist ein singulärer Grenzwert: Die Schallgeschwindigkeit $c \to \infty$ und der hyperbolische Charakter des kompressiblen Systems degeneriert zur elliptischen Druckgleichung der inkompressiblen Strömung. Die akustischen Moden werden unendlich schnell und entkoppeln sich von der wirbelbehafteten Dynamik.

Der Bereich $\text{Ma} < 0.3$ is an engineering heuristic. It reflects the empirical observation that relative density variation $\delta\rho / \rho \sim \text{Ma}^2 / 2$ stays below $\sim$5% in this range. The mathematical justification is the convergence theorem for the low-Mach limit, which requires well-prepared initial data and compatible boundary conditions.

Warum das Millennium-Problem den inkompressiblen Fall betrifft

Das Clay-Millennium-Problem stellt eine präzise Frage: Gegeben eine glatte, divergenzfreie Anfangsgeschwindigkeit auf $\mathbb{R}^3$, erzeugt das inkompressible Navier-Stokes-System dann stets eine glatte Lösung, die für alle Zeiten existiert?

Warum speziell inkompressibel? Drei Gründe.

Erstens ist es ohnehin schon schwierig genug. Die inkompressiblen 3D-Gleichungen widersetzen sich seit Lerays grundlegender Arbeit von 1934 einem Beweis der globalen Regularität. Variable Dichte, Thermodynamik und Stoßwellen hinzuzufügen, würde das Problem enorm erschweren, nicht leichter behandelbar machen.

Zweitens liegt die Schwierigkeit in der reinen Strömungsmechanik. Das inkompressible System isoliert die zentrale mathematische Herausforderung, das Zusammenspiel von nichtlinearer Advektion $(u \cdot \nabla)u$ und viskoser Dissipation $\nu \Delta u$, ohne thermodynamische oder akustische Komplikationen. Es ist die klarste Umgebung, um die Regularitätsfrage zu stellen.

Drittens ist die Physik klar. Die inkompressiblen Gleichungen modellieren die häufigsten Alltagsströmungen. Ob sie aus glatten Daten Singularitäten erzeugen können, ist eine grundlegende Frage nach der mathematischen Konsistenz der klassischen Strömungsmechanik.

Das kompressible System hat eigene tiefgehende offene Probleme (Existenz globaler Lösungen für große Daten, Bildung und Wechselwirkung von Stoßwellen), aber das sind andere Probleme mit anderen Strukturen. Der Clay-Preis zielt auf den inkompressiblen Fall, weil dies die spezifische Regularitätsfrage ist, die Fefferman für die 3D-Navier-Stokes-Gleichungen formuliert hat.

Die offizielle Clay-Formulierung (Fefferman, 2000) spezifiziert das inkompressible System auf $\mathbb{R}^3$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

mit $u_0 \in C^\infty(\mathbb{R}^3)$ divergenzfrei, und die Frage ist, ob $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ mit $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ beschränkt für alle $t \geq 0$.

Die Wahl des inkompressiblen Systems ist mathematisch motiviert. Das zentrale offene Problem, die Lücke zwischen schwachen Leray-Hopf-Lösungen (die global existieren, aber möglicherweise nicht eindeutig oder glatt sind) und klassischen glatten Lösungen (die lokal existieren, aber einen Blow-up entwickeln können), ist spezifisch für die inkompressiblen 3D-Gleichungen. In 2D ist globale Regularität für glatte inkompressible Navier-Stokes-Lösungen bekannt; der 3D-Fall bleibt offen.

Das kompressible System führt qualitativ andere Schwierigkeiten ein: Stoßbildung (die sogar für Euler-Gleichungen mit glatten Daten auftritt), Vakuumzustände ($\rho \to 0$), und die Kopplung zwischen Wirbelstärke und akustischen Moden. Das sind wichtige offene Probleme, aber sie unterscheiden sich strukturell von der inkompressiblen Regularitätsfrage.

Das inkompressible Problem isoliert die Konkurrenz zwischen der energie-superkritischen Nichtlinearität und viskoser Dissipation. Die natürliche Energieabschätzung liefert $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, was in 3D um eine halbe Ableitung hinter dem skalierungskritischen Raum $L^\infty_t \dot{H}^{1/2}_x$ zurückbleibt. Diese Lücke zu schließen oder zu beweisen, dass sie nicht geschlossen werden kann, ist der Kern des Millennium-Problems.

Was Sie als Nächstes lesen sollten

Beginnen Sie hier. Möchten Sie, dass jeder Term im inkompressiblen System einzeln zerlegt wird, wobei die physikalische Bedeutung und die mathematische Rolle jedes Bestandteils von Grund auf erklärt werden? Was sind die Navier-Stokes-Gleichungen?

Woher kommt dieses System? Herleitung der Navier-Stokes-Gleichungen.

Lässt man die Viskosität weg, erhält man die Euler-Gleichungen, die ein Jahrhundert älter sind, auf dem Papier einfacher aussehen und in mancher Hinsicht mathematisch noch schwerer zu verstehen sind, weil der glättende Effekt des Diffusionsterms verloren geht. Euler vs. Navier-Stokes.

Der Preis. Das Existenz- und Glattheitsproblem der Navier-Stokes-Gleichungen.

Empfohlene Wege von hier aus: