溫度對氣體標定裝置中噴嘴的特性影響

鄭 哈，朱 云，蔡 晴，鄭建英

(1.中國計量學院機電工程學院，浙江杭州 310018；2.浙江省計量科學研究院，浙江杭州 310018)

?

溫度對氣體標定裝置中噴嘴的特性影響

鄭 哈1，朱 云1，蔡 晴1，鄭建英2

(1.中國計量學院機電工程學院，浙江杭州 310018；2.浙江省計量科學研究院，浙江杭州 310018)

研究溫度對音速噴嘴內(nèi)部流場影響問題。當前研究對噴嘴性能產(chǎn)生影響的因素，研究濕度、入口段形狀等的較多，而研究溫度的較少。目前無法對噴嘴內(nèi)部流場進行檢測，文中提出一種計算流體力學方法，對不同溫度下噴嘴內(nèi)部流場進行數(shù)值模擬，研究不同溫度下噴嘴內(nèi)部流場情況。仿真結果表明噴嘴模型是可靠的，反應出溫度變化對噴嘴特性的影響。同時仿真結果對音速噴嘴標定裝置校準方法的研究有一定的參考價值。

音速噴嘴；膨脹系數(shù)；流出系數(shù)；流場仿真

0 引言

音速噴嘴又名臨界流文丘里噴嘴，主要應用于流量標準的傳遞、氣體流量測量和流量系統(tǒng)最大流量的限制。音速噴嘴是音速噴嘴氣體標定裝置[1]的主要組成部分，它選用符合國家標準ISO9300《臨界流文丘里噴嘴測量氣體流量》的臨界流文丘噴嘴作為標準。該裝置采用標準表法對多種流量計進行流量檢定和測試。音速噴嘴的性能直接影響到氣體標定裝置的檢定準確度，影響被檢流量計的準確度。

目前對噴嘴性能影響因素的研究比較多的有濕度對音速噴嘴性能影響的實驗研究及濕度修正方法以及入口段、擴散段形狀對流出系數(shù)影響的仿真研究，但是關于溫度對噴嘴性能影響的研究很少。音速噴嘴的檢定和喉部直徑測量是在恒溫室中進行的(溫度恒定在20℃)，裝置的檢定溫度會隨環(huán)境改變，溫度變化對噴嘴性能的影響直接影響裝置檢定準確度，因此研究溫度對噴嘴性能影響有重要意義。

噴嘴內(nèi)部的流場不能直接檢測，本文利用FLUENT流體仿真軟件仿真不同溫度下噴嘴內(nèi)部流場情況，對20℃噴嘴的速度場和壓力場進行分析，并與實測值進行比較，驗證了噴嘴模型的可靠性。對仿真結果進行比較和分析，得到溫度對噴嘴特性影響情況，仿真結果對對裝置校準方法的研究有一定的參考價值。

2 噴嘴原理

2.1 噴嘴輪廓

根據(jù)ISO9300，文丘里噴嘴的幾何形狀和使用條件是使其喉部流量達到臨界值的關鍵。本文所研究的為圓環(huán)形(也稱喇叭形)喉部文丘里噴嘴，其結構圖應符合圖1的規(guī)定。圖1為入口平面是入口直徑2.5d±0.1d處垂直于對稱軸線的平面，與曲率半徑為1.8d到有的結合處。3為擴散段，是一個具有2.5°至6°半角的平截頭圓錐體，擴散段的長度應不小于喉部直徑。

圖1 圓環(huán)形喉部噴嘴

上游滯止壓力為P0，噴嘴出口壓力為P2，當上游滯止壓力P0保持不變，逐漸減小噴嘴出口壓力P2，即減小背壓力比(P2/P0)，通過噴嘴的氣體流量先是不斷增加。當達到臨界壓力比時，氣體在喉部達到臨界速度即等于當?shù)匾羲?，流過噴嘴的氣體質量流量達到最大[2]。

2.2 噴嘴流場參數(shù)

2.2.1 流出系數(shù)

流出系數(shù)[3]是表明噴嘴特性的一個重要指標。它表示實際流量與相同滯止條件下，氣體一維、等熵地流過噴嘴的理想流量的比值。流出系數(shù)的計算公式如下：

Cd=qm/qmi

(1)

式中:qm為通過噴嘴的實際流量；qmi為相同滯止條件下，氣體一維、等熵地流過噴嘴的理想流量。

2.2.2 理想條件下的臨界流量

將氣體看成理想氣體，即忽略流體粘度的影響，可以導出臨界流噴嘴的流量為：

(2)

式中:A為噴嘴喉部內(nèi)截面積；R為氣體常數(shù)；M為氣體摩爾質量；P0、T0分別為噴嘴前氣體滯止壓力、滯止溫度；C*i為理想條件下的臨界流函數(shù)。

(3)

式中γ為絕熱指數(shù)，對理想氣體，它等于比熱比，取1.4。

2.2.3 實際條件下的臨界流量

實際真實氣體是有粘性的，通過噴嘴時會產(chǎn)生能量損失。所以將理想條件下的流量乘以一個系數(shù)C才能得到音速噴嘴的實際流量，也就是式(1)中所說的。

(4)

由于噴嘴的內(nèi)部流場不能直接測量得到，因此在建立數(shù)學模型的基礎上，對噴嘴內(nèi)部流場進行仿真，更加形象的展現(xiàn)內(nèi)部流場變化情況并得到噴嘴流場參數(shù)。

3 數(shù)學模型

噴嘴所選用的材料應有較高的硬度并有較好的耐腐蝕性，還應有穩(wěn)定的可復現(xiàn)的熱膨脹特性。本文選用名義流量為0.25 m3/h的圓環(huán)形喉部文丘里噴嘴為研究對象。所用材料為不銹鋼304L[4]，是一種通用性的不銹鋼材料，有良好的耐腐蝕性和成型性，本文主要考慮的是噴嘴的徑向線膨脹，線膨脹系數(shù)表示固體物質的溫度每改變1℃時，其長度的變化和它在0℃時長度之比。其計算公式如下

α=ΔL/(L·ΔT)

(5)

式中：ΔL為物體長度變化；ΔT為物體溫度變化；L為初始長度。膨脹系數(shù)見表1。

表1 不銹鋼304L線膨脹系數(shù)

由表1可以看出，在不同溫度范圍內(nèi)不銹鋼的線膨脹系數(shù)是不同的，本文所考慮的溫度范圍在0～100℃之內(nèi)，因此選取根據(jù)表1選取線膨脹系α=16.0×10-6/℃。

所選用圓環(huán)形噴嘴在20℃時的喉部直徑利用三坐標測量機測得為0.6 774 mm?；趪娮觳牧系臒崦浝淇s性能，溫度變化對噴嘴喉部尺寸產(chǎn)生影響。由于實驗室恒溫條件的限制，不能直接測量出不同溫度下的噴嘴直徑，因此根據(jù)式(5)通過計算得到不同溫度下噴嘴喉部尺寸如下表2所示。

表2 溫度與對應溫度下的噴嘴喉徑

由表2可得，喉部直徑與溫度成線性關系，且隨著溫度的升高而增大。

4 物理模型與仿真

4.1 物理模型的建立

本文選取名義流量為0.25 m3/h的圓環(huán)形臨界文丘里噴嘴作為研究對象。常溫下噴嘴喉部直徑為0.6 774 mm，入口段直徑為25 mm，出口段直徑為12 mm，擴散段角度為4°。根據(jù)圖1噴嘴的結構圖，分為建立不同溫度下噴嘴的物理模型，分別為模型1、模型2、模型3、模型4、模型5。

4.2 網(wǎng)格單元的劃分[5]

在二維(2D)空間中，可以使用的單元有三角形單元和四邊形單元。四邊形網(wǎng)格可以準確地處理邊界條件，計算精度高，但對復雜外形的網(wǎng)格較難生成，塊與塊的界面處理也十分復雜；三角形網(wǎng)格能夠方便地生成復雜外形的網(wǎng)格，但所需內(nèi)存更大、計算周期更長，并且會給計算精度帶來一定損失。本文在盡量保證計算精度的要求下，使用的網(wǎng)格格式單元是quad，指定的格式類型是pave，表明是拓撲結構的四邊形網(wǎng)格。為了保證計算精度又能提高計算效率，在噴嘴喉徑附近區(qū)域劃分相對較密的網(wǎng)格，而入口處和出口處區(qū)域的網(wǎng)格劃分相對較疏。

4.3 仿真模型的建立

為了能與實際測量的參數(shù)進行比較，驗證仿真模型的正確性，首先選擇20℃下的模型3進行仿真。

本文模擬的是理想氣體，求解器選擇為Pressure Based，計算方式選擇隱式?？諝庠谝羲賴娮熘袨橥牧髁鲃樱x擇Spalart-Allmaras湍流模型，該模型僅考慮了動量的傳遞方程。在氣體動力學中，對于有固動及壓力梯度較大的流動都可得到較好的結果。Spalart-Allmaras的模型方程為：

(6)

4.4 設定邊界條件

Fluent軟件提供了十余種類型的進、出口邊界條件，由于噴嘴入口處出口處的壓力為已知，因此入口邊界條件設置為壓力入口，出口邊界條件設置為壓力出口。影響噴嘴質量流量的為背壓比，只要正確控制背壓比即可。因此入口壓力設置為1個大氣壓(1 atm)，入口湍流邊界條件設置根據(jù)經(jīng)驗值湍流強度設置為1，水力直徑為25 mm，出口壓力設置為0.2 atm，湍流強度設置為5，水力直徑為12 mm。出入口的溫度均設為293.15 K即20℃。

5 仿真結果及分析

根據(jù)模型3仿真20℃時噴嘴內(nèi)部的流場，得到噴嘴前后部分管道及噴嘴內(nèi)部的速度矢量云圖如圖2。由圖2可以看出，理想氣體在噴嘴入口處加速，速度達到了200 m/s，在噴嘴喉徑處速度約338 m/s，即達到了當?shù)匾羲佟Ｔ跀U散段速度繼續(xù)增加，在接近噴嘴出口處，速度達到了400 m/s。由于回流作用，在其后區(qū)域內(nèi)速度降低。

圖2 噴嘴內(nèi)流速分布圖

從整個速度分布圖上看，擴散段速度分布區(qū)域化較為明顯，噴嘴喉部處的速度場區(qū)域化沒有擴散段明顯。但從速度矢量云圖中可明顯地看出其靠近噴嘴壁面的速度低于主流區(qū)域的速度，這是理想氣體的粘性作用導致的，氣體與壁間的黏滯作用使得靠近壁面流速大大低于靠近噴嘴中心流速。

根據(jù)模型3仿真20℃時噴嘴內(nèi)部的流場，得到噴嘴前后部分管道及噴嘴內(nèi)部的壓力云圖如圖3。

圖3 噴嘴內(nèi)部壓力分布圖

由圖3可以看出，噴嘴入口前連接在滯止容器上的管道的內(nèi)部壓力為1 atm，在噴嘴入口處噴嘴口徑驟縮而產(chǎn)生回流作用使其壓力發(fā)生改變，噴嘴出口后的管道壓力分布均勻，為0.2 atm。從入口處開始噴嘴內(nèi)部的壓力逐漸降低，最低的區(qū)域在擴散段為0.14 atm從管道入口到擴散段出現(xiàn)壓力最低的區(qū)域，其壓力損失接近86%。這是噴嘴出口處出現(xiàn)高速射流現(xiàn)象的根本原因。由圖可以看到在擴散段，壓力為0.14 atm的區(qū)域非常大，縮小該區(qū)域便可以降低壓力損失[6]，從而實現(xiàn)減小射流對噴嘴后管道流場的影響。

由于本文采用的模型是二維的，不能直接從軟件中讀出三維時的氣體質量流量?？捎萌缦滤闶睫D換

(7)

根據(jù)式(7)計算20℃下噴嘴喉徑口的平均質量流量為118.2 kg/s，其中理想氣體密度ρ=1.205 kg/m3，根據(jù)式(1)計算得Cd=0.94 238，再根據(jù)式(2)可以計算得到噴嘴的理論流量，20℃下噴嘴的理論流出系數(shù)、實測流出系數(shù)、仿真流出系數(shù)見表3。

表3 20℃時各流出系數(shù)比較

將實測流出系數(shù)、理論流出系數(shù)和仿真流出系數(shù)三者之間進行比較，由表4可以看出，流出系數(shù)之間的差值不大，仿真的流出系數(shù)比理論的流出系數(shù)的誤差要大一些，造成這個的原因主要有：仿真模型的精度，不能非常精確的反映噴嘴內(nèi)部的結構；計算模型的選擇，以及計算過程都會產(chǎn)生誤差。因此可以認為以上所建模型是正確的，在其他溫度條件下，可以應用以上模型進行仿真。

根據(jù)模型3的仿真步驟，分別對另外4個模型1、2、4、5進行仿真，仿真溫度依次設置為283.15 K，288.15 K，298.15 K，303.15 K。仿真后再根據(jù)式(7)進行計算得到噴嘴喉徑處平均質量流量和平均體積流量見表4。

由于噴嘴的標況流量會隨環(huán)境的溫度、壓力而改變，名義流量與標況流量允許的誤差為±5%，仿真所得與名義流量差值最大的標況流量為0.238 m3/h,與名義流量的誤差為-4.8%，在允許范圍內(nèi)。

表4 仿真噴嘴流量

噴嘴的理論流量可由式(2)計算得出，計算所得的理論流量如下表5所示。

表5 各溫度下的理論流量

再根據(jù)式(1)計算各溫度下的流出系數(shù)見表6。

表6 各溫度下流出系數(shù)

由表6可知，溫度低于標準溫度20℃時，流出系數(shù)變小，且溫度越低，流出系數(shù)越小，與常溫20℃下的誤差越大，溫度高于20℃時，流出系數(shù)偏高，且溫度越高，誤差越大，由此可知，不論溫度偏高還是偏低，都會對檢測的流出系數(shù)產(chǎn)生較大的影響。

6 結束語

本文利用FLUENT軟件對音速噴嘴進行仿真，由以上仿真得出的結果可知，20℃時仿真結果對比實測值說明仿真模型能夠準確反映內(nèi)部流場變化。在這個模型的基礎上，分別仿真了10℃、15℃、25℃、30℃時的不同模型，仿真結果說明溫度變化導致噴嘴的喉徑變化，最終導致了噴嘴內(nèi)部質量流量的變化，喉徑越大質量流量越大，得到的流出系數(shù)越大，溫度偏高或偏低都會對流出系數(shù)產(chǎn)生影響，為了保證噴嘴檢定過程的準確性，實際應用過程中，要盡量保持溫度恒定或者對噴嘴喉徑進行適當?shù)男拚齕7]。

[1] 王順利，李曉宇，鄭云林，等.負壓法臨界流文丘里噴嘴氣體流量標準裝置的研究及應用.工業(yè)計量，2013，23(5):44-47,53.

[2] 吳明清，劉慶華.臨界流文丘里噴嘴的最大允許背壓比.現(xiàn)代計量儀器與技術，2012(6):77-78.

[3] 宋德星.臨界流文丘里噴嘴測量不確定度分析.工業(yè)控制計算機，2013，26(12):59-60.

[4] 馬坤明，羅宇，王江超.304L不銹鋼結構焊接變形分析.焊接學報，2010，31(1):55-58.

[5] 陳斌，郭烈錦.非結構化網(wǎng)格快速成型技術.西安交通大學學報，2000，34(1):18-21,29.

[6] 楊友勝，馮輔周，朱玉泉，等.噴嘴內(nèi)壁壓力分布的研究.機械工程學報，2011，47(16):182-186.

[7] 呂運朋，禹金標，薛好.基于多功能氣體檢測儀的溫度補償及數(shù)據(jù)融合.儀表技術與傳感器，2013(11):31-34.

Temperature Effect on Sonic Nozzle Properties of Gas Calibration Device

ZHENG Ha1,ZHU Yun1,CAI Qing1,ZHENG Jian-ying2

(1.College of Mechanical and Electrical Engineering,China Jiliang University,Hangzhou 310018,China; 2.Zhejiang Province Institute of Metrology,Hangzhou 310018,China)

This paper studied the temperature effect on the flow field of the sonic nozzle.The current researches study the factor that affect the performance of the nozzle,that are mostly about humidity and entrance shape but less about temperature.Currently the flow field inside the nozzle is unable to detect,so this paper presented a method for computational fluid dynamics,simulated the flow field inside the nozzle under different temperatures,and researched on the flow field inside the nozzle under different temperature conditions.Simulation results show that the nozzle model is reliable,reflecting the effect of temperature changes on the nozzle characteristics.The simulation results also have some reference value for sonic nozzle calibration device calibration method.

sonic nozzle; expansion coefficient; discharge coefficient; flow field simulation

2014-09-22 收修改稿日期：2015-02-10

TP391.9

A

1002-1841(2015)07-0101-03

鄭哈(1991—)，碩士研究生，主要研究領域為蒸汽管網(wǎng)，流量檢測。E-mail:635429246@qq.com 朱云(1960—)，教授，主要研究領域為蒸汽管網(wǎng)，傳感器及信號處理，流量檢測。E-mail: yun7882@sohu.com