摘要：以某型號洗路車掃路車風(fēng)機內(nèi)部渦流和壓力脈動特性為研究對象，建立適宜于掃路車用風(fēng)機內(nèi)部湍流流動的C F D數(shù)值模擬方法，仿真與試驗結(jié)果的對比驗證了所提方法的合理性和可靠性.分析表明：掃路車用風(fēng)機靠近蝸舌區(qū)域的葉輪流道會產(chǎn)生強烈的漩渦流動乃至回流，該流動特征是風(fēng)機氣動噪聲的主要來源；蝸舌靠近出口始終存在一個低速、渦量值較小的漩渦區(qū)，蝸舌結(jié)構(gòu)形式及葉片與蝸舌間距是影響風(fēng)機氣動噪聲特性的關(guān)鍵參數(shù)；葉片壓力面比吸力面上的脈動周期性特征更強烈，正對葉輪流道的蝸舌表面壓力脈動很強.
關(guān)鍵詞：洗路車風(fēng)機；掃路車風(fēng)機；垃圾車吸風(fēng)機；蝸舌；數(shù)值模擬
隨著國家環(huán)境戰(zhàn)略的強勢推進，作為城鎮(zhèn)環(huán)衛(wèi)作業(yè)主力軍的掃路車的高噪高能耗特性與日益提高的低噪節(jié)能環(huán)保要求之間的矛盾日益突出.掃路車用風(fēng)機是掃路車氣力系統(tǒng)的關(guān)鍵零部件、動力源和主要噪聲源［1］，其降噪需求也日益迫切.
目前，針對洗路車掃路車風(fēng)機氣動噪聲的發(fā)聲機理、聲源特點以及聲傳播特性等研究仍顯不足.因此，開展掃路車用風(fēng)機氣動噪聲分析是掃路車技術(shù)領(lǐng)域中亟待探索的關(guān)鍵研究方向.
本文以某型號掃路車用風(fēng)機內(nèi)部渦流和壓力脈動特性為研究對象，對掃路車用風(fēng)機在額定工況下的內(nèi)部流動進行三維數(shù)值仿真.首先對風(fēng)機內(nèi)部渦流特性進行分析，然后結(jié)合時域分析方法對流場內(nèi)部壓力脈動進行分析，以揭示掃路車用風(fēng)機氣動噪聲產(chǎn)生的機理、氣動噪聲源的位置及噪聲類型，為掃路車用風(fēng)機氣動降噪提供有益的信息.
1 幾何模型及觀察點選取
本文數(shù)值模擬的掃路車用風(fēng)機主要由進口段、出口段、葉輪以及蝸殼構(gòu)成，具體結(jié)構(gòu)如圖1所示，其主要參數(shù)詳見表1.
圖1 風(fēng)機結(jié)構(gòu)示意圖
Fig.1 Schematic diagramofcentrifugal fan
表1 風(fēng)機主要結(jié)構(gòu)參數(shù)
Tab.1 Main structure parameters of centrifugal fan
進口直徑D1/mm葉輪直徑D2/mm葉輪寬度b/mm蝸殼寬度B/mm進口安裝角β1/（°）出口安裝角β2/（°）388 700 104.4 161 40 151.5
在蝸殼圓周上每隔30°取一個觀察點Pi，總共選取11個觀察點；在蝸舌上選取3個測試點A （P1）、B（P2）、C（P3）；在靠近蝸舌附近的3個葉片尾端上取6個觀察點，壓力面上的點為Di，吸力面上的點為di.所有點均處在離心風(fēng)機葉輪中截面上，部分觀察點具體位置如圖2所示.
圖2 部分測試點位置示意圖
Fig.2 Location drawing of some viewpoints
2 數(shù)值計算方法
2.1控制方程及湍流模型
掃路車用風(fēng)機的內(nèi)部流動是三維黏性不可壓縮流動，遵循物理守恒定律，控制方程包括連續(xù)性方程、動量守恒方程以及湍流輸運方程.
連續(xù)性方程：
動量守恒方程：
上述式中：ρ為密度；p為壓強；ν為運動黏度；fi為體積力；μe為湍流黏性系數(shù)；μ為分子黏性系數(shù)，μt為渦流黏性系數(shù).
由于掃路車用風(fēng)機葉輪高速旋轉(zhuǎn)和強曲率效應(yīng)的影響，其內(nèi)部流場結(jié)構(gòu)極其復(fù)雜，包含有葉片尾緣的分離尾跡流動、漩渦脫落以及泄漏流動等復(fù)雜流動現(xiàn)象，因此本文采用工程上應(yīng)用廣泛的R NGk-ε湍流模型.該模型通過對湍流黏性進行修正，考慮了旋轉(zhuǎn)和曲率效應(yīng)對流動的影響，與標準k-ε模型相比可以較好地處理高應(yīng)變率及流線彎曲程度較大的流動［2］.在R NGk-ε模型中，k方程和ε方程分別為：
式中：k為湍動能；ε為湍流耗散率；Gk為湍流的產(chǎn)生項；模型常數(shù)分別為Cμ=0.09，C1ε=1.42，C2ε=1.68，σk=1.0，σε=1.3.
2.2 網(wǎng)格生成
采用三維建模軟件進行實體建模，運用專業(yè)網(wǎng)格劃分軟件生成網(wǎng)格.為了使掃路車用風(fēng)機內(nèi)部流動充分發(fā)展，在風(fēng)機進、出風(fēng)口處各添加了一段延長管道.在劃分網(wǎng)格時，考慮到掃路車用風(fēng)機結(jié)構(gòu)的復(fù)雜性以及網(wǎng)格的生成質(zhì)量，采用分塊劃分網(wǎng)格，各個區(qū)域單獨生成合適的網(wǎng)格，相鄰的區(qū)域共用一個面.劃分網(wǎng)格時首*行了網(wǎng)格無關(guān)性計算，確保網(wǎng)格數(shù)量對計算精度不造成影響，確立了網(wǎng)格數(shù)約為1.14×106，2.4×106，6.33×106的3種網(wǎng)格.采用不同網(wǎng)格數(shù)進行計算所對應(yīng)的性能曲線與試驗性能曲線對比如圖3所示，考慮到數(shù)值預(yù)估的精度，本文選取網(wǎng)格數(shù)量約為6.33×106.蝸殼和葉輪采用適應(yīng)性較強的非結(jié)構(gòu)化四面體網(wǎng)格，并對曲率變化較大的面進行了網(wǎng)格加密［3］（如圖4所示）.模型各部分的網(wǎng)格數(shù)量分別為：進口段814339，葉輪區(qū)域2717480，蝸殼區(qū)域2067721，出口段729813，風(fēng)機網(wǎng)格總數(shù)為6329353.圖4葉片和蝸舌表面網(wǎng)格
圖3 全壓與流量關(guān)系曲線
Fig.3 Relation of total pressure and flow
圖4 葉片和蝸舌表面網(wǎng)格
Fig.4 Grid distribution of blade and scroll
2.3 近壁面處理方法
近壁面區(qū)域的流動由于其自身復(fù)雜特性及重要作用，一直是湍流模擬的熱點和難點，近壁面的處理對于湍流計算的精度也具有重要影響［4］.目前，應(yīng)用較廣的湍流數(shù)值計算方法是R AN S方法，它需要引入湍流模型才能使雷諾方程封閉可解.常用的R NGk-ε湍流模型是高雷諾數(shù)湍流模型，這些模型均是對于充分發(fā)展的湍流才有效.而近壁面附近區(qū)域的流動雷諾數(shù)較低，湍流的脈動影響不如分子黏性影響大，湍流發(fā)展并不充分，故高雷諾數(shù)湍流模型不能直接適用于近壁面附近區(qū)域的流場計算.為了使這些高雷諾數(shù)模型計算能夠延伸到壁面，可使用兩種方法為近壁區(qū)建模：標準壁面函數(shù)法和增強壁面處理法.
采用不同的壁面處理，對網(wǎng)格有不同的劃分要求.考慮到掃路車用風(fēng)機強旋轉(zhuǎn)和強曲率效應(yīng)，結(jié)合標準壁面函數(shù)法要求，本文y+的范圍確定為11. 25≤y+≤300［5］.
本文近壁面流動的處理方法采用標準壁面函數(shù)法，葉片和蝸舌表面的y+分布如圖5所示，從圖中可以看出本文數(shù)值計算的網(wǎng)格劃分能滿足上述要求.
圖5 葉片和蝸舌表面y+分布
Fig.5 y+distribution of blade and scroll
2.4 方程離散與求解
計算時采用“多重坐標系”（MR F）模型耦合動靜部分區(qū)域，即葉輪區(qū)域為旋轉(zhuǎn)區(qū)域，采用旋轉(zhuǎn)坐標系，給定相應(yīng)的旋轉(zhuǎn)速度；其余區(qū)域為靜止區(qū)域，采用靜止坐標系.離散方法采用有限體積法，壓力速度耦合關(guān)系采用S I MP L E算法；湍動能、耗散率、動量方程的離散采用二階迎風(fēng)格式，使用速度入口及壓力出口邊界條件，其中進口速度為24.76m/s，出口壓力為1個大氣壓.
葉片在葉輪內(nèi)沿圓周均勻分布，且轉(zhuǎn)速恒定，計算時間步長為6.41×10-5s，這個時間步長的選擇與葉輪的轉(zhuǎn)速有關(guān)，即每經(jīng)過360個時間步長，葉輪旋轉(zhuǎn)1周（時間為T）.這個時間很短，能夠獲取葉輪所需要的旋轉(zhuǎn)角度，每個時間步長迭代次數(shù)為20，葉輪轉(zhuǎn)速為2600r/m i m.
對掃路車用風(fēng)機內(nèi)流場先作定常計算，取殘差小于1×10-4時計算收斂.然后將定常計算的結(jié)果作為非定常計算的初始化數(shù)據(jù)，進行非定常計算.
3 計算結(jié)果及其分析
3.1方法驗證
本文對在不同工況下的流動情況進行了預(yù)測，并與試驗結(jié)果進行對比分析.按照G B/T 1236—1985《通風(fēng)機空氣動力性能試驗方法》的要求對掃路車用風(fēng)機進行試驗，試驗裝置如圖6所示.
圖6 試驗裝置
Fig.6 Testing device
本文采用無量綱流量系數(shù)和全壓系數(shù)對結(jié)果進行處理［6］，各工況下效率和全壓系數(shù)的仿真與測試結(jié)果的對比如圖7，8所示.從圖中可以看出，數(shù)值模擬結(jié)果的與試驗結(jié)果變化趨勢*，吻合良好，并且全壓內(nèi)效率曲線和全壓系數(shù)曲線的很大誤差在4%以內(nèi)，準確預(yù)測了掃路車用風(fēng)機的性能變化情況，為后續(xù)掃路車用風(fēng)機內(nèi)流場特性和壓力脈動分析奠定了基礎(chǔ).
圖7 試驗和模擬的效率曲線
Fig.7 Efficiency curve of experiment and simulation
圖8 試驗和模擬的全壓系數(shù)曲線
Fig.8 Total pressure coefficient curve of experimentand simulation
3.2 葉輪內(nèi)部渦流特征分析
由于流體的黏性作用以及旋轉(zhuǎn)部件與靜止部件的動靜干涉，掃路車用風(fēng)機內(nèi)部流場呈現(xiàn)出全三維、非定常的流動特征.這種流動特征引起壓力場的壓力脈動，導(dǎo)致掃路車用風(fēng)機氣動噪聲的產(chǎn)生.
圖9為掃路車用離心葉輪1個葉片流道內(nèi)同一截面、不同時刻的瞬態(tài)流線圖.從圖中可以看出，在長葉片進口吸力面處，由于進口沖角過大，使得長葉片吸力面上產(chǎn)生流動分離，整個流道流動情況惡化，該過程一直延續(xù)到葉輪出口，并在葉輪流道出口靠近蝸舌附近形成一個回流區(qū)，該回流區(qū)不斷地與蝸舌表面相互作用，使得蝸舌表面的壓力產(chǎn)生脈動，形成掃路車用風(fēng)機氣動噪聲的一個來源.同時，在短葉片吸力面處，流道中一直存在一個大的漩渦區(qū)，幾乎堵塞整個流道，該漩渦區(qū)域與回流區(qū)的流動相互干涉，進行能量交換，使得漩渦在1個葉輪旋轉(zhuǎn)周期內(nèi)，周期性地經(jīng)歷產(chǎn)生、發(fā)展、耗散過程，該過程加劇了掃路車用風(fēng)機內(nèi)部壓力的脈動，成為掃路車用風(fēng)機內(nèi)部渦流噪聲的主要來源［7］.上述分析結(jié)論表明，在掃路車用風(fēng)機內(nèi)部，由于蝸殼結(jié)構(gòu)的非對稱性作用，靠近蝸舌區(qū)域的葉輪流道會產(chǎn)生強烈的漩渦流動乃至回流，該非定常流動增加了葉輪氣流作用在蝸舌表面的脈動作用力，增大了掃路車用風(fēng)機的氣動噪聲.
圖9 葉輪與蝸舌相互作用瞬態(tài)流線圖
Fig.9 Transient streamlines of the impellerand volute tongue interaction
如圖10所示，在蝸舌區(qū)域，始終存在一個低速漩渦區(qū)域，該漩渦周期性地經(jīng)歷產(chǎn)生、發(fā)展、壯大、脫落過程.漩渦的產(chǎn)生過程使得蝸舌表面受到非定常脈動力的作用，從而產(chǎn)生氣動噪聲.在葉片的尾部，由于漩渦的脫落，使得該處的速度場極其復(fù)雜，流體微團受到拉伸、擠壓和扭曲變形作用［8-9］，從而產(chǎn)生較大的渦量值.湍流渦聲理論［10］認為，氣動噪聲與漩渦的拉壓和變形密切相關(guān)，渦量是關(guān)聯(lián)氣動噪聲源的關(guān)鍵湍流量［11］.上述分析表明在掃路車用風(fēng)機靠近出口的蝸舌區(qū)域，始終存在一個低速、渦量值較小的漩渦區(qū)，而葉片尾部較大渦量值的存在加劇了風(fēng)機氣動噪聲的產(chǎn)生.因此，蝸舌結(jié)構(gòu)是影響氣動噪聲特性的關(guān)鍵參數(shù)，在其設(shè)計中需要著重考慮才能使得其氣動噪聲性能達到很佳.
圖10 蝸舌附近速度大小隨時間變化
Fig.10Changes of volocity near volutetongue with time
3.3 計算結(jié)果的時域分析
由于實際采集得到的壓力信號值較大，不利于分析觀察，因此本文將采集到的360個壓力脈動信號進行去平均處理，即用各時間點壓力脈動值減去采集時間內(nèi)壓力平均值［12］：
式中：pi為原始采集壓力信號每一時間點壓力脈動值；Pi為去平均處理后每一時間點壓力脈動值；Pj為原始采集壓力信號每一時間點壓力脈動值.
如圖11所示，在蝸殼周向表面，蝸舌附近的壓力脈動很強烈，且B點處的壓力脈動很強，這是由于葉輪流道流出的流體垂直沖擊到蝸舌表面，使得此處的壓力脈動很強，因此是氣動噪聲的主要來源.圖12中葉片與蝸舌表面噪聲源分布（脈動壓力大小分布）也說明了這點.從圖12中可以看出，蝸舌表面和靠近蝸舌附近的葉輪流道噪聲源分布值很大，其他位置處噪聲源值要小得多.這主要是由于葉輪和蝸舌的動靜干涉，使得該流道附近存在漩渦流動，流體渦團周期性地從葉片尾端脫落，因此導(dǎo)致此處壓力脈動增強，從而產(chǎn)生較大的噪聲源.上述分析說明在掃路車用風(fēng)機內(nèi)部，正是由于葉輪與蝸舌非定常作用產(chǎn)生的漩渦，導(dǎo)致該區(qū)域的噪聲源值很大，成為風(fēng)機氣動噪聲的主要來源.
圖11 蝸殼周向壓力脈動情況
Fig.11 Pressure fluctuation of thecircumferential volute
圖12 葉片與蝸舌表面噪聲源分布
Fig.12 Noise source distribution ofimpeller and volute tongue
從圖13可以看出，在1個周期內(nèi)，葉片壓力面上的脈動峰值有30個，這與葉片數(shù)量（15片長葉片+15片短葉片）對蝸舌表面的非定常作用相吻合，呈明顯的周期性.而吸力面上壓力脈動的周期性則要微弱得多，這主要是由于吸力面上容易產(chǎn)生流動分離，形成漩渦從而產(chǎn)生渦流噪聲.由于葉片與蝸舌的非定常作用存在周期性，說明可以通過增大葉片與蝸舌的間距來削弱它們之間的非定常作用，由此降低掃路車用風(fēng)機的離散噪聲.
從圖14可以看出，與葉片壓力面上的脈動值不同，蝸舌表面的壓力脈動并沒有產(chǎn)生周期性變化.蝸舌表面的3個觀察點中，B點的壓力脈動值很大，A、C兩點的壓力脈動相對較弱.這主要是由于B點正對著葉輪流道氣流的打擊作用，因此該點壓力脈動很強［13］.
圖13 葉片表面壓力脈動隨時間變化圖
Fig.13 Pressure fluctuation of the blade with time
圖14 蝸舌表面壓力脈動隨時間變化圖
Fig.14 Pressure fluctuation of the volute tongue with time
4 結(jié)語
（1）掃路車用風(fēng)機靠近蝸舌區(qū)域的葉輪流道會產(chǎn)生強烈的漩渦流動乃至回流，該非定常流動增加了葉輪氣流作用在蝸舌表面的非定常脈動作用力，使得該區(qū)域的噪聲源值很大，成為掃路車用風(fēng)機氣動噪聲的主要來源.
（2）掃路車用風(fēng)機蝸舌結(jié)構(gòu)形式及葉片與蝸舌間距對其氣動噪聲的產(chǎn)生具有重要影響，在設(shè)計時需要著重考量才能使得風(fēng)機氣動噪聲性能達到很佳.
（3）葉片壓力面比吸力面上的脈動周期性特征更強烈；正對葉輪流道的蝸舌區(qū)域壓力脈動很強.
