離心風(fēng)機(jī)流場(chǎng)大渦模擬與管中噪聲數(shù)值預(yù)報(bào)

李亞，劉忠族 ，許影博 ，張楠

1船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，江蘇無錫214082

2中國船舶科學(xué)研究中心，江蘇無錫214082

3中國船舶科學(xué)研究中心水動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室，江蘇無錫214082

0 引 言

離心風(fēng)機(jī)管道中的噪聲一般來源于葉輪輻射聲、蝸殼和管道壁面邊界層噪聲、尾流場(chǎng)的湍流噪聲、蝸殼和管道振動(dòng)引起的噪聲。噪聲在離心風(fēng)機(jī)管道中傳播時(shí)，還會(huì)受到管道噪聲模態(tài)的影響。另外，離心風(fēng)機(jī)運(yùn)行時(shí)，電機(jī)噪聲和環(huán)境噪聲也會(huì)傳入管道。葉輪輻射聲和壁面邊界層噪聲屬于偶極子類型，輻射效率高于四極子類型的湍流噪聲，離心風(fēng)機(jī)振動(dòng)噪聲相比氣動(dòng)聲小很多［1］。

將計(jì)算流體動(dòng)力學(xué)（CFD）與計(jì)算聲學(xué)相結(jié)合，可以進(jìn)行氣動(dòng)噪聲預(yù)報(bào)。Kato等［2］運(yùn)用大渦模擬（Large Eddy Simulation，LES）方法得到葉片表面的脈動(dòng)壓力，采用Lighthill聲類比理論計(jì)算遠(yuǎn)場(chǎng)輻射聲。Younsi等［3］采用滑移網(wǎng)格方式及SST k-ω湍流模型計(jì)算了離心風(fēng)機(jī)內(nèi)部流場(chǎng)，利用非定常流場(chǎng)參數(shù)作為FW-H方程［4］輸入項(xiàng)預(yù)報(bào)了輻射聲。陳敏等［5］運(yùn)用多塊混合網(wǎng)格和動(dòng)網(wǎng)格技術(shù)對(duì)螺旋槳的水動(dòng)力性能進(jìn)行大渦模擬，并結(jié)合FW-H方程預(yù)報(bào)了螺旋槳的低頻流噪聲。在上述噪聲預(yù)報(bào)中，采用的是自由空間格林函數(shù)，其必須假設(shè)離心風(fēng)機(jī)的幾何尺寸滿足聲學(xué)緊致結(jié)構(gòu)特征的要求。實(shí)際上，由于風(fēng)機(jī)蝸殼和葉輪壁面的存在，機(jī)殼內(nèi)部聲源激發(fā)的噪聲會(huì)在葉輪和蝸殼內(nèi)壁面上發(fā)生多次反射及散射，最終通過風(fēng)機(jī)進(jìn)、出口向外輻射。

早在 1974 年，Moreland［6］就通過實(shí)驗(yàn)證實(shí)了蝸殼腔體結(jié)構(gòu)對(duì)聲波反射及散射會(huì)導(dǎo)致部分頻段向外輻射的聲功率大大增強(qiáng)。劉厚林等［7］采用大渦模擬方法獲得了離心泵內(nèi)部瞬態(tài)流場(chǎng)，在計(jì)算葉片表面偶極子輻射聲場(chǎng)時(shí)，將葉片分為10個(gè)部分，分別積分得到3個(gè)方向的時(shí)域力，以此作為葉片旋轉(zhuǎn)偶極子聲源；以進(jìn)、出口為吸聲邊界，其他壁面作為全反射壁面條件，聲場(chǎng)計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合良好。

離心風(fēng)機(jī)形狀復(fù)雜，過去流場(chǎng)模擬多采用非結(jié)構(gòu)化網(wǎng)格，本文擬采用ICEM-CFD軟件繪制結(jié)構(gòu)化網(wǎng)格，使用Fluent軟件模擬流場(chǎng)，通過得到的聲源預(yù)報(bào)葉輪輻射聲和壁面邊界層這2種途徑下的管道噪聲，并將兩者疊加結(jié)果與實(shí)際結(jié)果進(jìn)行比較分析。

1 基于大渦模擬的風(fēng)機(jī)流場(chǎng)模擬

1.1 模擬對(duì)象

計(jì)算對(duì)象為某型離心風(fēng)機(jī)，如圖1所示。離心風(fēng)機(jī)重約117 kg，高約1.25 m；葉輪為后向葉輪，包含10個(gè)葉片，葉輪直徑為503 mm；轉(zhuǎn)速為1 470 r/min。采用管道進(jìn)口、管道出口方式［8］對(duì)離心風(fēng)機(jī)進(jìn)行實(shí)驗(yàn)。進(jìn)口管道與出口管道直徑均為267 mm，進(jìn)口管道長約12.3 m，出口管道長約11.6 m。測(cè)量情況如圖2所示，其中傳聲器安裝鼻錐，迎著來流，氣體從鼻錐光順地流過，可減少湍流脈動(dòng)的影響。

圖1 離心風(fēng)機(jī)與葉輪Fig.1 Centrifugal fan and impeller

圖2 離心風(fēng)機(jī)噪聲實(shí)驗(yàn)Fig.2 Noise measurement of the centrifugal fan

1.2 網(wǎng)格劃分

根據(jù)實(shí)驗(yàn)中離心風(fēng)機(jī)管道的連接情況，建立了整個(gè)流域的幾何模型。流域含旋轉(zhuǎn)葉輪，因此采用滑移網(wǎng)格方法模擬風(fēng)機(jī)葉片的旋轉(zhuǎn)運(yùn)動(dòng)：在風(fēng)機(jī)葉片所在區(qū)域劃分出一個(gè)小圓柱體區(qū)域（圖3），部分管道伸入葉輪內(nèi)部，導(dǎo)致交界面有階梯。小圓柱體內(nèi)包含葉片，此區(qū)域作為內(nèi)域，為轉(zhuǎn)動(dòng)部分；其他區(qū)域作為外域，為靜止部分。靜止部分和轉(zhuǎn)動(dòng)部分設(shè)置交界面并進(jìn)行數(shù)據(jù)交換，以此來保證各物理量守恒。

采用結(jié)構(gòu)化網(wǎng)格形式劃分網(wǎng)格，內(nèi)域網(wǎng)格如圖4所示；外域網(wǎng)格存在由方變圓段，無法保證結(jié)構(gòu)化網(wǎng)格質(zhì)量，因此也分為2部分，網(wǎng)格劃分如圖5所示。

圖3 交界面劃分示意圖Fig.3 Schematic diagram of interface division

圖4 風(fēng)機(jī)內(nèi)域結(jié)構(gòu)化網(wǎng)格Fig.4 Structured grid of internal zone of the centrifugal fan

圖5 風(fēng)機(jī)外域結(jié)構(gòu)化網(wǎng)格Fig.5 Structured grid of exterior zone of the centrifugal fan

在大渦模擬之前，采用劃分方式相同而網(wǎng)格數(shù)量不同的多套網(wǎng)格（第1層網(wǎng)格厚度保持不變）進(jìn)行穩(wěn)態(tài)流場(chǎng)模擬，以出口流量為判據(jù)進(jìn)行網(wǎng)格收斂性分析，考慮到大渦模擬網(wǎng)格數(shù)量要求以及單機(jī)計(jì)算能力，選取了其中1套網(wǎng)格：內(nèi)域網(wǎng)格數(shù)為 5.19×106，蝸殼與管道網(wǎng)格數(shù)為 2.34×106，出口管道網(wǎng)格數(shù)為2.7×105；蝸殼第1層網(wǎng)格尺度設(shè)置為0.6 mm，葉片為0.3 mm，出口段為0.3 mm。在運(yùn)算穩(wěn)定后，無因次壁面法向高度 y+的范圍為5～35，滿足計(jì)算要求。

1.3 大渦模擬計(jì)算風(fēng)機(jī)內(nèi)部流場(chǎng)

湍流數(shù)值模擬方法可以分為直接數(shù)值模擬（DNS）方法和非直接數(shù)值模擬方法。大渦模擬既有直接模擬又有非直接模擬特征。它將流動(dòng)分解為大、小2種尺度的湍流，對(duì)大尺度湍流進(jìn)行直接模擬，而對(duì)小尺度湍流做近似處理。由于小尺度渦不直接依賴邊界條件，而且一般具有各向同性性質(zhì)，可以采用合適的湍流模型模擬。

大渦模擬計(jì)算對(duì)初始條件要求較高。首先采用RNG k-ε湍流模式下的結(jié)果作為初始條件，然后采用壁面自適應(yīng)局部渦粘性（Wall-Adapting Local Eddy-viscosity，WALE）亞格子模型，進(jìn)行大渦模擬計(jì)算，其中內(nèi)域采用運(yùn)動(dòng)參考系。在計(jì)算穩(wěn)定后，內(nèi)域采用網(wǎng)格運(yùn)動(dòng)形式進(jìn)行模擬，在運(yùn)行一段時(shí)間后進(jìn)行數(shù)據(jù)采集。數(shù)據(jù)采集的起始步數(shù)為第4 200步，對(duì)應(yīng)時(shí)刻為0.105 s，采樣間隔時(shí)間為2.5×10-5s。

葉輪旋轉(zhuǎn)1/10周（即一個(gè)葉片旋轉(zhuǎn)到上一個(gè)葉片位置）需要163步，1/10周再分為3個(gè)部分，分別繪制截面的壓力分布云圖，如圖6所示。從圖中可以看出，在計(jì)算過程中，葉片在不斷旋轉(zhuǎn)。在第4 212步與第4 374步，下一個(gè)葉片來到了同樣位置，但葉輪附近的壓力分布并不一致，這說明受渦流影響，葉輪所受壓力并不是完全重復(fù)的。在第4 212步時(shí)，葉片離渦舌部位很近，氣體受到擠壓導(dǎo)致壓強(qiáng)增大，葉片轉(zhuǎn)過之后，壓力下降。

圖6 葉輪旋轉(zhuǎn)1/10周時(shí)流場(chǎng)絕對(duì)壓力的變化Fig.6 Absolute pressure changes of flow field in the 1/10 cycle of impeller

2 離心風(fēng)機(jī)管中噪聲預(yù)報(bào)

離心風(fēng)機(jī)管中噪聲主要來源于旋轉(zhuǎn)葉片的氣動(dòng)噪聲，另外蝸殼表面壓力脈動(dòng)也會(huì)輻射噪聲，這2個(gè)激勵(lì)源均由大渦模擬計(jì)算得到。噪聲在管中傳播時(shí)還應(yīng)考慮管道的影響。

2.1 葉輪輻射聲計(jì)算

2.1.1 旋轉(zhuǎn)噪聲理論

首先，分析旋轉(zhuǎn)機(jī)械在自由場(chǎng)的輻射聲。葉輪工作時(shí)，葉片面元上會(huì)產(chǎn)生脈動(dòng)壓力，特別是在葉輪出口附近存在渦脫落時(shí)，壓力脈動(dòng)更明顯，這些壓力脈動(dòng)會(huì)產(chǎn)生噪聲。而且葉輪在不斷旋轉(zhuǎn)，因此葉輪噪聲相當(dāng)于旋轉(zhuǎn)偶極子輻射聲。

對(duì)于其中一塊面元，F(xiàn)={Fx，F(xiàn)y，F(xiàn)z} ，可以得到各分力產(chǎn)生的總聲壓［9］：

2.1.2 聲學(xué)有限元計(jì)算原理

葉輪在蝸殼中旋轉(zhuǎn)時(shí)，其噪聲傳播會(huì)受到管道聲模態(tài)的影響。對(duì)于有限空間中的輻射聲問題，一般采用邊界元法或有限元法。邊界元法盡管可以將維數(shù)降低，但求解速度慢［10］。聲學(xué)有限元法的一般思路是：聲場(chǎng)滿足微分方程（含邊界條件），可以構(gòu)建與之等價(jià)的泛函，整個(gè)聲場(chǎng)泛函可以用網(wǎng)格單元的形函數(shù)表達(dá)，對(duì)其求極值，即其變分為0，可以構(gòu)建整體的矩陣方程。

在均勻介質(zhì)、非粘性和絕熱狀態(tài)下，流體內(nèi)的聲學(xué)波動(dòng)方程為

式中：ρ0為靜態(tài)密度；c=ω/k。令聲壓p'=p(x，y，z)ejωt，聲源q'=q(x，y，z)ejωt，代入式（2），得到Helmholtz方程：

采用加權(quán)余量法，可以得到聲場(chǎng)滿足的矩陣方程［11］：

式中：K為聲學(xué)剛度矩陣；C為聲學(xué)阻尼矩陣；Μ為質(zhì)量矩陣；p為聲壓向量；Q為聲源向量；V為速度向量；P為輸入聲壓向量。這些量均由計(jì)算域單元的形函數(shù)表達(dá)。

2.1.3 葉輪輻射聲預(yù)報(bào)

采用大渦模擬計(jì)算葉輪旋轉(zhuǎn)1周（從4200步到5831步）時(shí)的流場(chǎng)，輸出了葉輪表面在每一步的受力情況，文件格式為計(jì)算流體力學(xué)通用符號(hào)系統(tǒng)數(shù)據(jù)（CFD General Notation System Data，CGNS）。在Virtual.Lab軟件中導(dǎo)入相應(yīng)文件，由于采用有限體積法計(jì)算得到的脈動(dòng)壓力結(jié)果存儲(chǔ)在網(wǎng)格形心上，需要將形心結(jié)果傳遞到網(wǎng)格各節(jié)點(diǎn)［10］。經(jīng)過傳遞，在某一時(shí)刻的壓力云圖如圖7所示。

圖7 第4 265步時(shí)葉片上的壓力云圖Fig.7 Pressure contours on impeller at the 4 265thstep

在Virtual.Lab軟件中進(jìn)行扇聲源計(jì)算時(shí)，扇聲源被當(dāng)作一類專門的計(jì)算對(duì)象。將其葉輪等效為一個(gè)個(gè)旋轉(zhuǎn)偶極子，即將葉片表面進(jìn)行劃分，得到的每個(gè)小區(qū)域相當(dāng)于1個(gè)聲源。聲源數(shù)與計(jì)算的最高頻率有關(guān)。圖8（a）所示的彩色色塊為在設(shè)置計(jì)算頻率后葉輪表面區(qū)域自動(dòng)劃分的結(jié)果。然后，添加扇聲源的邊界條件，旋轉(zhuǎn)速度設(shè)置為1 470 r/min，傅里葉變換加窗選擇漢寧（Hanning）窗。圖8（b）為葉輪旋轉(zhuǎn)參數(shù)設(shè)置結(jié)果，由圖可以看出旋轉(zhuǎn)方向以及各個(gè)面積塊的離散力分布。

圖8 在Virtual.Lab軟件中的葉輪設(shè)置Fig.8 Impeller setting in software Virtual.Lab

通常在劃分聲學(xué)有限元網(wǎng)格時(shí)，1個(gè)波長不少于6個(gè)單元，但由于外域尺度大，網(wǎng)格數(shù)量較多，計(jì)算時(shí)間長且對(duì)計(jì)算機(jī)硬件要求高，所以采用自適應(yīng)階次有限元（FEM Adaptive Order，F(xiàn)EMAO）算法，其主要原理是采用高階多項(xiàng)式表達(dá)單元內(nèi)的參數(shù)，比較粗糙的網(wǎng)格也能得到同樣精度的計(jì)算結(jié)果［12］。

原來蝸殼入口處有部分伸到葉輪中，如圖9（a）所示，這是由于蝸殼厚度為2 mm，而且有卷邊。單個(gè)網(wǎng)格尺寸比較大，難以擬合其邊界，導(dǎo)致網(wǎng)格生成質(zhì)量不高，所以將卷起部分進(jìn)行了簡化，如圖 9（b）所示。圖9（c）為劃分好的網(wǎng)格，其中網(wǎng)格均為四面體，單元最大邊長為64 mm，網(wǎng)格數(shù)為12.9×104。

圖9 聲場(chǎng)區(qū)域網(wǎng)格劃分過程Fig.9 Meshing process of sound field

進(jìn)風(fēng)管口與出風(fēng)管口暴露在空氣中，所以管口可作為無反射邊界條件。在進(jìn)、出風(fēng)管口處的單元組上定義聲阻抗來模擬無反射邊界，將阻抗設(shè)置為416.5 N·s/m3。最大計(jì)算頻率為6 kHz，頻率間隔為24.5 Hz。求解場(chǎng)點(diǎn)坐標(biāo)分別為：進(jìn)風(fēng)管（0 mm，0 mm，8 860 mm），出風(fēng)管（-8 150 mm，368.7 mm，175.5 mm），這2點(diǎn)與實(shí)際測(cè)量點(diǎn)位置一樣。圖10為1/3倍頻程噪聲預(yù)報(bào)結(jié)果。由圖可看出，噪聲最大值為85 dB，進(jìn)、出風(fēng)管測(cè)點(diǎn)的噪聲在高頻時(shí)比較接近。

圖10 葉輪輻射聲預(yù)報(bào)結(jié)果Fig.10 Prediction results of impeller noise

2.2 蝸殼脈動(dòng)壓力輻射聲計(jì)算

2.2.1 有固體壁面時(shí)湍流輻射聲

當(dāng)湍流區(qū)內(nèi)存在剛性表面時(shí)，湍流區(qū)內(nèi)的密度變化量可用Curle方程表示［13］：

2.2.2 蝸殼壁面脈動(dòng)壓力輻射聲預(yù)報(bào)

在大渦模擬時(shí)，蝸殼表面力無法同時(shí)輸出，可以采用另外一種軟件操作方法：即輸出計(jì)算過程中整體域的解，文件格式仍為CGNS。由于輸出文件包含了整個(gè)流體域的結(jié)果，自然也包括蝸殼與管道界面的結(jié)果。采用該方法后，每個(gè)文件均較大，因此步長增加1倍，為5×10-5s。數(shù)值計(jì)算從第4 200步開始，共計(jì)算816步，即葉輪旋轉(zhuǎn)1周。

管道內(nèi)表面是光順的，湍流強(qiáng)度低，僅產(chǎn)生邊界層的噪聲，且該噪聲的量級(jí)很低；而蝸殼內(nèi)表面有大量的渦分離，所以本次只對(duì)蝸殼內(nèi)表面部分進(jìn)行積分計(jì)算。

數(shù)據(jù)存儲(chǔ)于體單元中，實(shí)際計(jì)算只需要蝸殼表面網(wǎng)格上的數(shù)據(jù)，所以要進(jìn)行數(shù)據(jù)傳遞。流場(chǎng)計(jì)算時(shí)的網(wǎng)格較密，而殼體單元組的網(wǎng)格較稀疏，本文采用無能量損失算法（Conservative Maximum Distance）將力從流體網(wǎng)格映射到蝸殼殼體網(wǎng)格。在CFD時(shí)域計(jì)算時(shí)，考慮到聲學(xué)分析需要，對(duì)其進(jìn)行快速傅里葉變換，圖11為150 Hz時(shí)的脈動(dòng)壓力云圖。

圖11 150 Hz時(shí)的蝸殼壁面脈動(dòng)壓力云圖Fig.11 Pressure fluctuation contours of volute wall at 150 Hz

同樣，在入口單元和出口單元組分別定義無反射邊界條件，將完成數(shù)據(jù)傳遞后的結(jié)果作為邊界條件，物理類型為表面偶極子（Surface Dipole）。最大計(jì)算頻率6 kHz，頻率間隔為25 Hz。1/3倍頻程噪聲計(jì)算結(jié)果如圖12所示。由圖12的噪聲預(yù)報(bào)結(jié)果可以看出，不同位置處的輻射聲在低頻時(shí)有差別、在高頻時(shí)趨于一致。

圖12 蝸殼壁面脈動(dòng)壓力輻射聲預(yù)報(bào)結(jié)果Fig.12 Prediction results of volute wall pressure fluctuation noise

3 風(fēng)機(jī)管中噪聲計(jì)算結(jié)果

3.1 進(jìn)風(fēng)管中葉輪輻射聲與蝸殼壁面脈動(dòng)壓力輻射聲對(duì)比

圖13為前述進(jìn)風(fēng)管中葉輪輻射聲與蝸殼壁面脈動(dòng)壓力輻射聲的對(duì)比結(jié)果。由圖可以看出，進(jìn)風(fēng)管內(nèi)蝸殼流噪聲在低頻部分高于葉輪輻射聲，而在高頻時(shí)又低于葉輪輻射聲。這表明葉輪和蝸殼壁面的非定常力是離心風(fēng)機(jī)最主要的氣動(dòng)噪聲源，證實(shí)了文獻(xiàn)［14］中的結(jié)論。

圖13 進(jìn)風(fēng)管中葉輪輻射聲與蝸殼脈動(dòng)壓力輻射聲的對(duì)比結(jié)果Fig.13 Contrast between impeller noise and volute pressure fluctuation noise

一般葉片噪聲較大，但本例中離心風(fēng)機(jī)葉輪與蝸殼距離非常近，在蝸舌部位距離僅有28.5 mm，表明分離的渦（脫離渦將造成較大的脈動(dòng)，從而產(chǎn)生較大輻射聲）尾流還會(huì)直接作用于蝸殼表面。另外，從聲源作用面積來講，10個(gè)葉片的總面積（含葉面和葉背2面）為0.36 m2，而蝸殼內(nèi)表面面積為2.2 m2，后者遠(yuǎn)大于前者。這些因素導(dǎo)致蝸殼內(nèi)流噪聲較大。但高頻部分葉輪輻射聲又增加了，這部分跟直接渦發(fā)放有關(guān)。在大于臨界雷諾數(shù)（Re＞2×106）情況下，葉片渦的發(fā)放頻率滿足以下公式［15］：

式中：Yf為靠近葉片出氣邊處尾流的剪切層厚度，實(shí)際相當(dāng)于渦街的橫向?qū)挾龋籙s為出氣邊處的切向速度，St為斯特勞哈爾數(shù)。由于葉片為厚2.5 mm的直平頭，Yf厚度為2.25 mm；葉輪出氣邊速度由流場(chǎng)得到，約為38 m/s，由此得到渦發(fā)放頻率為3 378 Hz。一般在渦發(fā)放頻率附近的噪聲會(huì)明顯增高，由圖13可知，在3 150～4 000 Hz中心頻段處噪聲級(jí)較大，跟3 378 Hz這一頻率吻合。

3.2 計(jì)算結(jié)果與實(shí)驗(yàn)驗(yàn)證

葉輪與蝸殼表面脈動(dòng)力雖然從同一時(shí)刻開始計(jì)算，但步長不一樣，所以將他們當(dāng)作無規(guī)相位聲波進(jìn)行疊加，并對(duì)比疊加結(jié)果與實(shí)驗(yàn)結(jié)果（圖14），由圖可以看出，與試驗(yàn)結(jié)果相比，數(shù)值計(jì)算結(jié)果在中低頻處偏低，在高頻處有點(diǎn)偏高，最大偏差約14 dB，兩者整體趨勢(shì)一致。

圖14 實(shí)驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果對(duì)比Fig.14 Comparison between experimental results and predicted results

數(shù)值計(jì)算結(jié)果的精度與整個(gè)預(yù)報(bào)過程的每個(gè)環(huán)節(jié)都有關(guān)系：在流域建模、流場(chǎng)預(yù)報(bào)、輻射聲預(yù)報(bào)都會(huì)存在影響。其中流場(chǎng)預(yù)報(bào)時(shí)與網(wǎng)格個(gè)數(shù)、交界面設(shè)置（由于內(nèi)域旋轉(zhuǎn)需設(shè)置此項(xiàng)）、大渦模擬計(jì)算本身精度有關(guān)；聲學(xué)計(jì)算只考慮力源項(xiàng)的結(jié)果，采用軟件Virtual.Lab求解時(shí)，將葉輪表面劃分為上千個(gè)集中面源，這也降低了精度。另外，由圖14還可以看出，預(yù)報(bào)結(jié)果出現(xiàn)了較多的尖峰，這與采樣時(shí)間短、在計(jì)算時(shí)沒有將采樣結(jié)果當(dāng)作隨機(jī)信號(hào)有關(guān)。而且葉輪對(duì)聲場(chǎng)的反射無法計(jì)算，流速對(duì)聲傳播的影響也未考慮。

4 結(jié) 論

本文采用結(jié)構(gòu)化網(wǎng)格，使用大渦模擬方法獲得了離心風(fēng)機(jī)流場(chǎng)多個(gè)時(shí)刻的壓力分布，分別計(jì)算了風(fēng)機(jī)管中葉輪輻射聲和蝸殼脈動(dòng)壓力輻射聲，得出以下結(jié)論：

1）葉輪輻射聲和蝸殼脈動(dòng)壓力輻射噪聲相比，在低頻時(shí)前者明顯低于后者，這與蝸殼輻射聲面積較大有關(guān)；在高頻時(shí)前者又會(huì)高于后者，這與葉片較強(qiáng)的渦發(fā)放有關(guān)。

2）與噪聲實(shí)驗(yàn)結(jié)果相比，噪聲數(shù)值計(jì)算結(jié)果在低頻時(shí)偏低，在高頻時(shí)偏高，兩者整體趨勢(shì)基本一致，能夠滿足工程分析需要。

3）數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果較吻合，驗(yàn)證了葉輪輻射聲和蝸殼壓力脈動(dòng)輻射噪聲預(yù)報(bào)方法的有效性，也從側(cè)面反映出兩者是管中噪聲主要的噪聲源。