李 亮，張 斌，付 玲，肖慶麟，滕新科，張 勁（.中聯(lián)重科股份有限公司，湖南長(zhǎng)沙4003；2.建設(shè)機(jī)械關(guān)鍵技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖南長(zhǎng)沙4003）

掃路車(chē)用風(fēng)機(jī)內(nèi)部渦流特性與壓力脈動(dòng)分析

李 亮1，2，張 斌1，2，付 玲1，2，肖慶麟1，滕新科1，2，張 勁1，2
（1.中聯(lián)重科股份有限公司，湖南長(zhǎng)沙410013；2.建設(shè)機(jī)械關(guān)鍵技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖南長(zhǎng)沙410013）

以某型號(hào)掃路車(chē)用風(fēng)機(jī)內(nèi)部渦流和壓力脈動(dòng)特性為研究對(duì)象，建立適宜于掃路車(chē)用風(fēng)機(jī)內(nèi)部湍流流動(dòng)的C F D數(shù)值模擬方法，仿真與試驗(yàn)結(jié)果的對(duì)比驗(yàn)證了所提方法的合理性和可靠性.分析表明：掃路車(chē)用風(fēng)機(jī)靠近蝸舌區(qū)域的葉輪流道會(huì)產(chǎn)生強(qiáng)烈的漩渦流動(dòng)乃至回流，該流動(dòng)特征是風(fēng)機(jī)氣動(dòng)噪聲的主要來(lái)源；蝸舌靠近出口始終存在一個(gè)低速、渦量值較小的漩渦區(qū)，蝸舌結(jié)構(gòu)形式及葉片與蝸舌間距是影響風(fēng)機(jī)氣動(dòng)噪聲特性的關(guān)鍵參數(shù)；葉片壓力面比吸力面上的脈動(dòng)周期性特征更強(qiáng)烈，正對(duì)葉輪流道的蝸舌表面壓力脈動(dòng)最強(qiáng).

聲學(xué)；掃路車(chē)風(fēng)機(jī)；渦流；蝸舌；數(shù)值模擬

隨著國(guó)家環(huán)境戰(zhàn)略的強(qiáng)勢(shì)推進(jìn)，作為城鎮(zhèn)環(huán)衛(wèi)作業(yè)主力軍的掃路車(chē)的高噪高能耗特性與日益提高的低噪節(jié)能環(huán)保要求之間的矛盾日益突出.掃路車(chē)用風(fēng)機(jī)是掃路車(chē)氣力系統(tǒng)的關(guān)鍵零部件、動(dòng)力源和主要噪聲源［1］，其降噪需求也日益迫切.

目前，針對(duì)掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲的發(fā)聲機(jī)理、聲源特點(diǎn)以及聲傳播特性等研究仍顯不足.因此，開(kāi)展掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲分析是掃路車(chē)技術(shù)領(lǐng)域中亟待探索的關(guān)鍵研究方向.

本文以某型號(hào)掃路車(chē)用風(fēng)機(jī)內(nèi)部渦流和壓力脈動(dòng)特性為研究對(duì)象，對(duì)掃路車(chē)用風(fēng)機(jī)在額定工況下的內(nèi)部流動(dòng)進(jìn)行三維數(shù)值仿真.首先對(duì)風(fēng)機(jī)內(nèi)部渦流特性進(jìn)行分析，然后結(jié)合時(shí)域分析方法對(duì)流場(chǎng)內(nèi)部壓力脈動(dòng)進(jìn)行分析，以揭示掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲產(chǎn)生的機(jī)理、氣動(dòng)噪聲源的位置及噪聲類(lèi)型，為掃路車(chē)用風(fēng)機(jī)氣動(dòng)降噪提供有益的信息.

1　幾何模型及觀(guān)察點(diǎn)選取

本文數(shù)值模擬的掃路車(chē)用風(fēng)機(jī)主要由進(jìn)口段、出口段、葉輪以及蝸殼構(gòu)成，具體結(jié)構(gòu)如圖1所示，其主要參數(shù)詳見(jiàn)表1.

圖1　風(fēng)機(jī)結(jié)構(gòu)示意圖Fig.1 Schematic diagramofcentrifugal fan

表1　風(fēng)機(jī)主要結(jié)構(gòu)參數(shù)Tab.1 Main structure parameters of centrifugal fan

在蝸殼圓周上每隔30°取一個(gè)觀(guān)察點(diǎn)Pi，總共選取11個(gè)觀(guān)察點(diǎn)；在蝸舌上選取3個(gè)測(cè)試點(diǎn)A （P1）、B（P2）、C（P3）；在靠近蝸舌附近的3個(gè)葉片尾端上取6個(gè)觀(guān)察點(diǎn)，壓力面上的點(diǎn)為Di，吸力面上的點(diǎn)為di.所有點(diǎn)均處在離心風(fēng)機(jī)葉輪中截面上，部分觀(guān)察點(diǎn)具體位置如圖2所示.

圖2　部分測(cè)試點(diǎn)位置示意圖Fig.2 Location drawing of some viewpoints

2　數(shù)值計(jì)算方法

2.1控制方程及湍流模型

掃路車(chē)用風(fēng)機(jī)的內(nèi)部流動(dòng)是三維黏性不可壓縮流動(dòng)，遵循物理守恒定律，控制方程包括連續(xù)性方程、動(dòng)量守恒方程以及湍流輸運(yùn)方程.

連續(xù)性方程：

動(dòng)量守恒方程：

上述式中：ρ為密度；p為壓強(qiáng)；ν為運(yùn)動(dòng)黏度；fi為體積力；μe為湍流黏性系數(shù)；μ為分子黏性系數(shù)，μt為渦流黏性系數(shù).

由于掃路車(chē)用風(fēng)機(jī)葉輪高速旋轉(zhuǎn)和強(qiáng)曲率效應(yīng)的影響，其內(nèi)部流場(chǎng)結(jié)構(gòu)極其復(fù)雜，包含有葉片尾緣的分離尾跡流動(dòng)、漩渦脫落以及泄漏流動(dòng)等復(fù)雜流動(dòng)現(xiàn)象，因此本文采用工程上應(yīng)用廣泛的R NGk-ε湍流模型.該模型通過(guò)對(duì)湍流黏性進(jìn)行修正，考慮了旋轉(zhuǎn)和曲率效應(yīng)對(duì)流動(dòng)的影響，與標(biāo)準(zhǔn)k-ε模型相比可以較好地處理高應(yīng)變率及流線(xiàn)彎曲程度較大的流動(dòng)［2］.在R NGk-ε模型中，k方程和ε方程分別為：

2.2網(wǎng)格生成

采用三維建模軟件進(jìn)行實(shí)體建模，運(yùn)用專(zhuān)業(yè)網(wǎng)格劃分軟件生成網(wǎng)格.為了使掃路車(chē)用風(fēng)機(jī)內(nèi)部流動(dòng)充分發(fā)展，在風(fēng)機(jī)進(jìn)、出風(fēng)口處各添加了一段延長(zhǎng)管道.在劃分網(wǎng)格時(shí)，考慮到掃路車(chē)用風(fēng)機(jī)結(jié)構(gòu)的復(fù)雜性以及網(wǎng)格的生成質(zhì)量，采用分塊劃分網(wǎng)格，各個(gè)區(qū)域單獨(dú)生成合適的網(wǎng)格，相鄰的區(qū)域共用一個(gè)面.劃分網(wǎng)格時(shí)首先進(jìn)行了網(wǎng)格無(wú)關(guān)性計(jì)算，確保網(wǎng)格數(shù)量對(duì)計(jì)算精度不造成影響，確立了網(wǎng)格數(shù)約為1.14×106，2.4×106，6.33×106的3種網(wǎng)格.采用不同網(wǎng)格數(shù)進(jìn)行計(jì)算所對(duì)應(yīng)的性能曲線(xiàn)與試驗(yàn)性能曲線(xiàn)對(duì)比如圖3所示，考慮到數(shù)值預(yù)估的精度，本文選取網(wǎng)格數(shù)量約為6.33×106.蝸殼和葉輪采用適應(yīng)性較強(qiáng)的非結(jié)構(gòu)化四面體網(wǎng)格，并對(duì)曲率變化較大的面進(jìn)行了網(wǎng)格加密［3］（如圖4所示）.模型各部分的網(wǎng)格數(shù)量分別為：進(jìn)口段814339，葉輪區(qū)域2717480，蝸殼區(qū)域2067721，出口段729813，風(fēng)機(jī)網(wǎng)格總數(shù)為6329353.圖4葉片和蝸舌表面網(wǎng)格

圖3　全壓與流量關(guān)系曲線(xiàn)Fig.3 Relation of total pressure and flow

圖4　葉片和蝸舌表面網(wǎng)格Fig.4 Grid distribution of blade and scroll

2.3近壁面處理方法

近壁面區(qū)域的流動(dòng)由于其自身復(fù)雜特性及重要作用，一直是湍流模擬的熱點(diǎn)和難點(diǎn)，近壁面的處理對(duì)于湍流計(jì)算的精度也具有重要影響［4］.目前，應(yīng)用較廣的湍流數(shù)值計(jì)算方法是R AN S方法，它需要引入湍流模型才能使雷諾方程封閉可解.常用的R NGk-ε湍流模型是高雷諾數(shù)湍流模型，這些模型均是對(duì)于充分發(fā)展的湍流才有效.而近壁面附近區(qū)域的流動(dòng)雷諾數(shù)較低，湍流的脈動(dòng)影響不如分子黏性影響大，湍流發(fā)展并不充分，故高雷諾數(shù)湍流模型不能直接適用于近壁面附近區(qū)域的流場(chǎng)計(jì)算.為了使這些高雷諾數(shù)模型計(jì)算能夠延伸到壁面，可使用兩種方法為近壁區(qū)建模：標(biāo)準(zhǔn)壁面函數(shù)法和增強(qiáng)壁面處理法.

采用不同的壁面處理，對(duì)網(wǎng)格有不同的劃分要求.考慮到掃路車(chē)用風(fēng)機(jī)強(qiáng)旋轉(zhuǎn)和強(qiáng)曲率效應(yīng)，結(jié)合標(biāo)準(zhǔn)壁面函數(shù)法要求，本文y+的范圍確定為11. 25≤y+≤300［5］.

本文近壁面流動(dòng)的處理方法采用標(biāo)準(zhǔn)壁面函數(shù)法，葉片和蝸舌表面的y+分布如圖5所示，從圖中可以看出本文數(shù)值計(jì)算的網(wǎng)格劃分能滿(mǎn)足上述要求.

圖5　葉片和蝸舌表面y+分布Fig.5 y+distribution of blade and scroll

2.4方程離散與求解

計(jì)算時(shí)采用“多重坐標(biāo)系”（MR F）模型耦合動(dòng)靜部分區(qū)域，即葉輪區(qū)域?yàn)樾D(zhuǎn)區(qū)域，采用旋轉(zhuǎn)坐標(biāo)系，給定相應(yīng)的旋轉(zhuǎn)速度；其余區(qū)域?yàn)殪o止區(qū)域，采用靜止坐標(biāo)系.離散方法采用有限體積法，壓力速度耦合關(guān)系采用S I MP L E算法；湍動(dòng)能、耗散率、動(dòng)量方程的離散采用二階迎風(fēng)格式，使用速度入口及壓力出口邊界條件，其中進(jìn)口速度為24.76m/s，出口壓力為1個(gè)大氣壓.

葉片在葉輪內(nèi)沿圓周均勻分布，且轉(zhuǎn)速恒定，計(jì)算時(shí)間步長(zhǎng)為6.41×10-5s，這個(gè)時(shí)間步長(zhǎng)的選擇與葉輪的轉(zhuǎn)速有關(guān)，即每經(jīng)過(guò)360個(gè)時(shí)間步長(zhǎng)，葉輪旋轉(zhuǎn)1周（時(shí)間為T(mén)）.這個(gè)時(shí)間很短，能夠獲取葉輪所需要的旋轉(zhuǎn)角度，每個(gè)時(shí)間步長(zhǎng)迭代次數(shù)為20，葉輪轉(zhuǎn)速為2600r/m i m.

對(duì)掃路車(chē)用風(fēng)機(jī)內(nèi)流場(chǎng)先作定常計(jì)算，取殘差小于1×10-4時(shí)計(jì)算收斂.然后將定常計(jì)算的結(jié)果作為非定常計(jì)算的初始化數(shù)據(jù)，進(jìn)行非定常計(jì)算.

3　計(jì)算結(jié)果及其分析

3.1方法驗(yàn)證

本文對(duì)掃路車(chē)用風(fēng)機(jī)在不同工況下的流動(dòng)情況進(jìn)行了預(yù)測(cè)，并與試驗(yàn)結(jié)果進(jìn)行對(duì)比分析.按照G B/T 1236—1985《通風(fēng)機(jī)空氣動(dòng)力性能試驗(yàn)方法》的要求對(duì)掃路車(chē)用風(fēng)機(jī)進(jìn)行試驗(yàn)，試驗(yàn)裝置如圖6所示.

圖6　試驗(yàn)裝置Fig.6 Testing device

本文采用無(wú)量綱流量系數(shù)和全壓系數(shù)對(duì)結(jié)果進(jìn)行處理［6］，各工況下效率和全壓系數(shù)的仿真與測(cè)試結(jié)果的對(duì)比如圖7，8所示.從圖中可以看出，數(shù)值模擬結(jié)果的與試驗(yàn)結(jié)果變化趨勢(shì)一致，吻合良好，并且全壓內(nèi)效率曲線(xiàn)和全壓系數(shù)曲線(xiàn)的最大誤差在4%以?xún)?nèi)，準(zhǔn)確預(yù)測(cè)了掃路車(chē)用風(fēng)機(jī)的性能變化情況，為后續(xù)掃路車(chē)用風(fēng)機(jī)內(nèi)流場(chǎng)特性和壓力脈動(dòng)分析奠定了基礎(chǔ).

圖7　試驗(yàn)和模擬的效率曲線(xiàn)Fig.7 Efficiency curve of experiment and simulation

圖8　試驗(yàn)和模擬的全壓系數(shù)曲線(xiàn)Fig.8 Total pressure coefficient curve of experimentand simulation

3.2葉輪內(nèi)部渦流特征分析

由于流體的黏性作用以及旋轉(zhuǎn)部件與靜止部件的動(dòng)靜干涉，掃路車(chē)用風(fēng)機(jī)內(nèi)部流場(chǎng)呈現(xiàn)出全三維、非定常的流動(dòng)特征.這種流動(dòng)特征引起壓力場(chǎng)的壓力脈動(dòng)，導(dǎo)致掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲的產(chǎn)生.

圖9為掃路車(chē)用離心葉輪1個(gè)葉片流道內(nèi)同一截面、不同時(shí)刻的瞬態(tài)流線(xiàn)圖.從圖中可以看出，在長(zhǎng)葉片進(jìn)口吸力面處，由于進(jìn)口沖角過(guò)大，使得長(zhǎng)葉片吸力面上產(chǎn)生流動(dòng)分離，整個(gè)流道流動(dòng)情況惡化，該過(guò)程一直延續(xù)到葉輪出口，并在葉輪流道出口靠近蝸舌附近形成一個(gè)回流區(qū)，該回流區(qū)不斷地與蝸舌表面相互作用，使得蝸舌表面的壓力產(chǎn)生脈動(dòng)，形成掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲的一個(gè)來(lái)源.同時(shí)，在短葉片吸力面處，流道中一直存在一個(gè)大的漩渦區(qū)，幾乎堵塞整個(gè)流道，該漩渦區(qū)域與回流區(qū)的流動(dòng)相互干涉，進(jìn)行能量交換，使得漩渦在1個(gè)葉輪旋轉(zhuǎn)周期內(nèi)，周期性地經(jīng)歷產(chǎn)生、發(fā)展、耗散過(guò)程，該過(guò)程加劇了掃路車(chē)用風(fēng)機(jī)內(nèi)部壓力的脈動(dòng)，成為掃路車(chē)用風(fēng)機(jī)內(nèi)部渦流噪聲的主要來(lái)源［7］.上述分析結(jié)論表明，在掃路車(chē)用風(fēng)機(jī)內(nèi)部，由于蝸殼結(jié)構(gòu)的非對(duì)稱(chēng)性作用，靠近蝸舌區(qū)域的葉輪流道會(huì)產(chǎn)生強(qiáng)烈的漩渦流動(dòng)乃至回流，該非定常流動(dòng)增加了葉輪氣流作用在蝸舌表面的脈動(dòng)作用力，增大了掃路車(chē)用風(fēng)機(jī)的氣動(dòng)噪聲.

圖9　葉輪與蝸舌相互作用瞬態(tài)流線(xiàn)圖Fig.9 Transient streamlines of the impellerand volute tongue interaction

如圖10所示，在蝸舌區(qū)域，始終存在一個(gè)低速漩渦區(qū)域，該漩渦周期性地經(jīng)歷產(chǎn)生、發(fā)展、壯大、脫落過(guò)程.漩渦的產(chǎn)生過(guò)程使得蝸舌表面受到非定常脈動(dòng)力的作用，從而產(chǎn)生氣動(dòng)噪聲.在葉片的尾部，由于漩渦的脫落，使得該處的速度場(chǎng)極其復(fù)雜，流體微團(tuán)受到拉伸、擠壓和扭曲變形作用［8-9］，從而產(chǎn)生較大的渦量值.湍流渦聲理論［10］認(rèn)為，氣動(dòng)噪聲與漩渦的拉壓和變形密切相關(guān)，渦量是關(guān)聯(lián)氣動(dòng)噪聲源的關(guān)鍵湍流量［11］.上述分析表明在掃路車(chē)用風(fēng)機(jī)靠近出口的蝸舌區(qū)域，始終存在一個(gè)低速、渦量值較小的漩渦區(qū)，而葉片尾部較大渦量值的存在加劇了風(fēng)機(jī)氣動(dòng)噪聲的產(chǎn)生.因此，蝸舌結(jié)構(gòu)是影響氣動(dòng)噪聲特性的關(guān)鍵參數(shù)，在其設(shè)計(jì)中需要著重考慮才能使得其氣動(dòng)噪聲性能達(dá)到最佳.

圖10　蝸舌附近速度大小隨時(shí)間變化Fig.10Changes of volocity near volutetongue with time

3.3計(jì)算結(jié)果的時(shí)域分析

由于實(shí)際采集得到的壓力信號(hào)值較大，不利于分析觀(guān)察，因此本文將采集到的360個(gè)壓力脈動(dòng)信號(hào)進(jìn)行去平均處理，即用各時(shí)間點(diǎn)壓力脈動(dòng)值減去采集時(shí)間內(nèi)壓力平均值［12］：

式中：pi為原始采集壓力信號(hào)每一時(shí)間點(diǎn)壓力脈動(dòng)值；Pi為去平均處理后每一時(shí)間點(diǎn)壓力脈動(dòng)值；Pj為原始采集壓力信號(hào)每一時(shí)間點(diǎn)壓力脈動(dòng)值.

如圖11所示，在蝸殼周向表面，蝸舌附近的壓力脈動(dòng)最強(qiáng)烈，且B點(diǎn)處的壓力脈動(dòng)最強(qiáng)，這是由于葉輪流道流出的流體垂直沖擊到蝸舌表面，使得此處的壓力脈動(dòng)最強(qiáng)，因此是氣動(dòng)噪聲的主要來(lái)源.圖12中葉片與蝸舌表面噪聲源分布（脈動(dòng)壓力大小分布）也說(shuō)明了這點(diǎn).從圖12中可以看出，蝸舌表面和靠近蝸舌附近的葉輪流道噪聲源分布值最大，其他位置處噪聲源值要小得多.這主要是由于葉輪和蝸舌的動(dòng)靜干涉，使得該流道附近存在漩渦流動(dòng)，流體渦團(tuán)周期性地從葉片尾端脫落，因此導(dǎo)致此處壓力脈動(dòng)增強(qiáng)，從而產(chǎn)生較大的噪聲源.上述分析說(shuō)明在掃路車(chē)用風(fēng)機(jī)內(nèi)部，正是由于葉輪與蝸舌非定常作用產(chǎn)生的漩渦，導(dǎo)致該區(qū)域的噪聲源值最大，成為風(fēng)機(jī)氣動(dòng)噪聲的主要來(lái)源.

圖11　蝸殼周向壓力脈動(dòng)情況Fig.11 Pressure fluctuation of thecircumferential volute

圖12　葉片與蝸舌表面噪聲源分布Fig.12 Noise source distribution ofimpeller and volute tongue

從圖13可以看出，在1個(gè)周期內(nèi)，葉片壓力面上的脈動(dòng)峰值有30個(gè)，這與葉片數(shù)量（15片長(zhǎng)葉片+15片短葉片）對(duì)蝸舌表面的非定常作用相吻合，呈明顯的周期性.而吸力面上壓力脈動(dòng)的周期性則要微弱得多，這主要是由于吸力面上容易產(chǎn)生流動(dòng)分離，形成漩渦從而產(chǎn)生渦流噪聲.由于葉片與蝸舌的非定常作用存在周期性，說(shuō)明可以通過(guò)增大葉片與蝸舌的間距來(lái)削弱它們之間的非定常作用，由此降低掃路車(chē)用風(fēng)機(jī)的離散噪聲.

從圖14可以看出，與葉片壓力面上的脈動(dòng)值不同，蝸舌表面的壓力脈動(dòng)并沒(méi)有產(chǎn)生周期性變化.蝸舌表面的3個(gè)觀(guān)察點(diǎn)中，B點(diǎn)的壓力脈動(dòng)值最大，A、C兩點(diǎn)的壓力脈動(dòng)相對(duì)較弱.這主要是由于B點(diǎn)正對(duì)著葉輪流道氣流的打擊作用，因此該點(diǎn)壓力脈動(dòng)最強(qiáng)［13］.

圖13　葉片表面壓力脈動(dòng)隨時(shí)間變化圖Fig.13 Pressure fluctuation of the blade with time

圖14　蝸舌表面壓力脈動(dòng)隨時(shí)間變化圖Fig.14 Pressure fluctuation of the volute tongue with time

4　結(jié)語(yǔ)

（1）掃路車(chē)用風(fēng)機(jī)靠近蝸舌區(qū)域的葉輪流道會(huì)產(chǎn)生強(qiáng)烈的漩渦流動(dòng)乃至回流，該非定常流動(dòng)增加了葉輪氣流作用在蝸舌表面的非定常脈動(dòng)作用力，使得該區(qū)域的噪聲源值最大，成為掃路車(chē)用風(fēng)機(jī)氣動(dòng)噪聲的主要來(lái)源.

（2）掃路車(chē)用風(fēng)機(jī)蝸舌結(jié)構(gòu)形式及葉片與蝸舌間距對(duì)其氣動(dòng)噪聲的產(chǎn)生具有重要影響，在設(shè)計(jì)時(shí)需要著重考量才能使得風(fēng)機(jī)氣動(dòng)噪聲性能達(dá)到最佳.

（3）葉片壓力面比吸力面上的脈動(dòng)周期性特征更強(qiáng)烈；正對(duì)葉輪流道的蝸舌區(qū)域壓力脈動(dòng)最強(qiáng).

［1］ 張斌，李亮，付玲，等.清掃車(chē)專(zhuān)用風(fēng)機(jī)內(nèi)流仿真分析與優(yōu)化設(shè)計(jì)［J］.機(jī)械設(shè)計(jì)與制造，2014，9（9）：246-248.ZHANG Bin，LI Liang，F(xiàn)U Ling，et al.Analysis of theoptimized design and numerical simulation on the centrifugalfan of sweepers［J］.Machinery Design & Manufacture，2014，9（9）：246-248.

［2］ Ansys Inc.Ansys fluent theory guide［M］.Canonsburg：Ansys Inc.，2011.

［3］ BALLESTEROS-TAJADURA R， VELARDE-SUARREZ.Numerical calculation of pressure fluctuations in the volute of acentrifugal fan［J］.Journal of Fluids Engineering，2006，128：359-369.

［4］ 張建.采用修正來(lái)流條件和粗糙壁面處理方法的繞流問(wèn)題研究［D］.北京：北京交通大學(xué)，2011.

ZHANG Jian.Analysis of flow around bluff body based on modified inlet profiles and rough wall treatments［D］.Beijing：Beijing Jiaotong University，2011.

［5］ 張斌.大渦模擬濾波網(wǎng)格分析及網(wǎng)格自適應(yīng)控制研究與應(yīng)用［D］.上海：上海交通大學(xué)，2011.

ZHANG Bin.Research and application of filtering grid scale and meshing adaptive-control strategy for large eddy simulation［D］.Shanghai：Shanghai Jiaotong University，2011.

［6］ 李亮，王瑞，金晗輝.基于非線(xiàn)性湍流模型的多翼離心風(fēng)機(jī)內(nèi)部流場(chǎng)三維數(shù)值模擬及性能預(yù)測(cè)［J］.流體機(jī)械，2014，40（10）：41-45.

LI Liang，WANG Rui，JIN Hanhui.Three dimensional numerical simulation of internal flow field in multi-blade centrifugal fan based on the non-linear turbulent model［J］.Fluid Machinery，2014，40（10）：41-45.

［7］ 付雙成，劉雪冬，鄒鑫，等.傾斜蝸舌對(duì)多翼離心風(fēng)機(jī)流場(chǎng)及噪聲的影響［J］.噪聲與振動(dòng)控制，2013，33（3）：87-89.

FU Shuangcheng，LIU Xuedong，ZOU Xin，et al.Numerical simulation of the influence of inclined volute tongue on flow field and noise of multi-blade centrifugal fans［J］.Noise and Vibration Control，2013，33（3）：87-89.

［8］ 陳慶光，房菲，張永超，等.對(duì)旋風(fēng)機(jī)級(jí)間流場(chǎng)渦量和噪聲分布規(guī)律的分析［J］.噪聲與振動(dòng)控制，2012，2（1）：89-93.

CHEN Qingguang，F(xiàn)ANG Fei，ZHANG Yongchao，et al.Vortex and noise characteristics analysis of the flow field between the two impellers of a contra rotating axial fan［J］. Noise and Vibration Control，2012，2（1）：89-93.

［9］JAATINEN A，GRONMAN A，TURUNEN-SAARESTI T.Effect of vaneless diffuser wid th on the overall performance of a centrifugal compressor［J］.Journal of Power and Energy，2011，225（5）：665-673.

［10］ 張強(qiáng).氣動(dòng)聲學(xué)基礎(chǔ)［M］.北京：國(guó)防工業(yè)出版社，2012.

ZHANG Qiang.Aero-acoustic foundation ［M ］.Beijing：National Defense Industrial Press，2012.

［11］ZHANG Bin， WANG Tong， GU Chuanggang. Blade optimization design and performance investigations of anultralow specific speed centrifugal blowers［J］.Science inChina Series E：Technological Sciences，2011，54（1）：203-210.

［12］ 卜鋒斌，趙博，代學(xué)昌，等.船用離心風(fēng)機(jī)壁面周向壓力脈動(dòng)測(cè)試及分析［J］.船舶工程，2013，35（6）：40-43.

BU Fengbin，ZHAO Bo，DAI Xuechang，et al.Dynamic test and analysis of circumferential distribution of static pressureon marine fan volute wall［J］.Ship Engineering，2013，35（6）：40-43.

［13］ 王軍，冉苗苗，姚姍姍.多翼離心風(fēng)機(jī)的氣動(dòng)噪聲預(yù)測(cè)［J］.工程熱物理學(xué)報(bào)，2014，29（7）：1144-1146.

WANG Jun，RAN Miaomiao，YAO Shanshan.Aero-acoustics prediction of multi-blade centrifugal fan［J］.Journal of Engineering Thermophysics，2014，29（7）：1144-1146.

Analysis on internal vortex and pressure pulsationproperties of road sweeper fan

LI Liang1，2，ZHANG Bin1，2，F(xiàn)U Ling1，2，XIAO Qing-lin1，TENG Xin-ke1，ZHANG Jin1，2（1.Zoomlion Heavy Industry Science &Technology Co.，Ltd.，Changsha 410013，China；2.State Key Laboratory of Construction Machinery，Changsha 410013，China）

Pertaining to the internal vortex and pressure pulsation properties of a specific road sweeperfan，a CFD numerical simulation method is first established.Then，the rationality and reliability of theproposed method are verified in comparison with simulation and testing results.Finally，it is indicatedthat，due to the strong vortex flow even recirculation near the volute tongue region，the fluidity is themajor source of fan aerodynamic noise.Owing that a low-speed vortex zone with small vorticity valuesexists near volute tongue outlet，the structural form of volute tongue and distance between blade andvolute tongue significantly impact aerodynamic noise.In addition，the periodic property of blade pressuresurface is stronger than that of blade suction surface，whereas the pressure pulsation of volute tongueagainst impeller flow passage is strongest.

acoustics；road sweeper fan；vortex；volute tongue；numerical simulation

TH 442

A

1672-5581（2015）06-0486-07

湖南省自然科學(xué)基金資助項(xiàng)目（13J J 3130）

李 亮（1985-），男，工程師.E-m a i l：20924010@z j u.e d u.c n

張 斌（1983-），男，高級(jí)工程師.E-m a i l：s j t u 2009@163.c Om