開孔分布對空腔流動特性的分析

張 航，馮大奎，蔡 瀟，肖志堅

(華中科技大學 船舶與海洋工程學院，湖北 武漢 430074)

開孔分布對空腔流動特性的分析

張 航，馮大奎，蔡 瀟，肖志堅

(華中科技大學 船舶與海洋工程學院，湖北 武漢 430074)

開孔分布是影響空腔流動的一個重要因素。為了對開孔空腔流動有更深的認識，采用大渦模擬（LES）的方法，以Suboff艇體母線建立二維模型，研究 4 種開孔分布對空腔流動阻力，頻譜特性及內外流交換的影響。對計算結果的分析表明，由于艇體表面的壓力分布不同，孔附近產生縱向壓力差，促使空腔內外流動交換，增加主艇體首尾壓差阻力，進而使得總阻力增大。計算結果表明艇體阻力增加與內外流交換的密切相關，開孔位于中部時總阻力增量最小，內外流增量最小，開孔均勻分布時引起內外流流動交換劇烈，阻力增量最大，而且開孔引起總阻力波動幅值增加，頻率分布特性發生相應的改變，開孔使得大幅波動頻帶變寬，可以預測噪聲強度增加，頻帶變寬。

大渦模擬；空腔流；頻譜特性；阻力增值

0 引 言

近年來，空腔結構的阻力分布及空腔附近流交換問題成為研究的熱點問題，空腔在海洋航行器和航空飛行器上廣泛應用，如潛艇上的流水孔等。潛艇的流水孔是指在潛艇的非耐壓非水密結構上，用于潛艇上浮下潛時供液體自由進出的開口。當流體流經流水孔時，流場的非定常性會產生各種尺度的漩渦及湍流引發的脈動，作用于艇體外表面導致壓力波動，破壞表面湍流附面層內的流動進而形成流噪聲，其對潛艇航行隱身性能產生一定的影響[1]，很可能會暴露潛艇的位置，嚴重影響潛艇的作戰能力和隱身性能。

空腔流動包含多種復雜的物理現象，如剪切層不穩定性、邊界層分離、漩渦運動、激波邊界層互相干擾等。研究空腔流動的方法主要有實驗方法、數值分析等方法。賴煥新等[2]運用大渦模擬的數值分析方法對空腔流動及氣動噪聲進行研究；鄔明[3]以潛艇流水孔為對象研究了流水孔引起的阻力增量集中在導流板和流水孔內壁上；陳燦等[4]運用大渦模擬研究了不可壓縮空腔流的特性及噪聲；孟生等[5]以潛艇流水孔為研究對象，通過改變流水孔的模型參數，得出了流水孔的噪聲特性；張楠等[6]以潛艇流水孔為對象通過求解RANS方程，結合PISO算法精確仿真了流水孔內外流場的情況。Rossiter等[7]通過實驗得出了空腔內存在渦和噪聲的反饋循環，從而提出了空腔流動振蕩頻率的模型，并給出了預估振蕩頻率的半經驗公式；涂海文等[8]運用RNG，k-ε湍流模型對潛艇的阻力分布和流場分布進行了模擬，為進一步優化潛艇的艇型和分析潛艇的流噪聲打下了基礎；李霖等[9]運用大渦模擬對于圓柱繞流進行了模擬，顯示出這種數值模擬方法的優勢。

本文采用大渦模擬方法，以Suboff艇型為對象，研究開孔縱向分布對空腔流動特性的影響，包括內外流的交換，頻譜特性及其對阻力的影響。

1 幾何模型

取Suboff主艇體母線建立二維模型如圖1所示。

從左至右孔依次標記為①②③④⑤，蓋板如圖所示標記為①②③④，擋板與蓋板延長線的夾角為45o，船長為4.356 m。

分別對光體、孔集中分布在平行中體首部、集中在平行中體中間、集中在平行中體尾部和均勻分布在平行中體上5種工況進行研究。各工況具體如表1所示。

2 數值分析方法

2.1 大渦模擬控制方程

大渦模擬（LES）主要利用濾波的方法將渦分為大尺度渦和小尺度渦，然后分別進行數值模擬分析。大渦模擬控制方程為連續性方程以及N-S方程進行過濾之后的方程：

表1 工況Tab.1 Conditions

本文采用Smagorinsky渦粘模型來模化亞格子應力，其張量的偏量部分為：

式中：vi為渦黏系數Cs為Smagorinsky常數；?為濾波尺度，通常取0.1～0.23；Sij為經過濾波后的速度變形張量，定義如下：

2.2 網格及計算域

整個流域采用非均勻的結構化網格，對于靠近艇體和空腔的擋板和蓋板處進行了網格加密，邊界層厚度約為0.8 mm。 5 種工況的網格量約為1 650 000，y+最大值大約為0.6，計算域上游取3倍艇長，下游取9倍艇長，寬度為3倍艇長。

2.3 邊界條件

本文對空腔數值模擬的邊界條件如下：

1）入口條件：入口處設置為速度入口，速度大小均為4 m/s，速度方向垂直于入口邊界，為x軸正方向。

2）出口條件：出口處設置為壓力出口。

3）壁面條件：艇體主體、水槽、擋板、蓋板均設置為無滑移壁面條件。

4）對稱條件：流域頂部和流域下部均設置為對稱邊界條件。

5）內部條件：5 個孔和內部流場均設置為內部條件。

2.4 數值離散

采用大渦模擬計算方法，壓力項采用二階隱式格式求解，時間項采用二階隱式，動量項采用中心差分求解，時間步長為0.000 25，每一時間步最大迭代次數為15，迭代收斂標準為0.000 001。 使用收斂性較好的SIMPLEC求解壓力耦合方程。

3 計算結果

3.1 各工況下各部位阻力

將各工況各部位阻力求解結果如表2所示。

表2 各工況各部位阻力情況Tab.2 The resistance of different parts on each condition

由表2可知，工況1與開孔工況2～工況5對比發現，工況1的總阻力最小，而且阻力增量主要在艇體上；將4種開孔工況之間對比發現，工況3的總阻力相對最小，且艇體的阻力增量也最小，這跟此種開孔分布下的孔的內外流交換和壓力分布都有關，此種工況下5個孔的流量最小，說明內外流交換最小，由內外流引起的主艇體的縱向壓差最小，因而阻力增量最小。

3.2 各工況的壓力分布與流線圖

面對長形式圖書閱讀過程中的問題，圖書出版業結合互聯網時代閱讀的大環境，開始從供給側的角度做相應的對接。

各工況流線圖如圖3～圖7所示。

由圖3～圖7可看出，工況1艇體附近幾乎沒有渦，流線沿艇體流過，并且越靠近艇體，流線越密集；工況2流線圖在艇體首部出現渦，這些渦分布在開孔附近；工況3流線圖上艇體中間部位出現渦，并且這些渦出現在開孔附近；工況4流線圖上艇體尾部出現渦，并且這些渦出現在開孔附近；工況5流線圖上艇體空腔內各部位都出現渦，這些渦的位置與開孔位置吻合。由于渦分布在孔的附近，即渦的形成位置與孔的位置有關，說明由于孔附近的內外流交換產生了渦，當水從一個孔流入后，在空腔內流動，受到空腔壁的阻擋改變方向，做螺旋運動，最終由其他孔流出，進而形成渦。

各部位壓力分布如圖8所示。由圖可知，各工況在開孔部位壓力變化較為劇烈，并引起首尾部部分壓力分布發生改變。

3.3 各工況下5個開孔流量

流量定義流入為負，流出為正，各工況下5個開孔流量如表3所示。

表3 各工況各孔流交換情況Tab.3 Flow exchange of each hole on each condition

由表3可知，工況2，水從孔③④⑤流入，從孔①②流出；工況3，水從孔③④⑤流入，從孔①②流出，但相比工況2，此種工況下5個孔的流量都較小，說明孔附近流交換小；工況4，水從孔①②③流入，從孔④⑤流出；工況5，水從孔②③④流入，從①⑤流出。綜合這 4 種開孔分布工況可知，工況3的孔的流量最小，即此種工況下孔的內外流交換最小。

3.4 頻譜分析

各個工況的頻譜曲線如圖9～圖13所示。

4 結果分析

結合圖3～圖8及表2并將光體工況與開孔工況對比知，開孔后阻力增加，而且阻力增量主要在主艇體上，這與開孔后孔附近流動交換有關。開孔后，水會從一部分孔流入空腔，從其他孔流出，從而影響首尾部壓力分布，進而影響主艇體壓差阻力，使得開孔后總阻力比光體大。

結合表2表3及流線圖可知，4 種開孔工況下，孔集中分布在中部時，總阻力最小，而且相比其他工況主要是艇體上阻力較小。由于水按照順壓流動，孔開在首部時，水從后部孔流入空腔，從靠近首部孔流出，與首部來流交匯阻礙了首部流動速度，從而使首部壓力增大，艇體首尾壓差增大，艇體壓差阻力增大，而摩擦阻力幾乎不變，從而主艇體總阻力增大；當孔開在尾部時，距離尾部較近，水從遠離尾部的孔流入，從靠近尾部孔流出，從而使得尾部水流動加快，尾部壓力相應降低，而首部壓力基本不變，故首尾壓差增大，主艇體壓差阻力增大，總阻力增大；當孔開在中部時，由于距離首尾均較遠，對于首尾壓力影響較小，所以阻力增量也最小；當孔均勻分布時，水從中間幾個孔流入，從靠近首尾的 2 個孔流出，分別阻塞首部流動，加速尾部流動，使得首部壓力增大，尾部壓力減小，首尾壓差明顯增大，主艇體壓差阻力增大，總阻力增加最為顯著。

由圖9～圖13可知，工況1阻力波動幅度較大頻帶較窄，約為f=13.6 Hz，阻力幅值峰值為10.1 N；開孔以后，阻力波動幅度較大的頻帶明顯增加，波動幅值相比于光體增加明顯。主要因為開孔引起的空腔內外流交換明顯；工況2下阻力波動幅值較大頻率范圍10～25 Hz，阻力幅值峰值為15.6 N；工況3下阻力波動幅值較大頻率范圍15～25 Hz，阻力幅值峰值為11.3 N；工況4下峰值頻率為13～30 Hz，阻力幅值峰值為13.5 N，開孔后峰值處頻率向高頻率值移動，并且開孔后阻力幅值均比光體阻力幅值大，工況5工況阻力波動幅值較大頻率范圍10～30 Hz，在這個頻率區間內沒有明顯的峰值頻率，最大波動幅值約為15 N，主要因為開孔引起大范圍的內外流的交換。4 種開孔分布之間對比，發現開孔集中分布在平行中體中間時阻力幅值最小，而且波動最小，說明此種開孔分布對阻力頻譜特性影響較小。

5 結 語

本文采用大渦模擬的數值分析方法，以Suboff艇型為對象，研究平行中體上開孔分布對于空腔流動特性的影響，研究了光體、開孔集中分布在平行中體左邊、開孔集中分布在平行中體右邊、開孔集中分布在平行中體中間、開孔均勻分布在平行中體上 5 種工況下的壓力分布、阻力情況、流交換情況以及阻力的頻譜特性。得出如下結論：

1）開孔后阻力相較于光體總阻力增大，而且當開孔位于平行中體中間時總阻力的增量最大，而且阻力增量是由于開孔附近流交換引起的首尾壓差變化，進而使壓差阻力增大，對于真實的潛艇設計時，流水孔設計位置因與光體壓力分布結合起來，盡量避免內外流的交換。

2）開孔引起的阻力增加主要是主艇體阻力增加，而非開孔部位的阻力增加，因此潛艇設計時注意開孔與主艇體的配合。

3）開孔分布對艇體水動力噪聲有較大影響，不僅增加了噪聲強度，而且改變了峰值頻率的分布。

[1]曹志敏, 杜陽華. 潛艇武器系統綜合仿真系統及關鍵技術[J].計算機仿真, 2007, 24(6): 31–37.

[2]賴煥新, 周邵萍, 蘇永升, 等. 空腔流動的大渦模擬及氣動噪聲控制[J]. 工程熱物理學報, 2008(2): 228–232.

[3]鄔明. 帶流水孔潛艇流噪聲特性分析[J]. 四川兵工學報, 2015(10): 27–31.

[4]陳燦, 吳方良, 李環, 等. 不可壓縮空腔流振蕩模式和聲學特性研究[J]. 水動力學研究與進展(A輯), 2015(3): 272–278.

[5]孟生, 張宇文, 王艷峰. 潛艇流水孔流噪聲特性研究[J]. 計算機仿真, 2011(7): 57–62.

[6]張楠, 沈泓萃, 姚惠之, 等. 帶流水孔潛體流場數值模擬[J]. 船舶力學, 2004(1): 1–11.

[7]ROSSITER J E. Wind tunnel experiments on the flow over rectangular cavities at subsonic and transonic speeds[R]. London, UK: Aeronautical Research Council, 1966.

[8]涂海文, 孫江龍. 基于CFD的潛艇阻力及流場數值計算[J].艦船科學技術, 2012,34(3): 19–25.

[9]李霖, 張志國, 王先洲, 等. 低雷諾數圓柱繞流的大渦模擬分析[J]. 艦船科學技術, 2013, 35(1): 22–26.

Analysis of cavity flow characteristics by hole distribution

ZHANG Hang, FENG Da-kui, CAI Xiao, XIAO Zhi-jian
(School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China)

Cavity flow depends greatly on hole distribution.To make hole cavity flow more clear, a large eddy simulation (LES) method is used to establish a two-dimensional model of SUBOFF hull busbar to study the effect of four kinds of hole distribution on cavity flow resistance, spectral characteristics and internal and external flow exchange in this paper. The analysis of the results shows that, due to the different pressure distribution on the surface of the hull, the longitudinal pressure difference is generated in the vicinity of the hole, which promotes the exchange of the inside and the outside of the cavity and increases the head-end pressure difference resistance of the main body. The results show that the increase of hull resistance is closely related to the exchange of internal and external flow. When the opening is in the middle, the total resistance increase is the smallest, the increment of internal and external flow is the smallest, and the flow of the internal and external flow is intense and the resistance is the largest. The opening of the hole causes the total resistance fluctuation amplitude to increase, the frequency distribution characteristic changes accordingly, the opening makes the large fluctuation frequency band widens, can predict the noise intensity increases, the frequency band widens.

large eddy simulation (LES)；cavity flow；spectral characteristics；resistance increment

O352

A

1672 – 7619(2017)06 – 0048 – 05

10.3404/j.issn.1672 – 7619.2017.06.010

2017 – 02 – 27

張航(1996 – )，男，本科生，研究方向為艦船與海洋工程。