周 慧，郝志勇，劉聯(lián)鋆

(浙江大學(xué)能源工程學(xué)系，杭州 310027)

?

2015059

渦輪增壓器出口消聲器的仿真和改進(jìn)

周 慧，郝志勇，劉聯(lián)鋆

(浙江大學(xué)能源工程學(xué)系，杭州 310027)

應(yīng)用計(jì)算流體動(dòng)力學(xué)(CFD)方法，計(jì)算無(wú)流條件下渦輪增壓器出口阻抗復(fù)合型消聲器的傳遞損失，與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了CFD模型的可靠性，進(jìn)而計(jì)算有流條件下阻抗復(fù)合型消聲器的消聲性能。結(jié)果表明，添加多孔材料后，部分頻段消聲性能略有提高，特定頻段消聲性能變差。最后對(duì)消聲器結(jié)構(gòu)進(jìn)行改進(jìn)，并采用CFD法重新評(píng)估了改進(jìn)后消聲器的消聲性能。

消聲器；CFD法；多孔材料；消聲性能

前言

渦輪增壓器產(chǎn)生的噪聲分為渦輪噪聲和增壓器噪聲。渦輪產(chǎn)生的氣動(dòng)噪聲與發(fā)動(dòng)機(jī)的排氣噪聲一起進(jìn)入排氣系統(tǒng)，后處理裝置與排氣消聲器可以降低渦輪產(chǎn)生的氣動(dòng)噪聲和發(fā)動(dòng)機(jī)排氣噪聲向下游傳播[1]。增壓器產(chǎn)生的氣動(dòng)噪聲通過(guò)進(jìn)氣管和空氣濾清器向進(jìn)氣口傳播并輻射，典型的增壓器寬頻流體噪聲和1 階高速旋轉(zhuǎn)噪聲的頻率范圍為1.5～3.5kHz[2]。如果進(jìn)氣空濾器設(shè)計(jì)恰當(dāng)，或者選用吸聲能力較強(qiáng)的濾芯，增壓器噪聲對(duì)進(jìn)氣口噪聲的貢獻(xiàn)量較小，因此增壓器與中冷器連接的管路輻射成為增壓器噪聲的主要表現(xiàn)形式。設(shè)計(jì)能消除寬頻、高頻噪聲的增壓器出口消聲器對(duì)提高發(fā)動(dòng)機(jī)聲品質(zhì)很有必要。文獻(xiàn)[3]中對(duì)螺旋式消聲器性能進(jìn)行了分析和評(píng)估，在大于1kHz的特定頻段傳遞損失達(dá)到20dB以上。文獻(xiàn)[4]中也針對(duì)whoosh噪聲設(shè)計(jì)了一款多腔消聲器，在1.8～2.8kHz頻率范圍傳遞損失大于15dB。文獻(xiàn)[5]中對(duì)國(guó)外某款汽車(chē)發(fā)動(dòng)機(jī)使用的增壓器出口消聲器的性能進(jìn)行了預(yù)測(cè)和評(píng)估，在1.8～2.6kHz有明顯的消聲效果，但是在其消聲高值區(qū)域2.6～3.3kHz消聲性能不佳，因此需要對(duì)該消聲器進(jìn)行合理的改進(jìn)設(shè)計(jì)。

進(jìn)、排氣系統(tǒng)中常用多孔介質(zhì)作為阻性消聲元件，以獲取寬頻帶的消聲特性。聲波進(jìn)入多孔材料時(shí)，引起孔隙內(nèi)的空氣振動(dòng)，并與周?chē)睦w維摩擦，使聲能耗散成熱能，因此多孔材料具有較好的吸聲性能[6]。多孔材料的消聲特性由流阻率、孔隙率、纖維直徑和纖維排列方式等因素決定[7]，可通過(guò)實(shí)驗(yàn)與仿真擬合的方法獲取多孔材料的特性參數(shù)[8]。

針對(duì)內(nèi)部結(jié)構(gòu)復(fù)雜，且消聲頻率較高的消聲器，須用三維仿真方法進(jìn)行分析，三維CFD法在時(shí)域中解析全部流場(chǎng)方程，計(jì)算時(shí)可以考慮到介質(zhì)的黏性、平均流場(chǎng)和溫度場(chǎng)對(duì)其聲學(xué)性能的影響[9-10]，因此相比有限元法和邊界元法[11-12]有著特定的優(yōu)勢(shì)。

本文針對(duì)文獻(xiàn)[3]中的增壓器出口消聲器進(jìn)行了仿真。首先不改變消聲器結(jié)構(gòu)，對(duì)頭段腔體內(nèi)添加多孔材料，并測(cè)量多孔材料的聲學(xué)參數(shù)，進(jìn)而用CFD法計(jì)算了該阻抗復(fù)合型消聲器在無(wú)流條件下的傳遞損失，并與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證模型的準(zhǔn)確性；然后計(jì)算了有流條件下消聲器的消聲性能，為消聲器的改進(jìn)設(shè)計(jì)提供基礎(chǔ)；最后改進(jìn)消聲器結(jié)構(gòu)，并用CFD法計(jì)算結(jié)構(gòu)改進(jìn)前后無(wú)流條件下和全負(fù)荷工況下的消聲性能，初步提出工程中的改進(jìn)設(shè)計(jì)方案。

1 模型分析

該消聲器內(nèi)部結(jié)構(gòu)剖視圖如圖1所示。在不改變結(jié)構(gòu)的基礎(chǔ)上，將工程中常用的多孔材料—玻璃纖維添加在頭段Helmholtz共振腔內(nèi)。由于消聲器整體尺寸較小，總長(zhǎng)為175mm，吸聲材料可均勻填充至腔內(nèi)，無(wú)需額外固定裝置。

1.1 多孔材料參數(shù)測(cè)量實(shí)驗(yàn)

由于實(shí)驗(yàn)條件限制，可通過(guò)實(shí)驗(yàn)與仿真擬合的方法得到多孔材料的特性參數(shù)，即將多孔材料置于管道內(nèi)，測(cè)量管道在多孔材料兩端處的消聲量，圖2為多孔材料實(shí)驗(yàn)示意圖，實(shí)驗(yàn)測(cè)量現(xiàn)場(chǎng)如圖3所示。根據(jù)測(cè)量結(jié)果和經(jīng)驗(yàn)判斷比流阻率和孔隙率的范圍，初步設(shè)定3個(gè)參數(shù)的初始值。建立包含多孔區(qū)域管道的二維CFD軸對(duì)稱(chēng)模型如圖4所示，計(jì)算管道的消聲量。對(duì)比消聲量計(jì)算值和測(cè)量值，反復(fù)修改參數(shù)，使計(jì)算值和測(cè)量值最大程度地吻合后，即確認(rèn)得到近似參數(shù)[8]。消聲量NR可以表示為

(1)

式中:NR為消聲量；Lpi和Lpo為消聲器進(jìn)、出口的聲壓級(jí)；pi和po為進(jìn)、出口的聲壓。

測(cè)得的多孔材料消聲量曲線(xiàn)如圖5所示?？梢钥闯?，在700Hz以上有無(wú)多孔材料的消聲量曲線(xiàn)才看得出區(qū)別，說(shuō)明多孔材料在700Hz以上才起作用。圖5中的兩條NR曲線(xiàn)在200Hz內(nèi)有較大差別，可能是因?yàn)閮蓽y(cè)點(diǎn)相距較近(200mm)，多孔材料在該頻率范圍內(nèi)起了一定的隔聲作用?？傮w來(lái)說(shuō)，存在多孔材料時(shí)的消聲量要比沒(méi)有多孔材料高。

根據(jù)經(jīng)驗(yàn)和反復(fù)仿真判斷玻璃纖維σ=3×103Pa·s/m2，孔隙率q=0.91。最終結(jié)果如圖6所示，內(nèi)置玻璃纖維的直管內(nèi)NR測(cè)量值與CFD仿真結(jié)果吻合良好，因此該多孔材料的參數(shù)選取合理，為后續(xù)仿真計(jì)算的可靠性提供了保證。

1.2 傳遞損失測(cè)量實(shí)驗(yàn)

由于傳遞損失是消聲器固有屬性，與噪聲源特性及管道系統(tǒng)無(wú)關(guān)，常用于評(píng)價(jià)消聲器聲學(xué)性能。這里通過(guò)對(duì)比傳遞損失的實(shí)驗(yàn)測(cè)量值及仿真計(jì)算值來(lái)驗(yàn)證CFD模型的可靠性，傳遞損失TL可表示為

(2)

式中:Lwi為入射聲功率級(jí)；Lwt為透射聲功率級(jí)；Ai為消聲器進(jìn)口的截面面積；Ao為消聲器出口的截面面積；pi為消聲器進(jìn)口處的入射聲壓；pt為末端為無(wú)反射條件下消聲器出口處的透射聲壓。

圖7為兩負(fù)載法[13]測(cè)傳遞損失的實(shí)驗(yàn)測(cè)量現(xiàn)場(chǎng)。將聲源置于半消聲室外，用管道將聲源引入到消聲室內(nèi)，并在半消聲室內(nèi)連接消聲器的入口端，出口端連接負(fù)載，在消聲器的上游和下游分別布置兩個(gè)傳聲器。忽略介質(zhì)的黏性，通過(guò)改變末端負(fù)載得到兩組測(cè)量結(jié)果。為了測(cè)量準(zhǔn)確，在負(fù)載末端均貼有吸聲材料，制作特定木樁控制傳聲器位置，并用橡皮泥保證密封性。其中消聲器上游兩傳聲器和下游兩傳聲器距離均為30mm，且4個(gè)傳聲器高度保持一致。

1.3 CFD模型

對(duì)消聲器和填充多孔材料部分分別劃分網(wǎng)格，同時(shí)為了減少入口邊界的反射聲波對(duì)消聲器入口測(cè)點(diǎn)的干擾，使出口測(cè)點(diǎn)能完整捕捉透射信號(hào)，在消聲器入口和出口端分別增加與消聲器入口和出口直徑相同的延長(zhǎng)管，長(zhǎng)度為0.8和0.6m，并劃分網(wǎng)格。網(wǎng)格類(lèi)型包括四面體網(wǎng)格、菱形網(wǎng)格和六面體網(wǎng)格，網(wǎng)格大小為3.5 mm[10]。消聲器CFD模型如圖8所示。

采用分離式三維二階隱式非定常求解器，二階迎風(fēng)格式對(duì)壓力、密度、動(dòng)量、湍動(dòng)能、湍流耗散率和能量方程進(jìn)行離散化處理，時(shí)間步長(zhǎng)取為5μs[9-10]。采用理想氣體模擬空氣介質(zhì)，多孔材料參數(shù)根據(jù)上述的結(jié)果進(jìn)行設(shè)置，考慮傳熱，用標(biāo)準(zhǔn)k-ε湍流模型計(jì)算無(wú)平均流和有平均流條件下的消聲性能。

在無(wú)流條件下，設(shè)置入口和出口均為無(wú)反射邊界，入口測(cè)點(diǎn)位置為距消聲器入口截面0.4m處，出口測(cè)點(diǎn)位置為距消聲器出口截面0.2m處。通過(guò)設(shè)置隨時(shí)間變化的馬赫數(shù)來(lái)產(chǎn)生脈沖激勵(lì)，該變化馬赫數(shù)為持續(xù)16μs的半個(gè)正弦波，幅值為0.004。模型其它表面均設(shè)為無(wú)滑移壁面[5]。該消聲器的材料是增強(qiáng)玻纖，其結(jié)構(gòu)緊湊，剛度較高，計(jì)算時(shí)將消聲器殼體當(dāng)作剛性壁面處理。

有平均流時(shí)，先計(jì)算給定馬赫數(shù)下消聲器內(nèi)部穩(wěn)態(tài)流場(chǎng)，入口流動(dòng)方向?yàn)榇怪比肟诮孛嫦騼?nèi)，出口流動(dòng)方向?yàn)榇怪背隹诮孛嫦蛲猓R赫數(shù)大小由入口截面面積和流量決定。待穩(wěn)態(tài)結(jié)果收斂后，在入口施加激勵(lì)，即在穩(wěn)態(tài)馬赫數(shù)的基礎(chǔ)上，加載一個(gè)和無(wú)流條件下相同的隨時(shí)間變化的馬赫數(shù)，其余設(shè)置與無(wú)流條件下相同。

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

2.1 無(wú)流條件下傳遞損失的計(jì)算

用CFD法計(jì)算得到的聲衰減曲線(xiàn)如圖9所示。由圖可見(jiàn)，曲線(xiàn)整體趨勢(shì)上都與實(shí)驗(yàn)測(cè)量值相吻合，消聲峰值對(duì)應(yīng)的頻率為1.9kHz。CFD法計(jì)算的傳遞損失峰值和測(cè)量結(jié)果基本一致，但在頻率高于2.5kHz以后，計(jì)算的傳遞損失比測(cè)量值低，可能是因?yàn)闇y(cè)點(diǎn)受到進(jìn)、出口邊界的微小反射波的干擾和建模時(shí)對(duì)模型進(jìn)行了簡(jiǎn)化，與實(shí)際幾何形狀存在微小差別，另外四面體網(wǎng)格模型在高頻區(qū)域的計(jì)算準(zhǔn)確性不足可能也是一個(gè)原因。實(shí)驗(yàn)測(cè)量的TL曲線(xiàn)存在不規(guī)則波動(dòng)，可能是因?yàn)闇y(cè)量時(shí)沒(méi)有校核傳聲器的相位差，另外雙負(fù)載法中的兩個(gè)負(fù)載阻抗差別小也可能導(dǎo)致測(cè)量結(jié)果不準(zhǔn)確。

通過(guò)對(duì)比實(shí)驗(yàn)測(cè)量值和仿真計(jì)算值，驗(yàn)證了CFD模型的可靠性，可用來(lái)計(jì)算有流條件下阻抗復(fù)合型消聲器的傳遞損失。

2.2 有流條件下傳遞損失的計(jì)算

當(dāng)入口溫度為27℃，不同入口相對(duì)壓力及不同入流馬赫數(shù)時(shí)，用CFD法計(jì)算的消聲器的傳遞損失頻譜曲線(xiàn)如圖10所示。有平均流時(shí)，隨著馬赫數(shù)的增加傳遞損失曲線(xiàn)向高頻區(qū)域移動(dòng)，傳遞損失幅值呈減小趨勢(shì)；馬赫數(shù)相同時(shí)，入口相對(duì)壓力為30 kPa與入口相對(duì)壓力為0時(shí)傳遞損失曲線(xiàn)幾乎重合，說(shuō)明不同的平均壓力對(duì)計(jì)算結(jié)果影響較小。

3 消聲器的改進(jìn)

圖11為添加多空材料前后消聲器傳遞損失的CFD仿真值結(jié)果對(duì)比。由圖可見(jiàn)，未添加多孔材料時(shí)，消聲器在1.9和2.5kHz附近有兩個(gè)消聲峰值，左邊的消聲峰值是由膨脹腔-插入管結(jié)構(gòu)導(dǎo)致的，右邊的消聲峰值是由Helmholtz共振腔導(dǎo)致的。添加了多孔材料后，第二個(gè)消聲峰值消失，這是由于Helmholtz共振腔內(nèi)添加了多孔材料后，占據(jù)了腔體內(nèi)的原有的消聲容積，阻抗變化，消聲效果變差。但是消聲器在1.0～2.3kHz頻率段傳遞損失略有提高，這是因?yàn)槁暡ㄟM(jìn)入多孔材料時(shí)，引起孔隙內(nèi)的空氣振動(dòng)并與周?chē)睦w維摩擦，使聲能耗散成熱能。整體上看，由于安裝空間的限制，增壓器出口消聲器體積較小，且容易造成堵塞，添加多孔材料并不適宜。

為改善消聲器在2.5～3.2kHz的消聲高值區(qū)和3.2 kHz以上的次消聲高值區(qū)的性能，將幾乎不起消聲作用的末端膨脹槽結(jié)構(gòu)改為膨脹腔-插入管結(jié)構(gòu),且根據(jù)文獻(xiàn)[14]中的研究表明，小孔-空腔組合極有可能引發(fā)空氣動(dòng)力噪聲，因此將Helmholtz共振腔改為加工更簡(jiǎn)單的內(nèi)插管結(jié)構(gòu)，同時(shí)調(diào)整第二段的插入管長(zhǎng)度。改進(jìn)后消聲器CFD模型局部剖視圖如圖12所示，整體結(jié)構(gòu)簡(jiǎn)單，加工方便，符合工程應(yīng)用。

用CFD法計(jì)算得到的改進(jìn)前后消聲器傳遞損失頻譜曲線(xiàn)如圖13所示。由圖可見(jiàn)，改進(jìn)后傳遞損失的兩個(gè)峰值都向右移動(dòng)，即TL峰值時(shí)刻都提高了約300Hz，且末段膨脹腔-插入管結(jié)構(gòu)在4.0kHz導(dǎo)致了新的峰值，總的來(lái)說(shuō)滿(mǎn)足性能要求。

同時(shí),計(jì)算發(fā)動(dòng)機(jī)全負(fù)荷，轉(zhuǎn)速為5 500r/min時(shí)，壓氣機(jī)出口流量Vm= 0.142kg/s(消聲器入口Ma= 0.2)，T=98℃，p=86kPa工況對(duì)應(yīng)平均流作用下，改進(jìn)前后消聲器的傳遞損失，結(jié)果如圖14所示。由圖可見(jiàn)，有流條件下，頻率為2.3～4.2kHz時(shí)，改進(jìn)后消聲器TL的幅值比改進(jìn)前高。因此，可初步判斷該消聲器的設(shè)計(jì)是合理的。

4 結(jié)論

(1)應(yīng)用CFD法分析渦輪增壓器出口復(fù)雜多腔體消聲器在添加多孔材料后的聲學(xué)性能，并與實(shí)驗(yàn)結(jié)果對(duì)比發(fā)現(xiàn)消聲器在1.0～2.3kHz頻率段傳遞損失略有提高，但由于填充多孔材料，占據(jù)了腔體內(nèi)原有的消聲容積，使2.3～2.9kHz消聲效果變差，從工程應(yīng)用角度，由于安裝空間的限制，增壓器出口消聲器體積較小，添加多孔材料并不適宜。

(2)用CFD法計(jì)算了有平均流條件下，在相同溫度，不同入口壓力的情況下不同入流馬赫數(shù)時(shí)阻抗復(fù)合型消聲器的傳遞損失。發(fā)現(xiàn)有平均流時(shí)，隨著馬赫數(shù)的增加傳遞損失曲線(xiàn)向高頻區(qū)域移動(dòng)，傳遞損失幅值呈減小趨勢(shì)，且不同的平均壓力對(duì)計(jì)算結(jié)果影響較小。

(3)改進(jìn)消聲器首段Helmholtz共振腔和末段膨脹槽結(jié)構(gòu)為內(nèi)插管結(jié)構(gòu)，消聲性能有較大幅度提高，且制作更為簡(jiǎn)便，符合工程應(yīng)用。

[1]R?mmalH, ?bomM.ExperimentalDeterminationofSoundTransmissioninTurbo-Compressors[C].SAEPaper2009-01-2045.

[2]R?mmalH, ?bomM.AcousticsofTurbochargers[C].SAEPaper2007-01-2205.

[3]GraefensteinA,WenzelW. “Herschel-QuinckeSpiral”ANewInterferenceSilencer[C].SAEPaper2003-01-1722.

[4]LeeIJ,SelametA,KimH.DesignofaMulti-ChamberSilencerforTurbochargerNoise[C].SAEPaper2009-01-2048.

[5] 劉聯(lián)鋆，郝志勇，錢(qián)欣怡，等. 渦輪增壓器出口消聲器的性能預(yù)測(cè)和評(píng)估[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2013(2): 197-201.

[6] 方丹群. 空氣動(dòng)力性噪聲與消聲器[M]. 北京：科學(xué)出版社, 1978.

[7]BeranekLL.AcousticImpedanceofPorousMaterials[J].TheJournaloftheAcousticalSocietyofAmerica， 1942, 13(3): 248-260.

[8] 劉聯(lián)鋆，郝志勇，錢(qián)欣怡. 空濾器濾芯聲學(xué)特性的仿真方法[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版)， 2012(10): 1784-1789.

[9] 徐航手，季振林，康鐘緒. 抗性消聲器傳遞損失預(yù)測(cè)的三維時(shí)域計(jì)算方法[J]. 振動(dòng)與沖擊， 2010, 25(4): 39-41.

[10] 石巖，舒歌群，畢鳳榮. 基于計(jì)算流體動(dòng)力學(xué)的內(nèi)燃機(jī)排氣消聲器聲學(xué)特性仿真[J]. 振動(dòng)工程學(xué)報(bào)， 2011, 24(2): 205-209.

[11] 季振林. 直通穿孔管消聲器聲學(xué)性能計(jì)算及分析[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2005, 26(3):302-306.

[12] 金巖，郝志勇，劉永. 空氣濾清器聲學(xué)性能的改進(jìn)設(shè)計(jì)[J]. 內(nèi)燃機(jī)工程, 2007(6): 58-60.

[13]TaoZ,SeybertA.AReviewofCurrentTechniquesforMeasuringMufflerTransmissionLoss[C].SAEPaper2003-01-1653.

[14]AshcroftGB,TakedaK,ZhangX.ANumericalInvestigationoftheNoiseRadiatedbyaTurbulentFlowoveraCavity[J].JournalofSoundandVibration， 2003, 265(1): 43-60.

Simulation and Improvement of Silencer at Turbocharger Outlet

Zhou Hui, Hao Zhiyong & Liu Lianyun

DepartmentofPowerEngineering,ZhejiangUniversity,Hangzhou310027

By applying computational fluid dynamics (CFD) technique, the transmission loss (TL) of impedance compound silencer at turbocharger outlet under non-flow condition is calculated and compared with experiment results, verifying the reliability of CFD model, and then the noise attenuation performance of impedance compound silencer under mean flow condition is calculated. The results show that after the silencer is covered by porous stuff, its noise attenuation performance slightly enhances in some frequency range but becomes poor in specific frequency range. Finally the structure of silencer is modified with its noise attenuation performance reevaluated with CFD technique.

silencer; CFD technique; porous material; noise attenuation performance

原稿收到日期為2013年7月31日，修改稿收到日期為2013年9月29日。