Exakte Lösungen der Navier-Stokes-Gleichungen

Die Navier-Stokes-Gleichungen sind berüchtigt nichtlinear. Wie kann sie also jemand exakt lösen?

Symmetrie. Das ist der ganze Trick.

Wenn die Geometrie einer Strömung einfach genug ist (ein gerades Rohr, zwei ebene Platten, eine unendliche Ebene), kann die Geschwindigkeit nur in eine Richtung zeigen und entlang einer oder zweier Koordinaten variieren. In vielen dieser symmetrischen Anordnungen verschwindet der nichtlineare Term entweder oder vereinfacht sich so stark, dass die Gleichungen auf eine lineare partielle Differentialgleichung kollabieren. In stationären Fällen bleibt oft eine gewöhnliche Differentialgleichung übrig, die man mit Papier und Bleistift lösen kann.

Hier ist die Intuition. Die Navier-Stokes-Gleichungen beschreiben alle möglichen Fluidbewegungen. Wenn man das Fluid aber in eine sehr geordnete Situation zwingt, ohne Wirbeln, ohne Chaos, alles in eine Richtung marschierend, wird der größte Teil der Komplexität der Gleichung irrelevant. Der schwierige Teil von Navier-Stokes ist die Rückkopplungsschleife, in der das Fluid sich selbst umherschiebt. In diesen symmetrischen Strömungen gibt es nichts, wogegen es schieben könnte. Der nichtlineare Selbstadvektionsterm fällt einfach weg.

Diese exakten Lösungen sind keine Kuriositäten. Sie sind die Grundlage der Ausbildung in der Strömungsmechanik, die Benchmarks für numerische Codes und der Ausgangspunkt, um zu verstehen, wann und wie reale Strömungen aus dem Ruder laufen.

Eine exakte Lösung der Navier-Stokes-Gleichungen ist ein Geschwindigkeitsfeld $u$ und ein Druck $p$ die erfüllen

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0$$

in geschlossener Form, ohne numerische Approximation.

Der zentrale Mechanismus ist das Verschwinden oder die Vereinfachung des nichtlinearen Advektionsterms $(u \cdot \nabla)u$. Für Parallelströmungen (Strömungen, bei denen $u = (u(y,t),\, 0,\, 0)$ in kartesischen Koordinaten oder $u = (0,\, 0,\, u(r,t))$ in zylindrischen) ist das Geschwindigkeitsfeld automatisch divergenzfrei, und der Advektionsterm verschwindet identisch, weil die Geschwindigkeit keine Komponente in Strömungsrichtung hat und entlang dieser Richtung nicht variiert. Die Impulsgleichung reduziert sich dann auf eine lineare partielle Differentialgleichung oder, für stationäre Strömungen, auf eine gewöhnliche Differentialgleichung.

Wenn die Strömung außerdem stationär ($\partial_t u = 0$), bleibt übrig

$$\nu \Delta u = \nabla p,$$

eine Poisson-Gleichung, deren Lösungen die klassischen exakten Lösungen viskoser Strömung sind: Poiseuille-Strömung, Couette-Strömung und ihre Verwandten.

Für instationäre Parallelströmungen ($\partial_t u \neq 0$ aber die Advektion verschwindet weiterhin) wird die Gleichung zu einer Diffusionsgleichung $\partial_t u = \nu\, \partial_{yy} u$, die das erste und zweite Stokes-Problem liefert.

Poiseuille-Strömung (auch Hagen-Poiseuille-Strömung genannt) ist die wichtigste exakte Lösung der Navier-Stokes-Gleichungen, und sie ist diejenige, der die meisten Ingenieure zuerst begegnen.

Stellen Sie sich Wasser vor, das stationär durch ein langes, gerades Rohr strömt. Eine konstante Druckdifferenz zwischen den beiden Enden treibt die Strömung an. Die Rohrwände bewegen sich nicht, daher haftet das Fluid, das die Wand berührt, bei Geschwindigkeit null fest (das ist die Haftbedingung). Weiter von den Wänden entfernt wird das Fluid schneller. In der Mitte bewegt es sich am schnellsten.

Das Geschwindigkeitsprofil? Eine Parabel. Null an der Wand. Maximum in der Mitte. Dazwischen eine glatte Kurve. Schneiden Sie das Rohr auf und betrachten Sie den Querschnitt: Er sieht aus wie eine umgedrehte Schüssel.

Die zentrale quantitative Tatsache ist diese: Der gesamte Volumenstrom skaliert mit der vierten Potenz des Rohrradius. Vierte Potenz. Verdoppelt man den Radius, erhält man nicht den doppelten Durchfluss, nicht einmal den vierfachen. Man erhält den sechzehnfachen Durchfluss. Das ist das Hagen-Poiseuille-Gesetz, und es erklärt, warum schon eine winzige Verengung einer Arterie die Blutversorgung abwürgen kann.

Annahmen: Die Poiseuille-Strömung setzt voraus, dass das Fluid inkompressibel und newtonsch ist (konstante Viskosität), dass die Strömung stationär und voll entwickelt ist (also nicht noch vom Eintritt her beschleunigt) und dass die Strömung laminar ist. Glatt. Geordnet. In der Praxis geht Rohrströmung bei einer Reynolds-Zahl von ungefähr 2.300 in Turbulenz über; das ist eine empirische Beobachtung, die noch niemand aus der zugrunde liegenden Theorie herleiten konnte.

Betrachten wir eine stationäre, voll entwickelte, axialsymmetrische Strömung in einem Kreisrohr mit Radius $R$, angetrieben durch einen gleichförmigen Druckgradienten $dp/dx < 0$ in the axial direction $x$. In cylindrical coordinates $(r, \theta, x)$ the velocity field has the form $u = (0,\, 0,\, u(r))$.

Die Inkompressibilitätsbedingung $\nabla \cdot u = 0$ ist automatisch erfüllt. Der Advektionsterm verschwindet, weil $u$ nicht von $x$ abhängt. Die axiale Impulsgleichung reduziert sich auf

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

oder äquivalent

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

Da $dp/dx$ konstant ist, ist dies eine gewöhnliche Differentialgleichung in $r$. Zweifache Integration mit den Randbedingungen $u(R) = 0$ (Haftbedingung) und $du/dr|_{r=0} = 0$ (Symmetrie) ergibt das parabolische Geschwindigkeitsprofil:

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

Die maximale Geschwindigkeit auf der Mittellinie ist $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$, und die mittlere Geschwindigkeit ist $\bar{u} = u_{\max}/2$.

Integration über den Querschnitt liefert das Hagen-Poiseuille-Gesetz für den Volumenstrom:

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

wobei $\Delta p > 0$ der Druckabfall über die Rohrlänge $L$.

Die $R^4$-Abhängigkeit ist das bestimmende Merkmal: Kleine Änderungen des Radius bewirken große Änderungen der Durchflussrate. In der Hämodynamik erklärt dies, warum der Blutfluss so empfindlich auf arterielle Stenosen reagiert.

Gültigkeit: Diese Lösung gilt für inkompressible, newtonsche, stationäre, voll entwickelte, laminare Strömung. Der Übergang zur Turbulenz tritt experimentell bei $\text{Re} = \bar{u} D / \nu \approx 2{,}300$ auf, wobei $D = 2R$ der Rohrdurchmesser ist. Diese Schwelle ist eine empirische Beobachtung; es gibt keinen Satz, der sie aus den Navier-Stokes-Gleichungen vorhersagt. Eine Diskussion der Reynolds-Zahl und des laminar-turbulenten Übergangs findet sich unter Reynolds-Zahl, Turbulenz und warum kleine Skalen wichtig sind.

Couette-Strömung ist die exakte Lösung für eine Flüssigkeit, die zwischen zwei parallelen Platten eingeschlossen ist, wenn sich eine Platte bewegt und die andere ruht. Sie ist eine der einfachsten exakten Lösungen der Strömungsmechanik.

Stellen Sie sich einen Kartenstapel vor, der flach auf einem Tisch liegt. Ziehen Sie die oberste Karte seitwärts, und die darunterliegenden Karten verschieben sich ebenfalls, jede ein wenig weniger als die darüberliegende. Genau das ist es. Die Geschwindigkeit variiert linear von null an der unteren Platte bis zur Geschwindigkeit der oberen Platte, und mehr steckt nicht dahinter.

Ein geradliniges Geschwindigkeitsprofil. Untere Platte: stationär. Obere Platte: bewegt. Alles dazwischen interpoliert einfach linear, kein Druckgradient erforderlich, kein komplizierter Aufbau; die Bewegung wird rein durch den bewegten Rand angetrieben, der die Flüssigkeit mitzieht.

Interessanter wird es, wenn zusätzlich ein Druckgradient entlang des Kanals angelegt wird, denn dann kombiniert man schergesteuerte und druckgetriebene Strömung in demselben Spalt. Das Geschwindigkeitsprofil verformt sich zu einer Parabel, die dem linearen Profil überlagert ist, manchmal ebene Poiseuille-Couette-Strömung genannt, und je nach Stärke des Druckgradienten relativ zur Plattengeschwindigkeit kann sogar Rückströmung nahe einer Wand auftreten.

Entfernt man die bewegte Platte vollständig, lässt beide Wände stationär und überlässt einer Druckdifferenz die gesamte Arbeit, dann erhält man ebene Poiseuille-Strömung, das Flachplatten-Analogon zur Rohrströmung. Parabolisch. Am schnellsten in der Mitte. Null an beiden Wänden.

Betrachten Sie eine stationäre Strömung zwischen zwei unendlichen parallelen Platten, die durch einen Spalt $h$ getrennt sind. Sei $y$ die Koordinate senkrecht zu den Platten, mit der unteren Platte bei $y = 0$ und der oberen Platte bei $y = h$. Die Strömung ist parallel: $u = (u(y),\, 0,\, 0)$.

Einfache Couette-Strömung. Die obere Platte bewegt sich mit der Geschwindigkeit $U$ in $x$-Richtung; die untere Platte ist stationär; $dp/dx = 0$. Die Impulsgleichung reduziert sich auf $d^2u/dy^2 = 0$, woraus folgt

$$u(y) = U\frac{y}{h}.$$

Dies ist die einfachste nichttriviale exakte Lösung der Navier-Stokes-Gleichungen: ein lineares Geschwindigkeitsprofil, das vollständig durch die Randbedingung angetrieben wird.

Ebene Poiseuille-Strömung. Beide Platten sind stationär; ein konstanter Druckgradient $dp/dx < 0$ drives the flow. The momentum equation becomes

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

mit $u(0) = u(h) = 0$. Die Lösung ist

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

ein parabolisches Profil, symmetrisch zur Mittellinie $y = h/2$.

Allgemeine Couette-Poiseuille-Strömung. Die Kombination einer bewegten oberen Platte mit einem Druckgradienten ergibt die Überlagerung

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

Der erste Term ist der schergesteuerte Beitrag; der zweite ist der druckgetriebene Beitrag. Je nach Vorzeichen und Betrag von $dp/dx$ relativ zu $\mu U / h^2$ kann das Profil monoton sein, ein inneres Maximum besitzen oder sogar Rückströmung nahe einer Wand aufweisen.

Poiseuille- und Couette-Strömungen sind stationär. Nichts ändert sich mit der Zeit. Stokes-Probleme sind die einfachsten instationären exakten Lösungen, und sie zeigen etwas Schönes: Viskosität lässt Impuls durch eine Flüssigkeit diffundieren, so wie Wärme durch einen Festkörper diffundiert.

Stokes’ erstes Problem (auch Rayleigh-Problem genannt): Stellen Sie sich ein ausgedehntes, ruhendes Flüssigkeitsbecken über einer flachen Platte vor. Zum Zeitpunkt null beginnt die Platte plötzlich, mit konstanter Geschwindigkeit seitwärts zu gleiten. Die Flüssigkeit direkt neben der Platte wird sofort mitgezogen, aber weiter entfernte Flüssigkeit braucht Zeit, um dies zu bemerken. Eine glatte Grenzschicht wächst von der Platte nach außen und wird mit der Zeit dicker.

Die Geschwindigkeit in einer beliebigen Höhe über der Platte hängt vom Verhältnis dieser Höhe zu einer charakteristischen Diffusionslänge $\sqrt{\nu t}$ ab, wobei $\nu$ die Viskosität und $t$ die verstrichene Zeit ist. Viskosere Flüssigkeit? Die Bewegung breitet sich schneller nach oben aus.

Stokes' zweites Problem: derselbe Aufbau, aber nun oszilliert die Platte sinusförmig hin und her, statt sich mit konstanter Geschwindigkeit zu bewegen. Die Oszillation dringt nur eine endliche Strecke in das Fluid ein. Weiter oben bemerkt das Fluid sie kaum. Die Amplitude der Bewegung nimmt mit der Höhe exponentiell ab und erzeugt eine dünne oszillatorische Grenzschicht. Das ist der Mechanismus hinter oszillatorischen Grenzschichten: Eine oszillierende Platte versetzt das nahe Fluid in Bewegung, aber die Störung klingt mit dem Abstand von der Platte exponentiell ab.

Beide Stokes-Probleme betreffen ein halbunendliches Fluid ($y > 0$) über einer ebenen Platte bei $y = 0$. Die Strömung ist parallel: $u = (u(y,t),\, 0,\, 0)$. Der Advektionsterm verschwindet, und die maßgebliche Gleichung ist die eindimensionale Diffusionsgleichung:

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

Stokes' erstes Problem (Rayleigh-Problem). Anfangsbedingung: $u(y,0) = 0$ für $y > 0$. Randbedingung: $u(0,t) = U$ für $t > 0$. Fernfeld: $u \to 0$ für $y \to \infty$.

Die Ähnlichkeitsvariable $\eta = y / (2\sqrt{\nu t})$ reduziert die PDE auf eine ODE. Die Lösung ist

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

wobei $\operatorname{erfc}$ die komplementäre Fehlerfunktion ist. Die Grenzschichtdicke wächst wie $\delta \sim \sqrt{\nu t}$, das Kennzeichen diffusiver Ausbreitung.

Stokes' zweites Problem. Die Platte oszilliert: $u(0,t) = U\cos(\omega t)$. Sucht man eine Lösung der Form $u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$, so ergibt sich

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

also

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

Die Amplitude nimmt mit der Eindringtiefe $\delta_s$ exponentiell ab, und die Phase hinkt linear nach: eine Transversalwelle, deren Energie vollständig durch Viskosität dissipiert wird. Höhere Frequenzen dringen weniger tief ein.

Poiseuille-, Couette- und Stokes-Strömungen erhalten die meiste Aufmerksamkeit. Sie sind jedoch keineswegs die einzigen exakten Lösungen. Bei Weitem nicht.

• Taylor-Green-Wirbel: ein abklingendes Muster wirbelnder Vortices in zwei Dimensionen mit echter Wirbelstruktur. Einen CFD-Code getestet? Wahrscheinlich hast du ihn daran laufen lassen. Es ist der Benchmark, zu dem alle zuerst greifen, und das schon seit Jahrzehnten.
• Jeffery-Hamel-Strömung: Strömung in einem keilförmigen Kanal, der konvergiert oder divergiert. Sie erfasst, wie Fluid in einen enger werdenden Spalt beschleunigt oder in einen sich erweiternden abgebremst wird.
• Hiemenz-Staupunktströmung: Fluid, das frontal auf eine ebene Wand trifft, an der Oberfläche bis auf null abbremst und seitlich abgelenkt wird. Wind, der auf ein Gebäude trifft. Ein Strahl, der auf eine Platte prallt.

Symmetrie. Jede dieser Lösungen nutzt eine bestimmte geometrische Symmetrie aus, um die Gleichungen handhabbar zu machen, und sie sind in spezialisierten Kontexten wichtig; für die alltägliche einführende Strömungsmechanik leisten Poiseuille und Couette aber weiterhin die Schwerarbeit.

Über parallele Strömungen hinaus nutzen mehrere weitere Familien exakter Lösungen bestimmte Symmetrien oder selbstähnliche Strukturen aus:

• Taylor-Green-Wirbel. In 2D, $u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$ mit $aA = bB$ (Inkompressibilität) und exponentiellem zeitlichem Abklingen $e^{-\nu(a^2+b^2)t}$. Dies ist eine exakte Lösung der vollständigen 2D-Navier-Stokes-Gleichungen, einschließlich des nichtlinearen Terms (der sich als Gradient herausstellt und in den Druck absorbiert wird). Sie ist ein Standard-Validierungsfall für DNS-Codes.
• Jeffery-Hamel-Strömung. Stationäre, rein radiale Strömung $u_r(r,\theta)$ in einem Keil mit Halbwinkel $\alpha$. Die Stromfunktion $\psi = \nu f(\theta)$ erfüllt eine nichtlineare ODE dritter Ordnung in $\theta$. Lösungen existieren sowohl für konvergierende als auch für divergierende Kanäle, mit reicher Bifurkationsstruktur bei höheren Reynolds-Zahlen.
• Hiemenz-Staupunktströmung. Eine 2D-Ähnlichkeitslösung für eine Strömung, die auf eine ebene Wand trifft. Die Stromfunktion $\psi = \sqrt{a\nu}\, x\, f(\eta)$, $\eta = y\sqrt{a/\nu}$, reduziert Navier-Stokes auf die Hiemenz-ODE: $f''' + ff'' - f'^2 + 1 = 0$ mit $f(0)=f'(0)=0$, $f'(\infty)=1$.

Wir können die Navier-Stokes-Gleichungen in all diesen Fällen exakt lösen. Warum gibt es also immer noch ein offenes Millionen-Dollar-Problem?

Tricks. Jede einzelne Lösung beruht auf einem Trick. Die Geometrie ist so sorgfältig gewählt, dass der schwierigste Teil der Gleichung (der nichtlineare Term) entweder vollständig verschwindet oder sich auf etwas Handhabbares reduziert, und das Problem wird gerade deshalb lösbar, weil ihm alles entzogen wurde, was Navier-Stokes schwierig macht. Rohrströmung? Eindimensional. Couette? Eine Linie. Der Taylor-Green-Wirbel versteckt seine Nichtlinearität im Druck.

Das Millennium-Preisproblem fragt nach allgemein dreidimensionale Strömungen. Keine Symmetrietricks. Keine vereinfachende Geometrie. Glatte divergenzfreie Anfangsdaten, volle nichtlineare Wechselwirkung über alle Skalen hinweg. In diesem Rahmen hat niemand bewiesen, dass Lösungen immer glatt bleiben, und niemand hat bewiesen, dass sie einen Blow-up entwickeln. Wir wissen es wirklich nicht.

Diese exakten Lösungen sagen uns also etwas, aber bei Weitem nicht genug. Sie beweisen, dass die Gleichungen explizite glatte Lösungen erzeugen können, wenn man ihnen starke Symmetrie als Stütze gibt. Die Millionen-Dollar-Frage ist, ob Glattheit immer gilt, für beliebige glatte divergenzfreie Anfangsdaten in den standardmäßigen 3D-Formulierungen, oder ob irgendwo in der vollen Wucht der Turbulenz etwas katastrophal und irreversibel schiefläuft.

Jede oben diskutierte exakte Lösung erreicht Handhabbarkeit dadurch, dass sie den nichtlinearen Advektionsterm eliminiert oder trivialisiert $(u \cdot \nabla)u$. In Parallelströmungen verschwindet er identisch. Im Taylor-Green-Wirbel ist er ein Gradient, der in den Druck aufgenommen wird. In Ähnlichkeitslösungen reduziert er sich durch einen Variablenwechsel auf ein niederdimensionales Problem.

Das Clay-Millennium-Problem betrifft das Anfangswertproblem für die dreidimensionalen inkompressiblen Navier-Stokes-Gleichungen mit beliebigen glatten, rasch abfallenden Anfangsdaten, genau den Bereich, in dem keine dieser Vereinfachungen gilt.

Die Frage ist, ob Lösungen von

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

auf $\mathbb{R}^3$ in $C^\infty(\mathbb{R}^3 \times [0,\infty))$ für alle Zeiten bleiben, gegeben glatte, divergenzfreie, rasch abfallende Anfangsdaten. Exakte Lösungen behandeln dies nicht, weil sie in Unterräumen des Lösungsraums leben, in denen die volle nichtlineare Dynamik nie zum Tragen kommt.

Kurz gesagt: Exakte Lösungen liefern explizite glatte Lösungen in hochsymmetrischen Regimen. Das offene Problem ist, ob Wohgestelltheit sich auf uneingeschränkte 3D-Strömungen erstreckt. Zur präzisen Formulierung siehe Existenz und Glattheit der Navier-Stokes-Gleichungen.

Um die Gleichungen selbst zu verstehen und was jeder Term bedeutet, beginnen Sie mit Was sind die Navier-Stokes-Gleichungen?

Um zu sehen, wie diese Gleichungen aus ersten Prinzipien aufgebaut werden, lesen Sie Wie die Navier-Stokes-Gleichungen hergeleitet werden.

Um zu verstehen, wann laminare Strömungen wie Poiseuille- und Couette-Strömung in Turbulenz übergehen, lesen Sie Reynolds-Zahl, Turbulenz und warum kleine Skalen wichtig sind.

Um die Millionen-Dollar-Frage zu verstehen, die exakte Lösungen nicht beantworten können, lesen Sie Das Millennium-Problem: Existenz und Glattheit.

Empfohlene nächste Seiten auf dieser Website:

• Was sind die Navier-Stokes-Gleichungen? (das zugrunde liegende System und seine Terme)
• Wie die Navier-Stokes-Gleichungen hergeleitet werden (von der Impulsbilanz und dem newtonschen Stoffgesetz zur PDE)
• Reynolds-Zahl, Turbulenz und warum kleine Skalen wichtig sind (der Übergang von den hier beschriebenen laminaren Regimen zur vollständig turbulenten Strömung)
• Das Millennium-Problem: Existenz und Glattheit (die präzise Formulierung des offenen Problems, das exakte Lösungen nicht lösen)