一種水下航行體附加質量數值計算方法

周景軍, 李育英, 項慶睿

?

一種水下航行體附加質量數值計算方法

周景軍, 李育英, 項慶睿

(中國船舶重工集團公司 第705研究所, 陜西 西安, 710075)

基于相對運動思想, 通過求解3D N-S方程, 給出了一種水下航行體附加質量高精度、高效的計算方法,可方便給出任意幾何外形航行體的附加質量。該方法中, 水下航行體的附加慣性力(力矩)通過非定常運動與定常運動流體動力之差獲得, 進而獲得附加質量。通過在動量方程添加源項求解航行體非定常運動流場, 一方面保證能夠模擬真實流場結構, 另一方面避免采用動網格過程中由于復雜網格重生帶來的困難, 同時由于只采用了一套網格, 計算精度得到進一步保證。本文計算方法的精度通過橢球精確解得到有效驗證, 具有重要的工程應用價值。

水下航行體; 附加質量; 數值計算方法; 相對運動

0 引言

物體在靜止的流體中開始運動時, 必然會推動周圍的流體質點, 使其克服慣性后也開始運動, 物體本身同時又受到這些流體質點的反作用力。流體中的物體在改變運動狀態時, 為了克服流體慣性而受到的作用力稱為慣性阻力。附加質量就是物體在流場中運動時流場慣性的一種度量。水下航行體的附加質量嚴重影響其運動特性, 是航行體操穩性分析、控制規律設計以及彈道預報的重要參數。目前, 計算附加質量的主要方法為基于勢流理論的數值分析方法[1-2], 其關鍵是確定物體的速度勢, 至今只有簡單的球、橢球等形狀物體能夠獲得理論精確解。工程上計算魚雷等細長回轉體的附加質量常用二元切片理論, 將雷體的附加質量看作是各橫截剖面附加質量疊加之和, 并用修正系數考慮三元影響[3]。采用試驗方法無疑更具可靠性, 但相比計算而言明顯周期長、費用高, 對于某些復雜航行體而言, 同樣存在著極高的難度。隨著計算機能力的不斷提高, 采用計算流體力學(computational fluid dynamics, CFD)數值仿真方法計算附加質量逐漸成為一種新的趨勢, 國內已有多篇文章相繼發表[4-6], 研究方法均是采用了動網格技術。盡管動網格技術更加接近實際流動, 但是對于具有復雜外形的航行體(例如帶有泵噴、導管對轉槳的魚雷), 采用動網格技術難度較大, 難以保證網格重生的質量, 上述方法適用對象必然受限。

針對上述問題, 本文采用相對運動的思想, 通過施加合適的邊界條件及在動量方程右側添加源項的方法, 發展了一種新的附加質量CFD計算方法, 避免了動網格帶來的一系列問題。通過采用橢球理論解的驗證, 本文的數值方法計算效率高, 精度完全可以滿足工程需求, 對于任意外形水下航行體均適用, 同時還給出了計算過程中一些特殊情況的處理方法。

1 附加質量的求解思路

附加質量的大小一般僅取決于流場中物體的形狀及運動方向, 而與運動速度無關。附加慣性力的影響只有通過非定常運動才能有所表現。具體求解思路如下。

1.1 水下航行體非定常運動流場仿真

水下航行體加速直線運動、轉動等非定常運動過程的數值仿真最直接的方法是采用動網格技術, 對于帶有泵噴推進器、對轉槳的水下航行體, 采用動網格技術由于涉及到網格重生, 網格質量難以保證, 從而嚴重影響求解精度。針對上述問題, 本文提出的求解方法, 網格只需要一次生成, 通過在動量方程中添加源項的方法, 保證了整個流場壓力分布的真實性, 同時由于采用了相同的網格, 可以將網格的影響降到最低。

源項法以航行體縱向加速為例, 添加源項后方向的動量方程如下。

通過添加源項, 可以使流域前后壓力場與實際流動情況一致, 避免產生壓力差。

1.2 附加慣性力(力矩)計算

2 數值方法及計算模型

2.1 數值方法

本文數值方法采用基于求解3D N-S方程的有限體積法, 空間離散采用2階迎風格式, 時間離散采用2階歐拉后差格式, 采用多重網格技術加速收斂。湍流模型采用SST模型。

2.2 研究對象及網格驗證

為了充分說明本文數值方法的有效性, 本文的計算模型選擇了具有精確解的橢球體。如圖1所示。橢球長半軸長40 cm, 短半軸長10 cm。建立了2套網格, 分別為單元數200萬的結構化網格, 用于粘性流計算和單元數60萬的非結構化網格用于無粘流計算。

圖1 橢球體及其坐標定義

結構化網格計算結果如圖2所示。橢球阻力系數理論解層流狀態為0.20, 湍流狀態為0.06。可以看出, 本文基于粘性流的數值方法和計算網格在預測橢球阻力系數時, 具有很高的精度, 下面給出附加質量的計算結果, 計算網格保持不變。

2.3 邊界條件設置

圖2 不同雷諾數下的阻力系數

圖3 加速過程壓力場分布

圖4 旋轉域速度矢量和壓力分布

Fig. 4 Pressure and velocity vector distributions in rotating field

3 附加質量特殊處理方法

圖5 恒定速度下x方向阻力計算結果

另外, 對于形狀復雜的航行體, 采用結構化網格難度很大。由于航行體的附加質量與流體粘性關系不大, 因此, 可以不求解邊界層, 網格可直接采用非結構化網格, 即上述第2套網格。

4 計算結果

表1給出了橢球體附加質量理論解[8], 與表1相比, 表2所示的計算結果采用粘性流和理想流, 2種計算方法均可獲得較好的計算結果, 誤差均在5%以內, 說明附加質量基本和流體粘性無關。采用結構化網格(如圖7所示)和非結構化網格差別不大, 非結構化網格適用于所有幾何外形, 因此大大增加了本文計算方法的適用范圍。

表1 橢球附加質量理論解

表2 附加質量理論解及計算結果

5 工程應用

圖6 水下航行體幾何外形

圖7 計算網格

表3 附加質量計算結果

6 結束語

本文基于相對運動思想, 提出了一種附加質量的CFD快速計算方法, 避免了采用動網格帶來的一系列問題, 由于該方法可直接采用非結構化網格, 因此適用于任意復雜的幾何外形, 大大增大了該方法的適用范圍, 具有一定的工程應用價值。

[1] 劉丹, 王曉亮, 單雪雄. 平流層飛艇的附加質量及其對飛艇運動的影響[J]. 計算機仿真, 2006, 23(6): 52-56. Liu Dan, Wang Xiao-liang, Shan Xue-xiong. Added Mass to Stratospheric Airship and Its Effect on Motion[J]. Computer Simulation, 2006, 23(6): 52-56.

[2] 林超友, 朱軍. 潛艇近海底航行附加質量數值計算[J]. 船舶工程, 2003, 25(1): 26-29. Lin Chao-you, Zhu Jun. Numerical Computation of Added Mass of Submarine Maneuvering with Small Clearance to Sea-bottom[J]. Ship Engineering, 2003, 25(1): 26-29.

[3] 張向強, 姜魯華, 王生, 等. 飛艇動力學仿真[J]. 計算機仿真, 2008, 25(6): 79-82. Zhang Xiang-qiang, Jiang Lu-hua, Wang Sheng, et al. Dynamic Equations of Airship[J]. Computer Simulation, 2008, 25(6): 79-82.

[4] 馬燁, 單雪雄. 數值計算復雜外形物體附加質量的新方法[J]. 計算機仿真, 2007, 24(5): 75-78. Ma Ye, Shan Xue-xiong. A New Numerical Computation Method for Added Masses of Complicated Object[J]. Com- puter Simulation, 2007, 24(5): 75-78.

[5] 傅慧萍, 李杰. 附加質量CFD計算方法研究[J]. 哈爾濱工程大學學報, 2011, 32(2): 148-152. Fu Hui-ping, Li Jie. Numerical Studies of Added Mass Based on the CFD Method[J]. Journal of Harbin Engineering Uni- versity, 2011, 32(2): 148-152.

[6] 弓三偉, 陸宏志, 鄒正平, 等. 彈性體和帶空泡體的附加質量動態數值計算[J]. 計算機仿真, 2010, 27(4): 349-353. Gong San-wei, Lu Hong-zhi, Zou Zheng-ping, et al. Hydro- dynamic Numerical Computation of the Added Mass of Elastic Body & Cavitation Bubble[J]. Computer Simulation, 2010, 27(4): 349-353.

[7] 吳望一. 流體力學[M]. 北京: 北京大學出版社, 2000: 154-166.

[8] 張宇文. 魚雷外形設計[M]. 西安: 西北工業大學出版社, 1998: 713-716.

A Numerical Computation Method of Additional Mass for Underwater Vehicle

ZHOU Jing-jun, LI Yu-ying, XIANG Qing-rui

(The 705 Research Institute, China Shipbuilding Industry Corporation, Xi′an 710075, China )

A high-precision and high-efficiency method of computing additional mass for an underwater vehicle is proposed by solving the 3-D N-S equations based on the relative motion principle. This method is convenient to compute the additional mass for an underwater vehicle with complex shapes. The additional inertia force(moment) is achieved from the difference of hydrodynamics of fluids in unsteady and steady motions, then the additional mass is obtained. By adding the source term to the momentum equation, the unsteady flow field is solved to simulate the real structure of flow field, and the difficulty generated by complex grid rebirth in using the dynamic grid is avoided. The accuracy of the proposed method is also assured because of adopting the same grid. The validity of the proposed method is verified by the exact solution of ellipsoid.

underwater vehicle; additional mass; numerical computation method; relative motion

TJ630.2; TV131

A

1673-1948(2013)04-0246-04

2012-10-18;

2013-03-21.

中國博士后科學基金特別資助(2013T60895).

周景軍(1981-), 男, 博士, 主要研究方向為魚雷總體技術.

(責任編輯: 陳 曦)