潛艇近水面航行自由液面干擾效應(yīng)數(shù)值研究

2021-03-09 07:18:36周廣禮歐勇鵬

艦船科學(xué)技術(shù) 2021年1期

王陸，畢毅，周廣禮，向國，歐勇鵬

(1.海軍工程大學(xué) 艦船與海洋學(xué)院，湖北武漢 430033；2.中國人民解放軍 92578 部隊，北京 100161)

0 引言

作為現(xiàn)代海軍作戰(zhàn)力量的重要組成部分，潛艇在空間立體作戰(zhàn)及戰(zhàn)略威懾中均發(fā)揮著無可替代的作用[1]。當前按照動力來源，潛艇可分為核動力及常規(guī)動力兩大類，前者多在深海中隱蔽航行，其續(xù)航力及水下潛航時間均遠優(yōu)于常規(guī)動力潛艇，然而對于近海岸偵察、打擊及巡邏等軍事任務(wù)，各國則更多地應(yīng)用排水量相對較小的常規(guī)動力潛艇。受于作戰(zhàn)環(huán)境及戰(zhàn)術(shù)限制，近岸航行的常規(guī)動力潛艇往往需要完成近水面換氣發(fā)電、救援、發(fā)射導(dǎo)彈及潛望鏡深度瞭望等戰(zhàn)術(shù)動作，較于深水潛航狀態(tài)，近水面航行的潛艇受到文丘里效應(yīng)的影響，潛艇水動力隨潛深、航速均呈現(xiàn)較為明顯的變化。

在早期潛艇操縱過程中，國內(nèi)外學(xué)者已認識到潛艇近水面航行時水動力性能的變化[2]，并認為當潛深ζ≤3D（D為潛艇直徑）時，艇體受到較為明顯的吸力和力矩，且吸力隨下潛深度的增加呈指數(shù)規(guī)律衰減；隨著潛艇設(shè)計能力與技術(shù)的不斷提升，在較長的時間內(nèi)，很多學(xué)者采用切片理論及勢流方法來求解潛艇近水面狀態(tài)下的水動力性能[3-4]，目前仍有較多相關(guān)研究單位基于Havelock源[5]、Frank源[6]等的勢流方法及Ship-flow等勢流軟件[7]開展考慮波浪影響等復(fù)雜工況下的潛艇近水面水動力特性研究，但上述方法忽略了流體粘性的影響，適用范圍往往受到工況條件與計算經(jīng)驗的限制。

近年來，隨著粘流CFD技術(shù)的快速發(fā)展，基于Fluent，Star-CCM+及CFX等商用軟件平臺可快速獲取不同航態(tài)下潛艇的水動力性能，極大地提高了工程效率和研發(fā)周期，張楠等[8-9]較早地應(yīng)用粘流計算軟件開展了近水面航行潛艇數(shù)值計算，其重點分析了潛艇周圍流場特性及近水面對螺旋槳伴流的影響，并未對潛深及航速等因素展開分析；ZHANG[10]通過試驗及數(shù)值方法分析了AUV近水面航行時下潛困難的現(xiàn)象，分析了影響AUV下潛的主要因素；Wilson H S[11]則應(yīng)用CFX軟件以Suboff標模裸艇體為研究對象，開展了不同潛深及航速下潛艇水動力性能分析，但其選用模型長僅為1.556 m，其結(jié)論對于帶有附體及更大模型尺度潛艇是否適用還有待進一步研究。

1 數(shù)值方法

1.1 控制方程及湍流模型

文中以RANS方程作為求解潛艇近水面航行的基本方程，其具體形式如下：

式中： ρ為流體密度； μ為流體粘度；p為靜壓；fi為單位質(zhì)量的質(zhì)量力；ui和uj為速度分量。

湍流模式選取SSTk- ω模型，詳細推導(dǎo)過程參見文獻[12]，下面僅給出該模型的數(shù)學(xué)表達式：

式中： Γk和 Γω表示k和 ω 的有效擴散率；Gk表示由于平均速度梯度產(chǎn)生的湍流動能；Gω表示特殊湍流動能耗散 ω 的產(chǎn)生；Yk和Yω表示由于湍流k和 ω 的耗散；Sk和Sω為用戶自定義源項。

1.2 網(wǎng)格離散方法

Star CCM+商用軟件中提供了多種網(wǎng)格劃分方式，由于潛艇近水面航行計算的關(guān)鍵在于自由液面興波的準確捕捉，因此本文中選用六面體切割體網(wǎng)格劃分策略對計算域進行網(wǎng)格離散。

切割體網(wǎng)格技術(shù)的優(yōu)勢在于能夠確保流體域內(nèi)絕大部分網(wǎng)格為排列規(guī)則的六面體網(wǎng)格，同時通過對近艇體壁面區(qū)域進行網(wǎng)格加密設(shè)置可對潛艇周圍的湍流區(qū)域進行準確模擬，且對于過于復(fù)雜的幾何結(jié)構(gòu)，邊界層內(nèi)網(wǎng)格自適應(yīng)地采用其他形狀的網(wǎng)格結(jié)構(gòu)形式，以滿足物面形狀及網(wǎng)格質(zhì)量的要求。其生成原理及步驟如下：

步驟1首先根據(jù)艇周圍加密區(qū)設(shè)置及域邊界確定網(wǎng)格的生成邊界，如圖1（a）所示，邊界層外邊界與域邊界圈定的區(qū)域為中心網(wǎng)格區(qū)，艇體壁面與邊界層外邊界間的區(qū)域為邊界層區(qū)；

步驟2在整個計算域及周圍生成初始體網(wǎng)格（見圖1（a）），用第1步確定的邊界進行裁剪，裁剪后的網(wǎng)格如圖1（b）所示；

步驟3根據(jù)網(wǎng)格密度設(shè)置對中心網(wǎng)格區(qū)內(nèi)網(wǎng)格疏密程度進行控制，并進一步提高網(wǎng)格質(zhì)量；

步驟4在邊界層區(qū)域內(nèi)，根據(jù)物面形狀自適應(yīng)地生成網(wǎng)格。

圖1 網(wǎng)格生成示意圖Fig.1 Configuration of meshing strategy

2 研究對象及計算工況

以美國國防高級研究計劃局（DARPA）公布的Suboff全附體模型為研究對象，該模型由主艇體、指揮臺圍殼及十字型尾舵組成，模型主參數(shù)如圖2所示。將潛艇中心線距自由液面的距離與艇體直徑的比值定義為無因次潛深（見圖3），即H*=H/D；無因次航速用艇長傅氏數(shù)表示，即式中U為潛艇航速，L為艇長，g為重力加速度，取值為9.81 m/s-2。

圖2 Suboff潛艇模型主尺度參數(shù)Fig.2 Main particulars of Suboff

圖3 潛深無因次化定義Fig.3 Definition of submarine’s non-dimensional depth

本文開展傅氏數(shù)Fr=0.1～0.5、潛深H*=1.0～2.9D及無限潛深工況條件下潛艇近水面航行水動力數(shù)值計算。

3 計算域設(shè)置及網(wǎng)格劃分

由于Suboff潛艇模型及流場均關(guān)于艇體中縱剖面對稱，為減小計算量取一半潛艇開展計算，如圖4所示。計算域大取為長寬高為5.5L×1.5L×3.0L的方形域，邊界條件設(shè)置如下：

1）距離艇首1.5倍艇長，邊界條件設(shè)置為速度入口；

2）距離艇尾4倍船長，邊界條件設(shè)置為壓力出口，壓力為未擾動時邊界壓力；

3）流域的上、下及右邊界設(shè)置無滑移、不可穿透壁面；

4）流域左邊界設(shè)置為對稱面；

5）艇體表面定義為無滑移、不可穿透邊界條件，邊界條件設(shè)置為無滑移壁面。

圖4 計算域大小及邊界條件設(shè)置Fig.4 Computational domain and boundary conditions

在網(wǎng)格離散過程中，為準確地捕捉兩相流的水氣分界面，對自由液面附近區(qū)域進行網(wǎng)格加密。此外，在艇體近壁面設(shè)定了局部加密區(qū)，確保壁面y+～1，并對指揮臺及尾舵處的網(wǎng)格進行了細化處理，船體表面及域剖面網(wǎng)格如圖5所示，總網(wǎng)格數(shù)為215萬。

圖5 艇體表面及域剖面網(wǎng)格Fig.5 Mesh for submarine's surface and domain profile

4 計算方法可行性分析

為驗證上述計算方法的可行性，以Suboff模型裸艇體為計算對象，模型長1.556 m，試驗數(shù)據(jù)參見文獻[11]。計算域設(shè)置及模型無因次潛深的定義與前述模型相同。圖6為該模型無因次潛深H*=1.1，F(xiàn)r= 0.23～0.64工況下艇體阻力的計算結(jié)果與試驗值。對比可知，本文方法可較為準確地計算潛艇近自由液面的阻力性能。

圖6 阻力計算值與試驗結(jié)果對比（H* = 1.1）Fig.6 Comparison of calculation and experimental results for submarine’s resistance at H* = 1.1

由圖6中近水面裸艇體阻力隨傅氏數(shù)的變化可知，隨著航速的增大，潛艇近水面航行時艇體阻力呈現(xiàn)先增大后減小再增大的變化趨勢，這與深水航行存在較大差異。

5 計算工況與結(jié)果分析

5.1 力、力矩定義

為便于計算與分析，設(shè)定潛艇重心與浮心縱向位置相同，即距艇首2.009 5 m，并定義潛艇近水面航行時所受阻力、垂向力及對重心的縱傾力矩的正方向如圖7所示。其中，潛艇垂向力中不包含艇體所受浮力。力及力矩系數(shù)定義如下：

圖7 潛艇受力及力矩定義Fig.7 Definition of force and torque for submarine

式中：R′， Z′及M′分別為阻力系數(shù)、縱傾力矩系數(shù)，ρ為流體密度，L為艇長，U為潛艇航速。

5.2 航速的影響

圖8～圖10分別給出了不同潛深下潛艇阻力系數(shù)、垂向力系數(shù)及縱傾力矩系數(shù)隨傅氏數(shù)Fr的變化規(guī)律，分析可得：

1）潛艇無因次潛深H*≤1.6時，阻力系數(shù)隨Fr呈先增大后減小再增大的趨勢，而當潛深H*≥1.9時該變化規(guī)律并不明顯，在無限潛深下，潛艇阻力系數(shù)隨傅氏數(shù)逐漸減小并最終趨于穩(wěn)定值。

2）不同潛深下，潛艇垂向力系數(shù)隨傅氏數(shù)Fr的變化規(guī)律基本一致，除無限水深工況外，均在Fr=0.4時出現(xiàn)峰值狀態(tài)。由圖也可看出，無因次潛深H*=2.9時，自由液面對潛艇向上的吸力作用較弱，且隨航速變化垂向力系數(shù)變化幅度較小。

圖8 阻力系數(shù)隨傅氏數(shù) Fr 的變化Fig.8 Resistance coefficient of submarine at different Fr

圖9 垂向力系數(shù)隨傅氏數(shù) Fr的變化Fig.9 Vertical force coefficient of submarine at different Fr

圖10 縱傾力矩系數(shù)隨傅氏數(shù) Fr的變化Fig.10 Trimming moment coefficient of submarine at different Fr

3）在無限水深下，潛艇縱傾力矩系數(shù)均遠小于近水面航行工況，其值趨近于0；而在近水面工況下（H*≤2.9），當傅氏數(shù)Fr≤0.3時，潛艇所受縱傾力矩系數(shù)M′隨航速的變化規(guī)律并不明顯，而當Fr≥0.3時，縱傾力矩系數(shù)絕對值隨航速逐漸增大。

5.3 潛深的影響

通過對比分析不同潛深下潛艇水動力的差異可直觀分析潛艇水動力隨潛深的變化趨勢，潛艇阻力系數(shù)、垂向力系數(shù)及縱傾力矩系數(shù)隨無因次潛深的變化曲線如圖11～圖13所示。

由圖可知：

1）在低航速下（Fr≤0.3），潛艇阻力系數(shù)隨潛深波動變化，而在Fr=0.4與Fr=0.5工況下，潛艇阻力系數(shù)隨潛深增大而減小；在研究航速范圍內(nèi)，當潛深H*=2.9時，潛艇阻力系數(shù)與無限潛深下潛艇的阻力系數(shù)較為相近，即可認為當潛深H*≥2.9時，潛艇阻力性能受自由液面的影響較小。

2）相較于阻力性能，潛艇近水面航行時受到向上的垂向吸力，且一定航速下，隨著潛深的不斷增大，自由液面的吸力作用逐漸減弱，并在H*≥2.9后的較深水域趨于0。

圖11 阻力系數(shù)隨潛深的變化Fig.11 Resistance coefficient of submarine at different depths

圖12 垂向力系數(shù)隨潛深的變化Fig.12 Vertical force coefficient of submarine at different depths

圖13 縱傾力矩系數(shù)隨傅氏數(shù) Fr的變化Fig.13 Trimming moment coefficient of submarine at different depths

3）由圖13可知，在研究航速及潛深范圍內(nèi)，潛艇均受到埋首的縱傾力矩，而當Fr≤0.3時，潛艇所受縱傾力矩值較小，且隨潛深的變化規(guī)律并不明顯；在Fr=0.5工況下，潛艇縱傾力矩系數(shù)絕對值隨潛深逐漸下降，潛深H*=2.9工況下潛艇縱傾力矩值約為H*=1.0時的3倍。

5.4 流場分析

由潛艇所受力及力矩分析可知，不同潛深下潛艇近水面航行時的水動力隨航速的變化規(guī)律基本趨于一致，潛深越小，潛艇水動力受自由液面的影響也越大，因此，取定潛深H*=1.0下，不同航速下潛艇表面壓力分布進行對比分析，如圖14所示。在近自由液面工況下，低航速時（Fr=0.1）艇體表面壓力主要受到靜水壓力的影響，而隨著航速的增大，艇體上表面壓力沿船長呈先增大后減小的趨勢，在Fr=0.5航速下，艇體尾部上表面低壓區(qū)較為明顯。艇體前后的壓力差使得潛艇在近水面高速航行時受到較大的埋首力矩，此現(xiàn)象與前述數(shù)據(jù)分析所得結(jié)論較為吻合。

圖14 不同航速下主艇體表面壓力分布（H* = 1.0）Fig.14 Pressure distribution on submarine's hull at different advancing velocity（H* = 1.0）

同等航速下，潛深是影響潛艇水動力的重要因素，圖15給出了航速Fr=0.5工況下，不同潛深時自由液面波高圖。分析可知，在H*=1.0工況下，潛艇興波較為明顯，艇首波高為正，艇尾尾舵處存在較明顯的負波高區(qū)域，而在艇后又會出現(xiàn)扇形波。隨著潛深的不斷增大，艇體興波衰減較快，至H*=2.9潛深下，興波只呈現(xiàn)較小的波高增值區(qū)域。