Rankine-Hugoniot-Bedingung

Die Rankine-Hugoniot-Bedingung oder auch Rankine-Hugoniot-Gleichung (nach William John Macquorn Rankine und Pierre-Henri Hugoniot) beschreibt das Verhalten von Stoßwellen durch eine eindimensionale hyperbolische Erhaltungsgleichung:

mit

der Geschwindigkeit .

Gegeben zwei Zustände u_L und u_R links und rechts eines Stoßes, besagt die Bedingung, dass die Stoßgeschwindigkeit die Gleichung

$f(u_L) - f(u_R) = s \cdot (u_L - u_R)$

erfüllt. Im Falle einer skalaren Gleichung $\left( u \in \mathbb{R}^1 \right)$ liefert dies direkt die Stoßgeschwindigkeit

$\Leftrightarrow s = \frac{f(u_L) - f(u_R)}{u_L - u_R}$ .

Bei Systemen mit $u \in \mathbb{R}^n; \, n \geq 2$ ist die Situation schwieriger.

Im Falle einer linearen Gleichung $u_t + A\cdot u_x = 0$ ergibt sich die Bedingung, dass die Stoßgeschwindigkeit ein Eigenwert der Matrix sein muss und die Differenz der Zustände ein Eigenvektor von . Dies ist nicht immer möglich, was dann bedeutet, dass diese Zustände durch eine Sequenz von Unstetigkeiten verbunden sind.
Dies kann auch auf nichtlineare Gleichungen angewandt werden, wobei dann zu beachten ist, dass sich hier die Stoßgeschwindigkeiten mit der Zeit ändern.

Umgekehrt bezeichnet man bei Systemen die Menge der Zustände, die mit einem gegebenen festen Zustand durch einen einzigen Stoß verbunden werden können, als Hugoniot-Lokus.

Beispiele

Advektionsgleichung in 1D

Eine sehr einfache Erhaltungsgleichung ist gegeben durch den skalaren Fluss:

$f\colon \mathbb {R} \to \mathbb {R} ,x\mapsto ax\quad {\textrm {mit}}\;a\in \mathbb {R} .$

Die Sprungbedingung ergibt somit sofort: $s=a$ .

Burgersgleichung in 1D

Die Burgersgleichung ist definiert über den folgenden Fluss:

$f\colon \mathbb {R} \to \mathbb {R} ,x\mapsto {\frac {1}{2}}x^{2}.$

Die Sprungbedingung ergibt somit: $s={\frac {u_{L}+u_{R}}{2}}$ .

Euler-Gleichungen

Im Falle der Euler-Gleichungen ergeben sich spezielle Beziehungen. Elimination der Geschwindigkeit führt auf:

$2 \cdot \left( h_R - h_L \right) = \left( p_R - p_L \right) \cdot \left( \frac{1}{\rho_L} + \frac{1}{\rho_R} \right)$

mit

dem Druck
der Dichte $\rho$
der spezifischen Enthalpie . Das wird als hugoniotsche Adiabate bezeichnet (s.u.)
- der inneren Energie pro Masse (spezifische Größe)

Wird nun die Zustandsgleichung für das ideale Gas verwendet :

$p = \rho \cdot (\kappa - 1) \cdot e$

mit

dem Adiabatenexponenten $\kappa$ ,

so ergibt sich

$\Rightarrow \frac{p_L}{p_R} = \frac{(\kappa + 1) - (\kappa - 1) \cdot \frac{\rho_R}{\rho_L}} {(\kappa + 1) \cdot \frac{\rho_R}{\rho_L} - (\kappa-1)} \quad \Leftrightarrow \quad \frac{\rho_L}{\rho_R} = \frac{p_L \cdot (\kappa + 1) + p_R \cdot (\kappa - 1)} {p_L \cdot (\kappa - 1) + p_R \cdot (\kappa + 1)}$ .

Da die Drücke stets positiv sind, folgt daraus für das Dichteverhältnis:

$\Rightarrow \frac{\rho_L}{\rho_R} \leq \frac{\kappa + 1}{\kappa - 1}$

Für Luft mit $\kappa \approx 1,4$ beträgt das maximale Dichteverhältnis ungefähr 6. Dieses Ergebnis ist anschaulich nachvollziehbar, da eine Zunahme des Drucks auch zu einer Temperaturzunahme führt, die der Dichtezunahme teilweise entgegenwirkt. Während die Stoßstärke (der Überdruck) beliebig groß werden kann, erreicht das Dichteverhältnis also einen endlichen Grenzwert.

Allerdings kann hohe Temperatur bei starken Stößen zur Dissoziation oder sogar zur Ionisation und damit zur Zunahme der thermodynamischen Freiheitsgrade und damit wiederum zu einem kleineren Wert von $\kappa$ führen. Daher kann in realen Gasen die Obergrenze für das Dichteverhältnis wesentlich höher sein als in idealem Gas.

Die ersten beiden Erhaltungssätze folgen aus den Eulergleichungen bzw. führen zu diesen. Mit ihnen können die Sprungbedingungen für die Geschwindigkeit und die Dichte (bzw. den Druck) an der Stoßfront dargestellt werden. Die zentrale Idee von Rankine und Hugoniot war nun die Nutzung des dritten Erhaltungssatzes (der Energieerhaltung), um damit eine Sprungbedingung für die Entropie zu formulieren. Diese ist an der Stoßfront unstetig:

Daraus folgt, dass eine Stoßwelle kein adiabatischer (oder isentroper) Prozess mehr ist, sondern die Enthalpieänderung auch eine Entropiekomponente enthält (hugoniotsche Adiabate, auch als Stoßadiabate bekannt):

${\rm d}H = \int^{2}_{1} \frac{{\rm d}p}{\rho} + T \cdot {\rm d}S$

im Gegensatz zu

${\rm d}H = \int^{2}_{1}\frac{{\rm d}p}{\rho}$

für eine rein adiabatische Verdichtung.

Basierend auf einem Artikel in:

Wikipedia.de