Kelvin-Helmholtz-Instabilität

Als Kelvin-Helmholtz-Instabilität bezeichnet man das Anwachsen kleiner Störungen in der Scherschicht zweier Fluide mit unterschiedlichen Strömungsgeschwindigkeiten.

Kelvin-Helmholtz-Wirbel in der Atmosphäre des Saturn

Ein anschauliches Beispiel liefern Wellen auf einem See während eines Sturms oder der sich kräuselnde Rauch eines Räucherstäbchens in einem ansonst ruhigen Zimmer.

Wolkenbildung auf Grund einer Scherwelle

Kelvin-Helmholtz-Wirbel in der Atmosphäre hinter dem Monte Duval, Australien

Phänomenologie

Als Wetterphänomen kann es an seltsamen Wolken, die einzeln oder in gleich aussehenden Gruppen am Himmel zu sehen sind, erkennbar werden. Sie entstehen durch eine Verwirbelung von zwei übereinander liegenden Luftschichten, die sich mit unterschiedlicher Geschwindigkeit und/oder Richtung bewegen. Ähnlich wie wenn Wind über Wasser streicht, entstehen Wellen an der Grenzschicht, wobei Teile der meist feuchteren unteren Luftschicht so weit nach oben gewirbelt werden, dass ihr Taupunkt unterschritten wird und es zu Wolkenbildung kommt.

Physikalische Interpretation

Weit entfernt von der Grenzschicht sind die Strömungsgeschwindigkeiten konstant. Nahe der Grenzschicht muss sich aber ein Luftelement schneller über den Wellenbuckel bewegen als ein weiter entferntes (ähnlich wie bei einem Tragflügel). Nach der Bernoulli-Gleichung ist aber der Druck über der Welle infolge der höheren Windgeschwindigkeit kleiner als in der Umgebung. Infolgedessen gibt es eine Kraft, die den Wellenkamm nach oben zieht. Analog verhält es sich in einem Wellental: Die Luft fließt langsamer über die Oberfläche eines Wellentals als in der Umgebung, und darum ist der Umgebungsdruck lokal höher. Das Wellental wird nach unten gedrückt.

Theorie

Ein einfaches Modell für die Kelvin-Helmholtz-Instabilität erhält man, wenn man die folgende Frage zu beantworten sucht: Gegeben ist eine Strömung über eine Grenzschicht. Unter welchen Bedingungen ist die Grenzschicht stabil gegen kleine Störungen?

Störungsrechnung

Numerische Simulation der Kelvin-Helmholtz-Instabilität

Gegeben sei also eine Flüssigkeit mit der Dichte ρ ₋ , die sich horizontal mit der Geschwindigkeit V über eine Flüssigkeit mit der Dichte ρ ₊ bewegt. Bezeichne x eine Koordinate entlang der Scherschicht und y die Koordinate rechtwinklig dazu. Nun gibt man sich eine kleine Störung entlang der Scherschicht und bezeichnet sie mit ξ(x). Die dazu assoziierte Störung des Drucks P bezeichnen wir mit δP und die des Geschwindigkeitsfeldes $\vec v$ mit $\delta \vec v$.

Das Druckfeld lässt sich nun schreiben als

$P(\vec x, t) = P_0 + \delta P(\vec x, t)$

und das Geschwindigkeitsfeld als

$\vec v(\vec x, t) = V \Theta(y) \vec e_x + \delta v(\vec x, t) \vec e_y.$

wobei Θ(y) die Heaviside-Funktion bezeichnet. Diese zwei Störungen substituiert man nun in die einfachste Form fluiddynamischer Gleichungen, nämlich in die inkompressible Euler-Gleichungen: Die Inkompressibilitätsgleichung lautet

$\vec \nabla \cdot \vec v = 0$

und die Euler-Gleichung

$\frac{d \vec v}{d t} = \frac{- \vec \nabla P}{\rho}.$

Dort eingesetzt erhält man für die gestörten Größen

$\vec \nabla \cdot \delta v = 0$

und

$\frac{d \delta \vec v}{d t} = - \frac{\vec \nabla \delta P}{\rho}$

Diese zwei Gleichungen liefern für den gestörten Druck die Laplace-Gleichung

$\nabla^2 \delta P = 0$

Nun sucht man nach einer Wellenmode, die exponentiell mit dem Abstand y von der Grenzfläche abfällt. Aus der Laplace-Gleichung schließen wir, dass für den Druck gelten muss:

$\delta P = \delta P_0\, \exp\{-k |y| + i (kx - \omega t) \}$

Als Nächstes substituiert man dieses Resultat in die gestörten Euler-Gleichungen. Dabei erhält man

$\delta v_y = \frac{i k \delta P}{(\omega - k V) \rho_{+} }$ für $\quad y > 0$

und

$\delta v_y = \frac{- i k \delta P}{\omega \rho_{-} }$ für y < 0.

Nun müssen noch die Randbedingungen erfüllt sein: Zunächst muss die vertikale Komponente der Störung ξ(x) an der Scherschicht stetig sein. Weiter muss der Druck an der Grenzfläche stetig sein. Aus diesen beiden Bedingungen ergeben sich diese Bedingungen an die Störung:

$\xi = \frac{i \delta v_y}{\omega - k V}$

direkt über der Scherschicht (y = 0 ₊ ) und

$\xi = \frac{i \delta v_y}{\omega}$

direkt unter der Scherschicht (y = 0 ₋ ). Daraus lässt sich ein Zusammenhang zwischen der Dichte der Flüssigkeiten, den Wellenmoden und der Relativgeschwindigkeit der Flüssigkeiten herstellen:

$\rho_+ (\omega - k V)^2 \, + \, \rho_- \omega^2 = 0$

Löst man diese Gleichung nach ω auf, so erhält man ein Dispersionsgesetz für die linearen Kelvin-Helmholtz-Moden:

$\omega = k V ( \frac{\rho_{+} \pm i (\rho_{+} \rho_{-} )^{1/2} }{\rho_{+} + \rho_{-} } )$

Zeitliches Wachstum der Störung

Nun kann man das Wachstum der Störung berechnen. Bewegt man sich mit der Geschwindigkeit ρ ₊ V / (ρ ₊ + ρ ₋ ) entlang der Oberfläche, so ergibt sich für die Geschwindigkeit der oberen Flüssigkeit V' = ρ ₋ V / (ρ ₊ + ρ ₋ ). Die Störung entwickelt sich nun folgendermaßen:

$\delta P, \xi \propto \exp \{ k V' (\rho_{+} / \rho_{-} )^{1/2} t \} \cos (i k x' )$

Räumliches Wachstum der Störung

$k = \frac{\omega}{V} \left [ 1 \pm i (\frac{\rho_{+}}{\rho_{-}} )^{1/2} \right ]$

Siehe auch

• Scherschicht
• Karmansche Wirbelstraße
• Clear Air Turbulence

Weblinks

• Praktikumsunterlagen bei der gwdg (PDF-Datei; 675 kB)
• http://www.amath.washington.edu/~rjl/clawpack/euler/sester/KH.html - animiertes Beispiel
• http://www.math.lsa.umich.edu/~krasny/k-h.html - zeitliche Momentaufnahmen
• Wolkenfotos gib es: hier