王華偉,王 礫,劉漢陽,王 博,胡 溧
(1.武漢科技大學 汽車與交通工程學院,湖北 武漢 430070;2.中國船舶集團有限公司 第七二二研究所,湖北 武漢 430223)
隨著發動機性能的不斷提高與車用電器功率需求的日益增加,對發動機冷卻系統散熱性能的要求也越來越高。采用尺寸更大、轉速更高的冷卻風扇無疑是最直接的方法,但是其噪聲、成本也會相應增加。
隨著對汽車噪聲、振動與聲振粗糙度(noise、vibration、harshness,NVH)要求的不斷提高,冷卻風扇噪聲問題逐漸凸顯。因此,研究冷卻風扇的氣動及噪聲特性,探索風扇的降噪設計方法,對于性能優異的車用冷卻風扇設計具有重要意義。
國內外學者在軸流風機的噪聲預測與控制領域已經進行了大量的研究,并且做出了大量的貢獻。
LUO Bo等人[1]對低速軸流風機葉尖泄漏流和氣動聲學進行了分析,并基于試驗數據的對比,成功驗證了數值預測的可靠性和準確性,結果表明,當葉頂間隙增大時,葉尖渦強度的增強以及與葉片表面相互作用的增強是產生額外寬頻噪聲的主要原因。TOKUSHIGE T等人[2-3]對環形風機氣動噪聲進行了研究,結果表明,通過對開口風扇增加導流環,不僅可以改善氣動性能,還可以降低風扇噪聲,且葉片表面的離散頻率噪聲是由與葉片通過頻率同步的周期性壓力波動產生的。YE Xue-min等人[4]針對不同鋸齒后緣對單級變槳距軸流風扇氣動及噪聲性能的影響進行了數值模擬研究,并設計了一種具有新型鋸齒葉片的低噪聲軸流風扇。周敏等人[5]提出了一種在軸流式壓氣機葉片開孔的方法,使壓力面到吸力面的壓差產生射流,從而吹除葉片吸力面的附面層,增加附面層的湍流動能,并且延緩氣流分離。胡勇軍[6]設計了一種穿孔式雙轉子軸流風扇,并探究了不同穿孔參數對其性能的影響,其發現通過對葉片穿孔可有效降低雙轉子軸流風扇氣動噪聲,提升NVH性能。YANG Xing-lin等人[7]分析了葉片穿孔對大型船用軸流風機的性能影響,研究結果表明,軸流風機在葉尖區域進行開孔對離散噪聲具有抑制作用;但是該研究缺少對風扇性能曲線的分析。WANG Chen等人[8]基于小型反向旋轉風扇中的被動降噪技術,采用多孔后緣處理了由于兩個轉子之間相互作用產生的噪聲,并綜合考慮氣動性能和噪聲性能,發現應選擇盡可能小的孔徑和在遠離葉片后緣方向上遞減的穿孔率分布;但其未能精準仿真離散頻率噪聲。王善彬等人[9]針對后緣穿孔結構對葉片氣動特性和噪聲特性的影響,以NACA65019葉片為基礎,采用大渦模擬和FW-H方程,研究了孔型及穿孔傾角對葉片氣動特性及噪聲特性的影響機理。
上述研究結果表明,葉片穿孔結構已在部分旋轉機械或翼型結構改進中取得進展。但環形冷卻風扇導流環與葉尖為一體式結構,離散噪聲形成機理不同,穿孔降噪結構尚未應用在環形冷卻風扇葉片之上。
基于以上研究,筆者采用數值模擬的方法,探究葉片穿孔數對車用環形冷卻風扇氣動性能和噪聲性能的影響。
首先,建立準確的冷卻風扇計算模型;然后,設計具有不同穿孔數的葉片結構;最后,對原風扇模型和葉片穿孔模型的流場與聲場進行分析,得出其影響規律,以期為環形冷卻風扇噪聲控制技術提供參考。
FW-H聲學模型本質上是一個非均勻波動方程,可以從連續性方程和Navier-Stokes方程推導得出。FW-H聲學模型主要用于噪聲預測。
FW-H方程可寫成:
(1)
式中:u為流體速度;un為f=0的表面法向流體速度分量;v為表面速度;vn為表面速度法向分量;ρ為流體密度;δ(f)為狄拉克函數;H(f)為亥維塞函數;p為聲壓;p′為遠場聲壓,p′=p-p0;a0為遠場聲速;Tij為Lighthill應力張量;n為物體表面的外法矢量;下標i,j,k為xi、xj、xk方向的分量。
Tij的定義為:
(2)
式中:Pij為壓縮應力張量,自由流物理量由下標0表示;δij為克羅內克符號。
對于斯托克斯流體,其表達式如下:
(3)
在聲源和接收器之間沒有障礙物以及自由空間流的假設下,可以對式(1)進行綜合分析。
完整的解包括表面積分和體積積分。表面積分表示單極子和偶極子聲源,以及部分來自四極子聲源的貢獻,而體積積分表示聲源表面外部區域中的四極聲源。當聲源表面包圍源區域,且流動為低亞音速時,體積積分的貢獻變小。因此,式中忽略了體積積分。
遠場聲壓可表示為:
(4)
式中:t為到接收器時間。
厚度噪聲方程可表示為:
(5)
載荷噪聲方程可表示為:
(6)
葉片表面速度可表示為:
(7)
葉片表面載荷可表示為:
(8)
筆者以某6葉發動機環形冷卻風扇作為研究對象(風扇外徑為405 mm、輪轂直徑為165 mm、輪轂比為0.4、葉片厚度為2 mm),使用FLUENT軟件對原風扇模型進行數值模擬。
在保證計算結果準確的前提下,筆者首先對風架凸耳、加強筋和插接件等結構進行簡化,并且將風扇輪轂簡化為同尺寸的圓臺;然后,嚴格按照試驗條件建立計算域模型(計算域包括模擬風扇運動狀態的旋轉域和模擬風扇進出口通道的流體域),建立進口為邊長1 m的正方體,出口則采用長為2 m、寬和高為1 m的長方體作為流體域;最后,建立直徑比風扇外徑略大的圓柱作為旋轉域將風扇包絡,以模擬其近壁面氣流狀態。
風扇簡化模型和計算域如圖1所示。
圖1 數值模擬計算模型
由于冷卻風扇存在復雜曲面,所以筆者采用網格劃分的方式選擇數值模擬精度、自適應度高的非結構化網格。
筆者利用計算資源對網格無關性進行分析(風量不隨網格數量增加而變化),之后將流體域網格尺寸設置為30 mm,旋轉域和風架網格尺寸設置為5 mm,葉片表面網格設置為2 mm,葉片前、后緣以及其余小尺寸部位網格設置為1 mm。
由于葉片為主要噪聲源,需要在其表面建立邊界層網格。其中,第一層網格高度設定為0.01 mm,總層數為20層,增長率為1.2,網格總數量約為8.5×106。
邊界層局部網格示意圖如圖2所示。
圖2 葉片網格局部放大圖
筆者對照試驗環境設置仿真模型的邊界條件,進出口均設為壓力邊界,旋轉域采用多運動參考系模型(multiple reference frame,MRF)[10],風扇轉速按照實測工況下的轉速2 270 r/min進行設置。
在穩態計算中,湍流模型采用k-ε模型[11-12],壁面函數選擇兼容性更強的可放縮壁面函數(scalable)。
筆者運用SIMPLE算法求解壓力速度耦合方程,控制方程的離散格式運用二階迎風格式進行離散,收斂殘差設為1×10-4,迭代步數設置為2 000步,設置進出口流量和靜壓等監測項,進出口流量差小于0.5%則視為收斂。
噪聲計算的流體域是在網格尺寸不變的基礎上,改為用長4 m、直徑為2 m的圓柱將旋轉域包裹起來。
穩態計算方式與氣動性能計算保持一致,并且在穩態計算結果基礎上選擇大渦模擬(large eddy simulation,LES)[13]湍流模型,以便進行瞬態仿真計算,直到進出口流量計算穩定。
筆者設置與試驗情況相同的噪聲監測點,然后采用FW-H方程計算得到噪聲監測點的聲壓信號;最后,應用傅里葉變換對聲壓信號進行處理,以獲得噪聲的頻譜圖。
由于此處主要在20 Hz~5 000 Hz頻率范圍內進行噪聲研究,筆者結合奈奎斯特采樣定理[14],設置時間步長為0.000 1 s;并且遵循保證計算準確性且節約計算資源的原則,設置總時間為0.2 s。
按照《工業通風機用標準化風道性能試驗》(GB/T 1236-2017)[15]規定,筆者對研究對象進行氣動性能測試試驗。
風扇進口為管道式,出口為與外界大氣相連的自由式,風扇靜壓可由壓力傳感器測得(精度為±1%);采用超聲波流量計測定風量,即按照正態分布取95%的置信區間,包含因子k=2,風量擴展不確定度為:U(qv)=1.21%。
風量計算公式如下:
(9)
式中:D為管道直徑,m;v為流速,m/s;c為管道內的聲速,m/s;θ為流體靜止時超聲軌跡與管道軸線之間的夾角;Δt為順逆流時間差,s。
筆者將冷卻風扇布置于試驗臺架上,調節風扇工況,得到靜壓、風量、轉速、軸功率等數據,最終計算得出風扇性能曲線。
氣動性能試驗臺如圖3所示。
圖3 氣動性能測試臺架
筆者將試驗結果與仿真結果進行對比分析,結果顯示:將風量和靜壓效率的仿真與試驗對比,各工況誤差均在5%以內,滿足工程應用條件。
氣動性能試驗與仿真結果對比如圖4所示。
圖4 氣動性能試驗與仿真結果對比
由圖4可知:風量隨靜壓的升高而減小,靜壓效率與仿真誤差隨靜壓增加而升高;產生誤差的因素有:計算模型中對復雜結構的簡化、試驗儀器安裝誤差等。
根據《風機和羅茨鼓風機噪聲測量方法》(GB/T 2888-2008)[16],筆者將風扇置于半消聲室內(其中,風扇中心軸線距地面高度為1.5 m),進行噪聲性能試驗;同時試驗所在半消聲室也符合ISO3745標準。
筆者在風扇中心軸線上距離風扇中心點前后各1 m處設置兩個噪聲監測點,以此來采集聲音信號。
聲壓級計算公式如下:
(10)
式(10)中,筆者按照正態分布取95%的置信區間,包含因子k=2,噪聲擴展不確定度為U(SPL)=1.42%。
噪聲試驗測量環境和麥克風如圖5所示。
圖5 噪聲性能試驗
噪聲性能試驗設備參數如表1所示。
表1 設備名稱及參數
進氣側(監測點1)處的試驗與仿真的前三階離散噪聲聲壓級與總聲壓級對比如表2所示。
表2 離散噪聲聲壓級與總聲壓級對比
監測點1處試驗與仿真噪聲頻譜對比如圖6所示。
圖6 監測點1處噪聲頻譜圖對比
由表2與圖6可知:試驗結果與仿真計算結果整體上較為一致,冷卻風扇離散噪聲主要集中在低頻段,且最大峰值均出現在二階離散噪聲處,高頻段以渦流噪聲為主。
監測點處的A計權總聲壓級仿真結果較試驗值略低。筆者分析該誤差原因有:邊界層網格未精確模擬葉片表面邊界層運動情況;冷卻風扇旋轉時產生的結構振動噪聲與電機運轉產生的電磁噪聲使噪聲試驗值增加。
綜上可知,監測點1處的噪聲值與數值計算結果誤差在5%以內。由此可見,該數值模擬方法可用于進一步的風扇降噪研究。
當電機驅動環形冷卻風扇旋轉時,在摩擦阻力的[5]作用下,其葉片吸力面表面氣流改變了原有旋轉方向,沿徑向甩出后與導流環相互干涉形成回流;這些回流與葉片尾跡同后續葉片前緣相互作用,又形成離散頻率噪聲,葉片表面的漩渦脫落產生了渦流噪聲[17-19]。
葉尖附近是形成低壓回流的主要區域,也是離散頻率噪聲形成的主要部位,為此,筆者對葉尖部位進行穿孔設計,以達到降噪的目的。
穿孔定位方式如圖7所示。
圖7 穿孔定位方式
穿孔模型對比如圖8所示。
圖8 穿孔示意圖
將葉片尾緣繞輪轂圓心順時針分別旋轉3°、6°和9°,可以得到三條弧線L、M、N,作為孔的周向定位方式;然后,將葉片沿徑向進行15等分,取靠近導流環的兩條等分線O、P作為徑向定位方式,兩者交點為穿孔最終定位。
模型A以弧線M與等分線O、P相交,相交的兩個點作為圓心進行穿孔;模型B在模型A的基礎上增加圓心為弧線L與等分線O、P交點的兩個孔;模型C在模型B的基礎上增加圓心為弧線N與等分線O、P交點的兩個孔,穿孔半徑為2 mm。
三種穿孔模型與原模型的氣動性能對比結果如圖9所示。
圖9 氣動性能對比
根據圖9可知:與原模型相比,三種穿孔模型風量和靜壓效率隨穿孔數量的增加逐漸減小。這是由于葉片壓力面與吸力面壓差作用,使得部分氣體由壓力面流向吸力面,導致氣體泄漏,使得風量與靜壓效率下降。
葉片吸力面靜壓云圖如圖10所示。
圖10 吸力面靜壓云圖
葉尖速度矢量圖如圖11所示。
圖11 吸力面葉尖速度矢量圖
由圖10和圖11可知:原模型與三種優化模型葉片前緣靜壓達到最大,葉尖與導流環連接處出現低壓區[20]。這主要是由于吸力面表面的部分氣流在摩擦阻力的作用下,改變了原有旋轉方向沿徑向甩出,甩出后與導流環產生干涉,進而形成葉尖吸力面的低壓回流區。
由于穿孔設計,部分流體從壓力側流向吸入側形成射流,這使得葉片吸力面壓力和速度分布發生變化。隨著穿孔數的增加,葉尖穿孔部位低壓區域減小、壓力梯度減小,并且吸力面穿孔附近氣流速度增大,使得邊界層內的流體獲得新的動能以克服摩擦阻力,因此,有效地抑制了回流的產生。
t=0.1 s時刻的渦量云圖如圖12所示。
圖12 t=0.1 s時渦量云圖
由圖12可知:在導流環內圓以及葉片尾緣附近出現大量漩渦脫落現象。葉片穿孔引起的射流使邊界層低速氣流重新獲得動能,邊界層分離現象得到改善,尾緣渦量得到一定控制。但射流與吸力面融合時,穿孔位置附近也產生了部分漩渦。
模型A穿孔數量較少,尾緣附近渦量分布改善不明顯;增加模型B的穿孔后,尾緣渦量明顯減少;雖然模型C尾緣處渦量改善最明顯,但其穿孔數量過多且射流速度高,造成氣體泄漏過多,產生了額外的分離渦。
模型A、B、C與原模型噪聲頻譜圖對比如圖13所示。
圖13 原模型與其他三種模型噪聲頻譜圖對比
三種優化模型與原模型前三階離散頻率噪聲聲壓級與總聲壓級對比如表3所示。
表3 原模型與其他三種模型噪聲對比
由圖13與表3可知:三種優化模型對風扇噪聲性能的影響與壓力云圖展現情況一致,穿孔結構有效抑制了葉尖低壓回流區的產生,使得后續葉片切割的氣流更加穩定有序,從而降低了離散頻率噪聲(因葉片周期性切割氣流而產生)。
三分之一倍頻程圖如圖14所示。
圖14 三分之一倍頻程圖
由圖14可知:寬頻噪聲與渦量云圖情況展現一致。由此可見,穿孔設計不僅可以抑制回流的產生,并且對尾渦脫落的產生也有一定影響。
模型A由于穿孔數量較少,帶來的吸力面壓力和速度變化也較小,所以對高頻段渦流噪聲的影響同樣較小;模型B穿孔數量增加,穿孔射流使得邊界層內的流體獲得了新的動能,以克服摩擦阻力,在一定程度上抑制了尾渦脫落,減小了高頻段的渦流噪聲;模型C由于穿孔數量過多,而射流速度過大,因此產生了新的分離渦,致使高頻段渦流噪聲有所增加。
因此,數量適當的穿孔設計可以有效降低環形冷卻風扇的氣動噪聲。
為了降低車用環形冷卻風扇的氣動噪聲,筆者設計了一種葉片穿孔環形冷卻風扇,并采用數值模擬的方法,探究了葉片穿孔數對車用環形冷卻風扇氣動性能和噪聲性能的影響規律。
首先,建立了準確的冷卻風扇計算模型;然后,設計了具有不同穿孔數的葉片結構;最后,對原風扇模型和葉片穿孔模型的流場與聲場進行了分析,得出了其影響規律,以期為環形冷卻風扇噪聲控制技術提供參考。
研究結論如下:
1)葉片穿孔結構可以有效降低風扇的氣動噪聲,穿孔數量過少時降噪效果不明顯,但并不是數量越多越好,穿孔數量過多反而會產生多余的分離渦,致使噪聲增加。故其中模型B降噪效果最好,A計權總聲壓級較原模型風扇降低了4 dB;且隨著穿孔數量的增加,風量與靜壓效率相應降低;
2)合理的葉片穿孔結構產生的射流能夠改善葉片吸力面的壓力和速度分布情況,減少葉尖部位的回流以及尾渦脫落現象,在離散頻率噪聲與渦流噪聲方面取得一定的降噪效果;
3)葉片穿孔結構產生的射流造成了氣體泄漏,故風量和靜壓效率都有所降低,且隨穿孔數量的增加,風量和效率下降越明顯,因此,應合理選擇穿孔數目。
在接下來的工作中,筆者將研究其他穿孔參數對環形冷卻風扇性能的影響,并進行多目標參數優化,以便得到最優葉片穿孔結構。