流激孔腔噪聲特征及控制方法研究

劉璐璐，呂世金，劉進(jìn)

（中國(guó)船舶科學(xué)研究中心船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，江蘇無(wú)錫214082）

流激孔腔噪聲特征及控制方法研究

劉璐璐，呂世金，劉進(jìn)

（中國(guó)船舶科學(xué)研究中心船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，江蘇無(wú)錫214082）

文章利用大渦模擬方法，建立了流激孔腔自噪聲及輻射噪聲預(yù)報(bào)方法，分析了腔深、來(lái)流速度及孔腔流向尺寸對(duì)剪切振蕩特征的影響，并給出了流激孔腔近場(chǎng)輻射噪聲特性。在此基礎(chǔ)上探索了流激孔腔噪聲控制方法，研究了孔腔開口加格柵對(duì)流激孔腔渦流場(chǎng)和噪聲特征的影響。

大渦模擬；流激孔腔；剪切振蕩；格柵

0 引言

水下航行體由于潛浮的需要，開設(shè)了數(shù)量眾多、樣式各異的流水孔和壓載水艙開口以及其它大小切口。水下航行水流經(jīng)過(guò)開口表面時(shí)，由于表面的不連續(xù)性，流動(dòng)發(fā)生分離，沿流向跨越孔腔的內(nèi)、外流之間，存在不穩(wěn)定的剪切層波動(dòng)，一方面直接輻射噪聲，另一方面會(huì)在流水孔后緣局部區(qū)域造成湍流邊界層脈動(dòng)壓力增大，激勵(lì)航行體輕外殼振動(dòng)從而產(chǎn)生二次輻射噪聲。此外，對(duì)于采用雙層殼體的航行體，開口在輕外殼與耐壓殼之間形成孔腔，聲波在有限孔腔中傳播會(huì)形成駐波，其駐波共振頻率與孔腔尺寸有關(guān)，如果開孔的剪切振蕩頻率與孔腔駐波共振頻率相一致，即形成耦合，會(huì)產(chǎn)生強(qiáng)烈的聲輻射。這種聲輻射表現(xiàn)為若干個(gè)以耦合共振頻率為中心頻率的窄帶噪聲[1]。因而研究孔腔流動(dòng)與發(fā)聲機(jī)理并探索流激孔腔噪聲控制方案具有重要意義。早期的研究基于試驗(yàn)測(cè)試為主，隨著數(shù)值模擬方法尤其是大渦模擬方法的出現(xiàn)和發(fā)展，使得從數(shù)值模擬角度研究孔腔流激噪聲成為可能。

近年來(lái)，國(guó)內(nèi)外研究者針對(duì)流激孔腔噪聲控制方式進(jìn)行了大量的探索。被動(dòng)流動(dòng)控制技術(shù)由于結(jié)構(gòu)簡(jiǎn)單、應(yīng)用方便成為最早被研究的流動(dòng)控制方法。Kuo[2]采用激光多普勒測(cè)速儀及激光切片技術(shù)，在水筒中測(cè)量了傾斜腔底對(duì)空腔剪切層脈動(dòng)特征的影響。結(jié)果表明，低Re數(shù)情況下，在一定傾斜范圍內(nèi)，傾斜腔底對(duì)空腔剪切層振蕩產(chǎn)生調(diào)制作用，可明顯降低腔內(nèi)振蕩幅度，傾斜達(dá)到一定程度，自激振蕩完全抑制。Zhang和Rona[3]設(shè)計(jì)了一種在空腔前緣放置的斜坡塊，對(duì)空腔噪聲具有一定的抑制效果。但這種方法與來(lái)流的馬赫數(shù)有很大關(guān)聯(lián)，因此不同高度、角度的斜坡塊只能對(duì)特定流速的流體有控制效果。Tzborda[4]在空腔內(nèi)部添加不同高度的柵欄并進(jìn)行了大量的風(fēng)洞試驗(yàn)，試驗(yàn)結(jié)果表明這種控制方法能夠有效降低腔內(nèi)前后部的壓力差。Arunajatesan[5]采用數(shù)值模擬方法，計(jì)算分析了格柵肋條對(duì)剪切層振蕩的影響。結(jié)果表明孔腔開口加肋條后，開口仍然有剪切層振蕩現(xiàn)象出現(xiàn)，清晰可見類似于共振腔的大尺度旋渦，但是，格柵肋條對(duì)腔頸部位的流場(chǎng)有明顯的調(diào)制作用。朱習(xí)劍和何祚鏞[6]在水筒中進(jìn)行了變深度的、不帶導(dǎo)流板的突出矩形腔的振蕩和聲輻射研究。吳亞軍[7]在孔腔上游安裝了圓柱形、矩形、鋸齒和三角形四種形式的擾流體，控制孔腔的流激振蕩，對(duì)窄帶或?qū)拵г肼曈幸欢ǖ囊种菩Ч?/p>

本文利用大渦模擬方法，建立流激孔腔自噪聲及輻射噪聲預(yù)報(bào)方法，數(shù)值分析孔腔開口噪聲特征,考察腔深、來(lái)流速度、流向尺度等參數(shù)對(duì)流激孔腔剪切振蕩的影響以及流激孔腔近場(chǎng)輻射噪聲特征，并探索孔腔開口加格柵結(jié)構(gòu)對(duì)流激孔腔渦流場(chǎng)和噪聲特征的影響規(guī)律。

1 流激孔腔噪聲計(jì)算方法

1.1 計(jì)算方法

1.1.1 大渦模擬方法

流激孔腔流動(dòng)特征通過(guò)不可壓縮流動(dòng)方程求解。大渦模擬的理論建立在兩個(gè)基本假設(shè)之上：第一個(gè)是湍流的平均特性，主要由大尺度湍流運(yùn)動(dòng)來(lái)控制，幾乎不受小尺度湍流運(yùn)動(dòng)的影響。第二個(gè)假設(shè)是小尺度湍流，特別在高雷諾數(shù)下，表現(xiàn)出各向同性的特點(diǎn)。通過(guò)將非定常的N-S方程在波數(shù)空間或物理空間上進(jìn)行濾波，得到大渦模擬的控制方程。濾波過(guò)程有效的過(guò)濾掉了那些尺度小于濾波寬度（或網(wǎng)格尺度）的小渦。

濾波變量由下式定義：

其中：D為流體域，G為決定可求解渦尺度的濾波函數(shù)?？蛇x用不同的濾波函數(shù)，本文采用網(wǎng)格體積濾波。

濾波后的連續(xù)性方程和N-S方程可以表示為：

其中：σij為分子粘性引起的應(yīng)力張量，τij為亞格子應(yīng)力，需用亞格子渦模型進(jìn)行模擬：

Smagorinsky模型由Smagorinsky[8]于1963年提出。該模型基于小尺度渦是局部平衡狀態(tài)的假設(shè)，因此可用渦粘形式的亞格子雷諾應(yīng)力模型表示：

其中：Cs為Smagorinsky常數(shù)，動(dòng)態(tài)Smagorinsky模型采用多重網(wǎng)格濾波的方法將參數(shù)Cs與計(jì)算過(guò)程的湍流物理量的變化相關(guān)聯(lián)，通過(guò)大渦模擬本身的計(jì)算結(jié)果，可以在每個(gè)空間網(wǎng)格點(diǎn)、每個(gè)時(shí)間步計(jì)算模型參數(shù)Cs，完全擺脫了一般理論中的經(jīng)驗(yàn)參數(shù)。此時(shí)Cs是隨湍流物理量的時(shí)間和空間變化的函數(shù)。

1.1.2 聲學(xué)類比方程

流體的擾動(dòng)產(chǎn)生聲場(chǎng)，聲場(chǎng)模擬采用聲學(xué)類比方法。聲學(xué)類比方法就是對(duì)流動(dòng)和噪聲分別處理，將聲波的產(chǎn)生與傳播進(jìn)行解耦，即從流動(dòng)的數(shù)值模擬中計(jì)算出聲源項(xiàng)，忽略聲波的傳播與反饋對(duì)聲源的干擾，然后通過(guò)波動(dòng)方程的解析解亦即格林積分公式來(lái)預(yù)報(bào)遠(yuǎn)場(chǎng)噪聲。

Ffowcs Williams與Hawkings聲學(xué)類比方程，即FW-H方程[9]表達(dá)如下：

其中：ui為xi方向的流體速度分量；un為垂直于物體表面（f＝0）的流體速度分量；vi為xi方向的物體表面速度分量；vn為物體表面法向速度分量；δ（f）為Dirac delta函數(shù)；H（f）為Heaviside階躍函數(shù)。p′為遠(yuǎn)場(chǎng)聲壓。f=0表示物體表面，f＞0表示外部無(wú)界的自由空間。ni為物面外法線，指向流體內(nèi)部。c0為遠(yuǎn)場(chǎng)聲速，ρ0為遠(yuǎn)場(chǎng)密度，Tij為L(zhǎng)ighthill應(yīng)力張量。

利用自由空間格林函數(shù)δ（g）/4πr，結(jié)合Kirchhoff積分可求出FW-H方程的遠(yuǎn)場(chǎng)解。遠(yuǎn)場(chǎng)噪聲中包含單極子噪聲pT′（x,t），偶極子噪聲pL′（x,t）以及四極子噪聲pQ′（x,t）。

1.2 計(jì)算模型

本文在大渦模擬和聲學(xué)類比方法的理論框架下，結(jié)合動(dòng)態(tài)Smagorinsky亞格子應(yīng)力模型，計(jì)算了三維孔腔的流激噪聲，并與試驗(yàn)結(jié)果進(jìn)行了比較，驗(yàn)證了流激噪聲預(yù)報(bào)方法的可行性。在此基礎(chǔ)上探究了流激孔腔噪聲控制方案，并評(píng)價(jià)了噪聲控制效果。

圖1 開口空腔結(jié)構(gòu)示意圖Fig.1 Schematic of flow-induced cavity

圖2 計(jì)算模型網(wǎng)格劃分圖Fig.2 The computational domain and grid

研究的流激孔腔開口長(zhǎng)320 mm，寬194 mm，深15 mm，開口下方矩形腔體深度785 mm，長(zhǎng)360 mm，寬215 mm，自噪聲監(jiān)測(cè)點(diǎn)如圖1。數(shù)值模擬計(jì)算域長(zhǎng)和寬均為4 m，高2 m，進(jìn)口采用自由來(lái)流條件，來(lái)流速度為6 m/s，出口為壓力出口，其它均為無(wú)滑移壁面。網(wǎng)格劃分圖如圖2所示。亞格子應(yīng)力模型采用Dynamic Smagorinsk-Lilly模型，動(dòng)量方程的離散采用有限體積法，壓力速度耦合算法采用SIMPLE算法，時(shí)間項(xiàng)的離散采用二階隱式差分格式，而空間項(xiàng)的離散采用二階精度的限界中心差分格式。計(jì)算時(shí)間步長(zhǎng)為0.000 2 s，分析頻率為2 500 Hz。

1.3 計(jì)算方法驗(yàn)證

為了驗(yàn)證計(jì)算模型的可靠性，設(shè)計(jì)開口尺寸可變孔腔模型，可實(shí)現(xiàn)0～320 mm范圍任意選取。流激孔腔噪聲測(cè)試在中國(guó)船舶科學(xué)研究中心小型高速空泡水筒進(jìn)行。自噪聲測(cè)試?yán)梅叫嗡?25 mm× 225 mm×845 mm長(zhǎng)試驗(yàn)段，開口腔體換掉原有機(jī)玻璃觀察窗，開口寬度設(shè)計(jì)為194 mm，長(zhǎng)320 mm，開口下方腔體總深度800 mm，長(zhǎng)度360 mm，寬度215 mm，如圖3所示，在腔體內(nèi)與數(shù)值模擬監(jiān)測(cè)點(diǎn)相同位置處布置4個(gè)水聽器，研究開口附近及腔體內(nèi)不同位置的噪聲特征。

圖3 開口腔體實(shí)圖Fig.3 Experimental set-up of cavity noise measurement

圖4 流激孔腔噪聲窄帶頻譜預(yù)報(bào)與試驗(yàn)結(jié)果比較Fig.4 Comparison of computed and experimental frequency spectra for flow-induced cavity

Rossiter[10]給出了辨識(shí)流動(dòng)模態(tài)振蕩頻率的經(jīng)驗(yàn)公式為

其中：fn為模態(tài)振蕩頻率、n為模態(tài)數(shù)、U0為遠(yuǎn)方來(lái)流速度、M為馬赫數(shù)、L為孔腔開口長(zhǎng)度。

圖4給出了水速為6 m/s時(shí)同一監(jiān)測(cè)點(diǎn)處（P1）320 mm開口流激孔腔自噪聲計(jì)算與試驗(yàn)結(jié)果比較。表1給出了本文針對(duì)開口長(zhǎng)度為320 mm的孔腔剪切振蕩頻率的試驗(yàn)、計(jì)算和理論結(jié)果比較，可以看出計(jì)算和試驗(yàn)結(jié)果剪切振蕩頻率基本一致，總聲級(jí)偏差小于4 dB。這表明本文計(jì)算結(jié)果可信，可用來(lái)分析流激孔腔剪切振蕩特征。

表1 剪切振蕩頻率的計(jì)算值與試驗(yàn)、理論值比較Tab.1 Comparison of computed,experimental and theoretical shear oscillation frequency

2 流激孔腔噪聲特征分析

流激孔腔開口處流體的剪切振蕩產(chǎn)生孔腔線譜噪聲，因此考察孔腔流場(chǎng)特性是噪聲特征分析的第一步。流激孔腔中縱面展向渦量場(chǎng)如圖5所示。非定常繞流剪切層從孔腔開口上游角區(qū)開始分離并發(fā)展，剪切層中的渦發(fā)生對(duì)流傳遞，到達(dá)隨邊附近時(shí)擊打隨邊，使隨邊產(chǎn)生壓力脈沖。壓力脈沖以聲速向上游傳播反饋到導(dǎo)邊。滿足一定相位條件時(shí)，剪切層擾動(dòng)形成閉合反饋環(huán)，產(chǎn)生自持振蕩[1]。振蕩的流體作為一種等效聲源，在一定條件下有可能產(chǎn)生單調(diào)音或者寬帶噪聲。

圖5 320 mm開口渦量場(chǎng)Fig.5 Vorticity field of 320 mm cavity

2.1 流激孔腔自噪聲特征分析

為了分析腔深對(duì)孔腔自噪聲功率譜的影響，分別選取腔深為400 mm、800 mm及1 200 mm。針對(duì)開口長(zhǎng)度為50 mm的孔腔，來(lái)流速度為6 m/s時(shí)不同腔深自功率譜計(jì)算結(jié)果比較如圖6所示。腔深從400mm增加至800 mm，剪切振蕩頻率不變，一階剪切振蕩頻率峰值降低。若繼續(xù)增加腔深到1 200 mm，剪切振蕩頻率和峰值不再變化，自功率譜曲線重合。表明腔深大于800 mm時(shí)，流激孔腔自功率譜不再受腔深影響。故下文計(jì)算均取腔深800 mm。

圖6 不同腔深對(duì)孔腔自功率譜的影響Fig.6 Influence of cavity depth on flow-induced cavity noise

來(lái)流速度對(duì)流激孔腔自噪聲的影響如圖7，分別考察了4.5 m/s、6 m/s和9 m/s來(lái)流速度下，320 mm開口孔腔的自噪聲（P1）頻譜曲線。隨著流速增加，流激孔腔剪切振蕩峰值頻率向高頻移動(dòng)。來(lái)流速度從4.5 m/s增至9 m/s，對(duì)應(yīng)一階剪切振蕩頻率從8 Hz增加到16 Hz，二階剪切振蕩頻率從16 Hz向高頻偏移至32 Hz。此外，5～500 Hz頻段總聲級(jí)分別為166 dB、170 dB和175 dB，流速每增加一倍，總聲級(jí)增加約9 dB。

為了分析開口流向尺寸對(duì)流激孔腔噪聲的影響，計(jì)算了20 mm、50 mm、100 mm、150 mm、200 mm和320 mm等幾種開口附近自噪聲。圖8為6 m/s時(shí)不同流向尺寸開口附近（P1）自功率譜級(jí)比較曲線。隨著開口尺寸增加，流激孔腔自噪聲線譜峰值頻率向低頻偏移。開口尺寸從20 mm增加至320 mm，對(duì)應(yīng)的一階剪切振蕩頻率從140 Hz向低頻偏移至11 Hz，二階剪切振蕩頻率從281 Hz減小至21 Hz。由此可見開口流向尺寸越高，孔腔剪切振蕩頻率越低。

圖7 來(lái)流速度對(duì)流激孔腔自噪聲的影響Fig.7 Influence of inflow velocity on flow-induced cavity noise

圖8 6 m/s時(shí)不同開口尺寸對(duì)孔腔附近自噪聲的影響Fig.8 Influence of streamwise length on flow-induced cavity noise

圖9 U0/L為定值時(shí)流激孔腔自功率譜計(jì)算結(jié)果Fig.9 Comparison of auto-spectrum for flow-induced cavity when U0/L=constant

圖9給出了U0/L為定值情況下三個(gè)不同開口尺寸（150 mm、200 mm和320 mm）、對(duì)應(yīng)不同流速（2.1 m/s、2.8 m/s和4.5 m/s）的孔腔自功率譜曲線計(jì)算結(jié)果。可以看出剪切振蕩峰值頻率重合，由此可以推斷流激孔腔剪切振蕩頻率與U0/L成比例。該結(jié)論也與Rossiter[10]提出的流動(dòng)模態(tài)振蕩頻率經(jīng)驗(yàn)公式相吻合。在剪切振蕩峰值頻率一致的情況下，隨著流速增加，剪切振蕩峰值增加，自功率譜總聲級(jí)也不斷增加。

2.2 流激孔腔近場(chǎng)輻射噪聲特征分析

對(duì)流激孔腔的近場(chǎng)輻射噪聲計(jì)算監(jiān)測(cè)點(diǎn)如圖10所示，在孔腔開口中心正上方流向-法向平面R=1.7 m和R=1 m處分別布置五個(gè)聲壓接收器，考察6 m/s時(shí)流激孔腔輻射噪聲變化規(guī)律。

圖11-12分別給出了流向-法向平面距50 mm孔腔開口中心不同距離、不同方向處的近場(chǎng)輻射噪聲?？梢钥闯?，f＜100 Hz時(shí)流激孔腔輻射噪聲在不同方向上計(jì)算結(jié)果有顯著差異，表明流激孔腔輻射噪聲具有指向性，且主要集中在低頻區(qū)域。此外，近場(chǎng)輻射噪聲線譜峰值頻率與流激孔腔自噪聲一致，且線譜峰值強(qiáng)度較低，沒有出現(xiàn)強(qiáng)烈的振蕩，可以推斷腔口剪切振蕩和空腔聲模態(tài)沒有耦合共振。事實(shí)上本文數(shù)值模擬過(guò)程假定腔體壁面是剛性，此時(shí)由于聲波波長(zhǎng)較長(zhǎng)，水介質(zhì)空腔中的聲模態(tài)頻率遠(yuǎn)高于腔口流體動(dòng)力振蕩頻率，腔口剪切振蕩與腔體聲模態(tài)產(chǎn)生耦合的可能性很小，一般不會(huì)出現(xiàn)較強(qiáng)的低頻線譜噪聲。而實(shí)際腔體壁面是彈性的，彈性壁面提供的附加柔度會(huì)降低腔體的聲模態(tài)頻率，使得腔口剪切振蕩頻率與空腔的模態(tài)頻率耦合的可能性增大。

圖10 輻射噪聲計(jì)算監(jiān)測(cè)點(diǎn)（流向—法向平面）Fig.10 Monitoring points for radiated noise calculation (streamwise-wall-normal plane)

圖11 距50 mm孔腔開口中心R=1.7 m處中縱面不同方向測(cè)點(diǎn)的輻射噪聲Fig.11 Radiated noise of different directions in streamwise-wall-normal plane at R=1.7 m for 50 mm long cavity

圖12 距50 mm孔腔開口中心R=1 m處中縱面不同方向測(cè)點(diǎn)的輻射噪聲Fig.12 Radiated noise of different directions in streamwise-wall-normal plane at R=1 m for 50 mm long cavity

3 格柵對(duì)流激孔腔噪聲特征的影響

常用的流激孔腔剪切振蕩及聲輻射控制方式主要包括在腔口導(dǎo)邊布放擾流器及流動(dòng)擴(kuò)散器，或在隨邊布放流動(dòng)偏轉(zhuǎn)器。對(duì)流激孔腔進(jìn)行格柵處理可以改變孔腔流體的流動(dòng)特性和相位關(guān)系，并進(jìn)一步導(dǎo)致孔腔的聲場(chǎng)特性發(fā)生變化。本文探索了開口加格柵對(duì)流激孔腔剪切振蕩及聲輻射的控制效果。

3.1 格柵形式及參數(shù)

本文對(duì)320 mm開口孔腔進(jìn)行加格柵處理，分直立格柵和斜立格柵兩種方案，格柵間距均為50 mm。開口加直立格柵和斜格柵示意圖分別如圖13、圖14所示，模擬水速6 m/s時(shí)直立格柵和斜格柵對(duì)流激孔腔線譜噪聲的影響。

圖13 320 mm開口加直立格柵示意圖Fig.13 Schematic of 320 mm long cavity with upright grille

圖14 320 mm開口加斜格柵示意圖Fig.14 Schematic of 320 mm long cavity with oblique grille

3.2 格柵對(duì)流激孔腔渦流場(chǎng)的影響

要考察格柵對(duì)流激孔腔噪聲特征的影響，首先需要清楚格柵結(jié)構(gòu)對(duì)孔腔流場(chǎng)特性的調(diào)制作用。320 mm大開口孔腔加直立格柵和斜立格柵后的展向渦量場(chǎng)分別如圖15～16。格柵的存在使孔腔開口附近的大尺度渦分解為小尺度渦，且孔腔加斜格柵結(jié)構(gòu)時(shí)剪切層中的渦在下游壁面處破碎。如前文所述，流激孔腔剪切振蕩的成因在于孔腔內(nèi)部剪切層的不穩(wěn)定性，且這種不穩(wěn)定性會(huì)由于下游壁面的反饋機(jī)制被增強(qiáng)，下游角區(qū)的干擾以回饋環(huán)的形式反作用于上游角區(qū)邊界層分離點(diǎn)。斜立格柵相比于直立格柵角區(qū)平緩過(guò)渡，剪切層中的渦發(fā)生對(duì)流傳遞到達(dá)隨邊并擊打傾斜的下游壁面，壁面產(chǎn)生的壓力脈沖較小，下游壁面的反饋機(jī)制較弱，因而流激孔腔剪切振蕩強(qiáng)度較小。由此可以斷定，斜立格柵對(duì)流激孔腔渦流場(chǎng)的調(diào)制作用效果優(yōu)于直立格柵。

圖15 320 mm開口加直立格柵對(duì)應(yīng)的渦量場(chǎng)Fig.15 Vorticity field of 320 mm cavity with upright grille

圖16 320 mm開口加斜立格柵對(duì)應(yīng)的渦量場(chǎng)Fig.16 Vorticity field of 320 mm cavity with oblique grille

3.3 格柵對(duì)流激孔腔噪聲特征的影響

圖17給出了兩種格柵對(duì)320 mm開口孔腔結(jié)構(gòu)自功率譜的影響，320mm開口流激孔腔及加直立格柵和斜立格柵時(shí)自功率譜在5～500 Hz頻段總聲級(jí)分別為170 dB、161 dB和152 dB?？浊患痈駯攀沟迷趂＜50 Hz頻段范圍內(nèi)自噪聲功率譜曲線大幅下降。此外，加直立格柵后剪切振蕩峰值頻率向高頻偏移，一階峰值頻率從11 Hz向高頻偏移至66 Hz，二階峰值頻率從21 Hz增大至124 Hz，且峰值相比于未加格柵情況顯著下降。加斜立格柵情況下開口孔腔自功率譜曲線較平緩，沒有峰值出現(xiàn)，斜立格柵對(duì)開口孔腔自噪聲的調(diào)制效果較好，這與前文對(duì)渦流場(chǎng)的分析結(jié)果相吻合。

對(duì)320 mm大開口孔腔加均勻等間距格柵時(shí)自噪聲頻譜曲線與50 mm開口孔腔的比較如圖18所示。320 mm大開口孔腔加直立格柵對(duì)應(yīng)的自噪聲頻譜與50 mm開口孔腔幾乎一致，峰值頻率一一對(duì)應(yīng)，但峰值略有降低，而加斜格柵時(shí)自噪聲頻譜曲線峰值消失。說(shuō)明對(duì)于水下航行體上的大開口，格柵相當(dāng)于將大開口變?yōu)榈雀駯砰g距的多個(gè)小開口。因此在開口尺寸不可調(diào)的前提下，可通過(guò)加直立格柵或者斜格柵控制流激噪聲。

圖17 兩種格柵對(duì)320 mm開口自功率譜的影響Fig.17 Influence of grille on auto-power spectra of 320 mm long cavity

圖18 兩種格柵與50 mm開口自功率譜比較Fig.18 Comparison of auto-power spectra between grilled 320 mm long cavity and 50 mm long cavity

進(jìn)一步考察兩種格柵對(duì)320 mm開口孔腔結(jié)構(gòu)近場(chǎng)輻射噪聲的影響，如圖19所示。從圖中可以看出兩種格柵對(duì)流激孔腔輻射噪聲的控制效果無(wú)顯著差異，計(jì)算可知流激孔腔加格柵使近場(chǎng)輻射噪聲總聲級(jí)在5～500 Hz頻段下降約14 dB。

圖19 兩種格柵對(duì)320 mm開口輻射噪聲的影響Fig.19 Influence of grille on radiated noise power spectra of 320 mm long cavity

4 結(jié)論

本文利用大渦模擬數(shù)值方法，建立了流激孔腔自噪聲及輻射噪聲預(yù)報(bào)方法，并與試驗(yàn)結(jié)果吻合較好；分析腔深、來(lái)流速度和流向尺寸對(duì)流激孔腔自噪聲功率譜的影響以及近場(chǎng)輻射噪聲特性；研究了孔腔開口加不同格柵形式對(duì)孔腔渦流場(chǎng)和噪聲特征的影響，得到以下主要結(jié)論：

（1）流激孔腔剪切振蕩頻率與U0/L成比例。近場(chǎng)輻射噪聲線譜峰值頻率與流激孔腔自噪聲一致，且線譜峰值強(qiáng)度較低，沒有出現(xiàn)較強(qiáng)的低頻線譜噪聲。

（2）對(duì)大開口加直立格柵使剪切振蕩頻率向格柵間距開口剪切振蕩頻率偏移，且峰值下降；斜格柵對(duì)孔腔開口流場(chǎng)的調(diào)制作用效果較好，噪聲控制效果優(yōu)于直格柵，使流激孔腔自功率譜曲線較平緩，剪切振蕩峰值消失。流激孔腔加格柵使近場(chǎng)輻射噪聲總聲級(jí)在5～500 Hz頻段下降約14 dB。

本文對(duì)流激孔腔噪聲特征的數(shù)值模擬在假定腔體為剛性壁的前提下腔口剪切振蕩和腔體聲模態(tài)沒有耦合共振。未來(lái)進(jìn)一步的研究方向是重視腔壁彈性振動(dòng)與空腔聲共振和腔口剪切層振蕩之間的耦合作用，嘗試將腔壁彈性納入數(shù)值模擬過(guò)程，探索腔口剪切振蕩與腔體聲駐波的耦合共振，并通過(guò)調(diào)整水下空腔彈性壁參數(shù)來(lái)抑制腔口剪切振蕩與空腔聲模態(tài)水彈性耦合共振，多層面控制水下流激孔腔的線譜噪聲。

[1]俞孟薩,張錚錚,高巖.開口與空腔流激聲共振及聲輻射研究綜述[J].船舶力學(xué),2015,19(11)：1422-1430. Yü Mengsa,Zhang Zhengzheng,Gao Yan.A review of acoustic resonant and radiation of aperture and cavity by flow excitation[J].Journal of Ship Mechanics,2015,19(11)：1422-1430.

[2]Kuo C H,Huang S H.Influence of flow path modification on oscillation of cavity sheer layer[J].Experiments in Fluids, 2001,31：162-178.

[3]Zhang X,Rona A,Edwards J A．An observation of pressure waves around a shallow cavity[J]．Journal of Sound and Vibration,1998,214(4)：771-778．

[4]Taborda N,Bray D,Knowles K．Passive control of cavity resonances in tandem configurations[C].Proc 31st AIAA Fluid Dynamics Conference,2001．

[5]Arumajatesan S,Sinha N.Modeling approach for reducing Helmholtz resonance in submarine structures[R].AO-A458237, 2005.

[6]朱習(xí)劍,何祚鏞.水洞中突出矩形腔的流激駐波振蕩研究[J].哈爾濱船舶工程學(xué)院學(xué)報(bào),1993,14(4),41-52. Zhu Xijian,He Zuoyong.Study of flow-induced standing wave resonance of rectangular cavity in water tunnel[J].J of Harbin Shipbuilding Engineering Institute,1993,14(4)：41-52.

[7]吳亞軍,等.基于被動(dòng)控制的空腔脈動(dòng)實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào),2013,34(9)：1640-1644.Wu Yajun,Ouyang Hua,Huang You.Experimental investigation on the fluctuating pressure of cavity based on passive controls[J].J of Engineering Thermo Physics,2013,34(9)：1640-1644.

[8]Smagorinsky J.General circulation experiments with the primitive equations.I.The basic experiment[J].Month.Wea.Rev., 1963,91：91.

[9]Ffowcs Williams J E,Hawkings D L.Sound generation by turbulence and surfaces in arbitrary motion[C]//Proc.Roy.Soc.. London,264 A,1969：321-342.

[10]Rossiter J E.Wind tunnel experiments on the flow over rectangular cavities at subsonicand transonic speeds[J].Rep.Mem. 3438,Aeronautical Research Council,1964.

[11]Crighton D G,Dowling A P,Ffowcs Williams J E,Heckl M,Leppington F G.Modern methods in analytical acoustics[M]. Springer-Verlag London Limited,1992.

[12]張楠.孔腔流動(dòng)和流激噪聲機(jī)理及耦合計(jì)算方法研究[D].無(wú)錫：中國(guó)船舶科學(xué)研究中心,2010.

[13]劉進(jìn),呂世金,沈琪.流激開孔空腔噪聲數(shù)值預(yù)報(bào)方法研究[C].水下噪聲會(huì)議論文集,2015.

Characteristics and control of cavity noise induced by flow excitation

LIU Lu-lu,Lü Shi-jin,LIU Jin
(National Key Laboratory on Ship Vibration&Noise,China Ship Scientific Research Center,Wuxi 214082,China)

Using large eddy simulation method,a prediction method for flow induced cavity noise was established.The inlet velocity,cavity depth and streamwise length influence on shear oscillation characteristics were analyzed,and nearfield radiate noise characteristics were given.A research on noise control scheme was performed and the influence of grille on flow induced cavity noise characteristics was studied.

large eddy simulation;flow-induced cavity;shear oscillation;grille

O352

A

10.3969/j.issn.1007-7294.2017.04.014

1007-7294（2017）04-0493-10

2017-01-22

劉璐璐（1990-），女，工程師，E-mail：LL300100@aliyun.com；呂世金（1973-），男，研究員。