趙立杰,王新玲,署恒濤,段耀龍,竺菲菲

(1.沈陽航空航天大學(xué) 機(jī)電工程學(xué)院，沈陽110136；2.華南理工大學(xué) 機(jī)械與汽車工程學(xué)院，廣州 510641；3.寧波雪龍集團(tuán)股份有限公司 技術(shù)研發(fā)中心，浙江 寧波 315800)

軸流式發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能與氣動(dòng)噪聲數(shù)值分析

趙立杰1,王新玲1,署恒濤2,段耀龍3,竺菲菲3

(1.沈陽航空航天大學(xué) 機(jī)電工程學(xué)院，沈陽110136；2.華南理工大學(xué) 機(jī)械與汽車工程學(xué)院，廣州 510641；3.寧波雪龍集團(tuán)股份有限公司 技術(shù)研發(fā)中心，浙江 寧波 315800)

基于商用軟件Fluent，以《GB-T1236-2000工業(yè)通風(fēng)機(jī)用標(biāo)準(zhǔn)化風(fēng)道進(jìn)行性能試驗(yàn)》中C型裝置為標(biāo)準(zhǔn)建立了風(fēng)扇氣動(dòng)性能試驗(yàn)臺(tái)，介紹了風(fēng)扇氣動(dòng)性能的計(jì)算方法與相關(guān)模型參數(shù)的設(shè)置。以某款風(fēng)扇為例，計(jì)算了靜壓、功率和效率與流量的關(guān)系及噪聲水平。并與試驗(yàn)結(jié)果進(jìn)行對(duì)比，靜壓、功率、效率的誤差均在6%以內(nèi)；噪聲誤差雖然較大，但趨勢(shì)一致，可用于不同風(fēng)扇噪聲對(duì)比。通過研究風(fēng)扇內(nèi)部流場(chǎng)，分析了風(fēng)扇表面壓力與速度分布規(guī)律以及噪聲形成的主要原因，結(jié)果表明：葉尖間隙是風(fēng)扇效率低下的主要原因，同時(shí)葉尖渦是風(fēng)扇渦流噪音的主要貢獻(xiàn)者，如果能減小葉尖間隙不僅可提高效率，還可有效控制葉尖渦流，達(dá)到降噪的目的。

冷卻風(fēng)扇；計(jì)算流體力學(xué)；氣動(dòng)性能；氣動(dòng)噪聲

由于軸流式風(fēng)扇重量輕、成本低、結(jié)構(gòu)簡(jiǎn)單，在現(xiàn)代汽車中得到了普遍采用。作為發(fā)動(dòng)機(jī)冷卻系統(tǒng)的重要組成部分，其性能的好壞將直接影響發(fā)動(dòng)機(jī)的散熱、功耗與噪聲狀況[1-2]。冷卻風(fēng)扇的主要性能參數(shù)為流量、靜壓、功率、靜壓效率與風(fēng)扇轉(zhuǎn)速的關(guān)系，風(fēng)扇設(shè)計(jì)是在保證一定的靜壓下，盡可能降低功率和噪聲，提高風(fēng)扇效率[2-4]。

文獻(xiàn)[5]運(yùn)用CFD計(jì)算與粒子圖像測(cè)速技術(shù)(PIV)對(duì)風(fēng)扇流場(chǎng)進(jìn)行研究，研究表明，CFD可成功地預(yù)測(cè)風(fēng)扇實(shí)際流場(chǎng)細(xì)節(jié)問題。文獻(xiàn)[7]中運(yùn)用多重參考系方法(Multiple Reference Frame,MRF)，對(duì)冷卻風(fēng)扇及旋轉(zhuǎn)域大小進(jìn)行了研究。文獻(xiàn)[7]利用CFD技術(shù)對(duì)發(fā)動(dòng)機(jī)軸流冷卻風(fēng)扇的氣動(dòng)性能進(jìn)行了研究，并探討了計(jì)算聲學(xué)(CAA)在噪聲控制中的應(yīng)用。文獻(xiàn)[8]基于LES大渦模型與FW-H方程，用瞬態(tài)計(jì)算方法對(duì)風(fēng)扇進(jìn)行CAA數(shù)值模擬。文獻(xiàn)[9]用穩(wěn)態(tài)計(jì)算方法與FW-H模型對(duì)風(fēng)扇進(jìn)行氣動(dòng)噪聲的計(jì)算，驗(yàn)證了穩(wěn)態(tài)計(jì)算風(fēng)扇噪聲的方法可以用于風(fēng)扇設(shè)計(jì)初期噪聲的計(jì)算；文獻(xiàn)[11]用寬頻噪聲模型研究了葉尖間隙、入口形狀以及翼形形狀對(duì)噪聲的影響。

本文建立了發(fā)動(dòng)機(jī)冷卻風(fēng)扇性能計(jì)算的CFD模型，介紹了風(fēng)扇氣動(dòng)性能與氣動(dòng)噪聲的計(jì)算方法、求解技術(shù)以及模型參數(shù)的選擇，計(jì)算并分析了發(fā)動(dòng)機(jī)冷卻風(fēng)扇的靜壓、功率、效率與流量的關(guān)系以及噪聲水平。將計(jì)算結(jié)果與實(shí)驗(yàn)值進(jìn)行對(duì)比，驗(yàn)證了計(jì)算模型的正確性。同時(shí)對(duì)風(fēng)扇內(nèi)部壓力場(chǎng)、速度場(chǎng)以及噪音分布進(jìn)行分析，分析了風(fēng)扇效率低下的原因與影響風(fēng)扇渦流噪音的主要因素，并提出了改進(jìn)方法。

1 發(fā)動(dòng)機(jī)冷卻風(fēng)扇性能試驗(yàn)

為獲得發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能參數(shù)，實(shí)驗(yàn)采用以《GB-T1236-2000工業(yè)通風(fēng)機(jī)用標(biāo)準(zhǔn)化風(fēng)道進(jìn)行性能試驗(yàn)》為標(biāo)準(zhǔn)建立的風(fēng)扇性能試驗(yàn)臺(tái)。由于發(fā)動(dòng)機(jī)冷卻風(fēng)扇無前、后導(dǎo)葉，同時(shí)壓力較低，因此采用管道進(jìn)口-自由出口型測(cè)試裝置[11]，如圖1所示。

圖1 風(fēng)扇氣動(dòng)性能試驗(yàn)臺(tái)示意圖

該裝置主要由可變供氣系統(tǒng)、管道、驅(qū)動(dòng)電機(jī)等組件組成。其中：錐形進(jìn)口1和流量加載板2組成可變供氣系統(tǒng)來調(diào)控進(jìn)口流量q(單位：m3/s)；整流格柵3、5用于穩(wěn)定氣流；U型壓力計(jì)9、10測(cè)試管道進(jìn)口處及監(jiān)測(cè)面位置氣流壓力，進(jìn)而得到風(fēng)扇靜壓力psp(單位：Pa)；驅(qū)動(dòng)電機(jī)8調(diào)控風(fēng)扇轉(zhuǎn)速n(單位：r/min)；轉(zhuǎn)速、扭矩傳感器7測(cè)量驅(qū)動(dòng)電機(jī)扭矩值T(單位：N·m)；定義風(fēng)扇軸功率P(單位：kW)和靜壓效率η如式(1)和式(2)所示[13]。

(1)

(2)

噪聲試驗(yàn)在風(fēng)扇性能試驗(yàn)室內(nèi)進(jìn)行，聲壓計(jì)水平放置在風(fēng)扇出口中心45°方向lm處。由于風(fēng)扇性能試驗(yàn)室基礎(chǔ)噪聲較大，測(cè)量前需要分別對(duì)環(huán)境背景噪聲和測(cè)量轉(zhuǎn)速下電機(jī)空轉(zhuǎn)時(shí)的噪聲進(jìn)行測(cè)量。

風(fēng)扇的氣動(dòng)性能主要通過風(fēng)扇的3條性能曲線來評(píng)價(jià)，即：流量-靜壓曲線、流量-功率曲線、流量-效率曲線[13]。實(shí)驗(yàn)中，如圖2所示。通過改變調(diào)速電機(jī)8的轉(zhuǎn)速與節(jié)流加載板2的過流面積，就能得到風(fēng)扇在不同轉(zhuǎn)速下的性能曲線。

2 發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能與氣動(dòng)噪聲計(jì)算

2.1 風(fēng)扇性能計(jì)算模型的建立

本文以雪龍公司一款代號(hào)為Z760WH風(fēng)扇為例進(jìn)行研究，如圖3所示，這是一款輪轂收縮風(fēng)扇，為充分利用安裝空間其輪轂收縮為薄板。工作轉(zhuǎn)速為2100RPM，風(fēng)扇外徑為760 mm，葉片數(shù)為11，葉片等間距分布。利用三維造型軟件對(duì)實(shí)例風(fēng)扇進(jìn)行三維建模，取Z軸為風(fēng)扇的旋轉(zhuǎn)軸，風(fēng)扇的旋轉(zhuǎn)方向?yàn)檎D(zhuǎn)方向，Z軸負(fù)方向?yàn)閬砹鞣较颉?/p>

圖2 風(fēng)扇氣動(dòng)性能曲線示意圖

參照?qǐng)D1所示試驗(yàn)臺(tái)進(jìn)行建模，如圖4所示。根據(jù)管道進(jìn)口-自由出口的流場(chǎng)特性及風(fēng)扇附近區(qū)域的流線分布規(guī)律，所建立Fluent模型包括：等尺寸進(jìn)口管道1、過渡區(qū)4和出口區(qū)5。在Fluent管道模型中，質(zhì)量入口邊界條件可任意給定質(zhì)量流量，省略圖1中結(jié)構(gòu)1、2；根據(jù)管道模型入口流場(chǎng)均勻性，省略掉圖1中結(jié)構(gòu)3、5；管道模型壓力監(jiān)測(cè)面2的靜壓值可以直接獲取，省略掉圖1中結(jié)構(gòu)9、10。

圖3 風(fēng)扇三維模型

圖4 風(fēng)管試驗(yàn)臺(tái)仿真模型

因?yàn)槿~尖渦、分離流等大多出現(xiàn)在葉尖，為捕捉風(fēng)扇周圍復(fù)雜旋渦流動(dòng)，在區(qū)域4建立大小合適的旋轉(zhuǎn)區(qū)。為了更好地捕捉到葉尖復(fù)雜的流場(chǎng)信息，需將葉尖網(wǎng)格與旋轉(zhuǎn)區(qū)網(wǎng)格進(jìn)行細(xì)化[15-16]，如圖5所示；為保證計(jì)算精度并減少計(jì)算周期，根據(jù)流線分布規(guī)律，所建立模型中：區(qū)域1、5采用六面體網(wǎng)格，網(wǎng)格數(shù)約為80萬；區(qū)域5采用四面體網(wǎng)格數(shù)約為100萬。

圖5 風(fēng)管試驗(yàn)臺(tái)計(jì)算模型網(wǎng)格示意圖

2.2 模型參數(shù)及邊界條件設(shè)置

(1)模型參數(shù)設(shè)置

風(fēng)扇是典型的旋轉(zhuǎn)機(jī)械，是非定常流動(dòng)，MRF模型(多重參考系模型)可將非定常旋轉(zhuǎn)問題轉(zhuǎn)化為定常旋轉(zhuǎn)問題，即在風(fēng)扇周圍的旋轉(zhuǎn)區(qū)采用動(dòng)參考系，其他區(qū)域采用靜參考系，在保證計(jì)算精度的同時(shí)可大大節(jié)約計(jì)算資源與計(jì)算時(shí)間；湍流模型選用RNG湍流模型，在模擬復(fù)雜流動(dòng)時(shí)應(yīng)用較廣泛；選用無滑移的標(biāo)準(zhǔn)壁面函數(shù)；風(fēng)扇的流場(chǎng)屬于低馬赫數(shù)流動(dòng)，可以將空氣近似為粘性不可壓縮介質(zhì)，因此選擇壓力基求解器；采用SIMPLE算法求解壓力-速度耦合方程；采用標(biāo)準(zhǔn)格式離散求解連續(xù)性方程；采用一階迎風(fēng)格式離散求解動(dòng)量方程；采用二階迎風(fēng)格式離散求解湍流動(dòng)能與湍流耗散率。本模型的收斂準(zhǔn)則定為所有監(jiān)測(cè)項(xiàng)殘差小于1×10-3,且殘差隨著迭代步的推進(jìn)數(shù)值保持平穩(wěn)[2,17-18]。

寬帶噪聲模型也稱穩(wěn)態(tài)噪聲模，它是在穩(wěn)態(tài)結(jié)果的基礎(chǔ)上進(jìn)行模擬，只能計(jì)算出聲功率級(jí)，不能得到噪聲的頻譜，但消耗的計(jì)算資源僅有非穩(wěn)態(tài)模型的十分之一至百分之一，這是一個(gè)快速評(píng)估設(shè)計(jì)是否需要改進(jìn)的非常實(shí)用的模型[19]，對(duì)于工程問題非常有用。因此，本文的噪聲計(jì)算均采用穩(wěn)態(tài)噪聲模型，在噪聲模型中將Broadband Noise Resource 打開即可，噪聲源選擇風(fēng)扇葉片。

(2)邊界條件設(shè)置

入口選為質(zhì)量入口邊界，可根據(jù)需要設(shè)定不同入口流量，從而得到風(fēng)扇在不同流量下的性能曲線；出口選為壓力出口邊界，因?yàn)槌隹谂c大氣相連，故設(shè)定其靜壓值為0。風(fēng)扇表面選用旋轉(zhuǎn)壁面邊界并根據(jù)需要設(shè)定不同轉(zhuǎn)速，其余表面為靜止壁面邊界；六面體網(wǎng)格與四面體網(wǎng)格交界面設(shè)定為interface來實(shí)現(xiàn)數(shù)據(jù)傳遞[2,17-18]。

2.3 網(wǎng)格無關(guān)性檢驗(yàn)

為消除網(wǎng)格數(shù)量、質(zhì)量帶來的誤差，需進(jìn)行網(wǎng)格無關(guān)性檢驗(yàn)[18]。在保證相同網(wǎng)格質(zhì)量的前提下選取4種密度網(wǎng)格進(jìn)行仿真計(jì)算，結(jié)果如表1所示，當(dāng)網(wǎng)格數(shù)控制在一定數(shù)量時(shí)(約160萬)靜壓計(jì)算結(jié)果變化不大，僅相差0.46%≤2%，已滿足網(wǎng)格無關(guān)性要求，但噪聲值在網(wǎng)格達(dá)到210萬時(shí)才趨于穩(wěn)定。考慮計(jì)算時(shí)間成本(為節(jié)約計(jì)算資源)，選用210萬網(wǎng)格進(jìn)行計(jì)算。

表1 網(wǎng)格無關(guān)性檢驗(yàn)

3 計(jì)算結(jié)果

3.1 計(jì)算結(jié)果與試驗(yàn)結(jié)果對(duì)比

風(fēng)扇氣動(dòng)性能的計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比如圖6所示，在工作流量4 m3/s～6 m3/s區(qū)間，靜壓、功率與效率的最大誤差分別為5.57%、4.78%與5.69%，誤差均在6%以內(nèi)，驗(yàn)證了該模型與數(shù)值計(jì)算方法的正確性。選擇900～2 400 rpm中6個(gè)轉(zhuǎn)速進(jìn)行計(jì)算，噪聲計(jì)算值和實(shí)驗(yàn)值的對(duì)比見圖7。從圖7可以看到，計(jì)算值同實(shí)驗(yàn)值的趨勢(shì)較為一致，噪聲差值在16 dB左右，這是因?yàn)閷拵г肼暷Ｐ秃雎粤藢?duì)時(shí)間的導(dǎo)數(shù)項(xiàng)，只計(jì)算了穩(wěn)態(tài)壓力變化過程產(chǎn)生的氣動(dòng)噪聲。事實(shí)上，風(fēng)扇產(chǎn)生的氣動(dòng)噪聲還包括非穩(wěn)態(tài)壓力波動(dòng)所產(chǎn)生的噪聲。計(jì)算結(jié)果表明，穩(wěn)態(tài)噪聲模型能夠用于對(duì)不同風(fēng)扇噪聲的對(duì)比。

3.2 計(jì)算結(jié)果分析

冷卻風(fēng)扇表面靜壓力分布情況如圖8、圖9所示?？梢钥吹剑~片壓力面的靜壓分布整體較均勻，葉片中上部分壓力較大，表明這部分是風(fēng)扇做功的主要區(qū)域，靠近葉根部分壓力較小，因?yàn)檫@一區(qū)域的線速度較小。葉片頂端尤其是風(fēng)扇圓環(huán)處靜壓值較小，因?yàn)閴毫γ婧臀γ娴膲翰顚?dǎo)致氣體產(chǎn)生回流，降低了葉尖部分壓力面的靜壓，是風(fēng)扇效率損失的主要來源之一，減小葉尖間隙可降低葉尖回流損失，提高風(fēng)扇效率[1]。

圖6 風(fēng)扇性能曲線對(duì)比圖

圖7 風(fēng)扇噪聲對(duì)比圖

冷卻風(fēng)扇葉片表面動(dòng)壓分布規(guī)律圖10、圖11所示。風(fēng)扇葉片壓力面后緣靠近葉尖的部分，有小面積區(qū)域出現(xiàn)動(dòng)壓最大值，風(fēng)扇葉片吸力面后緣葉片上部有較大面積區(qū)域動(dòng)壓較大，在此區(qū)域內(nèi)出現(xiàn)動(dòng)壓的最大值。風(fēng)扇葉片吸力面前緣、壓力面前緣葉頂和葉根部分動(dòng)壓值較小。

圖8 冷卻風(fēng)扇壓力面靜壓分布圖

圖9 冷卻風(fēng)扇吸力面靜壓分布圖

圖10 冷卻風(fēng)扇壓力面動(dòng)壓分布圖

冷卻風(fēng)扇葉片表面總壓分布規(guī)律如圖12、圖13所示?？倝簽閯?dòng)壓和靜壓之和，所以綜合靜壓和動(dòng)壓的分布規(guī)律，可以得出總壓的分布規(guī)律。冷卻風(fēng)扇壓力面后緣葉片上部有較大面積區(qū)域總壓較大，靠近葉尖部分出現(xiàn)總壓最大值。冷卻風(fēng)扇吸力面前緣葉片上部總壓較小，出現(xiàn)負(fù)值，在吸力面后緣葉尖部分總壓值較大。

圖11 冷卻風(fēng)扇吸力面動(dòng)壓分布圖

圖12 冷卻風(fēng)扇壓力面全壓分布圖

圖13 冷卻風(fēng)扇吸力面全壓分布圖

冷卻風(fēng)扇葉片表面絕對(duì)速度分布和相對(duì)速度分布分別如圖14、圖15所示。由于動(dòng)壓為速度平方的函數(shù)，所以，冷卻風(fēng)扇表面速度大小分布與動(dòng)壓相同，但是速度是矢量，在壓力面，速度方向與來流方向夾角大于90度，而在吸力面，速度方向與來流方向夾角小于90度。

從圖16、圖17可以看到，風(fēng)扇表面的噪聲主要集中在吸力面的前緣以及葉尖前后緣處，尤其是葉尖處。其原因是流體粘性和壓力脈沖造成葉片邊緣形成分離流、葉尖渦以及脫落渦等，而這種渦流噪聲是風(fēng)扇氣動(dòng)噪聲的主要組成因素之一[18-19]。當(dāng)葉尖間隙很小時(shí)，其對(duì)葉尖渦流的控制是有利的，因?yàn)槿~尖間隙減小到一定程度時(shí)，抑制了葉尖二次渦流的形成。但過小的間隙無法滿足風(fēng)扇裝配與運(yùn)轉(zhuǎn)要求，葉尖間隙應(yīng)該控制在風(fēng)扇直徑的2%以內(nèi)。

圖14 冷卻風(fēng)扇表面絕對(duì)速度分布圖

圖15 冷卻風(fēng)扇表面相對(duì)速度分布圖

圖16 冷卻風(fēng)扇壓力面噪聲分布圖

圖17 冷卻風(fēng)扇吸力面噪聲分布圖

4 結(jié)論

(1)建立了與標(biāo)準(zhǔn)試驗(yàn)設(shè)備一致的風(fēng)管試驗(yàn)臺(tái)計(jì)算模型、介紹了風(fēng)扇試驗(yàn)臺(tái)的建模方法，計(jì)算方法以及相關(guān)模型參數(shù)的選擇，并將計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行了對(duì)比分析，驗(yàn)證了模型的正確性。所建模型能夠指導(dǎo)風(fēng)扇改進(jìn)設(shè)計(jì)。

(2)通過流場(chǎng)診斷研究分析了風(fēng)扇壓力面和吸力面的靜壓、動(dòng)壓、全壓與速度分布規(guī)律，分析了風(fēng)扇效率低下的主要原因，并提出應(yīng)盡量縮小葉尖間隙以降低葉尖的回流損失，提高風(fēng)扇效率。

(3)通過觀察風(fēng)扇表面的噪聲分布，分析噪聲形成的主要原因，并提出通過減小葉尖間隙的方法可有效控制葉尖渦流，達(dá)到降噪的目的。

[1]莫偉標(biāo).發(fā)動(dòng)機(jī)環(huán)形冷卻風(fēng)扇結(jié)構(gòu)與參數(shù)對(duì)其性能影響的研究[D].廣州:華南理工大學(xué),2015.

[2]上官文斌,吳敏.發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能的計(jì)算方法[J].汽車工程,2010,32(9):799-802.

[3]上官文斌,王益有.基于無量綱性能曲線的發(fā)動(dòng)機(jī)冷卻風(fēng)扇設(shè)計(jì)方法[J].汽車工程,2010,32(5):455-458.

[4]耿麗珍,袁兆成,李傳兵,等.轎車發(fā)動(dòng)機(jī)冷卻風(fēng)扇CFD仿真分析及降噪研究[J].汽車工程,2009,31(7):664-668.

[5]KOHRI I,KOBAYASHI Y,MATSUSHIMA Y.Prediction of the performance of the engine cooling fan with CFD simulation[J].Journal of Thermal Analysis & Calorimetry,2010,3(1):508-522.

[6]KOBAYASHI Y,KOHRI I,MATSUSHIMA Y.Study of influence of MRF method on the prediction of the engine cooling fan performance[J].British Journal of Dermatology,2011,163(4):670-671.

[7]張紅輝.發(fā)動(dòng)機(jī)軸流冷卻風(fēng)扇低噪聲氣動(dòng)性能分析與控制研究[D].重慶:重慶大學(xué),2002.

[8]唐永偉.汽車發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能和聲學(xué)性能研究[D].長(zhǎng)春:吉林大學(xué),2007.

[9]王益友.基于CFD的發(fā)動(dòng)機(jī)冷卻風(fēng)扇噪聲計(jì)算[J].機(jī)械工程師,2009,2009(11):11-12

[10]LI HUIBIN.Turbocharger noise prediction using broadband noise source model[J].Journal of Beijing Institude of Technology,2010,9(3):312-317

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

[12]昌澤舟.軸流式通風(fēng)機(jī)實(shí)用技術(shù)[M].北京:機(jī)械工業(yè)出版社,2005.

[13]BLEIER F P.Fan Handbook:Selection,Application,and Design[M].New York:McGraw-Hill,1997.

[14]JIN-HYUK,KIM,KWANG-YONG,KIM.Axial-flow ventilation fan design through muti-objective optimization to enhance aerodynamic performance[J].Jouranl of Fluids Engineering,2011，133(10):1-12.

[15]SZU HSIEN LIU，RONG FUNG HUANG.Computational and experimental investigations of performance curve of an axial flow fan using downstream flow resistance method[J].Experimental Thermal and Fluid Science,2010,(34):827-837.

[16]王福軍.計(jì)算流體力學(xué)分析-CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004:113-125.

[17]AYHAN NAZMI LIKAN,ERKAN AYDER.Influence of dihedral stacking on the performance of an axial fan[J].Engineering Applications of Computational Fluid Mechanics,2015,8(4):518-529.

[18]NASHIMOTO A,AKUTO T,NAGASE Y,et al.Aerodynamic noise reduction by use of a cooling fan with winglets[C].SAE World Congress & Exhibition,Detroit,MI,USA:2003,112(6):622-630.

[19]NASHIMOTO A,FUJISAWA N,AKUTO T,et al.Measurements of aerodynamic noise and wake flow field in a cooling fan with winglets[J].Journal of Visualization,2004,7(1):85-92.

(責(zé)任編輯：劉劃 英文審校：趙歡)

Calculation and analysis on aerodynamic performance and aerodynamic noise of cooling fan in axial-flow engine

ZHAO Li-jie1,WANG Xin-ling1,SHU Heng-tao2,DUAN Yao-long3,ZHU Fei-fei3

(1.School of Mechanical and Electrical Engineering,Shenyang Aerospace University,Shenyang 110136，China;2.School of Mechanical and Automotive Engineering,South China University of Technology,Guangzhou 510641，China;3.Technology Research and Development Centre,Ningbo Xuelong Group Co.,Ltd.,Ningbo 315806，China)

A test-bench for aerodynamic performance of cooling fan in axial-flow engine was set up based on a type C device in GB-T1236-2000(Standardized duct for industrial ventilators conducts performance tests)using a commercial software(Fluent).A method calculating aerodynamic performance of the fan and setting relevant parameters of the model was introduced.Taking a certain fan as an example,noise level and relationships between flow rates and other parameters(i.e.static pressure,power,and efficiency)were calculated.The calculated results were compared with the measured data.The comparisons show that the calculation errors of static pressure,power and efficiency are less than 6%,Change trend of the calculated parameters is consistent with that of the measured values despite the high noise error.The calculation method can be used to compare noise of different cooling fans.The distribution of pressure and velocity on surface of the fan and main causes for noise were further analyzed by studying flow field inside the fan.The results indicate that tip clearance is the main cause for low efficiency of the fan,and tip vortex led to the vortex noise of fan.If tip clearance is reduced,cooling efficiency will be improved and the tip vortex will also be controlled effectively to reduce noise.

cooling fans;CFD;aerodynamic performance;aerodynamic noise

2016-11-18

趙立杰(1964-)，男，遼寧沈陽人，副教授，主要研究方向：智能材料與智能結(jié)構(gòu)，E-mail:zhaolj@sau.edu.cn。

2095-1248(2017)01-0050-07

TK414.2

A

10.3969/j.issn.2095-1248.2017.01.008