Logo PTB

Verwendbarkeit von Niederdruck-Kalibrierungen kritischer Düsen im Hochdruck-Gas

15.11.2018

2016 wurde von der PTB für kritisch betriebene Venturidüsen eine Beschreibung des Düsenbeiwertes cD(Re) in Abhängigkeit von der Reynoldszahl vorgestellt, die sowohl den Bereich laminarer als auch turbulenter Grenzschichten gleichzeitig abdeckt [1]. Mit den Kalibrierdaten von insgesamt 30 Düsen kann nunmehr gezeigt werden, dass damit auch eine zuverlässige Extrapolation einer Kalibrierung bei kleinen Reynoldszahlen (Niederdruck) in den Bereich hoher Reynoldszahlen (Hochdruck) möglich ist.

Die Gleichungen für cD(Re) – im laminaren Bereich cD = a - blamRe-0,5 und für den turbulenten Bereich cD = a - bturbRe-0,139 – sind über die Relation bturb = 0,003654 ⋅ b1,736lam miteinander verknüpft. Gleichzeitig hängen die beiden Parameter a und blam von dem lokalen Krümmungsradius RHals an der Stelle des engsten Düsenquerschnittes - auch Hals genannt - ab [2], [3].

Damit kann die Anpassung der cD(Re)-Kurve an Messwerte mit Hilfe eines einzigen Parameters RHals erfolgen. Als Beispiel sind in Bild 1 die Messwerte aus einem Ringvergleich [4] zwischen PTB, NIST und LNE gezeigt. Die Anpassung der cD(Re)-Kurve erfolgte über eine Least Square Anpassung ausschließlich an die Messpunkte der PTB bei den niedrigen Reynoldszahlen, bei denen laminare Grenzschichten das Düsenverhalten bestimmen. Bei der Berechnung von cD(Re) ist zu beachten, dass in die Bestimmung von a und blam nach [2] und [3] auch noch der Isentropen-Exponent und die Prandtl-Zahl mit eingehen, was zu einer geringfügigen Differenz der Kurven für Luft/Stickstoff und Erdgas bei gleichem RHals führt.

Bild 1:    Messwerte für cD,mess aus einem Ringvergleich PTB, NIST und LNE [4] und cD(Re) auf Basis einer Least Square Anpassung an die Messwerte der PTB bei kleinen Reynoldszahlen. Für die Bewertung der Extrapolation werden die Differenzen ΔcD = cD,mess cD(Re) für Re > 1 ⋅ 106 herangezogen.

Für die Demonstration der Zuverlässigkeit der Extrapolation wird hier der Bereich für Re > 1 ⋅ 106 herangezogen, ab dem turbulente Grenzschichten vorherrschen. Dazu werden die Residuen ΔcD = cD,mess cD(Re) gebildet. Dies entspricht der Abweichung zwischen gemessenen und berechneten Werten. Bild 2 zeigt alle Residuen für 30 Düsen mit Halsdurchmessern dHals von 2,2 mm bis 46 mm, für die PTB-Kalibrierdaten sowohl im Bereich kleiner als auch ausreichend hoher Reynoldszahlen vorliegen. Der Mittelwert ΔcD,mittel aller dieser Residuen liegt nahezu bei Null und die doppelte Standardabweichung bei ca. 0,0015, was in etwa einem relativen Wert von 0,15 % entspricht. Berücksichtigt man die relative Messunsicherheiten der jeweilige Messwerte cD,mess, die ebenfalls im Bereich 0,08 bzw. 0,25 % liegen (je nach verwendeter Referenz bei der Kalibrierung), so ergibt sich bereits eine gute Übereinstimmung zwischen Messung und Berechnung ohne die Berücksichtigung einer zusätzlichen Unsicherheit für die Extrapolation. Die für die Anpassungen ermittelten Krümmungsradien RHals reichen dabei von 0,68 ⋅ dHals bis 21,1 ⋅ dHals, was eine große Spannbreite an cD(Re)-Verläufen abdeckt.

Bild 2:    Darstellung der Residuen ΔcD = cD,mess cD(Re) für alle 30 Düsen der Datenbasis sowie der Mittelwert der Residuen ΔcD,mittel und dessen doppelte Standardabweichung.

Eine Abschätzung der Unsicherheit aus der Extrapolation erfordert noch eine genauere Analyse der Datensätze. Darüber hinaus gilt es, die Mindestanforderungen an die Düsen zu bestimmen, bei denen eine solche Extrapolation verlässlich möglich ist. Derzeit ist davon auszugehen, dass dieser Ansatz nur für hydraulisch glatte Düsen gültig ist und Abweichungen von der idealen Düsenform nach ISO9300 ein bestimmtes Maß nicht überschreiten dürfen.

Literatur:

[1] B. Mickan, J.-P. Vallet, C. Li and J. Wright; Extended data analysis of bilateral comparisons with air and natural gas up to 5 MPa; FLOMEKO 2016, Sydney, Australia, September 26-29, 2016, Opens external link in new windowLink
[2] I. M. Hall, Transonic Flow in Two-Dimensional and Axially-Symmetric Nozzles, Quarterly Journal of Mechanics and Applied Mathematics, Vol. XV, Pt. 4, pp. 487-508, 1962, Opens external link in new windowLink
[3] Geropp, D., Laminare Grenzschichten in Überschalldüsen, Deutsche Luft- und Raumfahrt, Forschungsbericht 01 TM 8603-AK/PA 1, 1987.
[4] B. Mickan, R. Kramer, D. Dopheide, H.-J. Hotze, H.-M. Hinze, A. Johnson, J. Wright, J.-P. Vallet; Comparisons by PTB, NIST, and LNE-LADG in air and natural gas with critical venturi nozzles agree within 0.05 %; 6th ISFFM 2006, Queretaro, Mexico, May 16-18, 2006, Opens external link in new windowLink

Ansprechpartner:

Dr. Bodo Mickan, FB 1.4, AG 1.43, E-Mail: Opens window for sending emailbodo.mickan(at)ptb.de