運輸型載人潛器的阻力數值預報

李 明，王奎民，譚 浩，劉 峰

（1.海軍駐錦州地區軍事代表室，遼寧錦州 121000；2.哈爾濱工程大學 船舶工程學院，哈爾濱 150000）

運輸型載人潛器的阻力數值預報

李 明1，王奎民1，譚 浩2，劉 峰2

（1.海軍駐錦州地區軍事代表室，遼寧錦州 121000；2.哈爾濱工程大學 船舶工程學院，哈爾濱 150000）

首先對 Myring型回轉體的阻力數值計算方法進行了討論。其中著重對回轉體周圍控制域的網格離散方法及湍流模型的選擇進行了討論。將回轉體數值計算得到的結果與實驗值進行對比，選出適合的湍流模型，同時驗證了計算方法的可行性。其次應用之前提出的數值計算方法，對某型載人潛器的模型進行阻力數值計算，得到相應的阻力系數。最后通過換算得到實艇阻力，并繪制了實艇的有效功率曲線。研究表明，相比于k-ε模型，k-ω模型能更好地計算Myring型回轉體的阻力，并且在高雷諾數下有較高的精度，更適宜目標潛器的阻力預報。從目標潛器模型的直航阻力入手，對實艇的阻力進行預報，最后得到其有效功率，為該載人潛器的動力、控制等系統的設計提供了參考。

載人潛器；回轉體艇型；數值計算；阻力

0 引言

水動力是載人潛器的重要指標，對于其總體性能有著重要的影響。然而，載人潛器的外形復雜、附體眾多，這些均導致其水動力預報復雜且困難。模型實驗法是當前解決水動力實際問題較為有效的方法，但成本極高。隨著計算機技術的發展，數值模擬以其周期短、成本低等優點，逐漸成為研究潛器水動力性能的有效手段。目前，以“運動方程為基礎，拘束模型試驗為輔”[1]，從而獲得水動力系數的方法已成為潛器水動力預報的規范模式，并且得到了大多數研究人員的認可與應用。

李佳[2]研究了載人潛器的阻力數值計算方法，并與模型艇拖曳實驗進行了對比，證明了數值模擬的準確性。龐永杰等[3]提出了新的多面體網格運動求解器，并模擬了橢球體的平面運動。李剛[4]對穿梭艇進行了一系列的數值計算，最終提出了適用的平方項水動力模型。Fuglestad等[5]針對HUGIN 3000原型的阻力進行數值計算，通過計算值與模型試驗數據的比對，驗證了數值計算的有效性和高雷諾數阻力預報的可能性。涂海文等[6]運用雷諾平均 N-S方程，采用RNGk-ε湍流模型，對四種潛艇模型的阻力進行數值模擬。趙金鑫[7]針對某型帶吊艙潛器進行數值計算，通過比對試驗結果驗證了計算方法的有效性，并建立了潛器的操縱運動仿真系統。邱遼原[8]對潛艇粘性繞流場進行了數值模擬，同時檢驗了數值計算結果，創建了基于該數學模型的型線自動生成系統。

本文首先以Myring型回轉體為基礎，針對其數值計算方法進行了研究與討論，并將計算得到的阻力系數與實驗值進行對比，以檢驗數值計算的精度和驗證計算方法的可行性。之后根據回轉體的數值計算方法對一型載人潛器的阻力進行數值模擬，得到潛器的阻力系數。最后，繪制實艇有效功率曲線。

1 計算方法的驗證

1.1 Myring型回轉體

Myring型曲線是常用的AUV艏艉形狀，具體形狀見圖1。

本文對此艇型進行數值模擬，并將計算值與實驗值對比，以證明計算方法的可行性。

Myring型曲線的艏部形狀方程為：

艉部形狀方程為：

式中，a為艏部長度，b為平行中體長度，c為艉部長度，d為中體直徑，x是長軸上點到艏部頂點的距離，r為x點處的半徑，n和θ分別是控制艏艉曲線飽和程度的參數。

本文Myring型回轉體模型參數為，平行中體直徑d=2 800 mm，長度b=724 mm。艏部參數分別為，a=280 mm，n=2；兩組艉部參數分別為，c=784 mm，θ=27°。

1.2 Myring型回轉體模型的數值計算

1.2.1 計算域的建立及網格劃分

由于計算模型左右對稱，為減少網格數量節約計算成本，本次數值計算僅針對半個回轉體，采用笛卡爾直角坐標系，坐標系原點位于回轉體頂部。計算域為半圓柱形，尺寸為：半徑為1.5L，來流域、去流域分別為1.5L，5L（L為回轉體長度）。

回轉體表面網格的離散方式采用混合網格劃分法。在首尾曲率變化較大的區域采用非結構化網格，平行中體采用結構化網格。為了在保證計算精度的同時盡量降低計算成本，將流體域分為靠近回轉體的內域和遠離回轉體的外域。內域采用細密的非結構化網格對回轉體周圍的流體進行捕捉，外域采用粗糙的結構化網格，在不影響精度的同時減少網格數量。由于選擇了合適的網格離散方法和網格尺度，最終取得了較高的計算精度并使用了較少的計算時間。

1.2.2 數值計算參數的設定與選取

求解器采用分離式求解器，基于有限體積法來離散控制方程，求解穩態三維RANS方程，單相流模型。邊界條件：入口條件為速度進口，出口條件為outflow，物體壁面條件選為無滑移條件，入口邊界條件的湍流參數項為湍流強度好湍流粘度比。

湍流強度和湍動粘度比的設置參考相關文獻[1]設為I=0.5%，μr/μ。壓力速度耦合方式選擇SIMPLEC，壓力插值選擇PRESTO。對流項及擴散項的選擇根據文獻[2]“應用二階迎風格式其結果比一階迎風格式更接近試驗值，并且這種結論不隨網格數目的增減而變化”，采用二階迎風格式離散。

湍流模型的選擇需通過對比試驗進行定量分析。在其余設置相同的情況下，分別采用Realizablek-ε湍流模型和SSTk-ω湍流模型對回轉體進行數值模擬。之后比對實驗值選擇適宜的湍流模型。

1.2.3 湍流模型

1）Realizablek-ε模型

標準k-ε模型對時均變率特別大的情形，有可能導致負的正壓力[2]。為使流動符合湍流的物理定律，需要對正應力進行某種數學約束。湍動粘度計算式中的系數Cμ不應該是常數，而應與應變率聯系起來。從而提出了Realizablek-ε模型。在Realizablek-ε模型中，關于k和ε輸運方程如下：

其中，

式中，k為湍動能；ε為湍流耗散率；ρ為密度；t為時間；ui(i=1,2,3)對應x、y、z三個方向；μ為動力粘度；Gk為湍動能k的產生項；σk和σε為Prandtl數，分別與湍動能k和耗散率ε對應；μt為各向異性的張量；C1和C2為經驗值。

2）SSTk-ω方程模型

SSTk-ω模型由Menter發展而來，以便其在廣泛的領域中可以優于k-ε模型，并且在近壁面自由流中k-ω模型有廣泛的應用范圍和精度。為了達到此目的，k-ε模型被修改為k-ω模型。SSTk-ω模型和標準k-ω模型相似，但SSTk-ω模型比標準k-ε模型在廣泛的流動領域有更高的精度和可信度。

湍動能k的輸運方程為：

湍流比耗散率ω的輸運方程為：

其中，湍流粘度為：

Gk表示由于平均速度梯度引起的湍動能的產生項：

Gω表示ω的產生項：

SSTk-ω湍流模型考慮了湍流剪應力的輸運特性，能夠準確地預報由于逆壓梯度導致的流動分離點和分離區域。

1.2.4 計算結果與實驗值對比分析

圖5為實驗值與不同湍流模型所測得的數值對比圖。

由對比圖可看出，當流速較低時k-ε模型與實驗值更為接近，最小誤差甚至達到了0.7%。雖然低速時k-ω模型比k-ε模型精度低，但是最大誤差為6.6%、最小誤差為4.8%，這種誤差工程上可以接受。

但隨著速度的增加，k-ε模型所得到的結果與實驗值的誤差越來越大，最大可達20.03%。而k-ω模型在高速時仍能較好的模擬流場，最小誤差達到2.8%。綜上，k-ω模型更適于模擬這種流動，并且在低、高速均具有較為滿意的精度，所以之后載人潛器的數值計算采用k-ω湍流模型。

2 目標載人潛器阻力預報

2.1 潛器模型的網格劃分及設置

針對目標艇的網格劃分，艇體面網格劃分方法采用非結構化網格；體網格內域采用 Tet/Hybrid，外域采用Hex/Map網格類型。

數值計算的初始設置與 Myring型回轉體的設置相同。在湍流模型的選擇方面，采用之前得到的結論，使用SSTk-ω湍流模型。

2.2 潛器模型阻力預報

載人潛器模型直航阻力數值計算結果見表1。

表1 載人潛器模型直航阻力計算

2.3 潛器實艇有效功率預報

由于本次研究的潛器作業地點遠離水面，所以總阻力只包含粘性阻力。因此模型與實體之間的阻力轉換采用雷諾相似。因為摩擦阻力系數是Re的函數，故相同雷諾數下潛器模型與實體的阻力系數相同。

再根據相同外形潛器的粘壓阻力系數恒定的假定可知，相同雷諾數下潛器模型與實體的阻力系數相同。根據雷諾數相等條件（Rem=Rep）可知：

λ=為縮尺比。（m為模型，p為實艇）當ν相同時，由Re數相似準則，得出

由此推斷出，若要進行拖曳驗證實驗，拖車速度需要達到實艇速度的6倍。這種要求明顯不能實現。考慮到該潛器主體部分為半球體與圓柱體的結合，在水下直航運動時。當Re數超過臨界Re數時，其阻力系數曲線將趨于水平。因此，可通過此法計得到實艇阻力。

根據圖7的數據，得出本潛器模型臨界雷諾數下的阻力系數，Cd=0.020 4。因此目標實艇在臨界雷諾數下的阻力系數

3 結論

本文通過對某運輸型載人潛器的阻力性能進行研究，得到的結論如下：

1）針對Myring型回轉體的數值計算，采用SSTk-ω湍流模型能更好的模擬流場情況，計算值與實驗值更為接近。

2）將Myring型回轉體的計算值與實驗值進行對比。由于計算誤差在工程允許范圍內，本次計算方法的可行性得以驗證。

3）在潛器直航阻力計算方面，通過對模型艇進行數值計算，得到了相應的阻力系數。同時采用“臨界雷諾數法”，得到實艇的阻力曲線及有效功率曲線。

[1]王小波.輕載HOV水動力系數數值預報[D].哈爾濱：哈爾濱工程大學, 2013.

[2]WANG Xiaobo Numerical Prediction for Hydrodynamic Coefficient of Light Load HOV[D].Harbin： Harbin Engineering University, 2013.

[3]李佳.載人潛器阻力性能的數值和試驗預報及外形優化研究[D].哈爾濱： 哈爾濱工程大學, 2009.

[4]LI Jia.Prediction with Numerical Calculation and Model Test for Drag Performance and Form Optimization for Drag Reduction of a Manned Submersible[D].Harbin Engineering University, 2009.

[5]龐永杰, 楊路春, 李宏偉, 等.潛體水動力導數的CFD計算方法研究[J].哈爾濱工程大學學報,2009,30(8)： 903-908.

[6]李剛.穿梭潛器水動力特性的數值模擬和試驗研究[D].哈爾濱： 哈爾濱工程大學, 2011.

[7]Fuglestad A L., Grahl Madsen M.Computational Fluid Dynamics Applied on an Autonomous Underwater Vehicle[C]// Proceedings of the 23rd International Conference on Offshore Mechanics and Arctic Engineering.Vancouver BC Cana-da, 2004, USA：American Society of Mechanical Engineers, 2004.

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

[9]趙金鑫.某潛器水動力性能計算及運動仿真[D].哈爾濱： 哈爾濱工程大學, 2011.

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

[11]楊卓懿, 宋磊.回轉體艇型阻力近似計算方法與試驗研究[J].船舶工程, 2015(6)： 18-21.

[12]劉峰, 韓端鋒, 王小波, 等.載人潛器水動力性能研究[J].科技導報, 2014,32(6)： 29-33.

Numerical Prediction on Resistance of Transportation Manned Submersible

LI Ming1, WANG Kuimin1, TAN Hao2, LIU Feng2
(1.Navy Military Representative Office, Liaoning Jinzhou 121000, China; 2.College of Shipbuilding Engineering,Harbin Engineering University, Harbin 150000, China)

Firstly, the numerical method of Myring type rotary body’s resistance is discussed.The method of grid discrete on control domain surrounding the axisymmetric body and the selection of turbulence model are emphatically analyzed.The calculation results are obtained and compared with the experimental data.Thus, an available turbulence model would be selected.The calculation results are compared with the experimental data.The feasibility of the calculation method is verified.Hereafter, a series of corresponding resistance coefficients can be obtained by calculating the model of a specific type manned submarine’s resistance.At last, the resistance of real submarine is got by conversion.And the effective power curve of the real submarine is drew.The results show that, comparing k-ε model, the drag of Myring type rotary body can be calculated better by k-ω model, and more accurate data can be received under the high Reynolds number.So k-ω model is more plausible for the predication of the aim submersible vehicle’s drag.Starting with the direct navigation resistance of the aim submersible vehicle model, the drag of real marine is estimated.Finally the effective power is obtained and the paper provides references for the design of the power system and control system of manned submarines.

manned submersible; rotatory submersible; numerical calculation; resistance

U661.31

A

10.14141/j.31-1981.2017.06.015

李明（1971—），男，高級工程師，研究方向：艦艇總體設計與管理。