消渦整流片對潛艇馬蹄渦的控制及其與輔翼效果的比較

劉志華，熊 鷹

（海軍工程大學船舶與動力學院，武漢 430033）

1 引 言

潛艇是一種具有良好隱身性和強大突襲能力的海戰裝備。自第二次世界大戰以來，隨著新型軍事裝備技術的全面快速發展，現代化海戰對潛艇綜合性能的要求越來越高。潛艇的隱身性能和水動力性能作為保障潛艇完成其作戰使命的基本性能，日趨受到各海軍強國的重視。進入二十一世紀以來，盡管探潛技術日趨向多樣化發展，但潛艇噪聲仍然是目前潛艇探測所依靠的主要信號特征，繼續降低潛艇噪聲仍然是高技術潛艇的重點發展方向之一。螺旋槳噪聲是潛艇的三大噪聲源之一。螺旋槳工作于潛艇的尾流區中，潛艇尾流品質與螺旋槳的噪聲有直接的關系。由于潛艇尾流高度的周向不均勻性，使得該處工作的螺旋槳槳葉剖面在旋轉過程中受到的來流攻角和速率劇烈變化，導致葉剖面上產生的推力和扭矩也隨之發生變化，從而產生了螺旋槳的轉動噪聲和艇尾殼體激勵振動噪聲，同時進流條件的高度不均勻性還會對螺旋槳空泡噪聲特性有明顯的影響，因此減小潛艇尾流的周向不均勻性是降低潛艇螺旋槳噪聲的重要途徑。

潛艇尾流高度的周向不均勻性是由于其幾何外形特點決定的。典型的潛艇都是由一主艇體與指揮臺圍殼、尾翼等附體組成，在主艇體與指揮臺圍殼、尾翼等附體的接合部，由于形狀的突變，潛艇表面湍流在生成與發展的過程中受到干擾，形成特殊的馬蹄渦結構，馬蹄渦強度高，耗散弱，當馬蹄渦傳播至螺旋槳盤面處時，與主艇體尾流、附體尾流發生相互作用，使得潛艇尾流成為以湍流脈動、粘性效應和漩渦運動為特征的復雜流場區域，導致潛艇尾流嚴重的不均勻性。國際上在開展潛艇尾流形態特征研究的同時，也對潛艇尾流的控制方法進行了大規模的研究，并得到了一些潛艇尾流控制的具體技術措施。總的指導思想是在不改變潛艇艇體和附體形狀的條件下，通過設置附加裝置來改善潛艇尾流形態，使得潛艇槳盤面處流場的均勻性提高，從而降低螺旋槳噪聲。

俄羅斯海軍曾經在潛艇上使用過一種特殊的整流裝置，該裝置長為2m，高為0.5m，整個裝置為流線型設計，像蓋子一樣安裝在艇體尾部離螺旋槳20%～30%的艇長處，有分析認為這種裝置具有抑制潛艇槳盤面流場不均勻性的功能[1]。

李新汶[2]、張楠[3]基于主艇體與附體結合部設置輔翼的整流方法，采用數值模擬進一步研究了不同的輔翼形狀對尾流的整流效果；張楠[4]還對指揮臺圍殼、尾翼前緣根部區域內的渦形態進行了顯示。

本文分析了潛艇主艇體與附體結合部馬蹄渦的特性，闡述了可抑制潛艇馬蹄渦的新型消渦整流方法的作用原理；運用分離渦（DES）計算方法，對全附體潛艇流場進行了數值模擬計算，比較了消渦整流片與輔翼對潛艇馬蹄渦的控制效果，展現了消渦整流片良好的降低潛艇尾流周向不均勻性、提高潛艇螺旋槳入流品質的功能，有利于控制潛艇螺旋槳的噪聲。

2 馬蹄渦的形態及消渦整流片原理

主附體結合部馬蹄渦的產生是由于附體迎流面前方區域產生了較強的逆壓力梯度，上游流體在逆壓力梯度作用下形成橫向渦流，橫向渦流受上游來流的沖擊要向下游傳播，但同時又受到附體的阻礙，所以橫向渦流的形狀在附體首部發生彎曲，其旋轉方向向縱向偏轉，產生縱向渦ωx。在這種情況下，又因為附體兩側的形狀類似于機翼的葉背，具有加速流體流動的作用，使縱向渦進一步拉伸，從而形成了馬蹄渦，馬蹄渦形成后促使結合部發生三維邊界層分離，而三維邊界層分離又進一步增強馬蹄渦強度。圖1為潛艇指揮臺圍殼與主艇體結合部馬蹄渦形態。潛艇尾翼與艇體的連接也是一種主附體連接形式，同樣有馬蹄渦產生，指揮臺圍殼處的馬蹄渦與尾翼處的馬蹄渦向下游傳播，對潛艇尾流有明顯的影響。

根據主附體結合部馬蹄渦的特點，本文設計了一種新型的消渦整流片，該消渦整流片安裝在馬蹄渦出生區域，利用附體兩側繞流產生的橫向速度分量以產生與馬蹄渦旋轉方向相反的“附渦”，這兩種漩渦在向下游的發展過程中就會相互削弱，從而實現對馬蹄渦的主動控制（如圖2所示）。在文獻[5]中，作者分別進行了該消渦整流片效果的試驗驗證和數值模擬，試驗和數值模擬所用的模型都為潛艇主艇體帶指揮臺圍殼的簡易模型，沒有包含潛艇尾翼，且所進行的研究為針對這種簡易模型在小尺度（模型總長1.0m）、低雷諾數（1.55×106）時的原理性研究。在此基礎上，還需以其實用化為目標，開展全附體、高雷諾數情況下的效果研究。

圖1 指揮臺圍殼結合部馬蹄渦的形成Fig.1 Horseshoe vortex at sail-body junction

圖2 消渦整流片產生的附渦（視圖從尾向首）Fig.2 The generation of attached vortex by vortex control baffler(viewed from afterbody to bow)

3 DES方法對潛艇馬蹄渦流場的模擬

3.1 DES方法

潛艇流場數值模擬是 CFD 方法的重點應用方向。 Bong Rhee（2003）[6]，吳寶山（2005）[7]，洪方文、常煜[8]等人分別采用RANS方法，結合不同的湍流模型對類似潛艇主艇體的簡單回轉體和帶附體潛艇流場進行了數值模擬計算，研究了帶附體潛艇的表面壓力分布、局部摩擦阻力系數分布以及尾流伴流分布。黃振宇[9]基于對大渦模擬方法優勢的分析，開展了大渦模擬在潛艇粘性流場計算中的初步應用研究。大渦模擬方法將湍流流場中湍流脈動分為強烈依賴于流動邊界、且對湍流脈動的動量和能量輸運過程起支配作用的大渦和主要起湍流耗散作用的小渦。在計算中只需要求解大渦的運動而將小渦對大渦的影響通過亞格子雷諾應力來表達。大渦模擬方法具有直接求解大于網格尺度的湍流脈動的能力，使其在湍流流場計算方面比求解雷諾平均N-S方程有明顯的優勢，但消耗的計算資源要大得多。

為了保證足夠的計算精度并且降低計算資源的消耗，人們提出了一種綜合應用RANS方法和LES方法的計算方法--分離渦方法（DES）。分離渦方法在壁面邊界層中應用RANS方法，而在邊界層以外流域應用LES方法，這樣較好地發揮了這兩種方法的優點。

進行模糊隸屬度矩陣，得到各分區一級評判結果。為了定量反映各級別隸屬度對目標的影響程度，對評價標準 V={Ⅴ，Ⅳ，Ⅲ，Ⅱ，Ⅰ}在 0～l之間進行離散。令 a=(0.9，0.7，0.5，0.3，0.1)，級別越高對可持續發展的貢獻就越高，評分值也就越高，表明可持續發展能力越大。綜合評價取上述計算所得到的aj和bj值進行計算。

潛艇流場中的馬蹄渦是一種強度高、持續范圍廣和對流場影響明顯的漩渦，采用DES方法可以對這種漩渦運動進行較為準確的模擬。

DES方法是LES方法和RANS方法的組合，即RANS在湍流邊界層內起作用，LES在外域流場起作用，這樣就需要一種轉移函數來控制湍流流動方程在兩種方法之間的轉換。

Δ是計算網格單元在x、y、z方向上最大的距離。

Eric[10]和Michael[11]對運用DES方法進行流場模擬的細節進行了詳細的描述。

各流動控制方程和湍流方程采用有限體積法進行數值離散。方程中的對流項采用二階迎風差分格式離散；擴散項采用中心差分格式離散；時間項采用二階隱式差分格式離散；壓力速度耦合方程采用PISO法求解。數值離散后的代數方程組用Gauss-Seidel迭代法求解，并以多重網格技術加速迭代收斂。

圖3 全附體SUBOFF模型流場計算網格Fig.3 Computational grid for full appended SUBOFF model

3.2 全附體SUBOFF計算模型與計算網格

SUBOFF模型為美國大衛·泰勒研究中心（DTRC）為了獲得有效的潛艇流場模擬仿真方法而建造的一條潛艇試驗模型，由于該模型進行了大量的流場與阻力測量試驗，所以成為國際上開展潛艇流場計算研究的標準模型。SUBOFF模型總長4.356m，其艏部長為1.016m，平行中體長為2.229m，艉部長為1.111m，中部最大直徑為0.508m，指揮臺圍殼高0.206m，長0.368m，最大厚度0.066m。

為了與試驗結果進行比較，本文在數值模擬計算中，以模型總長為特征長度的雷諾數設置為1.2×107。網格總數為9 130 425；時間步長為2.5×10-4s；為了縮短計算時間，首先運用RANS方程進行定常流場的計算，達到穩定后，將計算結果作為DES方法非定常計算t=0s時刻的初始值。

3.3 計算結果與分析

圖4給出了雷諾數為1.2×107時SUBOFF模型距艇艏部端點x/L=0.978的槳盤面處無因次軸向速度ux/U0的分布。圖5給出了r/R=0.25處無因次軸向速度隨周向角變化曲線的計算結果與風洞試驗結果的比較。比較結果表明DES數值模擬計算結果與試驗吻合良好。

圖4 槳盤面無因次軸向速度等值線的分布（計算結果）Fig.4 Iso-contours of non-dimensional axial velocity at propeller disc(calculated results)

圖5 軸向速度計算結果與試驗值的比較（r/R=0.25）Fig.5 Comparison on the axial velocity between the calculated results and tested results(r/R=0.25)

從圖4、5可以看出，在指揮臺圍殼、尾翼與主艇體結合部對應的區域，速度等值線有明顯的隆起，這反映了馬蹄渦在該處對流場的影響。即馬蹄渦將外半徑速度較高的流體帶入到內半徑，同時將內半徑速度較低的流體帶入到外半徑，使得速度等值線發生彎曲，造成了槳盤面流場明顯的周向不均勻性。

4 消渦整流片對馬蹄渦的控制及其與輔翼效果的比較

從圖2可以看出，要產生有利的附渦，消渦整流片需要安裝在附體與馬蹄渦渦核位置之間，且能充分吸收附體首部繞流中的橫向速度分量。根據指揮臺圍殼、尾翼的形狀與尺寸，分別設計了相應的消渦整流片，SUBOFF模型指揮臺圍殼處與尾翼處的消渦整流片尺寸與安裝位置如圖6、7所示，設置消渦整流片后的SUBOFF模型如圖8所示。

圖6 指揮臺圍殼處消渦整流片Fig.6 The vortex control bafflers set at the sail-body junction

圖7 尾翼處消渦整流片Fig.7 The vortex control bafflers set at the stern foil-afterbody junctions

圖9 圍殼處輔翼Fig.9 The fillet of the sail

圖10 尾翼處輔翼Fig.10 The fillets of the stern foil

同時，為了比較消渦整流片與輔翼的效果，開展了帶輔翼的SUBOFF流場計算，文獻[2]結果表明，輔翼的前伸長度較長時能較好地提高螺旋槳盤面流場均勻性，本文的計算中，指揮臺圍殼處的輔翼前伸長度為圍殼最大厚度的3.5倍，尾翼處的輔翼前伸長度為尾翼最大厚度的3.4倍，表面網格如圖9、10所示。

運用DES方法，進行了雷諾數1.2×107時各狀態全附體SUBOFF模型的流場。圖11、12給出了3種狀態下SUBOFF槳盤面r/R=0.25、r/R=0.3、r/R=0.4處流體的軸向和周向速度分量隨周向角的變化曲線，由于潛艇螺旋槳的直徑一般都不會超過主艇體最大直徑的40%，故以上3種半徑涵蓋了對潛艇螺旋槳性能有意義的半徑范圍。

圖11 軸向速度變化曲線的比較Fig.11 Comparison on the axial velocity curves

圖12 周向速度變化曲線的比較Fig.12 Comparison on the tangential velocity curves

從圖11、12可以看出，消渦整流片具有優良的減弱馬蹄渦對尾部流場影響、并提高槳盤面周向均勻性的效果。主附體結合部輔翼的設置也可控制馬蹄渦的強度，對槳盤面流場具有整流作用，這與李新汶[2]、張楠[3]等人的研究結論相符合。兩種方法效果的比較表明，消渦整流片可達到更好的改善潛艇螺旋槳盤面處流場品質的作用。

為了定量描述消渦整流片的效果，根據潛艇的尾流場特點，提出了流場速度分量的周向不均勻度系數的表達式為：uiΔ=（uimax-uimin）/U0， 式中，uimax為流體速度分量峰值，uimin為流體速度分量谷值，U0為無擾動時來流速度，i分別代表x、θ。各狀態全附體SUBOFF模型槳盤面流場的周向不均勻度的比較如表1所示。

表1 潛艇槳盤面速度分量不均勻度的比較Tab.1 Comparison on the circumferential non-uniformity of velocity components at propeller disc

表1說明消渦整流片的應用使得全附體SUBOFF模型槳盤面流場的周向不均勻性大幅度降低，且在對螺旋槳有意義的半徑范圍內，半徑越大效果越好。這顯示了消渦整流片明顯提高了潛艇螺旋槳的進流品質，會起到減小螺旋槳槳葉的非定常力變化范圍、改善螺旋槳振動與噪聲性能的作用。

5 雷諾數的變化對消渦整流片效果的影響

對于實尺度潛艇而言，在水下航行時，以艇總長為特征長度的雷諾數可達到108至109量級，本文關于運用消渦整流片技術控制潛艇主附體接合部馬蹄渦、改善槳盤面流場品質的研究中，雷諾數為107量級，因此超高雷諾數下整流片效果的研究是十分有必要的。并且實際情況下，潛艇經常需要改變航行速度，其雷諾數是不斷變化的。基于上述原因，運用數值模擬方法進行了增大雷諾數對潛艇整流片效果的影響研究。

根據DES方法對計算網格的要求和計算條件的限制，以全附體SUBOFF模型為對象，進行了來流速度8.304m/s、雷諾數3.6×107時的建模和計算，計算網格達到26 769 358，計算時間步長為1.0×10-4s，計算在HP工作站上進行。

圖13 雷諾數對uxΔ降幅的影響Fig.13 Influence of Reynolds number on uxΔ decline

圖14 雷諾數對uθΔ降幅的影響Fig.14 Influence of Reynolds number on uθΔ decline

整流片的尺寸和安裝位置參數與雷諾數1.2×107時完全相同。

將雷諾數為3.6×107時、設置整流片前后槳盤面流場的軸向、周向速度分量的周向不均勻性進行對比，并按照第4節的公式計算了軸向和周向速度分量的不均勻度系數。圖13、14為雷諾數3.6×107與雷諾數1.2×107下，SUBOFF模型槳盤面處軸向、周向速度分量的周向不均勻度系數下降幅度的對比。

從圖13、14可以看出，以全附體SUBOFF主艇體長度為特征長度的雷諾數從1.2×107變化到3.6×107，整流片作用下槳盤面流場速度分量的周向不均勻性下降幅度變化很小，說明了新型整流片的效果相對雷諾數的穩定性，整流片在實艇雷諾數時也可望獲得較好的整流效果。

6 結 論

本文為了抑制全附體潛艇槳盤面處流體速度分量的周向不均勻性，設計了可削弱主附體結合部馬蹄渦特征的新型消渦整流裝置，以數值模擬為基本方法，對整流片的效果進行了模擬計算，并與國內外廣泛研究的輔翼進行了對比研究，結果顯示了消渦整流片具有一定的優勢，可望成為降低潛艇螺旋槳振動與噪聲的一種新方法。

為了推進該方法的工程應用，還需結合實艇幾何特征，開展消渦整流片的測量試驗研究和設計研究，為工程應用提供有意義的技術參考。

[1]蔡 明.俄羅斯現代潛艇設計管窺[J].艦船科學技術,2002,24(3):18-21.

[2]李新汶,陳 源,王文琦.CFD在潛艇艉附體與艇體連接形式研究中的應用[J].船舶力學,2003,7(5):28-32.

[3]張 楠,沈泓萃,姚惠之.潛艇阻力與流場的數值模擬與驗證及艇型的數值優化研究[J].船舶力學,2005,9(1):1-13.

[4]張 楠,沈泓萃,姚惠之.用雷諾應力模型預報不同雷諾數下的潛艇繞流[J].船舶力學,2009,13(5):688-696.

[5]Liu Zhihua,Xiong Ying,Wang Zhanzhi,et al.Numerical simulation and experimental study on the new method of horseshoe vortex control[J].Journal of Hydrodynamics,2010,22(4):572-581.

[6]Bong Rhee,Chao Ho Sung,In Yong Koh.Validation of the flow around an appended SUBOFF body using parallelization and a new wall method[C]//The 8th International Conference on Numerical Ship Hydrodynamics.Busan,Korea,2003:223-232.

[7]Wu Baoshan,Xing Fu,Kuang Xiaofeng,Miao Quanming.Investigation of hydrodynamic characteristics of submarine moving close to the sea bottom with CFD methods[J].Journal of Ship Mechanics,2005,9(3):19-28.

[8]洪方文,常 煜.均勻流中潛艇水下運動表面尾跡的數值模擬[J].船舶力學,2005,9(4):9-17.

[9]黃振宇,繆國平.大渦模擬在水下航行體黏性流場計算中的應用[J].水動力學研究與進展,2006,21(2):190-197.

[10]Paterson E G,Peltier L J,Leonard J P.Detached-eddy simulation of high-Reynolds-number beveled-trailing-Edge boundary layers and wakes[J].Journal of Fluids Engineering,2005,127(3):897-906.

[11]Kinzel M P,Lindau J W,Peltier L J,Kunz R F.Detached-eddy simulations for cavitating flows[C]//18th AIAA Computational Fluid Dynamics Conference.Miami,USA,2007.