Lavaldüse

Lavaldüse im Schnitt mit Flussrichtung des Mediums, Strömungsgeschwindigkeit (v), Druck (p) und Temperatur (T)
Schnittmodell eines RD-107-Raketentriebwerks

Die Lavaldüse (auch Expansionsdüse) ist eine von Ernst Körting 1878 für Dampfstrahlapparate und dem Schweden Carl Gustav Patrik de Laval 1883 für die Beaufschlagung von Dampfturbinen mit Wasserdampf unabhängig voneinander entwickelte Düse.

Eine Lavaldüse ist ein Strömungsorgan, bei dem sich der Querschnitt zunächst verengt und anschließend weitet, wobei der Übergang von einem zum anderen Teil stetig erfolgt. Die Querschnittsfläche ist üblicherweise an jeder Stelle kreis- oder ellipsenförmig.

Lavaldüsen werden bereits seit der V2 und auch heute bei Raketentriebwerken verwendet. Das Ziel ist, ein durchströmendes Fluid auf Überschallgeschwindigkeit zu beschleunigen, ohne dass es zu starken Verdichtungsstößen kommt. Die Schallgeschwindigkeit wird kurz nach dem engsten Querschnitt der Düse erreicht. Die Entspannung im divergenten Teil der Düse setzt Wärmeenergie in Bewegungsenergie um. Ferner sollen möglichst große Anteile des ausströmenden Fluids eine parallel zur Achse verlaufende Geschwindigkeit haben, um schubwirksamer zu sein.

Herleitung der Form

Die Eulersche Bewegungsgleichung:

c d c d x = 1 ρ d p d x = 1 ρ d p d ρ d ρ d x {\displaystyle {\begin{aligned}c\cdot {\frac {\mathrm {d} c}{\mathrm {d} x}}&=-{\frac {1}{\rho }}\cdot {\frac {\mathrm {d} p}{\mathrm {d} x}}\\&=-{\frac {1}{\rho }}\cdot {\frac {\mathrm {d} p}{\mathrm {d} \rho }}\cdot {\frac {\mathrm {d} \rho }{\mathrm {d} x}}\end{aligned}}}

mit

  • der Strömungsgeschwindigkeit c {\displaystyle c} des Fluids
  • der Koordinate x {\displaystyle x} in Strömungsrichtung
  • dem Druck p {\displaystyle p}
  • der Dichte ρ {\displaystyle \rho }

ergibt zusammen mit der von der Dichte abhängigen Schallgeschwindigkeit a = d p d ρ {\displaystyle a={\sqrt {\frac {\mathrm {d} p}{\mathrm {d} \rho }}}} :

c d c d x = a 2 ρ d ρ d x {\displaystyle c\cdot {\frac {\mathrm {d} c}{\mathrm {d} x}}=-{\frac {a^{2}}{\rho }}{\frac {\mathrm {d} \rho }{\mathrm {d} x}}}

Einsetzen der Mach-Zahl M a = c a {\displaystyle {\mathit {Ma}}={\frac {c}{a}}} liefert:

1 ρ d ρ d x = M a 2 1 c d c d x ( 1 ) {\displaystyle {\frac {1}{\rho }}\cdot {\frac {\mathrm {d} \rho }{\mathrm {d} x}}=-{\mathit {Ma}}^{2}\cdot {\frac {1}{c}}\cdot {\frac {\mathrm {d} c}{\mathrm {d} x}}\qquad (1)} ,

Diese Gleichung sagt aus, dass die relative Dichteänderung längs des Stromfadens x {\displaystyle x} proportional ist zur relativen Geschwindigkeitsänderung mit dem Proportionalitätsfaktor M a 2 . {\displaystyle {\mathit {Ma}}^{2}.} Aus dem quadratischen Proportionalitätsfaktor folgt, dass

  • bei einer Unterschallströmung ( M a < 1 {\displaystyle {\mathit {Ma}}<1} ) die relative Dichteänderung (wesentlich) kleiner als die relative Geschwindigkeitsänderung ist
  • bei einer Überschallströmung ( M a > 1 {\displaystyle {\mathit {Ma}}>1} ) die relative Dichteänderung (wesentlich) größer als die relative Geschwindigkeitsänderung ist.

Ferner muss noch die Kontinuitätsgleichung betrachtet werden:

ρ c A = m ˙ = konst ln ρ + ln c + ln A = ln ( m ˙ ) = ln ( konst ) d ρ ρ + d c c + d A A = 0 {\displaystyle {\begin{alignedat}{2}&\rho \cdot c\cdot A&&=\quad {\dot {m}}\quad &={\text{konst}}\\\Leftrightarrow \quad &\ln \rho +\ln c+\ln A&&=\ln({\dot {m}})&=\ln({\text{konst}})\\\Leftrightarrow \quad &{\frac {\mathrm {d} \rho }{\rho }}+{\frac {\mathrm {d} c}{c}}+{\frac {\mathrm {d} A}{A}}&&=0\end{alignedat}}}

mit

  • der Querschnittsfläche A {\displaystyle A}
  • dem Massenstrom m ˙ {\displaystyle {\dot {m}}}

Differenziert man längs des Stromfadens, so ergibt sich

1 ρ d ρ d x + 1 c d c d x + 1 A d A d x = 0 {\displaystyle \Rightarrow {\frac {1}{\rho }}\cdot {\frac {\mathrm {d} \rho }{\mathrm {d} x}}+{\frac {1}{c}}\cdot {\frac {\mathrm {d} c}{\mathrm {d} x}}+{\frac {1}{A}}\cdot {\frac {\mathrm {d} A}{\mathrm {d} x}}=0}

Unter Berücksichtigung von Gleichung (1) folgt:

1 c d c d x = 1 M a 2 1 1 A d A d x {\displaystyle {\frac {1}{c}}\cdot {\frac {\mathrm {d} c}{\mathrm {d} x}}={\frac {1}{{\mathit {Ma}}^{2}-1}}\cdot {\frac {1}{A}}\cdot {\frac {\mathrm {d} A}{\mathrm {d} x}}}
Lavaldüse

Nimmt man die Querschnittsfläche  A ( x ) {\displaystyle A(x)} als gegeben, c ( x ) {\displaystyle c(x)} und  M a ( x ) {\displaystyle {\mathit {Ma}}(x)} hingegen als unbekannt an, so ermöglicht die letzte Gleichung die folgende qualitative Diskussion der Strömung durch eine Düse.

Will man eine Strömung beschleunigen, also d c d x > 0 {\displaystyle {\frac {\mathrm {d} c}{\mathrm {d} x}}>0} , so folgt aus der letzten Gleichung die Form der Lavaldüse:

  • Einlauf mit Unterschallströmung ( M a < 1 {\displaystyle {\mathit {Ma}}<1} ): hier muss d A d x < 0 {\displaystyle {\frac {\mathrm {d} A}{\mathrm {d} x}}<0} sein, die Düse muss sich also verengen (konvergenter Teil)
  • kurz nach dem engsten Querschnitt wird M a = 1 {\displaystyle {\mathit {Ma}}=1} erreicht
  • weitere Beschleunigung auf M a > 1 {\displaystyle {\mathit {Ma}}>1} im Auslauf: hier muss d A d x > 0 {\displaystyle {\frac {\mathrm {d} A}{\mathrm {d} x}}>0} sein, die Düse muss sich also erweitern (divergenter Teil).

Literatur

  • Ernst Messerschmid, Stefanos Fasoulas: Raumfahrtsysteme: Eine Einführung mit Übungen und Lösungen. Hrsg.: Springer Vieweg. 5. Auflage. 2017, ISBN 978-3-662-49637-4, S. 199 ff. 
  • Erich Hahne: Technische Thermodynamik. Einführung und Anwendung, 5. Auflage, Oldenbourg Verlag, München 2010, ISBN 978-3-486-59231-3.
  • Herbert Oertel jr., Martin Böhle, Thomas Reviol: Strömungsmechanik für Ingenieure und Naturwissenschaftler. 7. Auflage, Springer Fachmedien, Wiesbaden 2015, ISBN 978-3-658-07785-3.
  • Heinz Schade, Ewald Kunz, Frank Kameier, Christian Oliver Paschereit: Strömungslehre. 4. Auflage, Walter de Gruyter GmbH, Berlin 2013, ISBN 978-3-11-029221-3.

Siehe auch

Weblinks

  • Anleitung Lavaldüse auf TU Dresden (abgerufen am 24. März 2016)
  • Die Lavaldüse (abgerufen am 24. März 2016)
  • Untersuchung der turbulenten Grenzschicht einer Lavaldüse mittels CFD-Modellierung (abgerufen am 24. März 2016)