燃料電池車用旋渦風機氣動噪聲試驗及仿真研究*

康 強,左曙光,韓惠君,韋開君

(同濟大學新能源汽車工程中心，上海 201804)

前言

目前，燃料電池汽車研發(fā)的重點是動力系統(tǒng)的集成和控制，以及動力系統(tǒng)關(guān)鍵部件的開發(fā)，對于其振動和噪聲的研究尚未引起足夠重視。實際上處于開發(fā)階段的燃料電池車的噪聲水平與商業(yè)化的汽油車和柴油車相比并不具有明顯優(yōu)勢，最主要的原因是燃料電池車上的空氣輔助系統(tǒng)(簡稱“空輔系統(tǒng)”)產(chǎn)生很大的噪聲[1]。

本文中研究的燃料電池轎車雖然消除了內(nèi)燃機這個噪聲源，但卻增加了電機、燃料電池系統(tǒng)和冷卻系統(tǒng)等噪聲源，使振動噪聲源發(fā)生了根本的改變?？蛰o系統(tǒng)為燃料電池提供壓縮空氣，其中的空氣壓縮機在工作時產(chǎn)生明顯噪聲，而空氣濾清器在大進氣量時，會產(chǎn)生嘯叫聲。此外，驅(qū)動電機-變速器總成、冷卻水泵、空調(diào)壓縮機和動力轉(zhuǎn)向泵等作為運動件也將對車內(nèi)噪聲產(chǎn)生影響。測試結(jié)果表明，空氣壓縮機產(chǎn)生的氣動噪聲是燃料電池車的主要噪聲源，其噪聲特性為高頻、尖銳，聲品質(zhì)差[1]。文獻[2]中對一輛重型載貨車的燃料電池空氣壓縮機系統(tǒng)的噪聲進行了試驗研究，并提出了一些被動噪聲控制方法。文獻[3]中介紹了對NuCellSys新一代空輔系統(tǒng)(采用渦旋式壓縮機)采取的降噪措施。國內(nèi)對燃料電池空輔系統(tǒng)用空氣壓縮機的研究仍處于初級階段。大部分研制的壓縮機在能夠滿足流量、壓力和效率的情況下并不能滿足振動和噪聲的要求，因此低噪聲和小型化技術(shù)是空氣壓縮機研究的一個關(guān)鍵方向。

本文中研究的燃料電池汽車采用的空氣壓縮機為旋渦風機，其高頻噪聲對車內(nèi)外噪聲貢獻很大。目前國內(nèi)外對旋渦風機的研究主要應(yīng)用于工農(nóng)業(yè)生產(chǎn)中，對于其具體的噪聲機理及其相應(yīng)的降噪方法的研究還不多見[4]。通常在加裝了消聲器后，旋渦風機的噪聲值為70～80dB(A)，對于壓縮機或鼓風機來說是較低的。但對燃料電池轎車來說，由于取消了內(nèi)燃機等傳統(tǒng)汽車上的噪聲源，旋渦風機的噪聲更顯突出。因此，對燃料電池轎車旋渦風機的噪聲進行深入的研究具有重要的工程意義。

1 旋渦風機振動噪聲臺架測試

為進一步研究旋渦風機自身的振動和噪聲特性，以便從源頭上進行控制，設(shè)計加工了旋渦風機振動噪聲測試臺架，可分別測試進出口管道和進出口消聲器對旋渦風機噪聲的影響。測試臺架主要由旋渦風機、驅(qū)動電機、變頻器、消聲器、蝶閥、進出口管路、渦街流量傳感器、壓力傳感器、鋼絲軟管和支架等組成。為減少電機噪聲對風機噪聲的影響，將電機用吸音材料包裹起來。

根據(jù)GB T2888—1991《風機和羅茨鼓風機噪聲測量方法》和ISO 3744，在旋渦風機周圍布置聲學麥克風測點位置，如圖1所示。

試驗測試了多種工況，本文中只研究風機進出口不帶消聲器并用管道將進出口氣流引向室外時風機所產(chǎn)生的噪聲。根據(jù)旋渦風機的轉(zhuǎn)速范圍測試了風機在穩(wěn)定轉(zhuǎn)速(1 000～3 800r/min、間隔400r/min)和加速工況下的振動與噪聲信號。

將半球面上5個測點的聲壓級取對數(shù)平均，得到旋渦風機平均聲壓級隨轉(zhuǎn)速的變化曲線，如圖2所示。它在1 000r/min時的聲壓級僅為55dB(A)，在3 800r/min時達到78dB(A)，跟轉(zhuǎn)速基本成線性關(guān)系。

加速工況下C2點的噪聲云圖見圖3。

由圖3可以看出，旋渦風機噪聲的主要峰值頻率為風機轉(zhuǎn)頻的55倍及其諧頻，此外，在相對較低的500～1 300Hz范圍內(nèi)也存在較高的峰值。其中55倍為葉片通過頻率，而500～1 300Hz為風機殼體和葉輪的前幾階模態(tài)頻率范圍，故在此頻率范圍內(nèi)輻射的噪聲較高，尤其是在高轉(zhuǎn)速時。為研究其噪聲產(chǎn)生的原因和規(guī)律，以數(shù)值仿真的方法計算旋渦風機內(nèi)流場特性及其遠場輻射噪聲，為以后的優(yōu)化設(shè)計提供參考。

2 旋渦風機內(nèi)流場的計算與分析

2.1 結(jié)構(gòu)參數(shù)

所研究的旋渦風機結(jié)構(gòu)如圖4和圖5所示，其相應(yīng)的結(jié)構(gòu)參數(shù)見表1。

2.2 網(wǎng)格劃分

為更準確地設(shè)定進出口邊界條件和加快計算的收斂，將進出口段分別沿各自方向延長一段。計算域由4部分組成：進口流域、葉片旋轉(zhuǎn)區(qū)域、蝸殼區(qū)域和出口流域。為適應(yīng)旋渦風機內(nèi)部結(jié)構(gòu)的復雜性，采用多塊網(wǎng)格生成方法構(gòu)建高質(zhì)量非結(jié)構(gòu)網(wǎng)格?？紤]到不同流動區(qū)域的重要性差異，對旋轉(zhuǎn)葉輪內(nèi)部、拐角和間隙處壁面附近的網(wǎng)格進行了加密，最終生成的網(wǎng)格有550萬個單元,如圖6所示。

表1 旋渦風機流道主要結(jié)構(gòu)參數(shù)

2.3 CFD模型設(shè)置

CFD計算包括三維黏性可壓縮流的定常和非定常計算。定常計算的湍流模型采用RNGk-ε模型，近壁面采用標準壁面函數(shù)，壓力-速度的耦合采用SIMPLE方法。將定常計算的結(jié)果作為非定常計算的初值。非定常計算采用離散渦模擬(DES)，其中小尺度渦的模擬采用基于混合長度理論的亞格子Smagorinsky模型。非定常計算的壓力-速度耦合采用PISO方法。壓力項、動量項、湍動能和湍流耗散率項的離散都采用2階迎風格式，時間項離散采用2階隱式格式。旋轉(zhuǎn)葉片與殼體之間的動靜交界面置于葉頂間隙的中間位置。對于定常計算，采用多參考坐標系MRF模型處理動靜交界面；對于非定常計算，采用滑移網(wǎng)格技術(shù)。CFD計算邊界條件如下。

進口邊界條件：水力直徑為53mm，給定進口壓力，取相對壓力為0。

出口邊界條件：水力直徑為53mm，給定出口壓力，取相對壓力為20kPa。

本文中共設(shè)定了風機的4個轉(zhuǎn)速進行計算：1 400、2 200、3 000和3 800r/min。每個工況設(shè)定葉輪轉(zhuǎn)動一圈為800個步長。相應(yīng)的計算參數(shù)呈現(xiàn)周期性變化后即認為收斂。

此外，在流場內(nèi)部設(shè)置一些監(jiān)測點和監(jiān)測面，可以反映葉輪旋轉(zhuǎn)過程中流場內(nèi)部的壓力波動情況。從入口到出口設(shè)定5個監(jiān)測點，以監(jiān)測葉片和隔舌上的壓力隨時間的變化歷程，如圖7所示。

2.4 計算結(jié)果

圖8為3 000r/min工況下旋渦風機內(nèi)流場壓力云圖?？梢钥闯?，其壓力從入口到出口均勻增大。渦聲理論把輻射的噪聲大小與渦量聯(lián)系起來，主要關(guān)心流場中渦量的大小、變化和運動情況，對流場中細節(jié)問題和復雜流動現(xiàn)象進行了簡化。圖9為瞬時渦量云圖，由圖可見，流道中前半部分渦量大于后半部分的渦量，表示流道內(nèi)的前半部分是主要的氣動噪聲源。

圖10為旋渦風機葉輪表面聲功率級的分布，由圖可見，噪聲最強的位置也是位于從入口至大約流道中部的葉輪葉片上，與渦量分析的結(jié)果一致。

圖11分別為不同轉(zhuǎn)速工況下得到的監(jiān)測點p2、葉片和出口隔舌的壓力脈動信號頻譜圖。從圖中可見：在靜止的p2點和隔舌上，壓力脈動均反映出明顯的55倍頻，并且隨著轉(zhuǎn)速的升高其峰值逐漸增大；而旋轉(zhuǎn)的葉片上的壓力脈動則反映出42倍頻，正好等于進口到出口流道無隔板區(qū)域內(nèi)的葉片數(shù)量，如圖4和圖7所示。這是由于旋轉(zhuǎn)葉片在經(jīng)過隔板時，由于隔板有一個傾斜角度，使經(jīng)過隔板的葉片周圍流場逐漸減弱，直到基本消失，故葉片上的壓力脈動表現(xiàn)為42倍，且頻譜峰值模糊，不如55倍頻清晰；此外，在3 000r/min以上時，葉片上壓力脈動頻率反映出500～1 000Hz的成分，對應(yīng)于旋渦風機在高轉(zhuǎn)速下試驗測試的噪聲信號頻譜?？梢?，旋渦風機噪聲主要由流場內(nèi)部壓力波動引起，其頻率具有一一對應(yīng)的關(guān)系。

3 旋渦風機氣動噪聲的計算與分析

目前工程上計算氣動噪聲學的方法一般采用Lighthill聲類比法。對Lighthill方程的求解大多數(shù)采用積分解法[5]，如Curle方程和FW-H方程。其中對于外流場問題可采用自由空間格林函數(shù)求解；然而對于幾何邊界復雜的內(nèi)流場問題，格林函數(shù)難以得到，只有對部分很規(guī)則的空間，如自由空間、封閉球空間和無限大平板半空間等，才能得到其解析表達式，故該方法的應(yīng)用受到很大局限，求解風機遠場聲輻射問題時不得不將問題簡化為自由空間問題，將其作為緊湊聲源，忽略蝸殼和進出口管道等對聲傳播的影響[6-7]。因此，聲傳播在復雜固壁內(nèi)的反射與折射等問題得不到體現(xiàn)，也不能包含結(jié)構(gòu)的彈性效應(yīng)，計算結(jié)果存在較大誤差，限制了其適用性。另一種求解方法是將Lighthill方程轉(zhuǎn)化為變分形式，采用有限元法求解，這樣可以將聲源進行空間離散以考慮結(jié)構(gòu)和聲的相互作用[8]。這種方法可適用于任何復雜的幾何結(jié)構(gòu)，有利于工程應(yīng)用，本文中主要基于此方法計算旋渦風機氣動噪聲。

Lighthill方程的形式為

(1)

式中：Tij=ρvivj+(p-c02ρ)δij-τij，稱為Lighthill應(yīng)力張量；ρ為流體密度；t為時間；c0為未受擾動流體的聲速；p為流體中的壓力；δij為克羅內(nèi)克符號；τij為黏性應(yīng)力張量；xi和xj為位移；vi和vj為流體速度。

3.1 Lighthill方程變分解法

Lighthill聲類比公式也可通過變分原理進行改寫，然后通過有限元法進行求解。Lighthill聲類比理論的強變分形式[9]為

(2)

式中：δρ是一個測試函數(shù)；Ω是計算域。運用格林定理對空間導數(shù)積分，獲得弱變分形式：

(3)

Γ=?Ω，若定義總應(yīng)力張量

Σij=ρvivj+(p-p0)δij-τij

(4)

則式(3)變?yōu)?/p>

(5)

這就是Lighthill方程的變分形式。在方程右端有兩個氣動聲源項：一個是體積源項，另一個是面源項。體積源導入到包含能量的渦結(jié)構(gòu)區(qū)域中，面源導入到計算域的控制界面上。在計算旋轉(zhuǎn)機械時，為了更精確地描述源項，選擇靠近轉(zhuǎn)子的面。

3.2 氣動聲學有限元計算模型

建立的旋渦風機聲學模型見圖12。進口和出口用管道模態(tài)模擬無限長管道，即不考慮旋渦風機進出口處的射流噪聲。將非定常流場計算結(jié)果插值到圖12所示的體聲源網(wǎng)格和面聲源網(wǎng)格，其中面聲源網(wǎng)格用于計算葉輪旋轉(zhuǎn)區(qū)域的聲源。體聲源網(wǎng)格外部建立旋渦風機殼體的實體網(wǎng)格，在殼體以外的區(qū)域建立一個包圍旋渦風機的聲傳播區(qū)域，最外層用無限元的面網(wǎng)格模擬自由空間的無反射邊界條件。

3.3 噪聲計算結(jié)果

3.3.1 聲壓級對比

在風機外部1m球面處分別沿x-y，y-z和x-z平面布置一圈場點，可計算出各點處的聲壓級頻譜圖，將所有點的聲壓級頻譜曲線取平均，可得旋渦風機平均聲壓級頻譜曲線。本文中僅計算旋渦風機轉(zhuǎn)速在3 000r/min下產(chǎn)生的噪聲，結(jié)果如圖13所示。

圖中將變分解法、積分解法計算得到的平均聲壓級與試驗測試得到的平均聲壓級進行了對比。可以看出，用變分方法時考慮到了殼體對聲波的作用，計算結(jié)果與試驗值較吻合，平均聲壓級僅相差2.7dB(A)。而用格林函數(shù)的積分解法進行計算時，平均聲壓級比試驗值高32.9dB(A)，可見對于內(nèi)流場的問題所做的簡化將導致結(jié)果產(chǎn)生很大的誤差。此外，從頻譜圖上看，數(shù)值計算結(jié)果均可反映出其最主要的峰值頻率2 750Hz(55倍頻)。在中低頻時，變分解法得到的噪聲頻譜與試驗結(jié)果比較吻合，但在高頻時頻譜相差很大，這主要受限于數(shù)值算法本身的精度，它的色散和耗散性使其在高頻時的計算誤差很大[10]。

3.3.2 聲場分布

圖14為有限元法計算的2 750Hz下旋渦風機流場內(nèi)、外部空間分布的聲壓級云圖。可以看出，由于殼體的存在，其聲壓級從內(nèi)部輻射到外部時顯著減小。

3.3.3 噪聲指向性

綜合各場點的聲壓級可得到旋渦風機的噪聲指向性圖，如圖15所示。

定義軸向為y方向，垂向為z方向，水平方向為x方向?？梢钥闯觯趚-y和y-z平面，不同方向聲壓級最大相差10dB(A)左右，x-z平面，即旋渦風機徑向平面，不同方向聲壓級最大相差4dB(A)左右。從圖15(a)和圖15(b)可見，沿旋渦風機軸向方向的聲壓級最低，為61dB(A)左右；從圖15(b)和圖15(c)可見，進口側(cè)聲壓級比出口側(cè)聲壓級大3dB(A)左右。據(jù)此，在進行整車布置時，應(yīng)使旋渦風機軸線方向或者出口端方向面向駕駛員位置，以減小駕駛員位置處的噪聲。此外，也可根據(jù)上述規(guī)律對旋渦風機采取相應(yīng)的降噪措施。

4 結(jié)論

(1) 試驗測試的燃料電池車用旋渦風機噪聲特性表明，其主要峰值頻率為55倍頻和在500～1 300Hz范圍內(nèi)，其平均聲壓級與轉(zhuǎn)速成線性關(guān)系。

(2) 利用離散渦模擬方法計算了旋渦風機內(nèi)部非定常流動，從渦量和表面聲功率級的對比說明，流道內(nèi)前半部分是氣動噪聲的主要聲源；而由壓力波動分析可知，流場內(nèi)壓力脈動頻譜基本上反映遠場噪聲的頻譜。

(3) 用有限元法計算遠場氣動噪聲，克服了積分方法不能求解復雜邊界格林函數(shù)的缺點，雖然建模較復雜，計算時間較長，但計算結(jié)果與實際更加吻合，非常適用于內(nèi)流場問題。這為進行噪聲優(yōu)化設(shè)計提供了重要的參考。

(4) 根據(jù)旋渦風機噪聲指向性圖，在進行整車布置時，應(yīng)使旋渦風機軸線方向或出口端方向面向駕駛員位置，以減小駕駛員位置處的噪聲。

[1] 申秀敏,左曙光,蔡建江,等.燃料電池車聲振測試及噪聲源識別[J].聲學技術(shù),2008(4):570-576.

[2] Palan V, Shepard Jr. W S, Lim T C. Case History: Noise Control Approaches for an Air-compressor in a Fuel-cell Auxiliary Power Unit[J]. Noise Control Engineering Journal,2004,52(5):197-209.

[3] Sang J, Venturi M, Bocksch R. Nvh-challenges of Air Supply Subsystems for Automotive Fuel Cell Applications[J]. SAE International Journal of Engines,2009,1(1):258-266.

[4] Kim Y H, Kang S H. Performance Evaluation of a Regenerative Blower for Hydrogen Recirculation Application in Fuel Cell Vehicles[G]. Miami, FL, United States:2006.

[5] Fedala D, Kouidri S, Rey R. Numerical Study of Time Domain Analogy Applied to Noise Prediction from Rotating Blades[J]. Journal of Sound and Vibration,2009,321(3-5):662-679.

[6] Liu Q H, Qi D, Tang H. Computation of Aerodynamic Noise of Centrifugal Fan Using Large Eddy Simulation Approach, Acoustic Analogy, and Vortex Sound Theory[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science,2007,221(11):1321-1332.

[7] 毛義軍,祁大同.葉輪機械氣動噪聲的研究進展[J].力學進展,2009(2):189-202.

[8] Kaltenbacher M, Escobar M, Becker S, et al. Numerical Simulation of Flow-induced Noise Using LES/SAS and Lighthill’s Acoustic Analogy[J]. International Journal for Numerical Methods in Fluids,2010,63(9):1103-1122.

[9] Oberai A A, Roknaldin F, Hughes T J R. Computational Procedures for Determining Structural-acoustic Response Due to Hydrodynamic Sources[J]. Computer Methods in Applied Mechanics and Engineering,2000,190(3-4):345-361.

[10] Astley R J, Sugimoto R, Mustafi P. Computational Aero-acoustics for Fan Duct Propagation and Radiation. Current Status and Application to Turbofan Liner Optimisation[J]. Journal of Sound and Vibration,2011,330(16):3832-3845.