張 磊，樊林旭,2，沈理姣，張俊杰，孫江龍,2，王赤忠

(1.華中科技大學船舶與海洋工程學院，湖北武漢 430074；2.船舶與海洋工程水動力學湖北省重點實驗室，湖北武漢 430074；3.中國艦船研究設計中心，湖北武漢 430064；4.浙江大學海洋學院，浙江舟山 316021)

水下航行器頭部流噪聲性能數(shù)值模擬研究

張 磊1，樊林旭1,2，沈理姣3，張俊杰3，孫江龍1,2，王赤忠4

(1.華中科技大學船舶與海洋工程學院，湖北武漢 430074；2.船舶與海洋工程水動力學湖北省重點實驗室，湖北武漢 430074；3.中國艦船研究設計中心，湖北武漢 430064；4.浙江大學海洋學院，浙江舟山 316021)

流噪聲的降低對水下航行器的隱蔽性提升、主動聲吶系統(tǒng)探測能力的提高具有重要意義。首部作為水下航行器主要構成部分，其線型設計的優(yōu)劣直接關系到流噪聲性能和聲吶系統(tǒng)工作環(huán)境的好壞?；诖鬁u模擬（LES）的流場計算結果，采用ACTRAN這一聲學軟件對2種不同頭部線型、不同速度、不同攻角、不同監(jiān)測點處的流噪聲性能進行研究。分析研究發(fā)現(xiàn)，水下航行器流噪聲能量主要集中在低頻段，隨著速度的增大而增大。在速度、攻角相同的情況下，頭部端面半徑小、頭部長度較長的模型具有較好的流噪聲性能。在同一模型下，頭部駐點處的流噪聲比頭部端面一定半徑處任何點的流噪聲要小。將數(shù)值模擬與水下航行器流噪聲特性水洞試驗研究進行對比，兩者得到的結論一致。

水下航行器；頭部線型；大渦模擬；流噪聲

0 引 言

噪聲級是衡量水下航行器（潛艇）的一項重要性能參數(shù)。潛艇的噪聲分為3類：機械振動噪聲、螺旋槳噪聲和流噪聲。水下航行器噪聲產(chǎn)生機理十分復雜，不同部位受到三大噪聲源的影響也不同，如流噪聲對首部聲吶處自噪聲的影響較大。對于潛艇裝置的聲吶，控制流噪聲造成的自噪聲意義重大。據(jù)調查，聲吶平臺自噪聲降低5 dB，本艇聲吶探測距離增加60%，探測目標的海區(qū)面積為原來的3倍左右[1]。雖然流噪聲的數(shù)值比機械振動噪聲要低很多，但是隨著潛艇速度的不斷提高，其在潛艇總噪聲中的比重也迅速增加，從潛艇航行隱蔽性而言，降低流噪聲也十分重要。

目前，工程上對于流噪聲的研究主要分為試驗測量和數(shù)值仿真。試驗研究具有花費大、對試驗儀器測量精度要求高、等比例模型試驗難以完成等缺點，極大限制了工程上廣泛借助試驗方法來對水下航行器的流噪聲進行研究。隨著CFD技術以及計算機資源的極大提升，數(shù)值模擬成為可能。數(shù)值模擬雖不能完全代替試驗測量，但是作為試驗測量的補充，在對比分析、了解噪聲產(chǎn)生機理、優(yōu)化選型、為試驗減少工況等方面起到了重要作用。數(shù)值模擬在一些方面也具有試驗不可比擬的優(yōu)點，例如在數(shù)值仿真過程中可以在任意場點上布置虛擬接收器，有效避免了試驗中水聽器安裝困難的問題，也避免了試驗中聲接收器與流場之間的相互干擾問題。此外，試驗中環(huán)境噪聲或水洞洞壁的聲透射等問題在數(shù)值模擬中也不存在。

本文首先通過與已有文獻計算結果進行對比，說明LES和ACTRAN聲學軟件聯(lián)合計算流噪聲這種方法滿足一般工程的精度要求。隨后對2種不同頭部線型的水下航行器在不同航速、不同攻角、不同測試點的流噪聲進行研究，旨在分析頭部線型、來流速度和測試點位置對流噪聲的影響，然后與黃橋高[2]水下航行器流噪聲特性水洞試驗結果進行對比，本文數(shù)值模擬得到的結論與水洞試驗相同，進一步說明本文所用的流噪聲計算方法可以用于流噪聲的對比分析，最后根據(jù)數(shù)分析給出關于潛艇首部的線型設計和聲吶安裝位置的一些建議。

1 計算模型

1.1 頭部線型說明

本文選擇雙參數(shù)橢圓曲線作為頭部線型，雙參數(shù)橢圓線型的表達式為：

根據(jù)式（1），水下航行器頭部線型方程為：

式中：x為航行器物理橫坐標；y為航行器物理縱坐標；D為航行器的最大直徑；Df為航行器前端面直徑；Lh為航行器頭部長度；m，n為頭部參數(shù)[2,3]。

本文選擇如圖1所示，2種頭部線型。

1.2 計算域與計算網(wǎng)格

本文雖然只研究頭部線型的噪聲性能，但是建模時不能做截斷處理只建立頭部的線型，為了避免截斷面后會出現(xiàn)回流的影響，這里將模型平滑延長，模型總長為700mm，如圖2所示。根據(jù)水下航行器的實際尺寸，綜合考慮計算量、以及所需的流場詳細信息的準確性。流場域的大小滿足：入口在航器首部1倍長度處，出口在航器尾部2倍長度處[4]，流場域選擇圓柱域，半徑為航行器半徑的10倍，見圖3。為了節(jié)省計算資源，網(wǎng)格劃分在航器周圍采用外O-BLOCK的技術形成加密區(qū)，以達到盡可能獲得多的航行器頭部流場流動信息的目的。整個模型采用結構化的六面體網(wǎng)格，在網(wǎng)格生成過程中，分布應該稀疏合理?？拷吔鐚?，網(wǎng)格太稀不利于捕捉到流動細節(jié)；物面附近參數(shù)變化梯度大，網(wǎng)格太稀則與艇體實際形狀誤差大；離艇體較遠處，參數(shù)變化小，網(wǎng)格太密則浪費計算資源。圖4為頭部的局部網(wǎng)格示意圖，可見航行器壁面以及頭部的網(wǎng)格尺寸很小，以便捕捉詳細的邊界層流場信息。計算域網(wǎng)格總數(shù)為453 040，節(jié)點總數(shù)為443 526，為了準確捕捉壁面信息y+取50。

2 計算原理

2.1 湍流模型

對于流噪聲的預報，其不可壓縮粘性湍流流動優(yōu)先采用大渦模擬（LES）。LES的基本思想可概括為：用瞬時的N-S方程直接模擬湍流中的大尺度渦，對于小尺度渦的處理采用近似的模型來考慮。相對于RANS方法對控制方程中的時間量做平均得到脈動量的時均值的方法，LES可以為噪聲計算提供更詳細的更全面的流場信息，從而使得LES更適合于聲學計算。要實現(xiàn)LES必須完成以下2個重要環(huán)節(jié)的工作：一是建立一種數(shù)學濾波函數(shù)；二是加入被濾掉小渦對大渦運動的影響，LES中通過加入附加應力項實現(xiàn)，該應力項稱為亞格子尺度應力。

空間中位于x0處的濾波后的變量定義為

式中：G為濾波函數(shù)。濾波后參量的結構和大小由濾波函數(shù)決定，而濾波函數(shù)依賴于以及濾波寬度，其中為第i個空間坐標中的濾波寬度。

亞格子尺度模型采用Smagorinsky模型，該模型是基于均衡假設，即小尺度參量將從大尺度參量吸收的全部能量完全地、瞬時地耗散干凈。根據(jù)Smagorinsky的基本SGS模型，SGS應力具有下面的形式：

式中：μt是亞格子尺度的湍動粘度；亞格子尺度應力τkk的各向同性部分既可以通過建模添加到濾波后的壓力函數(shù)中，也可以忽略不計[5]。

亞格子尺度的湍動粘度為

2.2 流噪聲方程

從聲學基本概念來講，聲學其實是流場的不平衡擾動。通過緊致聲源的假設，1952年Lighthill[6,7]提出聲學類比理論，使噪聲的求解成為可能。該理論將聲學計算和流場計算結合在一起，通過求解流場信息進而求解流噪聲。Lighthill方程是由N-S方程導出，得到密度波動方程：

式中：c0為流體介質中等熵條件下的聲速；，其中ρ和ρ0分別是擾動與未擾動時流體的密度；Tij為Lighthill應力張量，可表示為

本文所用的ACTRAN對流噪聲的預報是基于Ligthhill聲類比的方法，其對波動方程進行了改進如下：

Ligthhill應力張量則改寫為式：

Ligthhill聲類比方程的弱變分方程為：

其中：方程右端第1項表征體聲源對聲場的貢獻，第2項表示面聲源對聲場的貢獻。體聲源包括3個部分：對流輸運的散度、熵源項、粘性應力的散度。ACTRAN在計算時，忽略了熵源項和粘性應力。在ACTRAN計算中，壁面邊界默認為剛性壁面，此時面源項等于0，對應與式（11）中右端第2項為0。

3 流噪聲計算方法驗證

本文流噪聲的計算首先采用LES方法獲得三維流場的瞬態(tài)信息，然后基于Lighthill的聲類比理論對流噪聲進行模擬。本文對研究數(shù)據(jù)較豐富的“大青花魚”簡化潛艇模型的三維流場及流噪聲進行了研究，得到了x方向聲壓總級變化特性對比。“大青花魚”簡化模型尺度：艇長L=3.2m，型寬B=0.4m。流噪聲對比計算結果見表1，3種結果流噪聲的輻射特性基本一致，在沿長度方向上呈逐漸遞減的規(guī)律。對于離艇體較近的區(qū)域聲壓2種不同方法（BEM和Actran）計算結果存在一定差別，而當艇體距離不斷增大時，3個結果具有很高的吻合度。本文計算結果與孟堃宇[8]、劉繼明[9]研究數(shù)據(jù)進行對比可知，除了在艇體表面較近處，本文采用的計算結果滿足工程要求的精度（誤差少于3 dB）；而對于流噪聲的計算，ACTRAN采用的是有限元和無限元來模擬聲源，而Sysnoise是采用邊界元，從所用方法來說，ACTRAN得到的近場結果更加準確。綜上討論，本文選用對近場和遠場流噪聲模擬具有更高精度的ACTRAN。

表1 仿真結果對比Tab.1 Comparison of calculation results

4 計算流程及數(shù)值結果

4.1 流場計算

本文劃分好網(wǎng)格后，將其導入Fluent中進行流場計算。流場數(shù)值模擬的計算步驟為：首先給定來流速度，應用RNG k-ε進行定常計算。為提高計算精度，計算時對動量、湍動能的處理采用了二階迎風差分格式，壓力速度耦合采用了SIMPLE算法。待定常計算充分收斂之后，通過阻力性能來驗證流場模擬的精度；在比較滿意的條件下以此定常解作為初值進行LES非定常計算。LES采用Smagorinsky-lilly粘性模型，壓力速度耦合選擇SIMPLEC算法，壓力空間離散采用PRESTO!方法，動量的處理采用有界中心差分（Bounded Central Differencing）。下面僅對第 1 種頭部線型在8m/s時的流場計算進行分析。

根據(jù)桑海公式，對于Re在1.0E+06到1.0E+09之間的估算公式：

4.2 流噪聲計算

表2 水下航行器阻力計算結果Tab.2 The drag coefficients of underwater vehicle

在非定常穩(wěn)定后，保存5 000步流場結果，時間步長為1.0E–04 s，取該值的原因為：LES需要滿足CFL在1附近的條件，其值的大小由最小求解渦的時間尺度來決定，若時間步長取太短固然能使CFL數(shù)小于1但會造成迭代步數(shù)增多，消耗大量時間的同時也會積累大量的舍入誤差，綜合考慮計算資源，所以時間步長取為1.0E–04 s。根據(jù)采樣頻率及采樣定理，轉換得到的頻域流場信息分辨率為2 Hz，帶寬是5 kHz。由于聲學計算軟件ACTRAN基于有限元、無限元的方法，則離散分為2個部分，一個是遠場的無限元區(qū)域，一個是近場的有限元區(qū)域；無限元不需要劃分網(wǎng)格，在ACTRAN里直接設置其階數(shù)即可（一般對于近場網(wǎng)格比較密的情況下，選擇默認設置為5）。這里ICEM的網(wǎng)格劃分只針對有限元區(qū)域，聲場中網(wǎng)格大小的劃分和分析頻率密切相關。分析頻率越高，波長越短，按最高頻率計算波長，要求每個波長范圍內至少要有6～8個線型單元。其計算模型設置如圖6所示：在聲源域，將有限元網(wǎng)格灰色區(qū)域的外表面設置為無限元基面，潛艇表面采用默認設置（隔聲面）。布置2個接收點，A點位于頭部駐點處，B點位于頭部端面15mm半徑處。

利用Fluent進行三維流場分析，通過ACTRAN中的ICFD模塊，將流體計算結果直接導入到A C TRAN中求解流噪聲輻射問題。ICFD為ACTRAN的中將CFD結果轉化為聲源的模塊，首先將流體計算中速度、壓強等信息通過積分算法插值到聲學網(wǎng)格上，得到聲學模型的時間歷程載荷，再利用其內部的傅里葉變換算法將時間歷程載荷轉化為頻譜載荷，最后利用ACTRAN進行求解。

圖7給出了模型1和模型2當來流無攻角時，在2個不同速度下，2個不同接收點處的流噪聲聲壓級特性曲線。

表3是根據(jù)聲壓級特征曲線計算出的聲壓總級值（單位為dB）。計算公式如下：

1）速度對流噪聲的影響分析

圖7給出了模型1和模型2在來流無攻角時不同速度下2個測試點的流噪聲聲壓級頻譜特性曲線。從圖中可以看出，流噪聲聲壓級隨著速度的增大而增大。聲學能量主要集中低頻下，隨著頻率的增大，流噪聲的聲壓級具有下降的趨勢。分析原因普遍認為，轉捩區(qū)的聲輻射是構成水下航行器流噪聲的主要組成部分，隨著來流速度的增加，轉捩點位置前移是導致輻射到相應位置在噪聲增大的原因。根據(jù)表3，模型1和模型2不同攻角情況下，在10m/s時的聲壓總級比8m/s時的聲壓總級有不同程度的升高，最高不到8 dB。

2）頭部線型對流噪聲的影響分析

根據(jù)表3，對比模型1和模型2的聲壓總級可以發(fā)現(xiàn)，相同速度下模型2的噪聲總聲壓級要比模型1大，分析原因認為，首部流噪聲主要源于流場轉捩區(qū)的聲輻射，模型1的頭部端面直徑小、且頭部曲線段長度長，這種線型使得最小壓力點和邊界層轉捩點（湍流起始點）的位置后移，這樣輻射到模型頭部相應監(jiān)測點處的噪聲要小。

3）測試點位置對流噪聲的影響分析

表3給出了模型1在0°，2°，4°攻角和模型2在0°攻角時不同測試點在不同速度下的流噪聲總聲壓級。頭部端面駐點處的聲壓總級比15mm半徑處的聲壓總級大約小11～17 dB，這里與黃橋高試驗在差值大小上有一定的差別，這里分析原因為：黃橋高試驗統(tǒng)計頻率范圍為15～35 kHz，轉捩區(qū)的輻射聲波入射到模型頭部端面時，聲波頻率高的時候，聲波的波長比較短，則需要考慮聲波傳播的繞射損失。而本文由于考慮到數(shù)值模擬分析最高頻率與聲波波長有關，在生成聲學計算網(wǎng)格需要劃分很密，超出計算機的計算能力，本文分析的最高頻率選為5 kHz，這里轉捩區(qū)的聲傳播過程中繞射損失要小很多。相比頭部端面其他位置，駐點處的流噪聲聲壓總級最小。

表3 水下航行器聲壓總級Tab.3 Overall sound pressure level of underwater vehicle

5 結 語

本文在不同速度、不同攻角下，分別對2種不同頭部線型的水下航行器流噪聲進行了數(shù)值仿真，對計算結果分析可得出以下結論：

1）在不同速度下，聲學能量主要集中低頻下。隨著頻率的增大，流噪聲的聲壓級具有下降的趨勢。對于同一線型的頭部，流噪聲隨著速度的增大，呈現(xiàn)出升高的趨勢。不同線型的頭部在同一速度下的流噪聲大小不同，就本文分析的2種線型而言，線型1具有較好的噪聲性能，同時結合流噪聲聲源主要來自于轉捩區(qū)的聲輻射這一共識，選擇頭部端面面積小而頭部較長的線型具有更好的噪聲性能。

2）對比各個速度各個線型A接收點（駐點）的流噪聲明顯小于B接收點的流噪聲，那么在潛艇首部聲吶的布置上，從降低自噪聲的出發(fā)點來看，盡量把聲吶布置在駐點附近，有利于提升聲吶的探測效能。

3）將數(shù)值仿真結果與試驗結果進行對比，說明本文計算流噪聲所用的方法在精度上具有一定的說服力。將這種數(shù)值仿真方法用于試驗前的優(yōu)化設計可以減少試驗方案，從而使設計者跳出考慮試驗高花費的限制，為大膽嘗試各種不同的線型提供了便利。

[1]繆旭弘,王振全.艦艇水下噪聲控制技術現(xiàn)狀及發(fā)展對策[C]//第十屆船舶水下噪聲學術討論會,北京,2005.M IAO Xu-hong,WANG Zhen-quan.The current situation and development strategies of underwater noise control technology in naval ships[C]//The10-th Symposium ofMarine Underwater Noise,Beijin,2005.

[2]黃高橋,潘光.水下航行器流噪聲特性水洞試驗研究[J].西北工業(yè)大學學報,2015,33(1):141–146.HUANG Gao-qiao,PAN Guang.W ater-tunnel experiment about flow-noise characteristic of underwater vehicle[J].Journal of Northwestern Polytechnical University,2015,33(1):141–146.

[3]武昊,潘光,黃橋高,等.水下航行器頭部流噪聲數(shù)值仿真與試驗[J].火力與指揮控制,2012,37(8):133–136.WU Hao,PAN Guang,HUANG Qiao-gao,et al.Numerical simulation and experiment research on underwater vehicle head flow noise[J].Fire Control&Command Control,2012,37(8):133–136.

[4]盧云濤,張懷新,潘徐杰.四種湍流模型計算回轉體流噪聲的對比分析[J].水動力學研究與進展,A 輯,2008,23(3):348–355.LU Yun-tao,ZHANG Huai-xin,PAN Xu-jie.Comparison between the simulations of flow-noise of a submarine-like body with four different turbulent models[J].Chinese Journal of Hydrodynamics,2008,23(3):348–355.

[5]蔣濤,馬軍,張萌.基于大渦模擬的潛艇流噪聲預測技術[J].海軍工程大學學報,2013,25(6):64–68.JIANG Tao,MA Jun,ZHANG Meng.Prediction of submarine hydrodynamic noise by using large eddy simulation[J].Journal of Naval University of Engineering,2013,25(6)64–68.

[6]LIGHTHILL M J.On sound generated aerodynamically I.General theory.Proc Royal Soc London,Ser A.1952,200(1107):564–587.

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

[8]孟堃宇.基于大渦模擬的潛艇脈動壓力與流噪聲性能數(shù)值計算[D].上海:上海交通大學,2011.

[9]劉繼明.水下航行器水動力噪聲分離預報與研究[D].武漢:華中科技大學,2014.

Numerical simulation research on flow-induced noise of underwater vehicle head

ZHANG Lei1,FAN Lin-xu1,2,SHEN Li-jiao3,ZHANG Jun-jie3,SUN Jiang-long1,2,WANG Chi-zhong4
(1.School of Naval Architecture and Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074,China;2.HubeiKey Laboratory of Naval Architecture and Ocean Engineering Hydrodynamics(HUST),Wuhan 430074,China;3.China Ship Development and Design Center,Wuhan 430064,China;4.Ocean College,Zhejiang University,Zhoushan 316021,China)

The decrease of flow-induced noise plays an important role in acoustic stealth condition of underwater vehicle and active sonar detection performance.The design of Underwater vehicle head,the important organ of underwater vehicle,directly affects the flow-induced noise performance and working circumstance of sonar.Based on the result of Large Eddy Simulation(LES),two different head lines are selected to research their flow-induced noise at different velocities,different incident angles and different testing point positions,by using software ACTRAN.Several conclusions are drawn through analysis.Flow-induced noise energy mainly concentrates on low frequency and increases with the increasing of velocity of underwater vehicle.At the same velocity and angle,underwater vehicle head with smaller radius and longer length owns a better flow-induced noise performance.For the same model,compared with other positions of front-section,the flowinduced noise overall sound pressure level of stagnation point is the smallest.Comparing numerical simulations with the Water-Tunnel experiment about Flow-Noise,the same conclusion can be drawn.

underwater vehicle；headline；Large Eddy Simulation(LES)；flow-induced noise

TP391

A

1672–7649(2017)12–0014–06

10.3404/j.issn.1672–7649.2017.12.004

2017–07–24

船舶振動噪聲重點實驗室資助項目（9140C280105150C28001）

張磊(1988–)，男，博士研究生，研究方向為流體力學、水動力噪聲。