# Düsenströmung

Als Düsenströmung wird die Strömung eines Fluids, d. h. eines Gases oder einer Flüssigkeit, durch eine Düse bezeichnet. Dabei wird das Fluid beschleunigt, während der Druck abnimmt. Bei einer Düsenströmung wird potentielle in kinetische Energie umgewandelt.
Der durch die Düse fließende Massenstrom wird bestimmt durch den Vordruck $p_{i}$, den Gegendruck $p_{a}$ und den engsten Querschnitt $A$ der Düse: bei konstantem Vordruck nimmt der Massenstrom mit sinkendem Gegendruck im Ausflussraum zu, bis bei einem bestimmten kritischen Druckverhältnis die Geschwindigkeit im engsten Querschnitt gerade die Schallgeschwindigkeit des Fluids erreicht. Bei weiterer Reduzierung des Gegendrucks unter den kritischen Wert bleibt der Massenstrom konstant.

## Eindimensionales Berechnungsmodell

Das eindimensionale Berechnungsmodell (mit der Koordinate $x$ längs der Düsenachse in Strömungsrichtung) stützt sich auf die Erhaltungssätze von Masse, Impuls und Energie sowie auf Zustandsgleichungen des Fluids. An jeder Stelle $x$ wird lokales thermodynamisches Gleichgewicht vorausgesetzt.
Zudem wird angenommen, dass die Strömungsgeschwindigkeit $w$, der Druck $p$ und die Temperatur $T$ in jedem Querschnitt $(y,z)$ senkrecht zur Strömungsrichtung gleichförmig sind:
* $w(y,z)$ ist konstant
* $p(y,z)$ ist konstant
* $T(y,z)$ ist konstant
Strömungsgeschwindigkeit, Druck und Temperatur sind also nur von der ersten Raumdimension $x$ abhängig.
Das Modell ergibt Beziehungen zwischen integralen (globalen) Eingangs- und Ausgangsgrößen (Massenstrom, mittlere Geschwindigkeit, Temperatur).
Die Erhaltung der Masse wird gewährleistet durch die Kontinuitätsgleichung, d. h. der Massenstrom ${\dot {m}}$ durch jeden Querschnitt $A$ ist gleich:
    ${\dot {m}}(x)=A(x)\cdot w(x)\cdot \rho (x)$ ist konstant
Hierbei bezeichnet
* $\rho (x)$ die Dichte.
Die Erhaltung der Energie wird durch die Bernoulli-Gleichung gewährleistet:
    ${\frac {w^{2}}{2}}+h=h_{0}$ ist konstant
wobei
* $h$ die spezifische Enthalpie des Fluids im Strömungszustand mit der Geschwindigkeit $w$
* $h_{0}$ die spezifische Enthalpie des Fluids im Kesselzustand ($w=0$) bezeichnet.
Wenn die Strömung adiabatisch verläuft und Reibungsverluste vernachlässigt werden können, so bleibt die Entropie des Fluids während der Beschleunigung durch die Düse in erster Näherung konstant (isentrope Düsenströmung):
    $S$ ist konstant
Sind die Dichte $\rho (s,p)$ und die spezifische Enthalpie $h(s,p)$ in Abhängigkeit von der spezifischen Entropie $s$ und dem Druck $p$ gegeben (Zustandsgleichungen des Fluids), dann gilt:
    $\Rightarrow w(p)={\sqrt {2\cdot \left(h_{0}-h(s_{0},p)\right)}}$
wobei $s_{0}$ die spezifische Entropie im Kesselzustand ist.
Der Strömungsquerschnitt $A$ in Abhängigkeit vom Druck folgt aus:
    $A(p)={\frac {\dot {m}}{w(p)\cdot \rho (p)}}$
Der Druck $p=p(x)$ und damit alle anderen Größen sind Funktionen der Koordinate $x$.
Der Strömungsquerschnitt $A(p)$ hat ein Minimum bei dem Druck, bei dem die Strömungsgeschwindigkeit $w$ gleich der Schallgeschwindigkeit $c_{s}$ ist.
Die Schallgeschwindigkeit ist definiert durch:
    $c_{s}={\sqrt {\left({\frac {\partial p}{\partial \rho }}\right)_{s}}}$
Zudem gilt allgemein:
    $\left({\frac {\partial h}{\partial p}}\right)_{s}={\frac {1}{\rho (s,p)}}$
Damit ergibt sich:
    $\left({\frac {\partial A}{\partial p}}\right)_{s}={\dot {m}}\cdot {\frac {c_{s}^{2}-w^{2}}{c_{s}^{2}\cdot w^{3}\cdot \rho ^{2}}}$

## Düsenströmung eines idealen Gases

Bei der reversibel adiabaten Strömung eines idealen Gases gilt folgender Zusammenhang zwischen Dichte und Druck:
    $p(\rho )=p_{i}\cdot \left({\frac {\rho }{\rho _{i}}}\right)^{\kappa }$
mit
* $\kappa ={\frac {c_{p}}{c_{V}}}$: Adiabatenexponent (z. B. für Luft: $\kappa =1{,}4$)
* $p_{i}$: Druck im Vordruckbereich
* $\rho _{i}$: Dichte im Vordruckbereich.
Wird die Strömungsgeschwindigkeit in der Vorkammer vernachlässigt ($w_{i}\approx 0$), so ergibt sich bei adiabater Strömung folgender Zusammenhang zwischen Massenstrom und Druckverhältnis:
    ${\dot {m}}=A\cdot \mu \cdot \Psi (p_{a}/p_{i},\kappa )\cdot {\sqrt {2\cdot p_{i}\cdot \rho _{i}}}$
mit
* der Ausflusszahl $\mu$, welche die Verringerung des tatsächlichen Massenstroms durch Strahleinschnürung und Reibung beschreibt. Die Ausflusszahl hängt nur von der Düsengeometrie ab: für gut abgerundete Düsen ist $\mu \approx 1$, für scharfkantige Blenden kann der Wert bis auf $\mu =0{,}59$ verringert sein.
* der Ausflussfunktion $\Psi (x,\kappa )={\begin{cases}\left({\frac {2}{\kappa +1}}\right)^{\frac {1}{\kappa -1}}\cdot {\sqrt {\frac {\kappa }{(\kappa +1)}}}&x\leq \left({\frac {2}{\kappa +1}}\right)^{\frac {\kappa }{\kappa -1}}\\{\sqrt {{\frac {\kappa }{\kappa -1}}\cdot x^{1/\kappa }\cdot \left(x^{1/\kappa }-x\right)}}&\left({\frac {2}{\kappa +1}}\right)^{\frac {\kappa }{\kappa -1}}<x<1\end{cases}}$
    wobei hier statt $x$ das Druckverhältnis $p_{a}/p_{i}$ eingesetzt werden muss. Der Wert $\left({\frac {2}{\kappa +1}}\right)^{\frac {\kappa }{\kappa -1}}$ entspricht dem oben beschriebenen kritischen Druckverhältnis.

### Anwendungsbeispiel: Gasbrennerdüse

Der Massenstrom eines Brenngas-Luft-Gemischs durch die Düse eines Gasbrenners ergibt sich mit sehr guter Genauigkeit aus:
    ${\dot {m}}=A\cdot \mu \cdot \Psi \left({\frac {p_{0}}{p_{0}+p_{D}}},\kappa \right)\cdot \left(p_{0}+p_{D}\right)\cdot {\sqrt {\frac {2\cdot M}{R\cdot T}}}$
wobei
* $A$: kleinster Querschnitt der Düse
* $p_{0}$: Druck im Ausflussraum (Umgebungsdruck)
* $p_{D}$: Überdruck vor der Düse
* $M$: mittlere molare Masse des Brenngas-Luft-Gemischs
* $R=8{,}3145\,\mathrm {J/(mol\,K)}$: Gaskonstante
* $T$: absolute Temperatur des Brenngas-Luft-Gemischs vor der Düse.
Der effektive Düsenquerschnitt $A\cdot \mu$, der als einziger Parameter von der Bauform der Düse abhängt, kann aus Messwerten bestimmt werden.