傾斜式雙導葉擴散器氣動特性的數值模擬

（河南理工大學 機械與動力工程學院，焦作 454000）

0 引言

擴散器是指礦井主扇風機出口處外接的具有一定長度、斷面逐漸擴大的附屬裝置。它通過流道截面的逐漸擴大，減小出口處的氣流速度（出口動壓），將部分動壓轉變為通風機靜壓，回收風機出口的動能損失，減小能量損失，從而提高礦井通風機的運行效率。工程實踐表明，主扇通風機擴散器能夠回收風機出口動能的30%～80%[1]。

擴散器的通風阻力是影響其能量回收效率的主要因素，阻力越小，回收效率越高；而擴散器的通風阻力與其幾何結構等相關。為減小擴散器的通風阻力、提高其回收效率，國內外學者已經利用數值計算和實驗相結合的方法從擴散角、彎道形狀、斷面擴大系數等方面做了許多工作[2～7]。潘地林等利用數值模擬的方法研究了平面擴散器和錐形擴散器在不同工況下的擴散性能，給出了兩種結構型式擴散器的適用場合及結構參數的選擇原則[2]；修東亮利用fluent軟件模擬了軸流式通風機所用兩種結構型式擴散器的內部流動，分析了它們的運行效率[3]；戴巨川等采用粒子群優化算法對礦用對旋式軸流通風機擴散器進行了結構參數的多目標優化設計[4]；陳世強提出了雙切流線型擴散器，并數值計算方法和實驗方法對該型擴散器進行了結構參數的優化[5]。

本文采用計算流體力學方法對某型含雙導葉的60°傾斜式擴散器的內流場進行數值模擬研究，分析其流場結構和阻力特性，以期為擴散器幾何結構優化設計提供參考。

1 計算模型

圖1為某60°傾斜式擴散器數值計算模型的平面幾何示意圖。該擴散器模型高510mm，寬198mm，入口斷面高214mm，出口斷面高353mm，斷面擴大系數為1.65；導葉固定于擴散器彎道段，沿徑向等距分布，導葉沿擴散器中心線方向的長度為121mm，高為2.5mm，寬度與擴散器寬度相同。

圖1 擴散器計算模型示意圖

擴散器內流場利用ICEM CFD軟件生成六面體結構化網格，它具有占用內存少且易于計算邊界層的優點，網格總數約為86萬。

2 數值模擬

2.1 基本控制方程

假定流體為不可壓縮連續性介質，瞬時流動控制方程主要包括以下連續性方程和雷諾平均N-S方程：

式中：ui、uj、p均為湍流時均值，為雷諾應力張量；ρ、μ分別為流體的密度和流動動力粘性系數。

2.2 SSTk?ω湍流模型

擴散器內流場數值模擬計算湍流模型為剪切應力輸運k?ω模型（SSTk?ω），其湍流動能k方程以及湍流耗散率ω方程如下：

式中，vt可以用和γ等常數表示，常數的取值可參閱文獻[8]。

2.3 邊界條件及求解參數

數值計算使用有限體積法離散控制方程，離散格式采用高精度（High Resolution）格式，該方法使離散方程具有二階精度的截段誤差。求解參數的收斂標準為：最大殘差小于1×10-3；平均殘差小于1×10-4。為提高精度，計算中采用二階迎風格式。

數值模擬采用的流體介質為20℃的空氣，其密度為1.2255kg/m3，動力粘性系數為0.00001789Pa.s。根據工程實際情況，具體邊界條件設置如下：1）入口邊界條件：入口為速度入口，來流垂直于入口斷面，均勻分布，流速為v，壓強為101325Pa；2）出口邊界條件：出口邊界為壓力出口，相對靜壓為0Pa；3）壁面邊界條件：擴散器的壁面及導葉均設定為靜止無滑移壁面，壁面粗糙度為0，垂直于壁面的壓力梯度為0。

3 計算結果與分析

為驗證本文所用方法的精度和可靠性，采用本文的方法對某60°傾斜式擴散器[9]在給定來流速度13.56m/s的條件下的受力特性進行數值模擬分析，入口斷面橫截面面積為特征值，計算得到的擴散器阻力系數為1.165，此結果與其實驗值1.17的相對誤差為0.427%，相差很小，表明本文所用的數值模擬計算方法是可行的。

3.1 流場特性分析

擴散器內流場縱向中間截面上的詳細情況如圖2所示。

圖2 擴散器中心截面處的速度矢量分布圖

由圖2可知：在擴散器內流場的流動下游、擴散器彎道與出口之間、上壁面附近出現較大面積的漩渦區域，即吸風區；導葉的尾部附近也存在漩渦區，但是面積很??；出口斷面處，氣流的速度與出口截面并不垂直，而是存在一定的角度。出口處吸風區的存在將使得擴散器流場內的壓差阻力增大，并較大程度地影響擴散器出口的動能回收效率；導葉尾端漩渦區由于面積很小，對擴散器壓差阻力影響較小。

圖3所示為擴散器中心截面處的速度分布云圖。

圖3 擴散器中心截面處的速度分布云圖

由圖3可以看出，擴散器內流域的不同區域的速度分布不均勻。在擴散器內，由于流動方向的改變，在流體質點的離心慣性力作用下，使得擴散器彎道區域的流速值沿徑向由上壁面到下壁面逐漸減小，但下壁面附近沒有大面積的低速區存在，這主要是因為導葉的擾流作用，減小了速度的梯度；在擴散器的彎道部分上壁面和導葉1之間形成高速氣流分布區域；在內流場的流動下游、擴散器彎道與出口之間、上壁面附近出現較大面積的低速區域，主要是因為漩渦的存在使得流速值急劇減小，甚至出現逆流。

為清晰地表示出擴散器的內流場，給出擴散器內流場計算域流線分布圖如圖4所示。圖中顯示擴散器彎道與出口之間、上壁面附近出現較大面積漩渦區域，受壁面的影響，沿導葉展長方向漩渦的結構不同，因此利用二維來模擬三維流動時并不能真實地描述出擴散器內流場的具體情況。

圖4 擴散器內流場計算域流線圖

3.2 擴散器氣動特性分析

擴散器的局部通風阻力的大小可以用擴散器阻力系數來衡量，其計算公式如下：

式中：ξ表示擴散器的局部阻力系數；F表示擴散器的局部結構通風阻力，N；ρ表示流體介質密度，kg/m3；A表示特征面積，本文取為入口斷面面積，m2；V表示流速，m/s。

圖5為不同擴散器模型的阻力系數與來流速度之間的關系圖，其中，模型a表示無導葉擴散器模型、模型b表示單導葉擴散器模型、模型c表示雙導葉擴散器模型。

圖5 不同模型的阻力系數對比

從圖5可以看出：三種擴散器數值計算模型的阻力系數均隨流速的增大先逐漸減小，而后又增大；相對于無導葉型擴散器模型，在彎道處設置導葉可以顯著地減小擴散器阻力系數，改善擴散器的氣動性能；單導葉擴散器和雙導葉擴散器的局部阻力系數接近，大小與來流速度有關；流速較低時雙導葉擴散器的阻力系數略小一些，具有更好的氣動性能；流速較大時單導葉擴散器模型的阻力系數略小。

3.3 擴散效率分析

壓力恢復系數和擴散效率是衡量擴散器性能優劣的兩個重要參數，可以反映擴散器回收動壓能力的大小。

壓力恢復系數與進、出口斷面的面積比和擴散器局部阻力系數有關，計算公式如下：

式中，A1和A2分別為擴散器進口、出口橫截面積。

擴散器的實際壓力恢復系數與理論壓力恢復系數的比值稱為擴散效率，即：

式中，cpi表示理想狀態下，擴散器局部阻力系數為零時的理論壓力恢復系數。

根據式(6)和式(7)計算出雙導葉擴散器的氣動性能參數如表1所示。

表1 雙導葉擴散器的氣動性能參數

從表1可以看出，擴散器的氣動性能除了與其結構形式有關外，還受來流速度的影響。因此，對于某固定形式的擴散器來說，它們都有一個最優的來流速度使得擴散器的動能回收效率最高。

4 結論

1）數值模擬結果和試驗結果表明，利用CFD技術進行擴散器內流場數值模擬的計算方法是可行的，在新構型擴散器的開發、結構參數的最優化等工程應用中具有較大的指導作用；

2）傾斜式擴散器出口處速度矢量與擴散器出口斷面并不垂直，內流場中速度分布不均勻，在擴散器彎道之后、出口內側有渦流區域，使得出口處發生吸風現象，導致擴散器氣動特性降低；

3）在擴散器彎道處加設導葉可以改善內流場的速度分布均勻性，減小擴散器的局部阻力系數，提高擴散器出口動能損失的回收效率。

[1]王海橋,趙云勝,陳世強.礦井主要通風機擴散器結構形式及其性能評價[J].礦業安全與環保,2008,35(5):25-27.

[2]潘地林,程凱,楊春魚.礦用通風機直壁式擴散器參數計算分析[J].煤礦安全,2014,45(1):84-87.

[3]修東亮.軸流式通風機擴散器數值模擬與研究[J].煤礦機械,2013,34(3):57-58.

[4]戴巨川,劉德順,金永平,等.礦用對旋式軸流通風機擴散器結構多目標優化設計[J].中國機械工程,2013,24(3):386-391.

[5]陳世強.礦井通風系統擴散器性能與結構優化研究[D].湖南大學,2015.

[6]Effect of geometrical parameters on the performance of wide angle diffusers[J].International Journal of Innovative Research in Science, Engineering and Technology,2013,2(9):4178-4193.

[7]Majumder S, Roy D, Bhattacharjee S, et al. Experimental investigation of the turbulent fluid flow through a rectangular diffuser using two equal baffles[J]. Journal of The Institution of Engineers (India): Series C,2014,2(2):1-5.

[8]Menter F R. Zonal two equationk-ω turbulence mod-els for aerodynamic flows[J].AIAA Paper No.93-2906,1993.

[9]蔣軍成,張惠忱.礦井軸流式主扇擴散塔的實驗研究[J].煤炭科學技術,1995,23(8):10-13.