Düsenströmung: Unterschied zwischen den Versionen

Düsenströmung: Unterschied zwischen den Versionen

129.247.247.240 (Diskussion)
(Formel für Ausflussfunktion bei überkritischem Druckverhältnis war falsch im Exponenten. Das ergibt sich durch Stetigkeitsbedingung und Einsetzten des kritischen Druckverhältnisses in die Formel für unterkritich.)
 
imported>Morneo06
K (Änderungen von 2001:638:513:E80:2C6D:E68B:A100:D7E (Diskussion) rückgängig gemacht (HG) (3.4.10))
 
Zeile 1: Zeile 1:
Als '''Düsenströmung''' wird die [[Strömung]] eines Fluids, d. h. eines [[Gas]]es oder einer [[Flüssigkeit]], durch eine [[Düse]] bezeichnet. Dabei wird das Fluid [[Beschleunigung|beschleunigt]], während der [[Druck (Physik)|Druck]] abnimmt. Bei einer Düsenströmung wird [[potentielle Energie|potentielle]] in [[kinetische Energie]] umgewandelt.  
Als '''Düsenströmung''' wird die [[Strömung]] eines Fluids, d. h. eines [[Gas]]es oder einer [[Flüssigkeit]], durch eine [[Düse]] bezeichnet. Dabei wird das Fluid [[Beschleunigung|beschleunigt]], während der [[Druck (Physik)|Druck]] abnimmt. Bei einer Düsenströmung wird [[potentielle Energie|potentielle]] in [[kinetische Energie]] umgewandelt.


Der durch die Düse fließende [[Massenstrom]] wird bestimmt durch den Vordruck <math>p_i</math>, den Gegendruck <math>p_a</math> und den engsten Querschnitt <math>A</math> 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.
Der durch die Düse fließende [[Massenstrom]] wird bestimmt durch den Vordruck <math>p_i</math>, den Gegendruck <math>p_a</math> und den engsten Querschnitt <math>A</math> 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.
Zeile 8: Zeile 8:
Zudem wird angenommen, dass die [[Strömungsgeschwindigkeit]] <math>w</math>, der Druck <math>p</math> und die Temperatur <math>T</math> in jedem Querschnitt <math>(y, z)</math> senkrecht zur Strömungsrichtung gleichförmig sind:
Zudem wird angenommen, dass die [[Strömungsgeschwindigkeit]] <math>w</math>, der Druck <math>p</math> und die Temperatur <math>T</math> in jedem Querschnitt <math>(y, z)</math> senkrecht zur Strömungsrichtung gleichförmig sind:


:<math>w(x = \text{konst}; y, z) = \text{konst}</math>
* <math>w(y, z)</math> ist konstant
:<math>p(x = \text{konst}; y, z) = \text{konst}</math>
* <math>p(y, z)</math> ist konstant
:<math>T(x = \text{konst}; y, z) = \text{konst}</math>
* <math>T(y, z)</math> ist konstant


Das Modell ergibt Beziehungen zwischen integralen (globalen) Eingangs- und Ausgangsgrößen (Massenstrom, mittlere Geschwindigkeit, Temperatur).  
Strömungsgeschwindigkeit, Druck und Temperatur sind also nur von der ersten Raumdimension <math>x</math> 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.&nbsp;h. der Massenstrom <math>\dot{m}</math> durch jeden Querschnitt <math>A</math> ist gleich:
Die Erhaltung der Masse wird gewährleistet durch die [[Kontinuitätsgleichung]], d.&nbsp;h. der Massenstrom <math>\dot{m}</math> durch jeden Querschnitt <math>A</math> ist gleich:


:<math>\dot{m}(x) = A(x) \cdot w(x) \cdot \rho(x) = \text{konst}</math>
:<math>\dot{m}(x) = A(x) \cdot w(x) \cdot \rho(x)</math> ist konstant


Hierbei bezeichnet
Hierbei bezeichnet
Zeile 23: Zeile 25:
Die Erhaltung der Energie wird durch die [[Bernoulli-Gleichung]] gewährleistet:
Die Erhaltung der Energie wird durch die [[Bernoulli-Gleichung]] gewährleistet:


:<math>\frac{w^2}{2} + h = h_0 = \text{konst}</math>
:<math>\frac{w^2}{2} + h = h_0</math> ist konstant


wobei
wobei
* <math>h</math> die [[spezifische Größe|spezifische]] [[Enthalpie]] des Fluids im Strömungszustand mit der Geschwindigkeit <math>w</math>
* <math>h</math> die [[spezifische Größe|spezifische]] [[Enthalpie]] des Fluids im Strömungszustand mit der Geschwindigkeit <math>w</math>
* <math>h_0</math> die spezifische [[Totalenthalpie|Enthalpie des Fluids im Kesselzustand]] (<math>w = 0</math>) bezeichnet.  
* <math>h_0</math> die spezifische [[Totalenthalpie|Enthalpie des Fluids im Kesselzustand]] (<math>w = 0</math>) bezeichnet.


Wenn die Strömung [[Adiabatische Zustandsänderung|adiabatisch]] verläuft und [[Reibung]]s<nowiki/>verluste vernachlässigt werden können, so bleibt die [[Entropie (Thermodynamik)|Entropie]] des Fluids während der Beschleunigung durch die Düse in erster Näherung konstant ([[Isentrop|isentrope]] Düsenströmung):
Wenn die Strömung [[Adiabatische Zustandsänderung|adiabatisch]] verläuft und [[Reibung]]s<nowiki/>verluste vernachlässigt werden können, so bleibt die [[Entropie (Thermodynamik)|Entropie]] des Fluids während der Beschleunigung durch die Düse in erster Näherung konstant ([[isentrop]]e Düsenströmung):


:<math>S = \text{konst}</math>
:<math>S</math> ist konstant


Sind die Dichte <math>\rho(s,p)</math> und die spezifische Enthalpie <math>h(s,p)</math> in Abhängigkeit von der spezifischen Entropie <math>s</math> und dem Druck <math>p</math> gegeben (Zustandsgleichungen des Fluids), dann gilt:
Sind die Dichte <math>\rho(s,p)</math> und die spezifische Enthalpie <math>h(s,p)</math> in Abhängigkeit von der spezifischen Entropie <math>s</math> und dem Druck <math>p</math> gegeben (Zustandsgleichungen des Fluids), dann gilt:
Zeile 37: Zeile 39:
:<math>\Rightarrow w(p) = \sqrt{2 \cdot \left( h_0 - h(s_0,p) \right)}</math>
:<math>\Rightarrow w(p) = \sqrt{2 \cdot \left( h_0 - h(s_0,p) \right)}</math>


wobei <math>s_0</math> die spezifische Entropie im Kesselzustand ist.  
wobei <math>s_0</math> die spezifische Entropie im Kesselzustand ist.


Der Strömungsquerschnitt <math>A</math> in Abhängigkeit vom Druck folgt aus:
Der Strömungsquerschnitt <math>A</math> in Abhängigkeit vom Druck folgt aus:
Zeile 45: Zeile 47:
Der Druck <math>p = p(x)</math> und damit alle anderen Größen sind Funktionen der Koordinate <math>x</math>.
Der Druck <math>p = p(x)</math> und damit alle anderen Größen sind Funktionen der Koordinate <math>x</math>.


Der Strömungsquerschnitt <math>A(p)</math> hat ein [[Minimum]] bei dem Druck, bei dem die Strömungsgeschwindigkeit <math>w</math> gleich der Schallgeschwindigkeit <math>c_s</math> ist.  
Der Strömungsquerschnitt <math>A(p)</math> hat ein [[Extremwert|Minimum]] bei dem Druck, bei dem die Strömungsgeschwindigkeit <math>w</math> gleich der Schallgeschwindigkeit <math>c_s</math> ist.


Die Schallgeschwindigkeit ist definiert durch:  
Die Schallgeschwindigkeit ist definiert durch:


:<math>c_s=\sqrt{\left( \frac{\partial p}{\partial \rho}\right)_s}</math>
:<math>c_s=\sqrt{\left( \frac{\partial p}{\partial \rho}\right)_s}</math>

Aktuelle Version vom 22. April 2021, 14:29 Uhr

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

Die Erhaltung der Energie wird durch die Bernoulli-Gleichung gewährleistet:

$ {\frac {w^{2}}{2}}+h=h_{0} $ ist konstant

wobei

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 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.