彭 楓，李志鵬，張程鈔，廖志芳，王東福，謝 科

（1.長沙理工大學 能源與動力工程學院，長沙 410004；2.博納斯威閥門股份有限公司，天津 301802）

0 引言

空氣閥在輸送管線系統(tǒng)中重要的安全元件，在坡度起伏較大的輸送系統(tǒng)中，空氣閥的主要作用是在輸水管道的排水階段從外部吸入空氣，避免管道由于真空而造成坍塌，在輸水管道的充水階段排出管道中的空氣，防止管道中的空氣對供水造成阻礙。空氣閥的另一個重要功能是防止水錘，水錘會在斷電、泵停機等事故情況下導致水柱分離和水柱彌合［1-3］。

ROBERTO等［4］指出在空氣閥內部應該放置緩沖裝置或者節(jié)流孔板，可以用來減少管道中空氣閥突然關閥會引起內部壓力劇烈變化。ZLOCZOWER［5］認為空氣閥在運行時的進排氣都會對防水錘效果有較大的影響，從而指出空氣閥進排氣口徑應該有適中的規(guī)格，太大或者太小都會影響空氣閥水錘防護效果，進而造成水錘現象加劇。CARLOS等［6］建立了一套基于西城/代爾夫特液壓的完整試驗模型，用于解決空氣閥各個參數中對于傳統(tǒng)數值模型的不確定性。OSCAR等［7］就關于管路中水流放空時遇到的問題做了深入的研究，發(fā)現在管路中設置空氣閥就可以解決管路水流放空而導致的負壓。BALACCO等［8-9］研究發(fā)現當下游空氣閥閥門打開部分時，排氣孔直徑對壓力峰值影響較小，其最大壓力大概是作用壓頭的1.5倍。杜建軍［10］就空氣閥在長距離輸水管道中的設置原則做了歸納總結，并研究了空氣閥的特性參數、種類類型、技術規(guī)范以及設置原則作用。胡建永等［11］經過實際工程分析，結果表明對于空氣閥的設置，從理論分析得到的關于空氣閥進排氣孔口直徑和布置間距要比依據經驗得出的值小。高潔等［12-13］模擬計算了空氣閥的進排氣特性。結果表明，空氣閥的流通面積和進排氣流量系數之乘積隨壓差的變化規(guī)律；并且通過不同壓差的流通面積和進排氣流量系數的乘積值，對不同進氣流量系數下的進排氣流通面積進行了計算。王玲等［14］研究發(fā)現，當選用微量排氣和大口徑進氣的空氣閥，能夠防止空管在注水時產生的水柱分離甚至彌合高壓水錘。徐放等［15］通過實際工程計算驗證了，傳統(tǒng)高速進排氣閥內部結構優(yōu)化后的水錘防護效果明顯是要優(yōu)于傳統(tǒng)的高速進排氣閥。

由于空氣閥浮球的重力小于浮球的氣動力，從而被快速吹起導致排氣口被堵塞，空氣閥內部的氣體無法排出，造成空氣閥失效。為了解決浮球的吹堵問題，對空氣閥的結構進行改進，然后模擬分析空氣閥排氣特性，計算分析不同護筒高度下輸入不同壓力下空氣閥排氣量大小和浮球升力大小，為浮球式空氣閥的設計提供依據。

1 數學模型

空氣閥的排氣和吸氣是一個極其復雜的動態(tài)過程，這2個不同的階段過程中涉及到氣液兩相流。當前對空氣閥的數值模擬依舊采用WYLIE等［16］和 STREETER 等［17］所提出的數學模型，該模型需建立在下列4個假設條件上［18-22］：

（1）空氣等熵地流入流出閥門；

（2）通過空氣閥進入管道中的空氣一直處于等溫狀態(tài)，由于管道內空氣相對較少，同時與管壁和液體的接觸面積很大，因此管道中空氣和液體的溫度相同；

（3）通過空氣閥進入管道內的空氣停留在空氣閥位置的管道周圍；

（4）液體表面高度基本不變，空氣體積與管段內液體體積相比很小。

基于以上4個假設，流過閥的空氣質量取決于管外大氣的絕對壓力P0、絕對溫度T0以及管內的絕對壓力P和絕對溫度T，分成以下4種情況建立模型。

（1）空氣以亞音速流入（P0＞P＞0.528P0）：

（2）空氣以臨界流速流入（P≤0.528P0）：

（3）空氣以亞音速流出（P0/0.528＞P＞P0）：

（4）空氣以臨界速度流出［（P＞（P0/0.528）］：

（5）空氣閥流量系數c為實際流量大小與同孔徑噴管理論流量大小之比，其表達式如下：

C1——進氣時空氣閥的流量系數；

A1——進氣時空氣閥的流通面積；

ρ0——大氣密度；

P/P0——壓力比；

C2——排氣時空氣閥的流量系數；

A2——排氣時空氣閥的流通面積；

R ——氣體常數；

c ——進排氣質量流量系數；

m ——通過空氣閥實際質量流量，kg/s。

2 數值模型

2.1 幾何模型

本文研究的空氣閥結構圖1所示，該空氣閥是一種復合式高速進排氣閥，根據其結構特點和工作原理，進行空氣閥模型處理，由于只關心空氣閥流動特性，因此應先對模型的非關鍵部位進行結構簡化，確保流速變化劇烈區(qū)域的幾何特征不變，對流體域中圓角倒角進行填平處理，采用Solidworks三維軟件按照1：1大小繪制三維模型，繪制后的三維模型簡化了空氣閥的結構，分析空氣閥排氣特性，通過對空氣閥護筒高度進行調整，將分為3種模型，護筒高度分別為285，380，440 mm。

圖1 空氣閥護筒模型Fig.1 Model diagram of air valve protective sleeve

2.2 網格與邊界條件

對簡化后的空氣閥三維模型采用Fluent meshing繪制多面體網格，劃分網格后的模型如圖2，3所示，逐步增加網格數量，當網格劃分后網格數量達到240萬時，空氣閥排氣速度和浮球升力均不再明顯增加，此時說明再增大網格數量對計算的精度不會有提高，因此對該空氣閥模型，采用網格數量為240萬時，網格數量達到最佳值，最小網格尺寸0.2 mm，最大網格尺寸10 mm，網格質量良好，能夠保證模擬計算相對準確。

圖2 空氣閥外殼網格劃分Fig.2 Meshing diagram of air valve housing

圖3 浮球閥座網格劃分Fig.3 Mesh division diagram of float valve seat

空氣閥在進排氣過程中流體域的工作介質為空氣，空氣閥的流體域介質設為理想氣體，并對空氣閥壁面設置為固定的無滑移壁絕熱壁面，由于空氣閥在運行時流速較大，是屬于湍流模型，因此對空氣閥采用標準的湍流模型進行計算，通過Fluent軟件對能量、動量、質量守恒方程進行求解。對空氣閥邊界參數設置，空氣閥護筒高285，380，440 mm 3組模型分別采用相同的處理方式，計算模型進口邊界選擇壓力入口分別設置35，70，100 kPa，出口邊界選擇壓力出口，出口相對壓力設置為0 Pa，對空氣閥3種模型進行FLUENT數值模擬。

3 流場計算及分析

對3組護筒高度空氣閥模型分別設置入口壓力35，70，100 kPa計算，收斂后取3組護筒高度空氣閥不同入口壓力的流體域壓力場、速度場進行分析。根據流體域云圖，主要得出的數據為空氣閥在排氣過程中排氣量大小、速度大小以及吹動氣流對浮球的升力大小，并且分析空氣閥排氣量、速度以及浮球升力的變化規(guī)律。

3.1 壓力場分析

圖4～6分別示出空氣閥護筒高度為285，380，440 mm 時，入口壓力分別為 35，70，100 kPa的壓力云圖。

圖4 空氣閥護筒285 mm壓力云圖Fig.4 Pressure nephograms at the height of the protective sleeve of 285 mm

圖5 空氣閥護筒380 mm壓力云圖Fig.5 Pressure nephograms at the height of the protective sleeve of 380 mm

圖6 空氣閥護筒440 mm壓力云圖Fig.6 Pressure nephograms at the height of the protective sleeve of 440 mm

當壓力為35 kPa時，空氣閥護筒高度285，380，440 mm 的最大壓力分別為 36.542，36.461，36.216 kPa；在排氣過程中，氣體流經浮球產生的壓差力與浮球的重力，二者的合力向下。由于空氣流繞過浮球閥座，浮球底部表面壓力比頂部表面壓力大，受壓較大的主要集中在浮球底部；當壓力為70 kPa時，空氣閥護筒高度285，380，440 mm的最大壓力分別為 73.093，72.865，71.863 kPa，在相同的排氣壓差下，隨著護筒高度的增加，空氣閥表面壓力逐漸降低。

3.2 速度場分析

圖7～9為空氣閥護筒高度為285，380，440 mm時，入口壓力分別為35，70，100 kPa的速度云圖。

圖7 空氣閥護筒285 mm速度云圖Fig.7 Velocity nephograms at the height of the protective sleeve of 280 mm

圖8 空氣閥護筒380 mm速度云圖Fig.8 Velocity nephograms at the height of the protective sleeve of 380 mm

圖9 空氣閥護筒440 mm速度云圖Fig.9 Velocity nephograms at the height of the protective sleeve of 440 mm

結合圖7～9可以看出，當壓力為35 kPa時，護筒高度285，380，440 mm的最大速度分別為210.74，296.50，284.30 m/s；當壓力為 70 kPa 時，護筒高度285，380，440 mm的最大速度分別為302.29，348.12，337.03 m/s，在相同的排氣壓差下，隨著護筒高度的增加，空氣閥內流場最大速度先增加后微量的減小，說明隨著護筒高度的增加空氣閥節(jié)流通徑變小。

4 排氣特性分析

根據GB/T 36523—2018《供水管道復合式高速進排氣閥》要求排氣閥的排氣量應不小于表1所列的數值。其中DN300空氣閥排氣壓差為35 kPa，排氣量為38 000 m3/h，排氣壓差為70 kPa時，排氣量為49 400 m3/h。從表2排氣量計算結果可得，空氣閥排氣壓差為35 kPa時，護筒高度285，380，440 mm均大于排氣量標準規(guī)范，偏差量分別為4.3%，3.4%，0.4%；空氣閥排氣壓差為70 kPa時，護筒高度285，380，440 mm的偏差分別為14.4%，13.5%，8.1%，滿足規(guī)范要求。

表1 排氣閥排氣量Tab.1 Air displacement of exhaust valve

表2 空氣閥排氣量計算結果Tab.2 Calculation results of air displacement of the air valve

從圖10可以看出，當護筒高度從285 mm增加到380 mm時，雖然隨著護筒高度的增加，空氣閥排氣量在逐漸降低，但空氣閥排氣量較為接近，變化不大；而護筒高度從380 mm增加到440 mm時，空氣閥排氣量變化較為明顯，隨著護筒高度的增加，排氣量顯著降低。在相同的護筒高度下，隨著排氣壓差的增大，空氣閥排氣量在逐漸增大，成正比；在相同的排氣壓差下，隨護筒高度的增加，排氣量逐漸減小，因此護筒高度不宜過高。

圖10 空氣閥排氣量隨排氣壓差變化Fig.10 Variation curve of air displacement of the air valve with exhaust pressure difference

排氣流量計算結果見表3。

表3 排氣流量計算結果Tab.3 Calculation results of exhaust flow

由表3可知，隨著壓差的增大空氣閥排氣量逐漸增大，在相同護筒高度下，實際質量流量、理論質量流量以及排氣流量系數都隨著排氣壓差的增加而增加；在相同壓差下的壓力比相同，但實際質量流量、理論質量流量以及排氣流量系數都隨著護筒高度的增加而減小。

空氣閥浮球升力變化如圖11所示。由圖11可知，空氣閥浮球升力全都為負值，說明力是向下的，排氣壓差為70 kPa時，空氣閥護筒高度為285，380，440 mm時，向下的力分別為312.63，449.7，335.4 N，在相同的排氣壓差下，隨護筒高度的增加，空氣閥浮球向下的力是先增加后減少的；空氣閥護筒高度為285 mm時，在35，70，10 kPa的排氣壓差下，浮球向下的力分別為312.63，675.6，885.6 N，在相同的護筒高度下，隨著排氣壓差的增大，空氣閥浮球向下的力在逐漸增大，說明排氣過程中浮球不會上升堵塞排氣孔。在工程實際應用中，隨著氣體的排出，空氣閥內液面不斷上升，水對浮球產生的浮力遠大于浮球向下的壓力，浮球會上升堵住排氣孔，防止液體流出，保證空氣閥安全運行。

圖11 空氣閥浮球升力隨排氣壓差變化Fig.11 Variation of floating ball lift of the air valve with exhaust pressure difference

5 結論

（1）空氣閥護筒高度為285，380，440 mm時，空氣閥排氣量均大于規(guī)范要求的排氣量，達到規(guī)范使用標準。

（2）在相同護筒高度情況下，排氣流量系數隨著排氣壓力增大而增大，在258 mm護筒高度，進出口壓差100 kPa時，達到最大排氣流量系數0.980；在相同排氣壓差情況下，空氣閥隨護筒高度的增加，排氣量逐漸減小，在440 mm護筒高度，進出口壓差為35 kPa時，最小排氣流量系數為0.868，因此護筒高度不宜過高。

（3）空氣閥在排氣過程中，浮球升力為負值，方向是向下的，并且排氣差壓越大浮球向下的力越大，空氣閥護筒高度為285 mm時，浮球向下的力分別為312.63，675.6，885.6 N，不會由于浮球的氣動力上升，導致堵塞排氣孔，并且隨著氣體排出，水位升高使浮球堵住排氣孔，避免水體流出，確?？諝忾y正常運行。