風洞射流和汽車繞流對聲傳播的影響

王毅剛，李方旭，沈哲，楊志剛

?

風洞射流和汽車繞流對聲傳播的影響

王毅剛，李方旭，沈哲，楊志剛

(同濟大學上海地面交通工具風洞中心，上海201804)

汽車風洞氣動噪聲測試中，由于射流和汽車繞流場的存在，流場外測量不能獲取汽車聲源的準確信息。研究以兩種流場對聲傳播方向(聲源漂移量)的影響開展。通過在風洞射流核心區安放特定聲源和實車的試驗，分別得到了射流和汽車繞流場對聲傳播的影響。研究表明兩種流場對聲傳播都會產生較大影響，試驗分析應加以考慮；研究分別對風洞射流和汽車繞流進行了氣動數值仿真，據此擬合出剪切層和汽車繞流場速度分布線性化模型，并結合幾何聲學方法，建立了連續多層均勻運動介質聲傳播模型，得到了預測射流和汽車繞流產生聲漂移量的方法。研究表明該方法可以有效預測風洞聲學實驗中介質運動對聲傳播方向的影響，并可用于風洞聲學試驗中流動對傳播方向影響的修正。

汽車聲學風洞；幾何聲學；聲傳播；聲源漂移量

0 引言

汽車聲學風洞是汽車氣動噪聲研究的基礎測試設施。聲學風洞射流結構中汽車的氣動噪聲源及其聲傳播特性是風洞聲學實驗測量的重點內容。在實際試驗測試中，許多聲學測量必須在射流區域外進行，一方面射流結構的存在使得聲從射流核心區內的聲源出發向流場外測試區域傳播的過程中經歷了勢流核心區和剪切層，其傳播的路徑和大小發生改變；另一方面汽車繞流場的非定常流動區域也使汽車表面的噪聲源向外傳播聲時其傳播路徑和大小發生改變。這樣，射流結構外的測試儀器測量到的噪聲信號和實際信號就存在差異，因此，風洞射流結構和汽車繞流場的存在成為準確了解汽車氣動噪聲源的重要技術障礙，在風洞聲學測量中必須解決。

國外在上世紀50年代開始對射流中的聲傳播問題進行研究，最初的目的是用于了解噴氣式發動機噪聲向周圍環境的傳播[1-2]。其中較為成功的是利用幾何聲學原理研究流體對聲的折射影響。莫爾斯[3]探索了在兩種均勻介質交界面上出現的聲折射、反射和投射現象，奠定了兩種不同介質中的聲傳播分析的基本方法。文獻[4-6]指出，Aimet于1975年以兩相對運動流體中平面聲波傳播的分析方法為基礎，利用幾何聲學原理，研究了流體對聲的折射及聲擴散影響，創建了經典的Amiet 剪切層修正理論。該理論將剪切層簡化成一個無限薄的渦面，射流和周圍空氣在各自的區域是各向同性的。而實際剪切層有一定的厚度，且沿順流方向剪切層厚度增加，所以該方法對流場的簡化與真實情況差別較大，僅在噴口附近能較好預測剪切層對聲傳播的影響，遠離噴口位置誤差較大。近年來國內學者也對該理論進行了數值和實驗方面的驗證[7]。另外，關于汽車繞流場對聲傳播的影響還沒有相關文獻描述。

基于上述研究狀況，本文利用上海地面交通工具風洞中心汽車聲學風洞，分別開展了射流結構和汽車繞流結構對聲源傳播方向影響的試驗研究，在此基礎上針對一定厚度紊流結構中的聲傳播，利用數值計算手段得到流場定常流動特性，進一步利用幾何聲學方法建立聲傳播模型，分析和確定射流剪切層及汽車繞流場對聲傳播方向的影響，為風洞中氣動噪聲源定位測試提供修正方法。

1 射流結構和汽車繞流對聲傳播影響的風洞實驗

1.1 射流結構聲傳播影響實驗

射流結構對聲傳播方向的影響包括了勢流核心區和剪切層區。認為勢流核心區為均勻定常流動，通常試驗件放置于該區域。為了了解這兩個區域對聲傳播的影響，試驗將如圖1所示的一個特定聲源(外形為翼型設計，中間黑色圓形即為內置揚聲器位置，支撐件也為翼型設計)，置于風洞核心區一定的高度位置，其有良好的氣動外形，經數值仿真驗證，在聲源位置附近對流場幾乎無擾動，自身產生的氣動噪聲較低，通過外接功放和信號源組成可控聲源系統。

試驗在同濟大學-上海地面交通工具風洞中心整車氣動聲學風洞中完成。該風洞為3/4開口回流風洞，試驗段背景噪聲在國際同類風洞中較低(在160 km/h風速下低于61 dB(A))。聲源揚聲器中心距離地面1.1 m，該高度對應的剪切層位置流動認為無地面效應(近似和后續的聲傳播模型一致)。

圖1 對稱翼型聲源整流罩和安裝效果

實驗中采用的相位傳聲器陣列聲源識別系統具體參數如下：陣列尺寸為1.8 m×1.8 m，通道數120；傳感器采用1/4 inch陣列專用壓力型傳聲器；數據采集系統為120通道高精度數采，單通道最大采樣率為192 kHz；分析系統為基于波束成型技術的聲源識別軟件NoiseImage。以上硬件設施是當前聲學風洞中聲源定位的先進方法，聲源的定位誤差穩定在cm量級，可以保證實驗的精確度和穩定性。試驗段噴口、聲源和測試設備相對位置如圖2所示。

圖2 聲源和傳聲器陣列相對位置示意圖

由于采用的相位傳聲器陣列是二維陣列，也即在實驗結果中可以觀測到由于介質運動產生的聲源在流向和高度方向的漂移。由于風洞試驗段中的流動以水平方向流動為主，高度方向的速度分量基本可以忽略。實驗結果中高度方向聲漂移量普遍比橫向漂移量小1~2個數量級，故在此忽略不計是合理的。

實驗主要研究不同射流速度對聲傳播的影響，保持其它實驗設置不變，實驗中射流速度設置為以下工況：無射流(0 m/s)、15~36 m/s區間內，速度每間隔3 m/s設置一個實驗工況。暫不考慮聲頻率對實驗的影響，故聲源信號采用白噪聲信號。

1.2 汽車繞流場對聲傳播影響實驗

為了得到汽車繞流場對聲傳播的影響，將某款實車置于試驗段內，在和上述特定聲源試驗相同射流結構的作用下，在汽車周圍會形成繞汽車的流場。同時，對應上述特定聲源實驗中的揚聲器相同位置，在汽車表面安裝揚聲器聲源，它被安裝在流線型的整流罩中，位于汽車左側前風窗外側，在汽車的繞流場之中。采用與1.1節相同的試驗測試儀器、聲源信號、測量狀態和測量方法，兩者試驗結果的差異就是汽車繞流場帶來的差異。圖3為汽車在風洞試驗段的位置及揚聲器安裝位置。

圖3 汽車表面聲源安裝示意圖

1.3 實驗結果

射流結構和汽車繞流結構對聲傳播路徑的影響用漂移量來描述。實驗中用相位傳聲器陣列可以測量到真實的聲信號，但在據此計算聲源位置時，將流動空氣介質當成靜止介質處理，由此確定的聲源位置和實際聲源位置有一定的偏移量，在此定義為聲源漂移量。

圖4為兩次試驗得到的測量結果。從圖4可以看出，射流結構會改變聲傳播方向，隨著射流流速的增加，聲源漂移量呈線性增加關系。在增加汽車繞流的情況下，聲源漂移量進一步增加，但增加量較小，說明汽車繞流對聲傳播方向的影響較小。由此可見，射流結構和汽車繞流結構對聲源位置的識別有較大影響，在聲學測量時應進行修正。

圖4 聲漂移量實驗結果

2 聲源漂移量分析及預測

2.1 風洞射流和汽車繞流現象

如上所述，汽車在聲學風洞進行氣動噪聲測試時，汽車放置于試驗段勢流核心區。來流和汽車相互作用時，汽車表面的突出件、凹槽等部件附近會成為氣動噪聲源，其產生的噪聲向車外空間傳遞時，聲波先經過汽車周圍的繞流場，進入射流核心區，再傳遞到剪切層區，最后傳到流場外靜止的空氣中，被測量設備接收。圖5為上述聲傳播過程示意圖。和靜止介質中的聲傳播相比，這些空氣介質流動會對聲傳播產生影響。

圖5 汽車繞流和射流對聲傳播的影響示意圖

由于核心區為均勻流動區，聲傳播可以通過解析法求解。而汽車擾流區和剪切層區為非定常流動區域，難以通過解析法求解，對其流動特性的了解是建立這兩部分聲傳播的前提。所以，下面利用數值仿真手段對這兩部分的流動進行計算，為建立聲傳播模型奠定基礎。

2.2 射流流動結構和汽車繞流特性

數值仿真的對象分別是上海地面交通工具風洞中心3/4開口風洞全尺寸模型和與實驗同款的實車全尺寸模型。

2.2.1 風洞射流流動分析

在數值模擬中，選取與實驗相對應的風洞結構和試驗工況。其中幾何建模立足于氣動聲學風洞具體結構和尺寸，其計算區域由收縮段、實驗段、收集口、擴散段等組成。為了使下游流動不出現回流，人為增加長度為10 m的延伸段。計算區域使用商業軟件ICEM進行網格劃分，網格模型見圖6。

圖6 數值模型示意圖

計算區域創建能提高數值精度的六面體網格。整個計算區域網格總數約為1 025萬。采用基于有限體積法的商業軟件Fluent對非結構網格進行數值求解。在滿足設定殘差和監控物理量為常量后停止迭代，并以風洞射流軸線靜壓系數與實驗值吻合作為數值計算準確的依據[8]。

以風洞射流速度30 m/s(108 km/h)為例對風洞流場進行分析，以風洞實驗段地板高度=0 mm為坐標平面，建立=1100 mm高度的水平面。圖7是距風洞實驗段地面1100 mm高度處水平截面速度云圖，可見射流剪切層隨距離噴口的距離增大而逐漸增厚。

圖7 1 100 mm高度處水平截面速度云圖

根據已有的對射流結構的流動研究，在射流的主體段，各個斷面的縱向速度分量的分布具有相似性(也稱作自保性)[9]。由圖8可知，在風洞實驗區(距離噴口9 m以內)內部，射流核心區的速度與來流速度相等，但隨著距噴口距離增大，剪切層區域逐漸侵蝕流動核心區和剪切層外部區域，使得剪切層區域逐漸變厚。從圖8中可見，剪切層內的氣流速度分布在橫向，基本呈線性遞減分布，故在之后的聲漂移量計算模型中，剪切層內的速度分布按照線性函數進行擬合，取邊界層范圍為速度分布在來流速度5%~95%的區域。

2.2.2 汽車繞流流動分析

對試驗中實車的繞流場仿真，汽車模型如圖9所示。計算區域尺寸為10×5×5，、、分別為汽車的長寬高。

數值仿真以各項計算殘差穩定，升阻力系數均穩定且接近風洞實驗數值作為數值計算準確的依據；以來流風速30 m/s(108 km/h)為例，對汽車繞流場進行分析。如圖10所示，在距離地面1 100 mm高度的水平截面中，除了后視鏡后方的氣流所受的擾動較大，汽車兩側一定范圍內的氣流也因為受到汽車的阻塞而流速加快。所以，若車表面(如門窗縫隙、臺階等)存在噪聲源，其產生的噪聲向外傳播時，在這些擾動區域聲的傳播特性會發生變化，因此在分析汽車的流場速度對聲傳播的影響時應該合理簡化汽車繞流場，使其更加接近真實流動狀況。

(a) 射流核心區、剪切層及遠場整體速度分布

(b) 剪切層部分速度分布

圖8 距離噴口不同距離1 100 mm高度速度分布

Fig.8 Velocity distributions at different distances from the nozzle at height of 1 100 mm

圖9 全尺寸汽車模型示意圖

圖10 距離地面1 100 mm高度汽車流場云圖

2.3 射流剪切層和汽車繞流的速度擬合

要建立考慮速度場的聲從核心勢流區到遠場的傳播模型，必須對聲傳播路徑上的速度場進行擬合。實驗中聲陣列和聲源正對放置，即沒有氣流影響時，聲波陣面的法線和來流方向夾角為90°。

2.3.1 射流剪切層的速度擬合

根據上述氣動數值計算結果，剪切層區域定義為速度分布為來流速度5%~95%的區域，該區域剪切層厚度隨距噴口距離線性變化；剪切層在核心區一側的速度為射流速度，外側的速度近似為0。結合風洞特定聲源聲學實驗中聲源和聲陣列的位置關系，對30 m/s射流速度下，距離風洞噴口3.9 m、高度1.1 m位置橫向直線上的速度進行擬合，結果如圖11所示。可見，對剪切層的線性擬合是比較理想的。

圖11 風洞射流剪切層速度擬合示意圖

根據以上條件，如圖11所示將整個剪切層從內邊界到外邊界均分為層，將勢流核心區(來流速度為)作為第1層，靜止空氣區為第2層。在剪切層內，每一層內的速度變化很小，故假設在每一層內速度是均勻的，第層內的速度如式(1)所示：

2.3.2 汽車繞流的速度場擬合

取來流速度為30 m/s，汽車繞流場在高為1100mm、與聲源位置相重疊的橫向直線上的速度分布如圖12所示。

由圖12可見，側窗外側的氣流速度分布經歷了由邊界層速度為零逐漸增長到一定速度之后又下降再次增大的過程，并且與所有的鈍體繞流類似[10-11]，繞流場中的最大速度大于來流速度并隨著位置遠離汽車表面速度逐漸下降。

圖12 車表面聲源向遠場水平直線的速度分布

圖13 汽車表面聲源至遠場聲陣列的速度擬合

2.4 分層速度場聲傳播模型

根據Snell聲折射定理，聲經過兩種不同介質的邊界時，沿邊界方向聲傳播的速度是連續的[3]。沿著邊界方向的聲波速度由介質內的固有聲速和流體速度在邊界上的分量組成，兩者應當是相等的。由此可得

圖14 聲經過不同馬赫數流動介質的折射關系

在低馬赫數情況下，可以認為流體中的固有聲速與靜止空氣中的聲速相等，即12。將等式兩邊同時除以，用和1替代1和2，得到聲在兩個不同馬赫數下的傳播角度關系：

本文中要考慮的流場，其流動速度變化較為復雜，故考慮將速度場分布劃分成多層勻速流動的集合。圖15展示了多層速度分布情形下聲通過每一個分層邊界的過程。

3 結果分析比較

3.1 實驗結果和預測結果的對比

通過建立上述的分層模型，求解每一層中聲法線方向和來流速度的夾角，從而求得總的聲漂移量數據。

圖16和圖17分別表示特定聲源和汽車表面安裝聲源兩種情況下聲漂移量的實驗結果和預測結果對比。圖18為兩種情況下預測結果和實驗結果的誤差對比。從總體上看，兩種情況在大部分風速下，速度場分層模型對聲漂移量的預測和實驗結果相一致，在風速較高為33 m/s時，由于實驗聲源出現小的晃動，致使試驗出現誤差，引起此風速下偏差較大。所以，本預測方法可以用于實驗中的聲漂移量修正。

圖16 翼型導流罩聲源結果對比

圖17 汽車表面聲源結果對比

圖18 新模型對兩種實驗的預測誤差

3.2 對兩種工況下聲漂移量數值接近的解釋

設計特定聲源和汽車表面安放聲源兩種工況下的實驗目的，是為了研究汽車繞流場對聲傳播的影響，但是所測得不同風速條件下兩者的聲漂移量偏差均較小，結合本文對聲漂移量的預測方法對該實驗現象進行以下解釋：

(1) 兩種實驗工況的差別除了聲源至聲陣列的流場分布不同之外，聲源在風洞流場中的位置、實驗測試設備、測量方法等都是相同的，即影響實驗結果的唯一因素是汽車繞流場和射流核心區速度分布的差異。

(2) 根據幾何聲學理論，流動介質的流速越快，聲在介質中傳播所偏離路徑的角度就會越大。圖19描述了兩種工況下聲源至聲陣列的速度擬合曲線，雖然特定聲源工況下，核心勢流區域內的最大流速相對較小，但是在整個核心勢流區內部其速度分布維持在來流速度的水平；而汽車繞流雖然在大部分核心勢流區內流速較快，但是汽車表面的邊界層流速由0變化到最大，故兩者對聲漂移量的影響雖然在局部有差異，但是在整個核心勢流區總體上是接近的，所以，聲漂移大小沒有表現出明顯差異，但汽車繞流場對漂移量的影響更復雜一些，對漂移量有一定的影響。

圖19 兩種工況下聲源至聲陣列速度擬合曲線對比

4 結語

本文針對汽車聲學風洞中由于介質流動導致聲源的聲傳播方向變化(漂移量)問題開展了研究。首先，通過安裝在風洞射流核心區的特定聲源試驗和實車的繞流試驗，分別得到了核心區流動及剪切層對聲傳播的影響，以及汽車繞流場對聲傳播的影響。研究表明，射流結構和汽車繞流場對汽車表面聲源向外界的聲傳播都會產生影響，在風洞聲學試驗和聲源特性識別時應考慮其影響，并需要進行修正；

其次，為了分析和預測上述影響，分別對風洞射流結構和汽車繞流結構進行了氣動數值仿真，得到了剪切層和汽車繞流場的平均速度分布，擬合出速度分布線性化模型；

最后，在此基礎上，結合幾何聲學理論中流速對聲傳播方向的影響機理，將速度場等效成連續多層均勻運動介質，建立了多層速度場分布的聲傳播模型，得到了預測射流結構和汽車繞流結構產生聲漂移量的預測方法，并和實驗結果比較，驗證了該方法可以有效預測風洞聲學實驗中介質運動對聲傳播方向的影響。這種方法也可以用于風洞聲學試驗中對傳播方向影響的修正。

[1] L. M. B. C. CAMPOS. The spectral broadening of sound by turbulent shear layers[J]. Part 2. the spectral broadening of sound and aircraft noise, J. Fluid Mech. 1978, 89(4): 761-783..

[2] 喬渭陽. 航空發動機氣動聲學[M]. 北京: 北京航空航天大學出版社, 2010.

[3] 莫爾斯. 聲學理論(下)[M]. 北京: 科學出版社, 1984: 829-835. Morsse P. M. Theoretical Acoustics[M]. Beijing: Science Press, 1984: 829-835.

[4] Thomas J Mueller. Aeroacoustic Measurement[M]. Springer, 2002.

[5] Amiet R K. Correction of open jet wind tunnel measurements for shear layer refraction[C]//2ndAIAA Aeroacoustics Conference, AIAA Paper. 75-532.

[6] Schlinker R H, Amiet R K. Refraction and scattering of sound by a shear layer[J]. J Acoust Soc Am, 1981, 70(6): 1797-1799.

[7] 張雪, 陳寶, 盧清華. Amiet剪切層理論的角度折射驗證研究[J]. 應用聲學, 2014, 33(5): 433-438. ZHANG Xue, CHEN Bao, LU Qinghua. Verification of angle refraction correction based on Amiet's shear layer theory[J]. Journal of Applied Acoustics, 2014, 33(5): 433-438.

[8] 李啟良, 鄭志強, 賈青, 等. 兩種改善汽車風洞軸向靜壓系數的方法[J]. 同濟大學學報(自然科學版), 2010, 38(3): 422-426. LI Qiliang, ZHENG Zhiqiang, JIA Qing, et al. Two methods to improve the axial static pressure coefficient of automotive wind tunnel[J]. Journal of Tongji University: (Natural Science), 2010, 38(3): 422-426.

[9] 余常昭. 紊動射流[M]. 北京: 高等教育出版社, 1993.

[10] 傅立敏. 汽車空氣動力學[M]. 北京: 機械工業出版社, 1998.

[11] 李玲, 李玉梁. 應用基于RNG方法的湍流模型數值模擬鈍體繞流的湍流流動[J]. 水科學進展, 2000, 11(4): 357-361. LI Ling, LI Yuliang. Numerical simulat ion of turbulent flow around bluf f bodies using the rng k-ε turbulent model[J]. Advances in Water Science, 2000, 11(4): 357-361.

The influences of wind tunnel jet flow and airflow around vehicle on sound propagation

WANG Yi-gang, LI Fang-xu, SHEN Zhe, YANG Zhi-gang

(Shanghai Automotive Wind Tunnel Center, Tongji University. Shanghai 201804, China)

For aero-acoustic measurement in open-jet automotive wind tunnel, the sound source information of the test vehicle cannot be accurately obtained due to the presence of jet shear layer and the airflow around the test vehicle. Through wind tunnel tests for a particular sound source and a test vehicle both positioned in the jet core of test section, the influences of jet flow and flow around the vehicle on sound propagation are studied. Result shows that the two flow fields all lead to significant acoustic drift phenomenon which should be taken into consideration in test data analysis. The linear velocity distribution profiles of jet shear layer and flow around the vehicle are acquired by fitting algorithm, and combined with geometric acoustics, a continuous multi-layer sound propagation model is established to predict the sound source drift quantities of jet flow and flow around vehicle. Comparison with test data indicates that the sound propagation model proposed in this paper can effectively predict the influence of medium motion on sound transmission direction; therefore, it can have widespread application in test data correction for automotive aero-acoustic wind tunnel test.

aero-acoustic wind tunnel; geometric acoustics; sound propagation; sound source drift quantity

O429

A

1000-3630(2017)-01-0057-07

10.16300/j.cnki.1000-3630.2017.01.011

2016-08-14;

2016-10-26

國家自然科學基金(51375342)資助項目。

王毅剛(1964－), 男, 陜西西安人, 博士, 教授, 研究方向為車輛氣動噪聲。

李方旭, E-mail: 103561@tongji.edu.cn