王超,岳廷瑞,萬振華,孫德軍
1.中國科學(xué)技術(shù)大學(xué) 近代力學(xué)系,合肥 230027 2.中國空氣動力研究與發(fā)展中心,綿陽 621000
方腔流動包含流致振蕩、聲不穩(wěn)定性和波的相互作用等豐富的流體動力學(xué)特征,被進行了廣泛研究。方腔流動問題最早被關(guān)注,是由于內(nèi)埋武器彈艙產(chǎn)生的強烈噪聲影響了飛行器結(jié)構(gòu)安全和武器精準度[1-2]。隨著航空運輸業(yè)的發(fā)展,飛機起飛和著陸階段的重要噪聲源—低速方腔—也得到了廣泛關(guān)注[3-4]。Roshko[5]報道了不同幾何形狀方腔的平均流動結(jié)構(gòu),發(fā)現(xiàn)方腔內(nèi)均存在回流結(jié)構(gòu)。Tracy[6]和Basley[7]等研究了馬赫數(shù)0.20~0.95方腔的流動特性,認為長深比(L/D)小于6的方腔為開式方腔,其前緣剪切層橫跨整個方腔并與后緣撞擊;長深比大于13時為閉式方腔,剪切層從前緣分離并在方腔底部再附。前人在開式方腔流動中發(fā)現(xiàn)了強烈的壓力振蕩和噪聲輻射,認為這是一個復(fù)雜的反饋過程。流動的不穩(wěn)定性導(dǎo)致方腔前緣剪切層內(nèi)渦卷起,在向下游傳播的過程中,渦結(jié)構(gòu)在黏性效應(yīng)作用下快速增長,渦與方腔后緣撞擊后輻射噪聲并產(chǎn)生壓力擾動前傳,從而形成閉環(huán),這一過程被稱為“剪切層模態(tài)”[8]。在前人研究結(jié)果的基礎(chǔ)上,Rossiter[9]對方腔流動反饋的影響參數(shù)進行了探索,給出了一個預(yù)測振蕩主頻率的半經(jīng)驗公式,但該公式未考慮三維參數(shù)影響。
Gharib等[10]研究了水中的方腔流動,發(fā)現(xiàn)了一種基于薄邊界層的流動模態(tài),即“尾跡模態(tài)”。在尾跡模態(tài)下,方腔流動的阻力和噪聲強度增大,而振蕩主頻不隨速度的增大而改變。Colonius等[11]在二維直接數(shù)值模擬(DNS)中也發(fā)現(xiàn)了尾跡模態(tài)。Shieh等[12]數(shù)值模擬了具有相同參數(shù)的二維和三維方腔流動,在二維情況下觀察到了剪切層模態(tài)和尾跡模態(tài),而在三維情況下僅觀察到剪切層模態(tài)。Martin等[13]研究了低馬赫數(shù)下三維開式方腔的流動和噪聲輻射,發(fā)現(xiàn)在剪切層模態(tài)下流動表現(xiàn)出部分三維性,在尾跡模態(tài)下流動表現(xiàn)出二維性,兩種流動模態(tài)的噪聲輻射有很大不同。
如前所述,方腔剪切層模態(tài)的振蕩機制主要是一個準二維和自持的過程,因此大部分實驗、理論和數(shù)值模擬研究工作都集中于探究雷諾數(shù)、馬赫數(shù)和長深比對反饋機制的影響,以降低方腔流致噪聲[14]。Block[15]、Ahuja[16]等對方腔寬長比(W/L)的影響進行了研究,前者發(fā)現(xiàn)隨著方腔寬度減小,輻射噪聲強度增大;而后者觀察到隨著方腔寬度減小,流致振蕩強度降低。產(chǎn)生這種差異的原因目前尚不明確,可能是馬赫數(shù)、雷諾數(shù)或邊界層厚度對流動結(jié)構(gòu)產(chǎn)生的影響所致。近年來,研究者開始關(guān)注三維方腔[17-18],但目前大多數(shù)研究集中于無限寬度(infinite width)方腔和固定寬度方腔,對方腔寬度影響的研究相對較少,特別是寬長比小于0.5的情況,幾乎沒有被研究過。
本文采用麥克風(fēng)陣列、脈動壓力傳感器和TRPIV(Time-Resolved Particle Image Velocimetry)等測量設(shè)備,開展側(cè)壁約束效應(yīng)對三維方腔自持振蕩和噪聲輻射影響的實驗研究,對比不同寬長比方腔的輻射噪聲強度、頻譜特性和指向性,結(jié)合方腔內(nèi)壁面壓力分布及流場結(jié)構(gòu),綜合分析減小寬長比對方腔流致振蕩及噪聲輻射的抑制機理。
實驗在中國空氣動力研究與發(fā)展中心0.55 m×0.40 m聲學(xué)風(fēng)洞中開展,如圖1所示。該風(fēng)洞為回流式低速風(fēng)洞,配備矩形試驗段,40 m/s風(fēng)速下的來流湍流度小于0.05%,80 m/s風(fēng)速下的背景噪聲小于78 dB(A)。
圖1 聲學(xué)風(fēng)洞實驗實物圖Fig.1 CARDC low-speed acoustic wind tunnel
矩形方腔實驗?zāi)P蜑橛袡C玻璃材質(zhì),通過平板(有機玻璃材質(zhì))與風(fēng)洞試驗段出口相連,平板上表面與試驗段出口下表面齊平。方腔沿試驗段中心線水平安裝,方腔上部開口所在平面與平板上表面齊平,前緣距平板前緣200 mm。為便于闡述,以方腔上部開口前緣中心點為原點O建立坐標系,沿流向為x軸正向,y軸與x軸垂直,指向上方,z軸垂直于xOy平面,其正向根據(jù)右手定則確定(圖2)。
圖2 實驗方案示意圖Fig.2 Test configuration of experiments
實驗測量了5個寬度不同的矩形方腔,其x向長度、y向深度均相同(長度L=200 mm,深度D=100 mm)。5個方腔的z向?qū)挾确謩e為20、40、60、80和100 mm,方腔的寬長比(W/L)在0.1~0.5之間變化。來流馬赫數(shù)Ma=0.03~0.25。
1.2.1 噪聲測量
噪聲測量使用B&K 4955麥克風(fēng),分別在水平面和垂直平面內(nèi)測量噪聲強度和指向性。如圖3所示,在y=200 mm水平面內(nèi),以方腔上部開口后緣中心點為圓心、1.5 m為半徑,左右兩側(cè)各均勻布置5個麥克風(fēng),彼此夾角15°;在z=0 mm的垂直平面內(nèi),以方腔上部開口后緣中心點為圓心、1.0 m為半徑,均勻布置8個麥克風(fēng),彼此夾角15°。麥克風(fēng)頻率響應(yīng)范圍為20 Hz~40 kHz,聲壓級(LSP)測量精度優(yōu)于0.1 dB,聲壓級公式為:
圖3 水平面和垂直平面的麥克風(fēng)陣列布置示意圖Fig.3 The arrangement of microphone arrays in horizontal and vertical directions
式中:p 為波動壓力;pref為參考聲壓,pref=2×10–5Pa。
噪聲測量采樣頻率25.6 kHz,采樣時間30 s。對原始數(shù)據(jù)進行分段處理,使用平均周期圖譜法計算功率譜密度[19],對每段數(shù)據(jù)加漢明窗函數(shù)[20],每個窗口的樣本數(shù)為25600,重疊50%,頻率分辨率為1.56 Hz。由于溫度、濕度以及麥克風(fēng)支撐裝置的影響,多次測量時發(fā)現(xiàn):對于垂直平面,噪聲測量誤差不大于0.2 dB,對于三腳架支撐的水平面,噪聲測量誤差不大于0.5 dB。
1.2.2 動態(tài)壓力測量
使用13個壓阻式脈動壓力傳感器(以P1~P13表示)測量壁面脈動壓力。脈動壓力傳感器表面直徑3.8 mm,最大量程為1 psi(約6.895 kPa),測量精度0.1%。如圖3所示,P1安裝于方腔上游,距離方腔前緣160 mm。其他12個傳感器安裝于方腔內(nèi)壁面,如圖4所示,P2~P4安裝于左側(cè)壁(迎來流左側(cè)),P5~P7安裝于右側(cè)壁,P8~P13安裝于底壁。所有傳感器均沿其所在壁面水平中心線等距安裝,并確保傳感器表面與方腔壁面齊平。Roshko[5]的實驗結(jié)果表明,即使傳感器與腔體壁之間存在間隙,傳感器安裝于方腔內(nèi)部對方腔流動也基本無影響。
圖4 方腔內(nèi)壁面脈動壓力傳感器布置示意圖Fig.4 The arrangement of high frequency dynamic pressure sensors on the cavity walls
脈動壓力傳感器采樣頻率51.2 Hz,采樣周期30 s,數(shù)據(jù)處理方法與麥克風(fēng)的數(shù)據(jù)處理方法相同,脈動壓力的功率譜密度采用系綜平均法計算[19],脈動壓力測量結(jié)果可以給出方腔內(nèi)壁面噪聲強度和平均壓力。
1.2.3 流場測量
使用高頻PIV測量流場(見圖5)。相機型號為PHOTRON SA-Z,雙曝光模式,配備85 mm尼康鏡頭,變焦比2.2,放大比例3.4 pixel/mm,視場最小分辨率0.3 mm。激光器為Nd:YLF激光器(30 mJ,雙脈沖)。示蹤粒子為甘油,直徑1 μm。采用同步控制器實現(xiàn)激光器和相機同步采集,采樣頻率1.8 kHz。以PIVTEC–PIVview2C軟件進行圖像分析,圖像處理采用帶窗口的迭代多重網(wǎng)格技術(shù)。
圖5 PIV測量實物圖Fig.5 PIV system in the acoustic anechoic room
使用PIV測量了Ma=0.12~0.25范圍內(nèi)5個不同寬長比的方腔中心平面(z=0 mm)的流場。為了比較近壁面流動與中心平面流動的差異,還測量了寬長比為0.5的方腔近壁面(z=40 mm)處的流場。
2.1.1 方腔振蕩及輻射噪聲主頻率分析
通過測量寬長比為0.5的方腔遠場噪聲,分析方腔流致振蕩及噪聲輻射特性。圖6給出了馬赫數(shù)從0.03增大至0.25過程中的方腔噪聲頻譜特性變化曲線。由圖6(a)可知,來流馬赫數(shù)在0.03~0.09之間時,噪聲頻譜波動出現(xiàn)于550 Hz附近,但波動幅值較小,沒有明顯的頻率尖峰,這說明方腔流動在較低風(fēng)速時即開始表現(xiàn)出振蕩和不穩(wěn)定性特征,但尚未形成自持振蕩。隨著來流馬赫數(shù)繼續(xù)增大,遠場噪聲頻譜在550 Hz附近開始出現(xiàn)明顯頻率尖峰,Ma=0.15時出現(xiàn)3個頻率尖峰(400、558和660 Hz),Ma=0.175時出現(xiàn)2個頻率尖峰(484和611 Hz),Ma=0.20時,僅有一個550 Hz頻率尖峰。噪聲頻譜峰值反映了方腔流致振蕩過程中的脫落渦動力學(xué)行為,該過程可視為一個多渦競爭自持振蕩形成過程。
圖6 不同風(fēng)速下的方腔遠場噪聲頻譜特性(W/L=0.5)Fig.6 PSD distributions for different wind speeds at W/L=0.5
同時可以看到,在來流馬赫數(shù)0.03~0.20范圍內(nèi),頻率尖峰中心頻率并未隨馬赫數(shù)增大而變化,說明在較低來流馬赫數(shù)下,方腔流致振蕩斯特勞哈爾數(shù)(Strouhal number)隨來流馬赫數(shù)不同而不同,對應(yīng)的流體動力學(xué)特征區(qū)別明顯,這一現(xiàn)象在以往方腔研究中并未引起注意。根據(jù)前人研究結(jié)果[21]可知,方腔遠場噪聲頻譜是方腔流動渦結(jié)構(gòu)的表征,前緣剪切層不穩(wěn)定性導(dǎo)致的渦卷起,隨著流動向下游移動,并與方腔后緣撞擊,產(chǎn)生噪聲輻射。本文認為:在低速條件下,脫落渦與方腔后緣撞擊產(chǎn)生的壓力擾動不能形成強相干運動,擾動能量不足以建立反饋回路,進一步證明了方腔自持振蕩是一個隨來流速度逐漸發(fā)展的過程。
當(dāng)來流馬赫數(shù)為0.20時,方腔噪聲尖峰頻率為550 Hz;來流馬赫數(shù)為0.25時,尖峰頻率為650 Hz。由圖6(b)可以看出,當(dāng)來流馬赫數(shù)大于0.20時,方腔輻射噪聲尖峰頻率隨著來流馬赫數(shù)的增大而增大,符合斯特勞哈爾數(shù)相似定理,說明此時方腔能夠形成穩(wěn)定的自持振蕩。為預(yù)測方腔噪聲主頻率,Rossiter[9]提出了一個半經(jīng)驗公式:
式中,Srn為振蕩模態(tài)n對應(yīng)的斯特勞哈爾數(shù);fn、L和u∞為方腔振蕩頻率、方腔長度和自由來流速度;α為渦運動與聲傳播的遲滯時間,k為渦對流速度。
式(2)中的經(jīng)驗參數(shù)α和k與方腔內(nèi)流動結(jié)構(gòu)和長深比相關(guān),是通過實驗數(shù)據(jù)擬合確定的[22],k的取值范圍在0.50~0.75之間,α的取值由長深比決定。對于長深比為2的方腔,k=0.5,α=0.25。用式(2)預(yù)測與本文實驗相同參數(shù)的方腔噪聲主頻率,得到來流馬赫數(shù)0.20和0.25時的主頻率分別為555 Hz和660 Hz,與實驗結(jié)果吻合得非常好,此時的方腔振蕩模態(tài)為Rossiter三階模態(tài)。
2.1.2 側(cè)壁約束效應(yīng)對方腔噪聲的抑制
對于寬長比為0.5的方腔,在來流馬赫數(shù)為0.20時,方腔流動已經(jīng)形成自持振蕩并輻射出尖頻噪聲。為研究側(cè)壁約束效應(yīng)對方腔流致振蕩的抑制作用,實驗對比了寬長比為 0.1、0.2、0.3、0.4和0.5的方腔遠場噪聲。圖7(a)給出了來流馬赫數(shù)0.20時不同寬長比方腔的遠場噪聲頻譜對比,當(dāng)寬長比W/L≤0.3時,方腔輻射噪聲的頻率尖峰明顯減弱甚至消除。由圖7(b)可以看出,當(dāng)來流馬赫數(shù)提高至0.25時也出現(xiàn)了相同現(xiàn)象。
圖7 不同寬長比方腔噪聲頻譜對比Fig.7 PSD comparison for different W/L ratios
圖8給出了方腔寬長比為0.3和0.4時的垂直和水平方向的噪聲總聲壓級(OASPL)對比結(jié)果??梢钥闯?,減小方腔寬度并不會改變水平方向的噪聲總聲壓級,但會大幅降低垂直方向的方腔上游總聲壓級,表現(xiàn)出明顯的側(cè)壁約束效應(yīng);與寬長比0.4的方腔相比,寬長比0.3的方腔上游噪聲強度最高能夠降低3 dB。對于長深比L/D≥1的方腔,產(chǎn)生的噪聲主要是縱波,具有向上游輻射的強烈指向性[21],因此,減小寬長比能夠顯著降低方腔流致振蕩輻射噪聲總聲壓級。
圖8 寬長比0.3與0.4的方腔噪聲總聲壓級對比(Ma=0.20)Fig.8 OASPL comparison between cases with W/L=0.3 and W/L=0.4 at a wind speed of Ma 0.20
前人對方腔寬度的研究多集中于無限寬度和寬長比≥0.5的單一狀態(tài),上述側(cè)壁約束效應(yīng)尚未被報道。方腔輻射噪聲與內(nèi)部流動結(jié)構(gòu)密切相關(guān),本文認為:當(dāng)方腔寬長比減小時,由于側(cè)壁的黏性效應(yīng),方腔流動的三維性[23]加強,內(nèi)部渦結(jié)構(gòu)的強度降低,導(dǎo)致方腔的自持振蕩和輻射噪聲受到抑制。后文將通過壁面脈動壓力和PIV流場測量結(jié)果進一步分析噪聲抑制機理。
圖9給出了來流馬赫數(shù)0.20、寬長比0.5條件下,方腔前緣及內(nèi)壁面壓力頻譜與遠場噪聲頻譜的對比(綠線為前緣測量結(jié)果,紅線為方腔壁面測量結(jié)果,藍線為遠場測量結(jié)果)。壁面脈動壓力傳感器測得的壓力頻譜和遠場噪聲頻譜具有相同的頻域分布及變化趨勢。由圖9可知,方腔內(nèi)壁面壓力頻譜的尖峰強度明顯高于方腔前緣,分析認為方腔上游的壓力脈動是由方腔內(nèi)輻射的二次波引起的。方腔前緣壓力頻譜在高頻段能量較高,是由于來流邊界層的高頻分量注入所導(dǎo)致。方腔內(nèi)部的壓力頻譜在高頻段存在多個尖峰頻率,分析認為是由于內(nèi)部回流誘導(dǎo)渦對流產(chǎn)生的,在遠場傳播中快速衰減,不能產(chǎn)生有效的噪聲輻射。
圖9 方腔內(nèi)壁面脈動壓力頻譜與遠場噪聲頻譜的對比(Ma=0.20,W/L=0.5)Fig.9 PSD distributions based on surface pressure fluctuations at local pressure sensors 1 and 2 and at far-field microphone 12 with a wind speed of Ma 0.20 and W/L=0.5
對脈動壓力數(shù)據(jù)進行時間平均,得到各點的平均壓力pmean,進而計算方腔壁面的時均壓力系數(shù)Cp。圖10給出了來流馬赫數(shù)0.20、寬長比為0.1和0.5的方腔壁面時均壓力系數(shù)Cp分布對比。圖中各點為實驗測量值,采用三次曲線對其進行擬合。由圖10(a)可以看出,寬長比減小,雙側(cè)壁Cp值減小,但沿流向的分布趨勢基本不變。圖10(b)則顯示:底壁的Cp值隨寬長比的減小而增大,但沿流向分布發(fā)生了明顯變化。由于方腔壁面時均壓力系數(shù)分布與內(nèi)部流動(特別是主回流渦結(jié)構(gòu))具有強相關(guān)性,因此由壁面壓力分布可知,方腔側(cè)壁約束效應(yīng)改變了內(nèi)部的流動及渦結(jié)構(gòu)。
圖10 寬長比0.1與0.5的方腔壁面時均壓力系數(shù)對比(Ma=0.20)Fig.10 Comparison of mean pressure coefficients between the cases with W/L=0.5 and W/L=0.1 at Ma 0.20
為進一步研究方腔寬度對流動結(jié)構(gòu)的影響,采用TR-PIV測量了不同寬度方腔內(nèi)各截面的流場結(jié)構(gòu)。圖11為來流馬赫數(shù)0.20、寬長比0.5時中心平面(z=0 mm)一個流致振蕩周期(T0)內(nèi)具有代表性的瞬時流線及渦量分布圖:方腔前緣初始產(chǎn)生一個小渦結(jié)構(gòu),隨來流向下游移動并逐漸增大,與方腔內(nèi)的大渦結(jié)構(gòu)融合并注入能量,形成一個更大的回流結(jié)構(gòu),最后與方腔后緣撞擊并輻射噪聲,部分渦結(jié)構(gòu)由后緣擠出,同時在方腔內(nèi)產(chǎn)生誘導(dǎo)渦,完成一個振蕩循環(huán)。該結(jié)果與前人對方腔流動的假設(shè)吻合。
圖11 一個振蕩周期內(nèi)方腔瞬時流場變化(Ma=0.20,W/L=0.5)Fig.11 Six snapshots of flow pattern in a period T0 in the central plane at a wind speed of Ma 0.20 and W/L=0.5.The streamlines are shown, while the background denotes the field of spanwise vorticity
考慮到方腔流動具有三維特性[23],對近壁面流場結(jié)構(gòu)進行了測量。圖12為來流馬赫數(shù)0.20、寬長比0.5時中心平面(z=0 mm)和近壁面(z=40 mm)的時均流場對比:中心平面流線分布以二維流動為主,而近壁面流線有著明顯的源流動,呈現(xiàn)出較高的三維特性;同時,近壁面的渦量值小于中心平面,回流強度被大幅抑制。
圖12 方腔中心平面和近壁面時均流場分布對比(Ma=0.20,W/L=0.5)Fig.12 Mean flow pattern comparison between central plane and near wall plane at the wind speed of Ma 0.20 (W/L=0.5)
通過對比不同寬長比方腔中心平面的流場結(jié)構(gòu),結(jié)合方腔近壁面流動特征,分析方腔側(cè)壁約束效應(yīng)對輻射噪聲的抑制機理。圖13為來流馬赫數(shù)0.20,寬長比 0.1、0.2、0.3、0.4和 0.5的方腔中心平面(z=0 mm)流場結(jié)構(gòu)??梢钥闯觯弘S著寬長比減小,方腔內(nèi)的主回流結(jié)構(gòu)逐漸向上游移動,中心平面流動逐漸表現(xiàn)出類似近壁面的三維流場特征。
圖13 不同寬長比方腔中心平面時均流場對比(Ma=0.20)Fig.13 The comparison of mean flow patterns in the central plane for different W/L at the wind speed of Ma 0.20
圖14給出了寬長比為0.3和0.4的方腔沿流向各截面的無量綱化速度型曲線(圖中,umax為方腔內(nèi)x方向最大流速,ux為方腔內(nèi)當(dāng)?shù)豿方向的實際速度)??梢钥闯觯S著寬長比減小,方腔內(nèi)回流速度降低,流動渦量減小。
圖14 寬長比0.3和0.4的方腔中心平面流向時均速度型對比(Ma=0.20)Fig.14 The mean streamwise velocity profiles inside the cavity for W/L=0.4 and W/L=0.3 in the central plane at a wind speed of Ma 0.20
綜上,本文認為:三維方腔側(cè)壁約束效應(yīng)之所以能夠抑制振蕩及輻射噪聲,是由于方腔寬度減小時,回流渦結(jié)構(gòu)向上游移動、流動呈現(xiàn)三維特性、渦量值減小等多種因素共同作用導(dǎo)致的。此時,渦與方腔后緣撞擊強度降低,自持振蕩被抑制甚至消除。
1)由于前緣剪切層的不穩(wěn)定性,來流馬赫數(shù)大于0.03時,方腔就開始出現(xiàn)壓力波動。隨著來流馬赫數(shù)增大,方腔內(nèi)部通過多渦競爭機制逐漸形成自持振蕩。來流馬赫數(shù)大于0.20時,方腔出現(xiàn)Rossiter三階模態(tài)流致振蕩,并輻射出指向方腔上游的尖頻噪聲,噪聲頻率變化滿足斯特勞哈爾數(shù)相似準則。
2)對于自持振蕩方腔流動,減小方腔寬長比能夠大幅度抑制方腔振蕩及輻射的尖頻噪聲。來流馬赫數(shù)0.20時,將方腔寬長比由0.4降至0.3,方腔上游噪聲總聲壓級能夠降低3 dB以上。
3)方腔壁面壓力分布和PIV流場結(jié)果表明:減小方腔寬長比,加強側(cè)壁約束效應(yīng),能夠改變方腔流動結(jié)構(gòu),增強方腔流動的三維性,降低渦強度;同時,還能夠促使方腔內(nèi)大渦結(jié)構(gòu)向前緣移動,削弱渦與方腔后緣的撞擊程度,進而抑制方腔流致振蕩和尖頻噪聲輻射。