附體對潛艇阻力性能的影響研究

徐 妍，高霄鵬

(海軍工程大學 艦船與海洋學院，湖北 武漢 430033)

0 引 言

潛艇在航行時，會受到流體阻礙其運動的阻力，潛艇航行過程中的阻力性能會影響潛艇的快速性。隨著各種反潛設(shè)備的發(fā)展，潛艇的航行安全問題不容忽視，提高潛艇的快速性已經(jīng)成為各國重要的軍事研究課題。國內(nèi)外眾多專家通過各種方法對潛艇阻力性能進行了研究[1-3]。隨著計算機技術(shù)的快速發(fā)展，數(shù)值方法開始廣泛應(yīng)用，對潛艇阻力的預(yù)報精度也有了更高的要求，潛艇阻力的計算一直是研究熱點。

潛艇表面布置有眾多附體，其形狀比較復(fù)雜，附體的存在會影響潛艇周圍的流場，從而對水動力性能產(chǎn)生影響，因此分析不同附體對潛艇阻力性能的影響便十分重要。胡斌等[4]基于分塊網(wǎng)格技術(shù)對全附體SUBOFF 潛艇模型的粘性繞流場進行了研究，分析了附體對流場的干擾。涂海文等[5]對帶有不同附體的潛艇模型的摩擦阻力系數(shù)進行了預(yù)報，并將數(shù)值模擬結(jié)果與經(jīng)驗公式計算結(jié)果進行了對比分析。王慶云等[6]針對不同舵翼位置的系列舵翼潛艇，進行了阻力性能模擬，分析了隨著舵翼的后移，潛艇摩擦阻力、壓差阻力和總阻力的變化情況。

本文基于STAR-CCM+對帶不同附體潛艇直航運動進行了數(shù)值仿真，探究附體對潛艇阻力性能的影響。

1 數(shù)值理論

1.1 控制方程

對不可壓縮粘性流體進行CFD 數(shù)值計算時，需要構(gòu)建基本控制方程，包括連續(xù)性方程和RANS 方程，其表達式如下：

式中：uiuj為速度分量；Si為廣義源項；為雷諾應(yīng)力。

1.2 湍流模型

湍流模型的選擇直接影響數(shù)值仿真結(jié)果的準確性，本文采用SSTk-ω湍流模型對RANS 方程進行封閉，從而進行數(shù)值仿真。SSTk-ω模型屬于兩方程模型，結(jié)合了標準k-ε 模型和k-ω模型的優(yōu)點，對湍流輸運方程進行了修正，不僅適用于近壁面的計算，也適用于遠處流場。SSTk-ω湍流模型的輸運方程可表示為[7]：

式中，Γk，Γω分 別為k和 ω 的有效擴散系數(shù)；GkGω為k和 ω 產(chǎn) 生的湍流動能；YkYω為 擴散項；SkSω為自定義項；Dω為交叉擴散項。

2 計算模型

本文的研究對象以SUBOFF 全附體模型為基礎(chǔ)，通過加裝左右對稱的圍殼舵而成。該潛艇模型包括主艇體、指揮臺圍殼、圍殼舵及4 個完全相同的尾翼，模型總長4.356 m，寬0.508 m，模型的三維視圖如圖1 所示。潛艇模型主要參數(shù)見表1。

表1 潛艇模型主要參數(shù)Tab.1 Main parameters of the submarine model

圖1 潛艇模型三維視圖Fig.1 The 3D view of submarine model

為研究附體對潛艇阻力性能的影響，通過拆附體來構(gòu)建不同的潛艇模型，本文共有7 種潛艇模型方案。

表2 潛艇模型方案Tab.2 Example of submarine models

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

3.1 計算域設(shè)置

完成各SUBOFF 潛艇模型的幾何建模后，需要建立合適的計算域來對潛艇周圍流場進行模擬。在設(shè)置計算域時應(yīng)當注意控制計算域的大小。以艇長L=4.356 m為特征長度，潛艇上浮計算域的尺寸設(shè)置如圖2 所示。

圖2 計算域及邊界設(shè)置Fig.2 Computational domain and boundary

3.2 邊界條件

邊界條件設(shè)置如下：

1）潛艇表面設(shè)置為無滑移壁面；

2）潛艇尾部的計算域邊界設(shè)置為壓力出口，與無限遠處流場相對壓力為0；

3）其余計算域邊界設(shè)置為速度入口，入口設(shè)置均勻來流，來流速度u=U0，v=w=0，其中U0分別取2，3，4，5 和6 m/s。

3.3 網(wǎng)格劃分

采用切割體網(wǎng)格對計算域進行離散，采用棱柱層網(wǎng)格劃分潛艇表面邊界層的網(wǎng)格。由于潛艇指揮臺圍殼、尾舵附近流場變化較劇烈，因此對指揮臺圍殼及尾流場附近等進行網(wǎng)格加密，合理布置網(wǎng)格，提高計算結(jié)果的精度。艇體表面網(wǎng)格劃分如圖3 所示。

圖3 計算域網(wǎng)格設(shè)置Fig.3 Mesh for computational domain

4 計算結(jié)果及分析

潛艇在無限流場中作直航運動時，會受到流體阻礙其運動的阻力。由于本文是對無限流場中潛艇的阻力性能進行研究，因此不考慮興波阻力。按阻力的成因來分，潛艇所受到的阻力Rt包括摩擦阻力Rf和粘壓阻力Rpv，其關(guān)系如下：

在給定來流速度U0=2，3，4，5，6 m/s 下，通過數(shù)值仿真計算可得到直航運動時各模型在不同流速下的摩擦阻力、粘壓阻力和總阻力。各模型的摩擦阻力和粘壓阻力隨速度的變化曲線如圖4 和圖5 所示。

圖4 不同附體模型的摩擦阻力Fig.4 The friction resistance of model with difference appendages

圖5 不同附體模型的粘壓阻力Fig.5 The viscous pressure resistance of model with difference appendages

由圖4 可知，隨著速度的增加，各模型的摩擦阻力均有增加的趨勢，且不同模型摩擦阻力之間相差較小，說明附體對模型摩擦阻力的影響比較小。

隨速度的增加，各模型的粘壓阻力也呈現(xiàn)出增大的趨勢，粘壓阻力約與速度的平方成比例。從圖中可以看出，不同模型的粘壓阻力相差較大，說明附體顯著影響潛艇的粘壓阻力，隨著速度的增大，附體對粘壓阻力的影響也隨之增大。

4.1 附體對潛艇阻力系數(shù)的影響

為了便于對比分析，通過式（5）～式（7）對各模型的阻力系數(shù)進行求解。

式中：Cf為摩擦阻力系數(shù)；Cp v為粘壓阻力系數(shù)；Ct為總阻力系數(shù)；S為潛艇模型的濕表面積。

以U0=4 m/s 為例，各SUBOFF 潛艇模型的摩擦阻力系數(shù)、粘壓阻力系數(shù)和總阻力系數(shù)情況如圖6 所示。

圖6 各潛艇模型阻力系數(shù)柱狀圖Fig.6 The histogram of resistance coefficient of model

在來流速度一定時，不同模型的摩擦阻力系數(shù)相差較小，說明附體對各潛艇模型的摩擦阻力系數(shù)影響較小。對于粘壓阻力系數(shù)，附體的增減對其影響較大，按粘壓阻力系數(shù)大小排列：SUBOFF1 ＞ SUBOFF2 ＞SUBOFF4 ＞ SUBOFF3 ＞ SUBOFF6 ＞ SUBOFF5 ＞ SUBOFF7，根據(jù)不同附體布置對模型粘壓阻力系數(shù)的影響，按影響由大到小排列分別為：尾舵、指揮臺圍殼、尾方向舵、尾升降舵、圍殼舵，這說明不同附體布置對模型粘壓阻力系數(shù)影響不同。

在7 種潛艇模型中，SUBOFF1 潛艇模型的總阻力系數(shù)最大，是裸艇體SUBOFF7 模型的1.13 倍，分析可得附體的布置影響潛艇阻力性能，尾舵的布置對阻力性能的影響最大。

4.2 附體對潛艇表面壓力系數(shù)的影響

進一步研究不同附體布置對潛艇直航時表面壓力分布的影響，壓力系數(shù)Cp的計算公式如下：

通過對各模型進行直航數(shù)值仿真，監(jiān)測了艇體表面的壓力分布，得到各模型的表面壓力系數(shù)分布。選取具有代表性的壓力系數(shù)分布曲線如圖7～圖12 所示。其中，圖7～圖10 為中縱剖面縱向壓力系數(shù)分布曲線，圖11～圖12 為橫剖面縱向壓力系數(shù)分布曲線。計算表明：

圖7 SUBOFF1 中縱剖面縱向壓力系數(shù)分布Fig.7 The pressure coefficient of profile section of SUBOFF1

圖8 SUBOFF3 中縱剖面縱向壓力系數(shù)分布Fig.8 The pressure coefficient of profile section of SUBOFF3

圖9 SUBOFF6 中縱剖面縱向壓力系數(shù)分布Fig.9 The pressure coefficient of profile section of SUBOFF6

圖10 SUBOFF7 中縱剖面縱向壓力系數(shù)分布Fig.10 The pressure coefficient of profile section of SUBOFF7

圖11 SUBOFF1 橫剖面縱向壓力系數(shù)分布Fig.11 The pressure coefficient of cross section of SUBOFF7

1）潛艇首部壓力系數(shù)存在突變現(xiàn)象。潛艇首部壓力較大，這是由于潛艇首部為迎流面，首部流體流速幾乎為零。

2）指揮臺圍殼前壓力系數(shù)變化率較大，表明帶指揮臺圍殼的潛艇在圍殼前局部壓力較高，原因是指揮臺圍殼阻礙流體運動，流體的流速較小。

3）尾舵前壓力系數(shù)變化率較大，表明帶尾舵的潛艇在尾舵前局部壓力較高，原因是尾舵阻礙流體運動，流體的流速較小。

4）圍殼舵的設(shè)置對潛艇中縱剖面壓力系數(shù)影響較小。

為直觀表示潛艇首部、指揮臺圍殼前以及尾舵前壓力分布情況，通過數(shù)值仿真中的壓力云圖來表示。以SUBOFF1 為例，潛艇直航運動中艇體表面壓力分布情況如圖13 所示。

圖13 SUBOFF1 潛艇模型表面壓力分布云圖Fig.13 The surface pressure profile of SUBOFF1 submarine model

4.3 附體對潛艇周圍速度場的影響

為了更加直觀把握附體的布置對潛艇直航運動阻力性能的影響，考慮到潛艇模型左右對稱，選取潛艇模型的中縱剖面、橫剖面及圍殼舵處橫剖面對潛艇周圍速度場進行分析。圖14～圖17 為潛艇中縱剖面速度分布，圖18 和圖19 為橫剖面速度分布，圖20 和圖21為圍殼舵處橫剖面速度分布，結(jié)論如下：

圖14 SUBOFF1 中縱剖面速度云圖Fig.14 The velocity of middle longitudinal section of SUBOFF1

圖15 SUBOFF3 中縱剖面速度云圖Fig.15 The velocity of middle longitudinal section of SUBOFF3

圖16 SUBOFF6 中縱剖面速度云圖Fig.16 The velocity of middle longitudinal section of SUBOFF6

圖17 SUBOFF7 中縱剖面速度云圖Fig.17 The velocity of middle longitudinal section of SUBOFF7

圖18 SUBOFF1 橫剖面速度云圖Fig.18 The velocity of cross section of SUBOFF1

圖19 SUBOFF7 橫剖面速度云圖Fig.19 The velocity of cross section of SUBOFF7

圖20 SUBOFF1 圍殼處橫剖面速度云圖Fig.20 The velocity of cross section around rudder of SUBOFF1

圖21 SUBOFF2 圍殼處橫剖面速度云圖Fig.21 The velocity of cross section around rudder of SUBOFF2

1）在潛艇首部、指揮臺圍殼前部和尾舵前速度等值線分布比較密集，速度變化大，會引起潛艇表面壓力分布變化，出現(xiàn)突變現(xiàn)象，從而影響潛艇阻力性能。

2）隨著附體的增加，原本在艇體首部和尾部的速度等值區(qū)域逐漸增大，這說明附體的布置會對潛艇中縱剖面流場產(chǎn)生擾動現(xiàn)象。

3）指揮臺圍殼的布置使?jié)撏е車鲌龈淖儯辉傧鄬τ跈M剖面對稱，指揮臺圍殼前部和后部速度等值線比較密集，表明此處潛艇表面壓力變化大，引起潛艇所受的阻力增加。另外，指揮臺圍殼后速度等值區(qū)域面積增大，說明指揮臺圍殼的設(shè)置對指揮臺圍殼后的潛艇繞流場以及潛艇的尾流場有一定的影響。

4）尾舵的設(shè)置使得艇體尾部尾舵后方收縮急劇，沿縱向速度變化劇烈，尾部速度等值區(qū)域面積增大，艇體表面壓力變化較大，呈現(xiàn)突變現(xiàn)象。尾部壓力減小，首尾壓力差增大，從而使粘壓阻力增大。

5）布置圍殼舵會改變流場速度分布情況，進而引起潛艇表面壓力變化，改變潛艇阻力。圍殼處速度等值區(qū)域減小，尾舵上方的速度等值區(qū)域增大。圍殼舵影響潛艇首部和尾部的流速分布，首尾部速度等值區(qū)域變小。布置圍殼舵可在一定程度上增大潛艇的粘壓阻力。

5 結(jié) 語

本文通過對帶不同附體的SUBOFF 潛艇模型進行直航數(shù)值仿真，分析附體對潛艇阻力性能的影響，得到以下結(jié)論：

1）通過比較不同SUBOFF 潛艇模型的摩擦阻力系數(shù)、粘壓阻力系數(shù)和總阻力系數(shù)，發(fā)現(xiàn)附體的增減對各SUBOFF 潛艇模型的摩擦阻力系數(shù)影響較小，對粘壓阻力系數(shù)影響較大，按影響由大到小排列分別為：尾舵、指揮臺圍殼、尾方向舵、尾升降舵、圍殼舵。可知，布置各種附體后，摩擦阻力系數(shù)、粘壓阻力系數(shù)和總阻力系數(shù)均有所增大，說明附體的布置影響潛艇阻力性能。

2）通過數(shù)值仿真，對潛艇表面壓力分布進行研究，分析不同附體對表面壓力系數(shù)的影響，發(fā)現(xiàn)在潛艇首部、指揮臺圍殼處和尾舵處，潛艇表面壓力系數(shù)均有突變現(xiàn)象，說明布置附體會使?jié)撏П砻鎵毫Ψ植及l(fā)生變化，從而改變潛艇的阻力。

3）以中縱剖面、橫剖面和圍殼舵處橫剖面為研究對象，研究潛艇周圍流場的速度分布情況，研究發(fā)現(xiàn)在潛艇首部、指揮臺圍殼前部和尾舵前部速度等值線比較密集，壓力變化較大，影響潛艇阻力性能。