離心風(fēng)機集流器的響應(yīng)面優(yōu)化設(shè)計

趙燕杰,譚俊飛,黨飛龍,李景銀

(1.西安交通大學(xué)能源與動力工程學(xué)院,710049,西安;2.中冶賽迪工程技術(shù)股份有限公司,400013,重慶)

?

離心風(fēng)機集流器的響應(yīng)面優(yōu)化設(shè)計

趙燕杰1,譚俊飛2,黨飛龍1,李景銀1

(1.西安交通大學(xué)能源與動力工程學(xué)院,710049,西安;2.中冶賽迪工程技術(shù)股份有限公司,400013,重慶)

對某高效低噪音離心風(fēng)機的集流器進行了優(yōu)化設(shè)計研究,首先通過數(shù)值計算得到了在初始設(shè)計的集流器擴張角、集流器收縮角和擴張段長度等參數(shù)下的離心風(fēng)機氣動性能,并與實驗結(jié)果進行對比驗證,在此基礎(chǔ)上采用響應(yīng)面方法進行二次回歸擬合,得到了錐弧型集流器擴張角、收縮角和擴張段長度與風(fēng)機效率及聲功率級之間的函數(shù)關(guān)系并進行了優(yōu)化分析,由此得到使風(fēng)機性能最優(yōu)的集流器參數(shù)。分析結(jié)果表明,錐弧型集流器對離心風(fēng)機氣動性能的影響比較顯著,對聲功率級的影響不太顯著,錐弧型集流器的擴張角、收縮角對離心風(fēng)機氣動性能的影響明顯,而集流器擴張段長度對氣動性能的影響很小。優(yōu)化后的集流器使得離心風(fēng)機的效率提升了3.6%,靜壓提高了4.1%,聲功率級降低了1.7 dB。

離心風(fēng)機;集流器;響應(yīng)面法;聲功率級

集流器作為離心風(fēng)機三大部件之一,起著降低流動損失、將氣流平穩(wěn)均勻送入葉輪的作用[1]。氣流由集流器進入葉輪時會在葉輪進口附近形成一個渦流區(qū),從而影響到葉輪中前盤附近的氣體流動狀況,因此集流器結(jié)構(gòu)和安裝位置會對風(fēng)機性能有很大影響。

國內(nèi)外學(xué)者對于離心風(fēng)機內(nèi)部的流動狀況及機理進行了廣泛的研究[2-5]。Refeal等通過數(shù)值研究發(fā)現(xiàn),葉道內(nèi)部的三維非對稱流動與集流器結(jié)構(gòu)參數(shù)有著直接的關(guān)系[6]。Montazerin等通過改進集流器結(jié)構(gòu)使風(fēng)機內(nèi)部氣流流動得以改善[7-8]。溫選鋒等采用實驗方法研究了橢圓形集流器對多翼離心風(fēng)機氣動性能及噪音的影響[9]。徐文明等研究發(fā)現(xiàn),收斂型集流器進口傾斜角度對風(fēng)機的氣動性能會產(chǎn)生較大影響,而對噪音的影響則很小[10]。楊昕等以雙吸前向多翼風(fēng)機為例實驗研究了集流器的偏心安裝位置對風(fēng)機性能及噪聲的影響[11]。王嘉冰等采用數(shù)值模擬方法研究了集流器安裝位置對葉道內(nèi)氣流分布及風(fēng)機性能的影響[12]。

本文首先采用數(shù)值模擬方法研究了離心風(fēng)機中常用的錐弧型集流器結(jié)構(gòu)參數(shù)對其氣動性能及噪聲的影響,并與實驗結(jié)果進行了對比驗證。在此基礎(chǔ)上,結(jié)合響應(yīng)面方法擬合出錐弧型集流器擴張角、收縮角和擴張段長度這3個關(guān)鍵參數(shù)與風(fēng)機氣動性能和氣動噪聲的函數(shù)關(guān)系式,從而對集流器參數(shù)進行了優(yōu)化,以提高離心風(fēng)機的效率、降低噪聲。

1 物理模型及數(shù)值計算和驗證

1.1 離心風(fēng)機各部件幾何參數(shù)

離心風(fēng)機主要包括集流器、葉輪和蝸殼三部分,其中集流器的結(jié)構(gòu)如圖1所示。葉輪是由16個周向均布的后向葉片及輪盤、輪蓋組成的閉式葉輪。蝸殼型線主要由4段圓弧組成,各部件主要結(jié)構(gòu)參數(shù)如表1所示。風(fēng)機的設(shè)計流量qv=0.65 m3/s,設(shè)計轉(zhuǎn)速n=2 900 r/min。

圖1 集流器結(jié)構(gòu)及主要參數(shù)

參數(shù) 數(shù)值葉輪葉片出口直徑D2/mm481葉輪葉片進口直徑D1/mm220葉輪進口直徑D0/mm209葉輪進口寬度b1/mm63.64葉輪出口寬度b2/mm24.54葉輪葉片數(shù)Z16葉片進口安裝角β1A/(°)27.2葉片出口安裝角β2A/(°)82蝸殼寬度B/mm108蝸舌半徑R/mm28集流器總長度L/mm40擴張段長度l/mm10擴張角σ/(°)30收縮角?/(°)25

1.2 離心風(fēng)機內(nèi)部數(shù)值計算方法

本文采用商業(yè)軟件Ansys Fluent模擬了風(fēng)機內(nèi)部的流動,獲得了離心風(fēng)機的氣動性能與氣動噪聲,同時在建模過程中忽略了集流器與葉輪套接部分徑向間隙造成的內(nèi)泄漏流動。風(fēng)機葉輪網(wǎng)格在Turbo Grid中生成,集流器及蝸殼網(wǎng)格在ANSYS Icem中生成。為了提高近壁面區(qū)域流動計算的準(zhǔn)確性,對葉輪內(nèi)部、蝸殼蝸舌壁面附近的網(wǎng)格進行了加密控制和非等距處理,同時對計算網(wǎng)格進行了網(wǎng)格無關(guān)性驗證,由此得到集流器、葉輪和蝸殼的網(wǎng)格數(shù)分別約為16萬、96萬和87萬,整機網(wǎng)格數(shù)約為199萬,如圖2所示。

由于風(fēng)機內(nèi)部流動馬赫數(shù)小于0.3,屬于不可壓縮流動,所以對離心風(fēng)機的氣動性能采用定常計算。定常流動計算時采用了多參考系(MRF)方法,葉輪區(qū)域在相對坐標(biāo)系下求解,集流器和蝸殼區(qū)域在靜止坐標(biāo)系下求解,動靜交界面用來傳遞信息。風(fēng)機進、出口邊界條件分別為流量進口、靜壓出口,湍流模型為Realizablek-ε湍流模型,穩(wěn)態(tài)流場采用SIMPLEC算法計算。

發(fā)動機服務(wù)性資料模塊化數(shù)據(jù)配置管理的研究………………………… 閆玉，劉棟，李小燕，譚磊，傅連朋（1-273）

(a)葉輪網(wǎng)格

(b)集流器網(wǎng)格

(c)蝸殼網(wǎng)格圖2 離心風(fēng)機網(wǎng)格

采用上述的計算模型、網(wǎng)格和數(shù)值計算方法對原始設(shè)計的風(fēng)機模型進行了數(shù)值計算,由此得到設(shè)計流量下風(fēng)機的靜壓Pst、總壓Ptot和效率η,并與實驗進行了對比,結(jié)果如表2所示。由表中結(jié)果可以看出,在設(shè)計工況(流量為0.65 m3/s)下,η最大誤差約為2.1%,靜壓Pst最大誤差約為1.9%,表明了本文數(shù)值計算方法的準(zhǔn)確性。

表2 設(shè)計流量下實驗與計算結(jié)果比較

在模擬非定常流動時,采用滑移網(wǎng)格(SM)技術(shù)處理各部件之間(集流器-葉輪、葉輪-蝸殼交界面)的信息傳遞,將定常計算結(jié)果作為非定常的初值進行計算,收斂后激活Fluent聲學(xué)模塊,并利用FW-H方程預(yù)測離心風(fēng)機噪聲,由此獲取聲源信息。選擇PISO算法求解非定常流動,而非定常流動中時間步長計算式如下

(1)

式中:n為轉(zhuǎn)速;K為單流道計算步數(shù),本文中K=30;Z為葉片數(shù)。

使用FW-H方程預(yù)測噪聲時忽略了四極子源的影響,離心風(fēng)機中的所有墻(Wall)將作為聲源,對數(shù)值計算得到的關(guān)于時間的離散聲壓進行FFT變換,進而得到監(jiān)測點上的線性聲壓級。聲功率級是聲源在單位時間內(nèi)向外輻射出的總聲能,本文根據(jù)ISO3745—2012[13]、采用包絡(luò)面法在離心風(fēng)機的包絡(luò)面上布置了20個檢測點,通過測量檢測點上的聲壓級,再由以下公式合成聲功率級

(2)

(3)

式中:Lpi表示檢測點i的聲壓級;Lw表示聲功率級;S1表示包絡(luò)面的面積;S0表示參考面積。

2 集流器優(yōu)化結(jié)果與分析

2.1 響應(yīng)面方法

響應(yīng)面方法是用來對設(shè)計目標(biāo)(響應(yīng))受多個設(shè)計變量(因素)影響的問題進行建模和分析的方法,目的是優(yōu)化響應(yīng),包括實驗設(shè)計、方程回歸、回歸方程檢驗和最優(yōu)化求解等。

(4)

(5)

其中k=1,2,…,m-1(j≠k)。

通過實驗設(shè)計,如中心復(fù)合設(shè)計(CCD)或Box-Behnken設(shè)計等,可以獲得一系列樣本點,響應(yīng)面方程的建立就是對這些實驗點進行回歸的過程。對于求解模型中的各個待定系數(shù),通常采用最小二乘法進行。在確定回歸方程后還要進行顯著性檢驗,包括方程顯著性檢驗、偏回歸系數(shù)檢驗和失擬性檢驗,回歸方程顯著性、失擬不顯著性是回歸方程可以用于預(yù)測優(yōu)化的前提。首先通過實驗設(shè)計建立回歸方程,再經(jīng)過檢驗得到可以用于預(yù)測優(yōu)化的數(shù)學(xué)模型,然后根據(jù)極值存在的條件得到最佳值,也可以通過繪制響應(yīng)面或等高線圖來直觀顯示預(yù)測的最佳值。要保證響應(yīng)變量取得最佳值時的設(shè)計變量取值落在實驗范圍內(nèi),可能需要通過多次實驗來確定設(shè)計變量所在的范圍。

2.2 實驗設(shè)計及參數(shù)選擇

集流器的主要參數(shù)(如圖1所示)包括L、σ、φ、l。本文通過響應(yīng)面方法,研究集流器的σ、φ、l對風(fēng)機性能及氣動噪聲的影響。

L為40 mm、出口直徑Dout為207 mm、集流器特征參數(shù)γ=L/Dout=0.02,σ、φ、ε作為實驗因素,ε=l/L。初始選定的取值范圍為ε=0.2～0.35,σ=15°～40°,φ=20°～40°,響應(yīng)變量選擇效率與聲功率級。采用Box-Behnken方法進行實驗設(shè)計,Design-Expert軟件給出實驗設(shè)計方案,如表3所示。

表3 Box-Behnken實驗設(shè)計方案

各實驗方案計算結(jié)果如表4所示。分析表4中的數(shù)據(jù)可以發(fā)現(xiàn),在集流器各因素取值范圍內(nèi),集流器結(jié)構(gòu)參數(shù)對離心風(fēng)機的聲功率級影響很小,對氣動性能影響較大。因此,將聲功率級從響應(yīng)變量中去除,僅保留效率作為響應(yīng)變量,再通過回歸方程得到集流器結(jié)構(gòu)參數(shù)與效率之間的近似函數(shù)關(guān)系。

在Design-Expert軟件中對3個實驗變量的9種組合ε2、σ2、φ2、εσ、εφ、σφ、ε、σ、φ選擇不同的項進行方程回歸,再對擬合出來的方程進行方差分析,根據(jù)回歸方程的顯著性檢驗、偏回歸系數(shù)檢驗、失擬性檢驗選擇參與方程回歸的項,通過多次試驗便可得到較好的效率擬合方程

0.006118φ2-0.005836σ2)

(6)

分析發(fā)現(xiàn),回歸方程的顯著性明顯、失擬性不顯著,回歸方程在各個實驗點擬合得很好。方程(6)中沒有包含集流器的ε項,這是因為在加入該項后,通過偏回歸系數(shù)檢驗得到的相關(guān)項對離心風(fēng)機效率沒有顯著性影響。從回歸方程可以看出,在集流器參數(shù)的取值范圍內(nèi),σ、φ對離心風(fēng)機的效率有顯著性影響,ε則顯著性影響不明顯。

對滿足顯著性檢驗的方程進行預(yù)測分析,通過Design-Expert軟件繪制出該回歸方程的等高線圖,由此求出響應(yīng)變量在各因素取值范圍內(nèi)的最佳值,如圖3所示。

圖3 第一次優(yōu)化效率等高線圖

從圖3中可以看出,在等高線圖中沒有出現(xiàn)最優(yōu)點,說明選取的實驗因素范圍不合適,需要進行調(diào)整。進一步由圖中的趨勢可以判斷,為了保證出現(xiàn)最優(yōu)點,可以將σ調(diào)整為20°～50°、φ調(diào)整為0°～40°。由于ε的顯著性影響不明顯,因此第2次優(yōu)化中僅選取σ、φ作為實驗因素,同時采用CCD來安排實驗。CCD實驗設(shè)計各方案及數(shù)值計算結(jié)果如表5所示。

表5 CCD實驗點設(shè)計及計算結(jié)果

通過Design-Expert軟件對計算結(jié)果進行方程回歸分析,經(jīng)過不斷調(diào)整參與擬合的項,得到滿足顯著性檢驗的回歸方程

(7)

分析發(fā)現(xiàn),該回歸方程的顯著性明顯、失擬性不顯著,回歸方程在各個實驗點擬合得很好。根據(jù)回歸方程繪制等高線圖如圖4所示。

圖4 第2次優(yōu)化效率等高線圖

由圖4中可以看出,在等高線圖中存在最優(yōu)點,即響應(yīng)變量可以在實驗因素的范圍內(nèi)得到極大值。對回歸方程進行最優(yōu)化求解,即可得到最佳值。

2.3 結(jié)果驗證及分析

由回歸方程最優(yōu)化求解結(jié)果顯示,當(dāng)σ為42.8°、φ為23.8°時,回歸方程的響應(yīng)變量η最高。在集流器設(shè)計參數(shù)中,由于ε的影響不顯著,因此改進后的l仍然保持不變。對優(yōu)化后的風(fēng)機模型進行數(shù)值計算,得到設(shè)計流量(0.65 m3/s)下風(fēng)機的性能及聲功率級,并與模型優(yōu)化前的計算結(jié)果進行了對比,如表6所示。

表6 改進前后風(fēng)機設(shè)計工況下氣動性能及聲功率級

為了驗證優(yōu)化結(jié)果的可靠性,本文還對集流器優(yōu)化后的風(fēng)機性能進行了實驗測試,實驗遵循了GB/T1236—2000[14]、采用了進口側(cè)實驗風(fēng)筒裝置。改進前后風(fēng)機在設(shè)計工況(0.65 m3/s)下的性能實驗結(jié)果如表7所示。

表7 改進前后風(fēng)機設(shè)計工況下性能實驗結(jié)果對比

對比表6和表7可以發(fā)現(xiàn),集流器改進前后風(fēng)機性能的計算結(jié)果與實驗值都在允許的誤差范圍內(nèi),并且趨勢一致。由表7還可以看出,采用響應(yīng)面法對集流器結(jié)構(gòu)參數(shù)進行優(yōu)化后,離心風(fēng)機的靜壓升高了4.1%,總壓提升了3.9%,效率提高了3.6%。

集流器優(yōu)化前后風(fēng)機內(nèi)部流場的分布如圖5所示?？梢钥闯?集流器優(yōu)化后風(fēng)機主流區(qū)大部分流道的流動性能得到了改善,流動的非對稱性減弱,吸力面低速區(qū)減小,只有在靠近蝸舌附近的流道內(nèi)流動性能惡化。

(a)集流器優(yōu)化前 (b)集流器優(yōu)化后圖5 集流器優(yōu)化前后50%葉高處截面流線圖

3 結(jié) 論

本文采用數(shù)值模擬與響應(yīng)面相結(jié)合的方法,對某高效、低噪聲離心風(fēng)機集流器參數(shù)進行了優(yōu)化研究,通過對氣動性能及噪聲的分析,得到以下結(jié)論。

(1)將數(shù)值模擬方法與響應(yīng)面方法相結(jié)合,能夠應(yīng)用于離心風(fēng)機的改進和優(yōu)化。

(2)采用響應(yīng)面法對錐弧型集流器的優(yōu)化結(jié)果顯示,集流器對離心風(fēng)機的聲功率級的影響不明顯,對氣動性能的影響顯著;錐弧型集流器的擴張角、收縮角對離心風(fēng)機的氣動性能的影響明顯,擴張段長度的影響則很小。

(3)優(yōu)化后離心風(fēng)機的效率提高了3.6%,靜壓提高了4.1%,總壓提高了3.9%,聲功率級降低了1.7 dB。

[1] 李慶宜. 通風(fēng)機 [M]. 北京: 機械工業(yè)出版社, 1981: 61-64.

[2] KIND R J, TOBIN M G. Flow in a centrifugal fan of the squirrel-cage type [J]. ASME Journal of Turbomachinery, 1990, 112(1): 84-90.

[3] KIND R J. Prediction of flow behavior and performance of squirrel-cage centrifugal fans operating at medium and high flow rates [J]. ASME Journal of Fluids Engineering, 1997, 119(3): 639-646.

[4] MINER S M. Evaluation of blade passage analysis using coarse grids [J]. ASME Journal of Fluids Engineering, 2000, 122(2): 345-348.

[5] LIN S, HUANG C. An integrated experimental and numerical study of forward-curved centrifugal fan [J]. Experimental Thermal and Fluid Science, 2002, 26(5): 421-434.

[6] RAFAEL B, FRANCISCO G, SANDRA V. Numerical model for the unsteady flow features of a squirrel cage fan [C]∥Proceedings of the ASME Fluids Engineering Division Summer Meeting. New York, USA: ASME, 2009: 173-183.

[7] MONTAZERIN N, DAMANGIR A, MIRIAN S. A new concept for squirrel-cage fan inlet [J]. Proceedings of the Institution of Mechanical Engineers: Part A Journal of Power and Energy, 1998, 212(5): 343-349.

[8] MONTAZERIN N, DAMANGIR A, MIRIAN S. Inlet induced flow in squirrel-cage fans [J]. Proceedings of the Institution of Mechanical Engineers: Part A Journal of Power and Energy, 2000, 214(3): 243-253.

[9] 溫選鋒, 楊昕, 祁大同, 等. 橢圓形進口集流器對多翼離心風(fēng)機性能影響的實驗研究 [J]. 西安交通大學(xué)學(xué)報, 2011, 45(11): 45-51. WEN Xuanfeng, YANG Xi, QI Datong, et al. Experimental study on the effect of elliptical inlet nozzle on performance of multi-blade centrifugal fan [J]. Journal of Xi’an Jiaotong University, 2011, 45(11): 45-51.

[10]許文明, 鄭祖義, 萬忠民, 等. 空調(diào)柜機多翼離心風(fēng)機集流器的優(yōu)化研究 [J]. 制冷學(xué)報, 2012, 33(3): 58-61. XU Wenming, ZHENG Zuyi, WAN Zhongmin, et al. Optimization study of bell mouth on multi-blade centrifugal fan used in floor standing air-conditioner [J]. Journal of Refrigeration, 2012, 33(3): 58-61.

[11]楊昕, 溫選鋒, 袁民建, 等. 雙吸多翼離心風(fēng)機進口集流器偏心安裝的試驗研究 [J]. 流體機械, 2012, 40(2): 1-4. YANG Xi, WEN Xuanfeng, YUAN Minjian, et al. Experimental research on eccentric inlet of a multi-blade centrifugal fan with dual inlet [J]. Fluid Machinery, 2012, 40(2): 1-4.

[12]王嘉冰, 區(qū)穎達, 吳克啟. 空調(diào)風(fēng)機葉道內(nèi)旋渦流動分析及進氣口偏心的影響 [J]. 工程熱物理學(xué)報, 2005, 26(6): 53-55. WANG Jiabin, QU Yinda, WU Keqi. Analysis on the vortices flow in the blade passages of the fan for air-conditioner and the influence of the eccentric fan inlet [J]. Journal of Engineering Thermophysics, 2005, 26(6): 53-55.

[13]Technical Committee ISO/TC 43, Acoustics, Subcommittee SC 1, Noise. ISO3745-2012 Acoustics-determination of sound power levels and sound energy levels of noise sources using sound pressure: precision methods for anechoic rooms and hemi-anechoic rooms [S]. Geneva, Switzerland: ISO, 2012.

[14]國家技術(shù)監(jiān)督局. GB/T 1236-2000 工業(yè)通風(fēng)機用標(biāo)準(zhǔn)化風(fēng)道進行性能試驗 [S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2000.

(編輯 苗凌)

Optimization Design of Centrifugal Fan Inlet Collector by Response Surface Methodology

ZHAO Yanjie1,TAN Junfei2,DANG Feilong1,LI Jingyin1

(1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China;2. CISDI Engineering Co., Ltd., Chongqing 400013, China)

Optimization design method for the inlet collector of a centrifugal fan with high-efficiency and low-noise is investigated. The aerodynamic performances of the centrifugal fan are obtained by numerical simulation in the original parametric design, including the divergence angle and shrinking angle of the fan inlet collector, and the length of the divergent section. The numerical calculation is compared with the experimental result to validate the accuracy of the algorithm. To improve the aerodynamic and acoustic performances of the centrifugal fan, quadratic polynomial regression is considered to analyze the function relation of the centrifugal fan performance and the fan inlet collector parameters, and the optimal structure parameters of centrifugal fan inlet collector are obtained. The results indicate that the three fan inlet collector parameters exert remarkable effects on the aerodynamic performances, but slight effect on sound power level. The divergence angle and shrinking angle of the fan inlet collector obviously affect the aerodynamic performances, while divergent section length does not. After optimizing the centrifugal fan inlet collector, the efficiency and static pressure ratio increase by 3.6% and 4.1%, respectively, and the sound power level decreases by 1.7 dB.

centrifugal fan; inlet collector; response surface methodology; sound power level

2015-02-09。

趙燕杰(1986—),男,博士生;李景銀(通信作者),男,教授,博士生導(dǎo)師。

國家自然科學(xué)基金資助項目(51276137);中央高?；究蒲袠I(yè)務(wù)費專項資金資助項目。

時間:2015-08-13

10.7652/xjtuxb201511009

TH432

A

0253-987X(2015)11-0049-06

網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20150813.1014.006.html