潛艇艏舵繞流場的數值模擬

劉明靜 馬運義 吳 軍

中國艦船研究設計中心，湖北武漢430064

潛艇艏舵繞流場的數值模擬

劉明靜 馬運義 吳 軍

中國艦船研究設計中心，湖北武漢430064

為了研究艏水平舵位置對于艉部流場的影響，將艏水平舵布置在指揮室或主艇體上（分別稱為圍殼舵和艏舵）。通過對SUBOFF進行數值模擬計算，驗證網格劃分方式和數值計算方法的正確性，并將該方法應用于圍殼舵和艏舵模型。計算結果表明，在采用較好的流線型指揮室的情況下，采用艏舵形式有利于艉部伴流場的均勻性。

艏舵；圍殼舵；數值模擬；尾部伴流場

1 引言

隨著世界形勢的變化，人們對于海洋權益越來越重視，潛艇也隨之受到人們的重視。對于潛艇周圍粘性繞流場的研究也日益增多。在阻力方面，為了避免艇體與附體的相互干擾，需要充分考慮附體的布置。

潛艇的附體一般都是由3對舵組成，除尾鰭之外，還使用一對水平舵，這對舵布置在指揮室圍殼上（圍殼舵）或艏部（艏舵）。裝設艏水平舵的原因是能夠獨立控制縱傾和下潛深度。對于無縱傾的下潛和上浮的機動情況，單靠一對艉水平舵是無法保證操作的完成。此外，在低于逆速（1.5～2.0 kn）的較低航速上，垂直面內只靠艉舵操縱是很困難的。與高于逆速的速度相比，在逆速和低于逆速的速度上，轉舵的規律改變了，因此有必要布設兩對水平舵。

將艏舵布置在主艇體上時，為了避免潛艇橫靠碼頭時造成舵的損壞，可以設計成可收折式，以較大的展弦比布置在艏部載重水線以上的上層建筑空間內。航速較高時，為避免增大阻力、消除繞流艏舵時的水動力噪聲，也可將艏舵收回。由于舵體距離主艇體水動力中心較遠，故舵效較高，只需較小的舵角就可以產生較大的垂向操縱效果，特別是在發生尾舵卡舵的緊急情況，艏水平舵的這種優點更為明顯。但是艏舵也有明顯的缺陷，由于艏部有效空間較小，可收折式的艏舵在布置上與水聲設備容易產生矛盾。受舵機和水動力噪聲的影響，對水聲器材的使用也會產生影響。

相對于艏舵，圍殼舵一般不會超出潛艇的最大寬度，可以避免橫靠碼頭時的碰撞問題。同時由于圍殼舵遠離艏端，可以減少操縱機構產生的噪聲對聲吶的不利影響。不過由于圍殼舵的位置較高，當潛艇處于接近水面的深度上航行時，與裝設在艇體上的艏水平舵相比，圍殼舵受到海面波浪的影響相對更大一些，從而增加了保持潛望深度的難度。當潛艇從水面快速下潛時，由于圍殼舵位置高出水面很多，無法快速發揮其相應的作用。同時也由于圍殼舵靠近水動力中心，所以操舵效果不如艏舵明顯。

縱觀各國潛艇，俄羅斯的潛艇采用艏舵較多，而歐美等國則較多的采用圍殼舵。這兩種艏舵布置形式究竟哪種較好，需要根據具體的情況來分析，這里主要關注它們對于尾部伴流場的影響。近年來有很多學者應用CFD技術模擬潛艇周邊粘性繞流場［1］，但這些模擬大都是以SUBOFF為模型，沒有考慮艏水平舵對于尾部伴流場的影響。本文采用CFD方法，通過Fluent軟件模擬艏水平舵布置在主艇體和指揮臺上的兩種形式，研究這兩種水平舵形式對于尾部伴流場的影響。

2 計算模型

本文中所采用的計算模型主尺度取為SUBOFF標準模型數值［2］，即L×B×D=4.356×5.08× 5.08，其中前體Lf長1.016 m，平行中體Lp長2.229 m，后體La長1.111 m。

對于艏水平舵布置在指揮臺上的模型 （以下簡稱圍殼舵模型），指揮臺和尾翼布置及形式同SUBOFF，指揮臺長0.368 m，上部有一外凸的頂蓋。尾翼布置方式為翼型后緣，位于x=4.007 m處，4個尾翼剖面為NACA0020翼型，對稱布置。此外，在指揮臺垂向半高處加上圍殼舵，翼型剖面為NACA翼型，但厚度比略減。模型如圖1所示。

對于艏水平舵布置主艇體上的模型 （以下簡稱艏舵模型），主艇體及艉附體與SUBOFF相同。指揮臺為改型后的流線型指揮臺。艏舵布置在艏部上半部，翼型剖面為NACA翼型，厚度比與上述圍殼舵相同，模型如圖2所示。

圖1 圍殼舵模型

圖2 艏舵模型

3 潛艇三維流場計算

3.1 數值計算方法

對不可壓縮粘性湍流流動，采用RANS方程，同時采用RNG k-ε湍流模型進行封閉，其控制方程［3］中，不可壓縮流體連續性方程為：

RANS方程為：

湍動能方程為：

湍動能耗散率方程為：

采用有限體積法離散控制方程以及湍流模式［4-5］。壓力方程采用離散方程進行離散。對于動量方程、湍流方程和雷諾應力方程均采用二階迎風格式進行離散，壓力速度耦合迭代采用SIMPLEC算法［6-7］。

3.2 邊界條件及計算網格

計算過程中計算區域為與潛艇同軸的圓柱體，由入流邊界、出流邊界、壁面邊界及周向控制面組成［8］，如圖3所示。為保證Re≥1.2×107，故速度范圍為ν≥2.765 m/s。對邊界條件的說明如下：

1）入流邊界：位于艇前1倍艇長處，采用速度進口邊界條件，來流速度取為5 m/s；

2）出流邊界：位于艇后2倍艇長處，采用壓力出口邊界條件；

3）壁面邊界：采用無滑移壁面邊界條件；

4）周向控制面：直徑取為艇體直徑的10倍，采用速度進口邊界條件，來流速度與進口相同。

整個計算域采用結構化網格，在艇體周圍生成O型網格，并進行加密處理。網格總數為150萬左右。模型壁面網格如圖4所示。

圖4 模型壁面網格

4 計算結果

4.1 計算方法驗證

為驗證本文中網格劃分方式及計算方法的正確性，首先采用這種方法計算SUBOFF標準模型，并將該計算結果與試驗結果進行對比［9］。圖5所示為本文計算的槳盤面處伴流分數等高線圖，圖6所示為槳盤面處伴流分數實驗測量值等高線。由圖可見，兩者的趨勢是大體一致的。圖7所示為槳盤面上r=0.5R處尾流速度測量值與計算結果的比較，其中橫坐標為槳盤面處逆時針方向旋轉一周的角度值，縱坐標為x向的相對速度值。由圖可見，計算結果比較連續，均位于實驗測量值范圍內。這說明采用的網格劃分方式和數值計算方法能夠較好地模擬潛艇周圍的粘性流場。

4.2 計算結果

圖5 槳盤面處伴流分數等高線計算圖

圖6 實驗測量槳盤面處伴流分數等高線圖

圖7 槳盤面上0.5R處x向尾流速度測量值與計算值對比圖

采用上述方法，對兩種模型進行計算，所得的槳盤面處伴流分數等高線圖如圖8、圖9所示。由圖可見，相對于無艏舵的SUBOFF模型，受艏舵和圍殼舵影響，尾部伴流場發生了明顯變化。對于圍殼舵，在主艇體上部出現了明顯的圍殼舵影響區域。由于艏舵位于主艇體上，艏舵的影響區域與指揮室的影響區域混在一起。為了顯示出艏舵對于艉部伴流場的影響，圖10為沿艇長方向不同位置處伴流分數等高線圖。由該圖可見艏舵對于流場的影響隨著向艉部的發展，逐漸減小。由于本文艏舵模型中指揮室較寬，且與主艇體結合部位過渡不夠光順，以致在指揮室后方形成了較為明顯的渦。該漩渦在槳盤面處對伴流場的影響已明顯蓋過了艏舵，這也從側面顯示出艏舵經過沿艇向的發展，對艉部流場均勻性影響很小。

圖11為槳盤面上r=0.5R處尾流速度計算結果比較圖。圖中0°坐標取為垂向坐標正向，受圍殼舵和艏舵影響，在此處伴流場中有個速度的峰值［10］，相對于整個周向，在0°、360°處峰值最高。圖中實驗數據為全附體潛艇模型試驗結果。由該圖可見本文的計算結果大體趨勢是正確的。艏舵模型在槳盤面處伴流場的均勻度較圍殼舵稍優。但由于指揮室圍殼過寬，與主艇體結合部位不夠均勻，在指揮室后方形成了漩渦，使得上半部的速度等值線出現鼓起，削弱了艏舵的優勢。

圖8 圍殼舵模型槳盤面伴流分數等高線圖

圖9 艏舵模型槳盤面處伴流分數等高線圖

圖10 艏舵模型橫剖面速度等值線圖

圖11 槳盤面0.5R處x向尾流速度

5 結束語

本文采用結構化網格對于SUBOFF模型進行了數值模擬計算，通過將計算結果與試驗結果對比，驗證了該計算方法的可行性。并對帶有指揮臺圍殼舵和艏舵的潛艇模型進行了數值模擬，得出了尾部槳盤面處伴流分數等高線圖的大體趨勢。通過對于槳盤面處伴流分數等高線圖的分析，較低矮的指揮室和距離尾部槳盤面距離較遠的艏部水平舵對于尾部伴流場的影響較小。

但是指揮室的線型，尤其是指揮室的后半部的線型對于艉部伴流場的均勻性影響較大。為了減小阻力，采用艏部水平舵的時候必須采用線型優良的指揮室，這樣才能減少次生流的強度，從而削弱指揮室后方渦的影響。這些有待于后續開展進一步的計算研究。

［1］ KARLSSON A，FUREBY C.LES of the flow past a 6：1prolate spheroid ［C］//47th AIAA Aerospace Sciences Meeting including the new Horizons Forum and Aerospace Exposition，Florida，USA 2009.

［2］ AMIT T，DEBABRATA S.Calculation of transverse hydrodynamic coefficients using computational fluid dynamic approach［J］.Ocean Engineering，2006（33）：789-809.

［3］ QIU L Y，SHI Z K.Validation of numerical simulation of the flow over submarine geometries with full appendages［J］.Journal of ship Mechanics，2007（6）：341-351.

［4］ 馬崢，黃少峰，朱德祥.湍流模型在船舶計算流體力學中的適用性研究［J］.水動力學研究與進展，2009（3）：207-217.

［5］ 邱遼原.潛艇粘性流場的數值模擬及其阻力預報的方法研究［D］.武漢：華中科技大學，2006.

［6］ CHI Y.Prediction of flows over an axisymmetric body with appendages［C］//The 8th International Conference on Numerical Ship Hydrodynamics.Korea.2003

［7］ FUREBY C.Towards the use of large eddy simulation in engineering［C］//46th AIAA Aerospace Sciences Meeting and Exhibit.Nevada，2008.

［8］ 張懷新，潘雨村.CFD在潛艇外形方案比較中的應用［J］.船舶力學，2006，10（5）：1-9.

［9］ 張楠，沈泓萃，姚惠之.潛艇阻力與流場的數值模擬與驗證及艇型的數值優化研究［J］.船舶力學，2005，7（2）：1-12.

［10］ JEANS T L.High-resolution reynolds-averaged Navier-Stokes flow predictions over axisymmetric bodies with tapered tails［J］.AIAA JOURNAL，2009（1）：19-35.

Numerical Simulation on the Flow Around Fore Hydroplane of Submarine

Liu Ming-jing Ma Yun-yiWu Jun
China Ship Development and Design Center，Wuhan 430064，China

In order to investigate the impacts of the plane's position on the wake flow around the bow of submarine，two models for both the fore hydroplane and sailplane were established based on SUBOFF software，one positioned in the main hull and the other in the sail.Through a numerical simulation，SUBOFF model was computed to validate the method of grid plotting and numerical computation，and then applied this approach to the sailplane and fore hydroplane calculation respectively.The results show that the plane positioned in the main hull with streamlined sail will benefit the uniformity of wake flow.

fore hydroplane；sailplane；numerical simulation；wake flow

U661.2

A

1673-3185（2010）05-40-04

10.3969/j.issn.1673-3185.2010.05.008

2009-08-11

“十一五”國防科技預研課題（1010501××××）

劉明靜（1981-），女，博士研究生。研究方向：船舶與海洋結構物設計制造。E-mail：alice-lmj@163.com

馬運義（1942-），男，研究員，博士生導師。研究方向：船舶與海洋結構物設計制造