樊林旭，張 磊，陳德喜，劉國慶，王赤忠

(1. 華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074；2. 中國艦船研究設(shè)計(jì)中心，湖北 武漢 430064)

水下潛器首部聲基陣區(qū)流噪聲研究方法

樊林旭1，張 磊1，陳德喜2，劉國慶1，王赤忠1

(1. 華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074；2. 中國艦船研究設(shè)計(jì)中心，湖北 武漢 430064)

水下潛器在航行過程中，主要使用聲吶來探測(cè)敵方艦艇以及自身位置。主聲吶一般安放在潛艇首部位置，可分為主動(dòng)聲吶與被動(dòng)聲吶。在探測(cè)目標(biāo)時(shí)，噪聲是 2 種聲吶系統(tǒng)都必須克服的干擾因素。針對(duì)自噪聲中的流噪聲，首先使用LES模型對(duì)水下潛器的外流場(chǎng)進(jìn)行仿真計(jì)算。在獲得流場(chǎng)中的脈動(dòng)壓力分布后，將其導(dǎo)入基于Lighthill聲類比理論的聲學(xué)軟件ACTRAN中進(jìn)行聲場(chǎng)仿真計(jì)算，實(shí)現(xiàn)了對(duì)水下潛器首部聲基陣區(qū)流噪聲的數(shù)值預(yù)報(bào)。研究了航速、共形陣的安裝位置和基陣單元安裝面形狀對(duì)聲基陣區(qū)流噪聲傳播的影響。結(jié)果表明：航速越大，流噪聲越大；增大聲吶安裝面與導(dǎo)流罩的距離以及使用較光滑的安裝面，可以減小聲基陣區(qū)流噪聲的大小。

水下潛器；首部；流噪聲；LES；Lighthill聲類比

0 引 言

水下潛器在航行時(shí)產(chǎn)生的噪聲不僅會(huì)暴露自己的位置，還會(huì)影響聲吶系統(tǒng)的正常工作。一般認(rèn)為：水下潛器的噪聲主要由艇體的機(jī)械振動(dòng)、螺旋槳運(yùn)轉(zhuǎn)噪聲以及水動(dòng)力噪聲 3 部分構(gòu)成。而水動(dòng)力噪聲包括由艇體周圍湍流邊界層內(nèi)的擾動(dòng)及壁面上的脈動(dòng)壓力共同引起的噪聲，通常稱為流噪聲。在低航速下，機(jī)械噪聲是水下潛器噪聲的主要部分。而在中速下（10～20 kn），流噪聲是噪聲的主要因素[1]。而在高航速下，如果導(dǎo)流罩的材料不合適，可能會(huì)因?yàn)榭栈a(chǎn)生嚴(yán)重的噪聲。根據(jù)聲吶方程：（SL–2TL+TS）–（NL–DI）=DT（主動(dòng)聲吶）減小自噪聲可增大聲吶的探測(cè)距離。而在中高速下，流噪聲會(huì)成為自噪聲的主要成分。有研究表明：聲吶平臺(tái)區(qū)的自噪聲減小5 dB，本艇的聲吶探測(cè)距離可增加30%，探測(cè)目標(biāo)的海區(qū)可增大 2 倍[2]。因此，研究流噪聲具有很大的實(shí)際意義。

對(duì)于水下潛器這種有著固體邊界的流動(dòng)模型，存在著 2 種不同類型的聲源：壁面上的脈動(dòng)壓力所構(gòu)成的偶極子源與流場(chǎng)中的旋渦運(yùn)動(dòng)、動(dòng)量能量交換所產(chǎn)生的四極子源。其中壁面上的脈動(dòng)壓力是產(chǎn)生流噪聲的主要部分，而四極子聲源由湍流中的速度脈動(dòng)量等引起，在馬赫數(shù)較低的流動(dòng)情況中往往可以忽略。

水下潛器流場(chǎng)的數(shù)值模擬方法始于20世紀(jì)70年代，但受限于當(dāng)時(shí)的技術(shù)，模擬效果不佳。隨著計(jì)算機(jī)技術(shù)的蓬勃發(fā)展，80年代后，數(shù)值模擬逐漸成為水下潛器設(shè)計(jì)研究中的重要技術(shù)。張楠等[3]采用k-ε湍流模型計(jì)算SUBOFF不同附體形式的阻力和尾流場(chǎng)。涂海文[4]采用RNGκ-ε湍流模型研究 4 種不同形式的SUBOFF模型的外流場(chǎng)，結(jié)果很好地預(yù)報(bào)了水下潛器的摩擦阻力系數(shù)，證明了CFD方法研究水下潛器流場(chǎng)的可行性。張楠等[5]采用 5 種湍流模型預(yù)報(bào) 2 種不同附體形式水下潛器的外流場(chǎng)，并與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，比較了這些湍流模型各自的預(yù)報(bào)精度。趙鵬偉等[6]利用RNGκ-ε湍流模型研究不同艇體半徑，艇首長度，艇尾長度對(duì)流場(chǎng)的影響，結(jié)果表明縮短艇首長度可降低尾部的湍流強(qiáng)度。俞強(qiáng)等[7]用不同的湍流模型研究了流線型和鈍形水下潛器首部的外流場(chǎng)特性。

而流噪聲的研究從1952年Lighthill[8, 9]提出聲類比理論開始。Lighthill將流場(chǎng)劃分為近場(chǎng)和遠(yuǎn)場(chǎng) 2 個(gè)部分，近場(chǎng)描述聲波的非線性產(chǎn)生過程，為聲源域；遠(yuǎn)場(chǎng)描述聲波的線型傳播過程，為傳播域。俞孟薩[10]利用集成統(tǒng)計(jì)能量法計(jì)算了聲吶自噪聲的水動(dòng)力噪聲分量，結(jié)果表明利用集成SEA方法預(yù)報(bào)流噪聲在中高頻具有較高精度。盧云濤[11]研究了 4 種目前應(yīng)用廣泛的湍流模型RNGκ-ε，Realizableκ-ε，κ-ω，SSTκ-ω并結(jié)合RANS方程和FW-H聲學(xué)模型預(yù)報(bào)水下潛器的阻力與流噪聲，但是RANS方法丟掉了脈動(dòng)信息量，計(jì)算結(jié)果與實(shí)驗(yàn)值有所差距。Zhang等[12]采用LES湍流模型和FW-H聲學(xué)類比法模擬了孔腔流動(dòng)噪聲，并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比。劉明靜等[13]研究了水下潛器首部不同位置處流噪聲的聲壓分布并得出了分布規(guī)律。

本文首先在CFD前處理軟件ICEM CFD中進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分。劃分過程中，對(duì)壁面處進(jìn)行網(wǎng)格加密處理形成致密的邊界層。然后將網(wǎng)格導(dǎo)入Fluent進(jìn)行求解計(jì)算，計(jì)算的模型選擇LES模型，在獲得流場(chǎng)及壁面上的脈動(dòng)信息后將其導(dǎo)入ACTRAN中進(jìn)行聲學(xué)仿真獲得計(jì)算結(jié)果。

1 數(shù)學(xué)模型

1.1 LES模型

在湍流流場(chǎng)中，對(duì)整個(gè)流場(chǎng)影響最大的是大尺度渦。小尺度的漩渦主要是引起湍流動(dòng)量的擴(kuò)散。基于這種理論，LES（大渦模擬）提供了求解流場(chǎng)的一種新的有效的方法。LES的基本思想是：對(duì)大尺度渦直接求解非定常N-S方程得到精確解；對(duì)小尺度渦，則建立湍流模型模擬小尺度渦對(duì)大尺度渦的影響。LES方法是一種介于RANS（雷諾平均方法）和DNS（直接數(shù)值模擬）之間的方法，對(duì)計(jì)算機(jī)的要求較DNS低，同時(shí)計(jì)算結(jié)果又比RANS準(zhǔn)確。

LES的控制方程包括 2 個(gè)，一個(gè)是連續(xù)性方程，另一個(gè)是通過濾波函數(shù)處理的非定常N-S方程。

1.2 Lighthill聲類比

Lighthill聲類比理論是從納維-斯托克斯（Navier-Stokes）方程導(dǎo)出的。Lighthill的模型為：在無限大的均勻、靜態(tài)聲介質(zhì)中包含一個(gè)有限的湍流運(yùn)動(dòng)區(qū)域V。其中將有限的湍流流動(dòng)的區(qū)域設(shè)定為近場(chǎng)，即聲源區(qū)；將流動(dòng)區(qū)域以外的區(qū)域設(shè)定為傳播區(qū)，即輻射區(qū)。由此得到Lighthill聲類比方程：

1.3 評(píng)價(jià)指標(biāo)

聲學(xué)分析中聲壓級(jí)和聲壓總級(jí)的計(jì)算公式式（4）和式（5）。

其中：P（f）為f頻率下的聲壓大?。籔0為基準(zhǔn)聲壓，在水中，P0為 1×10–6Pa；?f為計(jì)算中頻率的分辨率?？紤]到現(xiàn)代水下艦艇安裝的低頻主被動(dòng)聲吶工作頻率的范圍，本文所計(jì)算的聲壓總級(jí)均從20 Hz開始統(tǒng)計(jì)。

2 聲基陣聲學(xué)模型

首部線型其原長度約為10 m，直徑約為8 m，本文對(duì)其進(jìn)行了1:5的縮尺。首部的聲吶導(dǎo)流罩在聲學(xué)計(jì)算中設(shè)置為透聲窗。其材料性質(zhì)為：厚度0.015 m，楊氏模量18.5 GPa，密度1 700 kg/m3。首部的其他部分在聲學(xué)計(jì)算中均設(shè)置為全反射壁面。

首部的內(nèi)部有聲基陣安裝面（見圖1），此為聲吶基陣單元的安裝位置。本文研究的是一個(gè)共形陣。所謂共形陣就聲吶單元附著與載體表面并與載體貼合，這樣既可以不破壞載體的外形結(jié)構(gòu)和流體性能又能相對(duì)擴(kuò)大安裝空間。在本文的計(jì)算中，將其簡化為從上至下橫截面逐漸縮小的柱面。在導(dǎo)流罩和聲基陣安裝面之間的介質(zhì)為水。為了研究不同位置處共形陣的聲壓大小分布，本文建立了 3 個(gè)共形陣面，每個(gè)陣面間距為0.1 m。按距離導(dǎo)流罩的遠(yuǎn)近，分別命名為“外陣面”，“中陣面”和“內(nèi)陣面”。每個(gè)陣面上均布了500個(gè)左右的水聽器以檢測(cè)不同位置處流噪聲大小。圖1中顯示了其中一個(gè)陣面上水聽器的分布情況。其中，首部聲基陣的隔聲障板圖1中黃色的部分為用來隔斷水下潛器尾部螺旋槳的噪聲，防止干擾聲基陣的工作。

考慮到首部較短，在計(jì)算中可能因回流產(chǎn)生負(fù)壓，本文將首部線型往后光滑延伸，得到流體計(jì)算的模型（見圖2）。整個(gè)計(jì)算域?yàn)樗颉Ｕ麄€(gè)流域約為首部艇長的4倍，流域的入口在艇上游1倍艇長處，出口在尾部2倍艇長處，流域的直徑為艇身直徑的3倍。

網(wǎng)格的劃分采用CFD前處理軟件ICEM進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分。計(jì)算域的網(wǎng)格總數(shù)為500萬，并且利用了ICEM中的prism技術(shù)在艇體表面生成了精細(xì)的邊界層網(wǎng)格，其中第一層厚度為0.2 mm，厚度遞增因子1.2，一共生成10層。

3 計(jì)算方法

本文的計(jì)算步驟可總結(jié)為以下幾步：

1）流場(chǎng)計(jì)算?；诓豢蓧嚎s流體的設(shè)定，將計(jì)算網(wǎng)格導(dǎo)入Fluent并采用LES模型湍流對(duì)水下潛器周圍的流場(chǎng)進(jìn)行定常計(jì)算。其中邊界條件的設(shè)置為：入口為速度入口；出口為壓力出口；水下潛器表面為不可滑移壁面；外域面選擇速度入口。定常計(jì)算穩(wěn)定后，以此定常解作為初始值進(jìn)行非定常LES計(jì)算，設(shè)定計(jì)算時(shí)間步長為2.5×10–4s。并通過對(duì)阻力系數(shù)的監(jiān)控來判斷流場(chǎng)是否計(jì)算穩(wěn)定。待流場(chǎng)穩(wěn)定后，將計(jì)算結(jié)果保存2 000步。根據(jù)采樣定理：計(jì)算結(jié)果的分辨率為1 Hz，計(jì)算最高頻率為2 000 Hz。這里的流場(chǎng)計(jì)算主要是為聲學(xué)計(jì)算提供聲源項(xiàng)（對(duì)于不可壓縮流體，主要是給出速度場(chǎng)的信息）。

2）聲場(chǎng)計(jì)算。通過ACTRAN中的ICFD模塊，將流體計(jì)算結(jié)果直接導(dǎo)入到ACTRAN中求解流噪聲輻射問題。具體求解過程為：首先通過ACTRAN中的ICFD模塊將流體計(jì)算結(jié)果中的速度、壓強(qiáng)等信息通過積分算法插值到聲學(xué)網(wǎng)格上，得到時(shí)間歷程載荷，然后通過快速傅里葉變換（FFT）將時(shí)間歷程載荷轉(zhuǎn)化為頻譜載荷，最后利用ACTRAN進(jìn)行求解。

聲學(xué)模型如圖3所示。為有限元域聲源域，需要?jiǎng)澐志W(wǎng)格；為無限元域傳播域，不需要?jiǎng)澐志W(wǎng)格。將有限元網(wǎng)格區(qū)域的外表面設(shè)置為無限元基面，水下潛器表面采用設(shè)置為隔聲面，另外噪聲透過導(dǎo)流罩向首部內(nèi)部傳播。

整個(gè)計(jì)算流程如圖4所示。

4 計(jì)算結(jié)果與分析

本文共模擬了6種不同工況下首部內(nèi)部的流噪聲特性：計(jì)算速度為6 kn，9 kn，12 kn，15 kn，18 kn，21 kn和24 kn，計(jì)算過程中設(shè)定聲基陣安裝面為全反射壁面。

4.1 流場(chǎng)結(jié)果分析

本文采用了大渦模擬方法對(duì)水下潛器首部外流場(chǎng)進(jìn)行仿真，下面給出6 kn航速下，首部外流場(chǎng)的特性。

由圖5可以定性而且直觀地看到流場(chǎng)模擬的情況。頭部是壓力較大的地方，中間部分逐漸趨于穩(wěn)定。當(dāng)流體流過艇體頭部時(shí)，在頭部附近受到阻滯，流速迅速降低，壓力隨即增大，形成所謂的全受阻壓力，艇體首部的局部壓力大于來流靜壓，而形成正壓駐點(diǎn)區(qū)。然后，流體沿著艇體表面，流動(dòng)速度逐漸增大，壓力迅速降低，并在平直段出現(xiàn)穩(wěn)定的壓力區(qū)。

其中艇體的首部是壓力梯度比較大的地方，這片區(qū)域也是流噪聲產(chǎn)生的主要區(qū)域。

4.2 聲場(chǎng)結(jié)果分析

以下數(shù)據(jù)都經(jīng)過無量綱化處理。

圖6中典型點(diǎn)位于外陣面上，其位置由圖1標(biāo)出。由圖可以看出，艇體的自流噪聲能量主要集中在低頻率范圍內(nèi)（20～500 Hz），在低頻率范圍內(nèi)聲壓級(jí)隨頻率的增大而迅速的衰減。高頻部分雖然有峰值的存在，但較低頻的聲壓級(jí)值較小，因此對(duì)聲壓總級(jí)影響較小。根據(jù)渦聲理論，流場(chǎng)中的聲由其中的渦產(chǎn)生，大渦頻率低，產(chǎn)生低頻噪聲，攜帶大量能量；小渦頻率高，產(chǎn)生高頻噪聲，攜帶較少能量。因此本文得到的頻譜曲線趨勢(shì)符合渦聲理論。另外該頻譜與s u boff流噪聲實(shí)驗(yàn)[14]的頻譜趨勢(shì)基本一致。

4.2.1 航速

對(duì)于水下潛器來說，航速是影響其性能的最重要的指標(biāo)，它對(duì)航行器的水動(dòng)力性能，聲學(xué)性能等都具有很大的影響。下面是計(jì)算得到的航速對(duì)內(nèi)部聲學(xué)性能的影響?？紤]到首部內(nèi)部不同位置處聲壓總級(jí)差別較大，本文采用統(tǒng)計(jì)學(xué)的方法，對(duì) 3 個(gè)陣面上的所有水聽器以及 3 個(gè)典型監(jiān)測(cè)點(diǎn)采集到的聲壓總級(jí)值進(jìn)行算術(shù)平均得到的值進(jìn)行分析。

由圖7可以看出，航速對(duì)聲壓總級(jí)值的影響顯著。無論對(duì)于典型檢測(cè)點(diǎn)還是 3 個(gè)陣面各監(jiān)測(cè)點(diǎn)的平均值來說，隨著航速的逐漸增大，聲壓總級(jí)值逐漸增大，但是其增大的幅值并不固定：在低速下，增長較快，高速下，增長緩慢，兩者近似呈現(xiàn)對(duì)數(shù)關(guān)系。

為了確定聲壓總級(jí)值和航速的定量關(guān)系，采用最小二乘法建立聲壓總級(jí)值的平均值與航速的回歸模型。假設(shè)

式中Lp為聲壓總級(jí)值，v為航速。通過回歸計(jì)算可得：

相關(guān)系數(shù)R2為0.997 2，因?yàn)镽2≈1，所以式（7）兩者關(guān)系擬合較好。由此也可以看到航速對(duì)首部自流噪聲聲壓總級(jí)的影響較大。

4.2.2 陣面位置

傳統(tǒng)的水下潛器受制于有限的內(nèi)部空間，首部的聲吶基陣一般布置成平面陣。盡管平面陣具有結(jié)構(gòu)簡單、信號(hào)處理方便、成本低廉等一系列優(yōu)點(diǎn)，但是其掃描角度有很大限制，且單個(gè)換能器尺寸較小。在無法增大首部內(nèi)部空間時(shí)，就限制了首部聲吶的探測(cè)性能。且平面陣聲吶在探測(cè)目標(biāo)時(shí)需要調(diào)整艦艇位置使平面陣正對(duì)著目標(biāo)，這限制了艦艇的機(jī)動(dòng)性。為了解決這個(gè)問題，共形陣就應(yīng)運(yùn)而生。所謂共形陣就是和安裝面表面共形，這樣不僅可以極大提高首部內(nèi)部空間的利用效率，也可以增大單個(gè)換能器的尺寸，增大聲吶探測(cè)距離，且不會(huì)影響艦艇的機(jī)動(dòng)性能。

盡管共形陣有諸多優(yōu)點(diǎn)，但是它也有自身的不足。如共形陣曲面過于復(fù)雜，對(duì)信號(hào)的處理技術(shù)要求很高，這就需要在艦艇上安裝高性能的計(jì)算機(jī)等硬件設(shè)備，且對(duì)艇員的素質(zhì)要求較高。此外，共形陣在首部的安裝位置也是一個(gè)很棘手的問題：太靠近導(dǎo)流罩，聲吶基陣受自身流噪聲影響較大；太遠(yuǎn)離導(dǎo)流罩，整個(gè)陣面的尺寸就受到限制。為了解決這個(gè)問題，本文在首部內(nèi)部建立了 3 個(gè)共形陣面，研究不同位置處自流噪聲大小與共形陣面的面積之間的關(guān)系

表1顯示了 6 個(gè)航速下，3 個(gè)陣面上水聽器接收到的流噪聲聲壓總級(jí)平均值，以及它們之間的差值。其中LP（I），LP（M）和LP（O）分別代表內(nèi)陣面，中陣面和外陣面上的聲壓總級(jí)平均值。由表可見在等間距的間隔內(nèi)，越往首部內(nèi)部，流噪聲衰減的幅度越小。但是在不同航速下，同 2 個(gè)陣面上的差值近乎相同。流體動(dòng)力聲源可分為單極子聲源、偶極子聲源和四極子聲源。單極子聲源在流體中不穩(wěn)定，一般不考慮。而四極子聲源與馬赫數(shù)四次方成正比，本文中最大的馬赫數(shù)為0.008，故四極子聲源也很小，亦不考慮。偶極子聲源聲壓的衰減成倒數(shù)關(guān)系，即距離偶極子聲源的等效中心R處聲壓為中心處聲壓的1/R倍。所以越遠(yuǎn)離導(dǎo)流罩，流噪聲減小的幅度越小。

實(shí)際工作中，為了求得最合適的陣面位置，可在同一坐標(biāo)系下做出每個(gè)陣面的面積和其上監(jiān)測(cè)點(diǎn)的聲壓總級(jí)平均值隨距離的變化關(guān)系。據(jù)此可以找出一個(gè)陣面，使得該聲吶的整體探測(cè)距離更遠(yuǎn)。

表 1 不同航速下各陣面檢測(cè)點(diǎn)聲壓總級(jí)平均值Tab. 1 The average sound pressure level of different array and speed

4.2.3 安裝面形狀位置

聲吶的聲基陣是以若干個(gè)換能器作為基陣單元，按照一定的幾何形狀和分布規(guī)律排列的陣列。聲吶基陣受到的噪聲越小噪聲分布越均勻，則越有利于提高聲吶的信噪比，從而增強(qiáng)聲吶的探測(cè)能力。圖8顯示了航速為12 kn，假設(shè)安裝面為全反射壁面情況時(shí)，內(nèi)陣面上的水聽器聲壓總級(jí)分布規(guī)律。

由圖8可以看到，同一個(gè)聲基陣上各水聽器接收到的流噪聲聲壓大小并不相同。聲壓總級(jí)較大的區(qū)域?yàn)轭^部以及舷側(cè)部分；而較小的區(qū)域?yàn)橄蟼?cè)與頭部中間的部位。整個(gè)陣面中流噪聲聲壓總級(jí)最大的點(diǎn)為圖中的4號(hào)點(diǎn)，此點(diǎn)對(duì)應(yīng)的艇體表面其壓力梯度較大（見圖5）。由圖中也可以看到 3 條帶狀的區(qū)域，其聲壓較周圍小。這是因?yàn)楹喕Ｐ蜁r(shí)，把安裝面簡化為一系列平面，平面交界的位置處為一條線（見圖1）。這樣當(dāng)流噪聲的聲波信號(hào)傳播過來的時(shí)候，該交界位置處對(duì)應(yīng)的區(qū)域只有由外部傳播過來的聲壓，沒有反射回去的聲壓疊加，故而聲壓總級(jí)較??；而這些區(qū)域以外的部分，由于疊加了安裝面反射回來的其他陣面的聲壓，因此很大。

由于聲吶基陣上各單元對(duì)整體信號(hào)探測(cè)的重要程度并不相同[1]，因此在設(shè)計(jì)的時(shí)候可以通過調(diào)整聲基陣安裝面的形狀來改變流噪聲在首部空間的傳播，以減少對(duì)整體信號(hào)較重要的單元接收到的流噪聲大小。

5 結(jié) 語

本文基于LES湍流模型和Lighthill聲類比理論對(duì)一種水下潛器的外流場(chǎng)和聲基陣區(qū)聲場(chǎng)進(jìn)行仿真計(jì)算。分析了不同情況下的流噪聲在首部內(nèi)部空間的傳播方式，對(duì)減小聲基陣區(qū)的流噪聲大小以及增大聲吶的整體工作性能提出了一些建議，得到以下結(jié)論：

1）隨著航速的增加，首部聲基陣區(qū)接收到的流噪聲的聲壓總級(jí)逐漸增大，其增大近乎成對(duì)數(shù)關(guān)系。

2）共形陣陣面與導(dǎo)流罩的距離決定了聲吶受到的流噪聲大小以及聲吶的有效探測(cè)距離。通過分析首部內(nèi)部空間流噪聲衰減以及陣面有效面積的變化，可選擇最合適的陣面位置。

3）可通過改變聲基陣安裝面的形狀來改變某些特定位置處基陣單元受到的流噪聲大小，以使聲吶的整體性能最優(yōu)。

[ 1 ]WAITE A D. 實(shí)用聲吶工程. 第三版[M]. 王德石譯. 北京: 電子工業(yè)出版社, 2004.

[ 2 ]繆旭弘, 王振全. 艦艇水下噪聲控制技術(shù)現(xiàn)狀及發(fā)展對(duì)策[C]//中國造船工程學(xué)會(huì)船舶力學(xué)學(xué)術(shù)委員會(huì)水下噪聲學(xué)組.第十屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集, 2005: 5.

[ 3 ]張楠, 沈泓萃, 姚惠之. 潛艇阻力與流場(chǎng)的數(shù)值模型與驗(yàn)證及艇型的數(shù)值優(yōu)化研究[J]. 船舶力學(xué), 2005, 9(1): 1–13.

[ 4 ]涂海文, 孫江龍. 基于CFD的潛艇阻力及流場(chǎng)數(shù)值計(jì)算[J].艦船科學(xué)技術(shù), 2012, 34(3): 1–13. TU Hai-wen, SUN Jiang-long, Numerical analysis of resistance and flow field of submarine based on CFD[J]. Ship Science and Technology, 2012, 34(3): 1–13.

[ 5 ]張楠, 沈泓萃, 姚惠之. 用雷諾應(yīng)力模型預(yù)報(bào)不同雷諾數(shù)下的潛艇繞流[J]. 船舶力學(xué), 2009, 13(5): 688–696.

[ 6 ]趙鵬偉, 盧曉平, 孫玉明. 基于潛艇模型尾流湍流強(qiáng)度和耗散率的CFD模擬[J]. 中國艦船研究, 2014, 9(3): 43–56.

[ 7 ]俞 強(qiáng), 魏子凡, 楊松林, 等. 基于CFD不同AUV艇體阻力性能研究[J]. 船海工程, 2014, 43(2): 177–181.

[ 8 ]LIGHTHILL M. J. On sound generated aerodynamically. General theory[J]. Proc Royal Soc London, Ser A, 1952, 211: 564–587.

[ 9 ]LIGHTHILL M. J. On sound generated aerodynamically. Turbulence as a source of sound[J]. Proc Royal Soc London, Ser A, 1954, 222: 1–25.

[10]俞孟薩, 朱正道. 集成統(tǒng)計(jì)能量法計(jì)算聲吶自噪聲水動(dòng)力噪聲分量[J]. 船舶力學(xué), 2007, 02: 273–283.

[11]盧云濤. 全附體潛艇的流場(chǎng)和流噪聲的數(shù)值研究[D]. 上海:上海交通大學(xué), 2008.

[12]ZHANG N, SHEN H C, YAO H Z. Numerical simulation of cavity flow induced noise by LES and FW-Hacoustic analogy [J]. Journal of Hydrodynamics(Ser. B), 2010, 22(Supp 1): 242–247.

[13]劉明靜, 馬運(yùn)義. 潛艇艏部聲吶流噪聲計(jì)算方法研究[J]. 船海工程, 2009, 38(5): 46–49.

[14]張平, 張宇文, 劉維, 等. 典型工況下回轉(zhuǎn)體流噪聲特性實(shí)驗(yàn)對(duì)比探究[J]. 計(jì)算機(jī)仿真, 2012, 08: 26–30. ZHANG Ping, ZHANG Yu-wen, LIU Wei, et al, Study on flow noise characteristics of symmetrical model under typical working conditions[J]. Computer Simulations, 2012, 08: 26–30.

Research method about flow-induced noise in the bow of underwater vehicles

FAN Lin-xu1, ZHANG Lei1, CHEN De-xi2, LIU Guo-qing1, WANG Chi-zhong1
(1. School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; 2. China Ship Development and Design Center, Wuhan 430064, China)

Underwater vehicle use sonar system to detect the enemy ship as well as its own position during its navigation. The main sonar system which can be divided into active sonar and passive sonar is usually placed in the bow. However, noise is an important interference factor that both systems must overcome. In order to research flow-induced noise in the bow, the LES model was applied to simulate flow field out of the underwater vehicle to get fluctuating pressure. Then the fluctuating pressure was import into ACTRAN- a software based on Lighthill's acoustic analogy theory- to simulate the sound field in the bow. Speed, position of the array and shape of the array carrier were studied to reveal their influences on the propagation of flow-induced noise. The results show thatthe higher speed, the bigger of flow-induced noise. Furthermore, flow-induced noise can be diminished by increasing the distance between sonar array and dome and using relative smooth installion surface.

underwater vehicle；bow；flow-induced noise；LES；Lighthill acoustic analogy theory

U661.39

A

1672 – 7649(2017)07 – 0048 – 06

10.3404/j.issn.1672 – 7649.2017.07.010

2016 – 10 – 25；

2016 – 12 – 07

樊林旭(1993-)，男，碩士研究生，研究方向?yàn)楣芟嫡駝?dòng)與流噪聲。