變面積臨界流文丘里噴嘴特性研究*

陳元杰， 葛　銳， 姚超玲， 嚴東旭, 沈文新

(1.浙江省計量科學研究院， 浙江 杭州 310000；2.浙江大學 海洋學院，浙江 杭州 310000)

?

變面積臨界流文丘里噴嘴特性研究*

陳元杰1， 葛銳1， 姚超玲2， 嚴東旭2, 沈文新1

(1.浙江省計量科學研究院， 浙江 杭州 310000；2.浙江大學 海洋學院，浙江 杭州 310000)

摘要:針對目前氣體流量計量領域臨界流文丘里噴嘴臨界流量不可調節或調節困難的情況，提出一種喉部過流面積可調節的變面積臨界流噴嘴。通過結構分析、建模計算、流體力學Fluent仿真等方式，對棒狀節流體、錐狀節流體、臺錐狀節流體三種不同調節方式進行研究分析。結果表明：采用三種不同節流體的變面積臨界流文丘里噴嘴均能實現臨界流流量的調節，在背壓比小于最大允許背壓比后，流量保持一穩定值，基本不變。建模計算得到的理論值與Fluent仿真到的數據較為符合。

關鍵詞:計量學； 臨界流； 變面積； Fluent仿真； 文丘里噴嘴

0引言

臨界流文丘里噴嘴作為流量計量裝置中的核心部件，具有良好的穩定性，原理簡單可靠。因此，其作為流量控制儀表或標準裝置中的標準表而得到廣泛應用。同時，其又作為量值傳遞標準器，在流量裝置和流量儀表的量值傳遞，特別是在天然氣計量工作中發揮著重要作用。

目前，臨界流文丘里噴嘴主要存在問題是流量測量范圍為“1”，即一個噴嘴只是一個流量，對流量的調節只能通過更換不同喉徑的噴嘴或者多個噴嘴并聯安裝的方式實現[1]，這在工業應用中往往較為困難，從而使其推廣應用受到極大的限制。

本文通過研究、分析單個臨界流文丘里噴嘴“多”流量調節的實現方法和可行性，設計一種喉部節流面積可調節的變面積臨界流文丘里噴嘴，并通過數學建模和計算流體動力學(computational fluid dynamics,CFD)流體仿真對其特性進行研究，以期代替一系列不同喉部面積的臨界流噴嘴的并聯安裝使用，從而大大提高量值傳遞和工業氣體流量計量的效率成本。

1結構分析

1.1傳統臨界流文丘里噴嘴結構

傳統臨界流文丘里噴嘴結構如圖1所示，其主要結構分為入口平面、入口收縮段、喉部、出口擴散段等幾個部分[2]。

圖1　傳統臨界流文丘里噴嘴結構原理圖Fig 1　Principle diagram of structure of traditional critical flow Venturi nozzle

工作時，通過減小出口壓力p2與入口滯止壓力p0之比p2/p0(背壓比)，使其減小到某一數值時，喉部流速達到一固定值(為當地音速)，當背壓比再減小時，喉部流速不再改變，氣體質量流量保持恒定，為臨界流狀態。臨界流狀態氣體質量流量可由式(1)計算得到[3]

(1)

式中Ant為喉部過流面積；Cd′和C*為流出系數和臨界流函數，可通過經驗公式計算、查表、實測等方法得到，在一定溫度和流體屬性下為常數；p0為上游滯止壓力；T0為上游滯止溫度，R為氣體常數，M為氣體摩爾質量。因此，當氣體屬性、上游滯止條件(T0，p0不變)一定的情況下，流量qm值與過流面積Ant呈正比，而常規的臨界流噴嘴喉部截面積為固定形式，無法進行調節。

1.2變面積臨界流文丘里噴嘴結構

根據上述問題，本文設計了三種可行的變面積臨界流文丘里噴嘴結構，通過在常規噴嘴喉部插入節流體的形式來調節喉部過流面積，從而實現噴嘴在臨界流狀態下流量調節。具體結構如圖2～圖4所示。

圖2為棒狀節流體結構，在常規噴嘴中插入一定直徑的圓柱棒來實現。根據所需臨界流流量的不同，更換不同直徑的圓柱棒來實現流量的調節[4,5]。優點是結構簡單，噴嘴內部流場相對于插入前變化不大；缺點是每次調節必須更換不同直徑的圓柱棒，較為麻煩。但是相對于更換不同喉部直徑的噴嘴來說，已經提高了效率，節約了成本。

圖3　文丘里噴嘴+錐狀節流體Fig 3　Venturi nozzle with cone-shaped throttle body

圖3為錐狀節流體結構，在常規噴嘴中插入一圓錐體來實現。通過改變錐狀節流體在喉部的位置來調節喉部過流面積的大小，從而實現臨界流流量的調節。其優點是調節較為方便，并可進行連續調節；缺點是由于節流體是圓錐狀，噴嘴內部流場變化較大，同時實際操作時外部驅動機構很難將圓錐體定位于同一位置，使其重復性較差。

圖4為臺錐狀節流體結構，在常規噴嘴中插入具有一定臺階數的臺錐體來實現。在結構上，其結合了棒狀節流體和錐狀節流體的特點，不僅調節方便，而且同一臺階插入到喉部時的臨界流流量都相同，提高了其定位精度和重復性。其缺點是:流量調節為有極調節，可調流量數目和臺階數目相同。

圖4　文丘里噴嘴+臺錐狀節流體Fig 4　Venturi nozzle with hybrid-shaped throttle body

2建模與計算

變面積臨界流文丘里噴嘴模型從結構上看，為在傳統臨界流文丘里噴嘴喉部插入一定結構的節流體，通過改變喉部過流面積來調節流量值。根據式(1)可得

(2)

式中AR為節流體在喉部的截面積。建模計算時，取噴嘴喉部直徑d=50mm，氣流為干燥空氣，其摩爾質量M=28.8g/mol，密度ρ=1.205kg/m3，大氣壓力p0=101 325Pa，溫度T0=293.15K，氣體常數R=8.314J/(mol·K)。臨界流函數C*和流出系數Cd′根據查表得到,分別為0.685 131 3和0.98。以下結合式(2)和上述參數值，分別計算三種節流體對噴嘴臨界流流量的調節。

2.1棒狀節流體模型

該模型較為簡單，建模時，插入喉部的圓棒半徑為rr=8.8，16.1，21.0mm。結合上述參數，根據式(2)計算可得該模型的相應喉部過流面積和臨界流流量，見表1。(選取合適的圓棒半徑rr，使計算出來的體積流量接近常用流量，下同)。

表1　文丘里噴嘴+棒狀節流體計算參數

2.2錐狀節流體模型

該模型的結構參數如圖5所示。

圖5　錐狀節流體模型結構參數(mm)Fig 5　Structure parameters of cone-shaped throttlebody model(mm)

圖中，錐體半角為6°，頭部半球半徑為7.5 mm。噴嘴出口擴散段半角為10°，入口收縮段圓環面截面半徑為100 mm。設噴嘴喉部位置為原點，節流體從右至左插入喉部距離為l。錐狀節流體插入后，噴嘴喉部位置(過流面積最小處)發生改變，從原喉部位置向擴散段方向偏移一小段距離。由圖5的幾何關系，可得沿x軸各處過流面積為(x取值為0到入口收縮段圓弧與擴散段切點)

(3)

喉部面積S(x)最小的位置出現在dS/dx=0處，對式(3)求導可得

(4)

在確定節流體位置l后，可以根據式(4)獲得喉部位置x0，代入式(3)后，可得新喉部過流面積S(x0)。

取l=10.6，49.4，78.5，102.8，124.1mm，經過式(2)～式(4)計算，可得出新喉部位置、喉部過流面積、臨界流流量值等參數，如表2所示。

表2　文丘里噴嘴+圓錐狀節流體計算參數

2.3臺錐狀節流體模型

該模型結構參數如圖6所示。

圖6　臺錐狀節流體模型結構參數(mm)Fig 6　Structure parameters of hybrid-shapedthrottle body model(mm)

圖中，各個臺階的半徑為Rs，臺階長度為ls，臺階間用6°錐面過渡，其他尺寸與錐狀節流體相同。由于節流體為臺錐狀，若臺階過短錐面過長，也會使過流面積最小的喉部向擴散段方向偏移，偏離原喉部位置。為保證最小過流面積在原喉部處，必須使臺階足夠長。根據圖6的結構，可得各個臺階后錐面處的過流面積為

(5)

求導后為

(6)

式中當x=ls時，能保證過流截面積在經過原喉部后，一直處于擴展狀態，使最小過流面積保持在原喉部。將x=ls，Rs數值(Rs=8.8，13，16.1，18.8，21mm)依次代入式(7)，可計算得到喉部不偏移的最小臺階長度ls=3.69，5.43，6.69，7.79，8.67mm。取最大值8.67mm(取整為9mm)作為臺階基礎長度，在此基礎上,再加3mm作為重復性定位誤差，因此，臺錐狀節流體每個臺階的長度為12mm。

由于此時喉部不偏移，為原喉部位置，根據臺階半徑Rs，喉部直徑d值，經式(2)計算可得表3相應參數。

表3　文丘里噴嘴+臺錐狀節流體計算參數

3CFD仿真研究

本文采用FLUENT作為CFD分析工具,使用的網格模型采用ICEM CFD軟件生成，仿真過程中選擇壓力基求解器，采用Spalart-Allmar湍流模型[6,7]。噴嘴入口平面上游設置長度、直徑分別為5d(d為噴嘴喉部直徑)的管道流場。入口邊界條件設置為為一個大氣壓，出口壓力根據臨界背壓比取多個點進行觀測。

該仿真流場為嚴格對稱的形狀，為了提高仿真效率，將流場沿中軸線對稱剖出1/4進行仿真，剖面設置為對稱面[8]。限于篇幅，本文只對質量流量隨背壓比的變化特性進行仿真分析。以下仿真分析中，水平直線為理論流量值(其值可參考表1～表3中的qm)，曲線為不同節流體在喉部不同位置時的仿真流量值。

3.1棒狀節流體仿真分析

將半徑為r=8.8，16.1，21 mm三種圓柱棒插入到常規噴嘴中進行仿真，背壓比在0.9～0.5間變化，仿真結果如圖7。

圖7　棒狀節流體模型流量—背壓比曲線圖Fig 7　Curve of flow as a function of back pressure ratio of rod-like throttle bady model

從圖中可以看出:流體在背壓比不斷減小的過程中，達到了臨界流體狀態,其最大允許背壓比在0.8～0.9之間(1和2約為0.8，3約為0.76)。此時，背壓比再減小，質量流量基本不變。

3.2錐狀節流體仿真分析

根據圖5的結構，將圓錐體插入噴嘴喉部不同深度進行仿真，分別為l=10.6，49.4，78.5，102.8，124.1 mm，并在不同的背壓比下，觀察其流量變化。仿真結果與理論計算結果對比如圖8所示。

從圖中可以看出:各個插入位置的氣體流量均能達到臨界流狀態，其最大允許背壓比在1，2，3位置約為0.84，第4個位置約為0.8，第5個位置約為0.76。此后，隨著背壓比的減小，流量值基本不變。

3.3臺錐狀節流體仿真分析

根據圖6的結構，將半徑r=8.8，13，16.1，18.8，21 mm的臺階依次插入噴嘴喉部，并在不同背壓比下，觀察其流量變化。仿真結果與理論計算結果對比如圖9所示。

從圖中可以看出:各個臺階位置的氣體流量均能達到臨界流狀態，其最大允許背壓比在1，2，3位置約為0.85，第4個位置約為0.80，第5個位置約為0.74。此后，隨著背壓比的減小，質量流量值基本不變。

圖8　錐狀節流體模型流量—背壓比曲線圖Fig 8　Curve of flow vs back pressure ratio of cone-shaped throttle body model

圖9　臺錐狀節流體模型流量—背壓比曲線圖Fig 9　Curve of flow vs back pressure ratio of hybrid-shapedthrottle body model

3.4現象分析

將氣體臨界流狀態下質量流量仿真值與理論計算值進行比較可知，兩者的大小基本一致，但隨著過流面積的增加(插入喉部的節流體直徑變小)，誤差有所增加。分析其原因，可能是在進行理論計算時采用的流出系數Cd′不夠準確。本文為計算方便，Cd′值統一為0.98，而實際上，該值隨著喉部雷諾數Re增大而增大，而Re會隨著過流面積的增加而增大。因此，隨著過流面積的增加，實際理論計算值應更大一些。Cd′與Re對應關系如表4所示。

表4　流出系數與喉部雷諾數對應關系

4結論

本文通過建模分析計算、CFD仿真的相互印證可知，變面積臨界流文丘里噴嘴在背壓比小于最大允許背壓比時，其流場能夠達到臨界流狀態。棒狀節流體流場簡單穩定，理論計算方便，但更換較為麻煩；錐狀節流體調節方便，其最小喉部位置較原喉部位置有所偏移，但重復性較差；臺錐狀節流體結合了棒狀節流體和錐狀節流體的特點，不但調節方便，而且具有較好的重復性。仿真雖與理論計算有一定誤差，但是通過查表對Cd′值的修正，或通過PVTt等設備的實流標定，得到準確的Cd′，能進一步提高計量的準確性。

參考文獻:

[1]ISO/TC 30/SC 2.ISO 9300—2005.Measurement of gas flow by means of critical flow Venturi nozzles[S].

[2]GB/T 21188—2007.用臨界流文丘里噴嘴測量氣體流量[S].

[3]梁國偉,周寧寧,李長武.臨界流文丘里噴嘴流量計的原理與應用[J].中國計量學院學報，2004,15(3):186-190.

[4]Kim J H,Kim H D,Park K A.Computational/experimental study of a variable critical nozzle flow[J].Flow Measurement and Instrumentation,2006,17:81-86.

[5]Kim H D,Kim J H,Park K A,et al.Computational study of the gas flow through a critical nozzle [J].Journal of Mechanical Engineering Science,2003,217(10): 1179-1189.

[6]王福軍.計算流體動力學分析—CFD軟件原理與應用[M].北京：清華大學出版社，2004:85-88.

[7]胡坤，李振北.ANSYS ICEM CFD工程實例詳解[M].北京：人民郵電出版社，2014:55-60.

[8]王麗辰,朱云,鄭哈,等.基于 Fluent的臨界流文丘里噴嘴的內部流場仿真分析[J].科學技術與工程,2013,13(34):10392-10396.

Research on characteristics of variable area critical flow Venturi nozzle*

CHEN Yuan-jie1， GE Rui1， YAO Chao-ling2， YAN Dong-xu2， SHEN Wen-xin1

(1.Zhejiang Institute of Metrology,Hangzhou 310000,China;2.Ocean College,Zhejiang University,Hangzhou 310000,China)

Abstract:In order to regulate critical flow rate of traditional critical flow Venturi nozzles,a new variable throat-area critical flow Venturi nozzle is proposed,throat area of which can be adjusted.Three different adjusting methods,including rod-style method,cone-style method and hybrid-style method,are analyzed by means of structure analysis,modeling calculation and fluid dynamics Fluent simulation.Results show that the critical flow of the new variable critical flow Venturi nozzles can be regulated well by using three different methods and when back-pressure ratio is smaller than the maximum allowable back pressure ratio,the flow rate of the new nozzle is almost kept constant,though the back-pressure ratio is still descending.And Fluent results coincides with the theoretical value calculated by modeling.

Key words:metrology； critical flow； variable area； Fluent simulation； Venturi nozzle

DOI:10.13873/J.1000—9787(2016)03—0044—05

收稿日期：2015—06—04

*基金項目：國家質檢總局一般項目(2013QK060)

中圖分類號:TH 814

文獻標識碼:A

文章編號:1000—9787(2016)03—0044—05

作者簡介:

陳元杰(1984-)，男，浙江慈溪人，碩士，工程師，主要從事氣體流量計性能研究。