吳思源，張懷新，姚慧嵐

（1.上海交通大學a.海洋工程國家重點實驗室，上海 200240；b.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240；2.中國海洋大學，山東青島 266100）

0 引 言

隨著船舶的大型化和高速化發(fā)展，螺旋槳的負荷和轉(zhuǎn)速相應增加，其空泡特性始終受到眾多研究者的關(guān)注?？张莸漠a(chǎn)生是造成螺旋槳的水動力性能惡化的重要原因，伴隨著空泡的產(chǎn)生，振動與噪聲的增加嚴重影響著船員的舒適性和船舶的隱蔽性。

隨著計算流體力學（Computational Fluid Dynamics，CFD）的不斷發(fā)展，數(shù)值計算方法在螺旋槳空泡現(xiàn)象模擬中得到了廣泛應用。朱志峰等以E779A螺旋槳為研究對象，重點探究了計算域大小和網(wǎng)格尺寸對螺旋槳空化流場預報，以及螺旋槳推力、軸向速度和片空化位置的影響。王順杰等采用動網(wǎng)格模型處理螺旋槳的周期性旋轉(zhuǎn)問題，對DTMB4381 螺旋槳在一定空泡數(shù)下的空泡形態(tài)進行了模擬，但缺乏對螺旋槳水動力性能的定量分析。蒲汲君等對來流速度和模型尺度變化對E799A 螺旋槳片產(chǎn)生空泡的影響進行了研究。劉登成等探究了在均勻流下，網(wǎng)格類型和湍流模型對PPTC（Potsdam Propeller Test Case）螺旋槳梢渦空泡模擬的影響，提出了一種在梢渦處進行網(wǎng)格加密的方法，并分析了梢渦和梢渦空泡的流動特性。馮學梅等通過計算給出了與試驗結(jié)果較吻合的螺旋槳片空泡形態(tài)，重點探究了不同不可冷凝質(zhì)量分數(shù)對空化流場模擬的影響。HELAL等以E779A螺旋槳為例，比較了低雷諾數(shù)下Transition k-kl-ω 轉(zhuǎn)捩模型與全湍流標準k-ε模型在螺旋槳空泡性能模擬方面的差異，前者對空泡位置分布和槳葉表面壓力分布的模擬結(jié)果更準確。SATO等利用CFX軟件對10 個不同商船螺旋槳的空泡形態(tài)和表面壓力進行了模擬，該方法對螺旋槳的片空泡分布和壓力脈動的一階葉頻分量的模擬效果較好。HUUVA 等利用OpenForm 平臺，以PPTC 螺旋槳為研究對象，驗證了數(shù)值模擬方法的準確性，并將計算程序應用于某一商業(yè)方位螺旋槳中，探究了其空泡性能。

本文基于德國波茨坦水池在2011 年組織的螺旋槳空泡研討會上選用的算例螺旋槳（PPTC 螺旋槳），利用StarCCM+軟件探究其在均勻來流下的敞水水動力和空泡性能，建立準確可靠的數(shù)值模擬計算方法，為相關(guān)研究提供參考。

1 數(shù)值方法

1.1 控制方程

連續(xù)性方程是質(zhì)量守恒定律在流體力學中的表現(xiàn)形式，當有空泡產(chǎn)生時，空間中存在著氣、液2 相物質(zhì)，因此基于均質(zhì)混合流的連續(xù)性方程為

式（1）中：ρ為混合密度，即ρ=ρ（1 -α）+ρα，ρ和ρ分別是液相和氣相的密度，α為氣相的體積分數(shù)。

動量方程為

式（2）中：μ為混合動力黏性系數(shù)；μ為湍流引起的混合渦流黏性系數(shù)，需通過求解湍流模型得到。

1.2 空泡模型

空泡模型選擇StarCCM+軟件中內(nèi)置的Schnerr & Sauer模型，該模型是基于質(zhì)量輸運方程得到的多相流模型，輸運方程的表達式為

在Schnerr & Sauer空泡模型中，R和R分別通過以下公式計算。

1）當p ＜p時，有

2）當p ＞p時，有

2 研究對象

本文選定的研究對象PPTC螺旋槳是VP1304 的可調(diào)螺距螺旋槳，在2011 年德國召開的船舶推進器國際研討會上，對該槳展開了豐富的試驗和數(shù)值計算研究，并公布了用以驗證數(shù)值計算方法的相關(guān)試驗結(jié)果。PPTC螺旋槳的主要參數(shù)見表1。

表1 PPTC螺旋槳的主要參數(shù)

對于均勻流場下螺旋槳的敞水和空化模擬，首先建立相應的計算域，整個計算域由旋轉(zhuǎn)域和靜止域2 部分組成。計算域流場示意圖見圖1。

圖1 計算域流場示意圖

為準確捕捉圍繞螺旋槳旋轉(zhuǎn)的流場特征，并準確監(jiān)測螺旋槳在工作時產(chǎn)生的推力和扭矩值，同時考慮到計算成本，采用切割體網(wǎng)格對計算域進行網(wǎng)格劃分，設(shè)置網(wǎng)格尺寸，保證槳葉壁面與第一層網(wǎng)格的Y +值大于30。試驗發(fā)現(xiàn)，螺旋槳除了在槳葉背產(chǎn)生明顯的片空泡，還在葉梢處產(chǎn)生脫離槳葉表面的梢渦空泡。為對梢渦空泡的范圍和位置分布進行準確模擬，需要更為精細的網(wǎng)格，因此進一步對葉梢部分網(wǎng)格進行加密。網(wǎng)格劃分和加密情況見圖2。

圖2 網(wǎng)格劃分和加密情況

3 無空化條件下的水動力性能計算

螺旋槳敞水水動力性能數(shù)值模擬的計算條件與試驗工況保持相同，對應的進速系數(shù)J =0.6，0.8，1.0，1.2，1.4，水的密度ρ=998.67 kg/m，運動黏性系數(shù)ν=1.070 ×10m/s，轉(zhuǎn)速n=15 r/s。

3.1 網(wǎng)格無關(guān)性驗證

為驗證數(shù)值計算結(jié)果與網(wǎng)格大小的關(guān)系，選取網(wǎng)格疏密程度不同的3 套網(wǎng)格，網(wǎng)格數(shù)量分別是181 萬個、286 萬個和601 萬個，3 套網(wǎng)格的計算結(jié)果見表2。其中，K為推力系數(shù)；K為扭矩系數(shù)；η為效率。

表2 網(wǎng)格無關(guān)性驗證結(jié)果

從表2 中可看出，網(wǎng)格1 在較低的進速系數(shù)下與網(wǎng)格2 和網(wǎng)格3 的計算結(jié)果相比，有一定的差別，而網(wǎng)格2 與網(wǎng)格3 的計算結(jié)果差別很小，綜合考慮計算的精度要求和計算成本，選取網(wǎng)格2 進行螺旋槳敞水水動力性能計算。

3.2 無空化水動力性能計算結(jié)果

模擬采用定常計算，湍流模型分別選擇Realizable k-ε 湍流模型和Standard k-ω湍流模型，比較不同湍流模型在計算時的差異。在此次模擬螺旋槳的敞水性能過程中，保持螺旋槳的轉(zhuǎn)速不變，改變來流速度和進速系數(shù)。

根據(jù)計算結(jié)果得到基于2 種不同湍流模型的PPTC 螺旋槳敞水性能曲線，并將其與試驗結(jié)果相對比，結(jié)果見圖3。

由圖3 可知，數(shù)值模擬曲線與試驗曲線的整體吻合程度較好，說明本文采用的計算模型和方法能有效預報螺旋槳的敞水水動力性能。2 種湍流模型的計算結(jié)果相差不大，考慮到k-ε湍流模型對螺旋槳推力的預報更準確，選用該模型作為空泡流場計算的基礎(chǔ)。

圖3 PPTC槳敞水性能曲線

4 空化條件下的水動力性能計算和空泡圖像

螺旋槳產(chǎn)生空泡可能造成一系列危害，例如槳葉表面材料剝蝕、螺旋槳水動力性能急劇下降、船舶噪聲和振動加劇等，因此對螺旋槳的空泡性能進行預報十分重要。

前述在德國漢堡召開的船舶推進器國際研討會對3 種工況下的PPTC 螺旋槳的空泡性能進行了試驗和觀測，分別命名為Test case 2.3.1、Test case 2.3.2 和Test case 2.3.3，螺旋槳轉(zhuǎn)速均為n =25 r/s。具體的空泡性能試驗工況參數(shù)見表3。

表3 空泡性能試驗工況參數(shù)

PPTC螺旋槳在Test Case 2.3.1、Test Case 2.3.2 和Test Case 2.3.3 等3 種工況下的空泡水動力性能計算結(jié)果及其與試驗結(jié)果的誤差見表4。從表4 中可看出，推力系數(shù)K和扭矩系數(shù)K的計算誤差均小于3%，與試驗數(shù)據(jù)很接近。

表4 PPTC螺旋槳空泡水動力性能計算結(jié)果和誤差

模擬計算得到的空泡在槳葉表面的位置分布云圖見圖4，以氣相的體積分數(shù)分布為標準，選擇氣相體積分數(shù)α=0.2 的等值面對試驗工況下的空泡范圍和位置分布情況予以說明。通過對計算結(jié)果與試驗結(jié)果進行對比可知：在Test Case 2.3.1 和Test Case 2.3.2 工況下，在螺旋槳吸力面的導邊，特別是在葉梢處出現(xiàn)了明顯的梢渦空泡，部分梢渦空泡已從葉梢曳出，脫離了螺旋槳葉面，而在葉根處出現(xiàn)了顯著的片空泡，在轂帽附近產(chǎn)生了轂渦空泡；在Test Case 2.3.3 工況下，在螺旋槳壓力面的導邊和葉根處出現(xiàn)了明顯的片空泡，但葉梢處的梢渦空泡沒有捕捉到，數(shù)值模擬計算結(jié)果與試驗結(jié)果整體吻合較好。

圖4 空泡形態(tài)及位置分布云圖

為進一步探究空泡對螺旋槳水動力性能的影響，通過改變環(huán)境壓力，進一步減少空泡數(shù)，分析在更少的空泡數(shù)下，螺旋槳水動力性能的變化規(guī)律。

圖5為螺旋槳空泡性能曲線。由圖5 可知：在一定的空泡數(shù)范圍內(nèi)，螺旋槳的水動力性能基本上保持不變；當空泡數(shù)減少到某一數(shù)值時，螺旋槳的推力和扭矩開始以近乎線性的趨勢不斷減小。

圖5 螺旋槳空泡性能曲線

為直觀地描述螺旋槳所處的旋轉(zhuǎn)流場，需定義并可視化漩渦結(jié)構(gòu)。本文選用Q準則對螺旋槳泄渦進行研究，通過Q等值面表示槳葉的三維泄渦結(jié)構(gòu)。Q準則由HUNT等提出，具體函數(shù)定義為

為便于實現(xiàn)渦量的可視化，取Q =800 s的等值面，并利用軸向速度大小對等值面進行上色，對應于Test Case 2.3.1 工況，結(jié)果見圖6。由圖6 可知，在該試驗工況下，螺旋槳尾流場的渦結(jié)構(gòu)由螺旋形的梢渦、轂渦和葉根渦等3 部分組成，而隨著與槳盤面軸向距離的增加，梢渦逐漸向內(nèi)收縮，轂渦的直徑略有增大，葉根渦在槳后傳遞過程中很快耗散并消失。

圖6 Q=800 s -2渦量等值面

對于梢渦的形成原因，參考胡健等的研究結(jié)論，根據(jù)PPTC螺旋槳在Test Case 2.3.1 工況下出現(xiàn)空泡時的表面壓力分布云圖（見圖7），通過分析可知：吸力面的壓力分布與空泡位置分布之間存在很強的關(guān)聯(lián)性，在出現(xiàn)空泡的導邊、葉梢和槳轂附近，同樣也是低壓區(qū)；對于壓力面而言，只有葉根處產(chǎn)生少量空泡，壓力較低。因此，葉梢處的壓力形成了從極大值到極小值的過渡，梢渦呈螺旋式向槳后流場泄出。

圖7 PPTC螺旋槳槳葉表面壓力分布

由上述分析可知，依據(jù)黏性流體多相流理論，采用的計算模型和方法能有效預報螺旋槳在出現(xiàn)空泡時的水動力性能。

5 空泡對泄渦頻率的影響

上文分析了螺旋槳尾流場的泄渦結(jié)構(gòu)，一般認為，當槳葉隨邊附近產(chǎn)生的旋渦頻率與螺旋槳的固有頻率接近時，螺旋槳會發(fā)出清脆的鳴音，稱為“諧鳴”現(xiàn)象，預報螺旋槳工作時旋渦產(chǎn)生的頻率，避免產(chǎn)生共鳴音，具有現(xiàn)實意義。下面以Test Case 2.3.1 為例，探究空泡對泄渦頻率的影響。

泄渦的產(chǎn)生與壓力的變化緊密相關(guān)，因此本文通過監(jiān)測壓力的變化實現(xiàn)對泄渦頻率的研究。在螺旋槳0.95R、0.90R和0.80R等3 個槳葉剖面與螺旋槳隨邊的交點處設(shè)置壓力監(jiān)測點，考慮到泄渦頻率可能出現(xiàn)在高頻區(qū)，將計算的時間步長設(shè)置為1 ×10s，此時可監(jiān)測到的最大頻率為5 000 Hz。計算完成之后，對監(jiān)測點的壓力數(shù)據(jù)作傅里葉變換，選擇有代表性的0.90R剖面上的監(jiān)測點，得到壓力的功率譜密度在頻域上的分布特征見圖8。

圖8 不同空泡數(shù)壓力監(jiān)測點頻譜圖

通過分析可知，在低頻區(qū)雖然也出現(xiàn)了峰值，但主要影響因素是螺旋槳的周期性運動造成的周期性脈動壓力，峰值主要出現(xiàn)在25 Hz和125 Hz左右，恰好是螺旋槳的軸頻和葉頻。因此，重點關(guān)注在高頻區(qū)出現(xiàn)的峰值。參考螺旋槳泄渦頻率附近存在多條線譜的結(jié)論，此時峰值的出現(xiàn)極有可能與梢渦的脫落有關(guān)。在不同空泡數(shù)工況下，高頻區(qū)的峰值分布不盡相同，基本上隨著空泡數(shù)的減少，峰值對應的頻率呈現(xiàn)先略有上升，再不斷下降的趨勢，這與空泡的產(chǎn)生改變了速度場和壓力場的特性有關(guān)。圖9 為梢渦的泄渦頻率隨空泡數(shù)的變化情況。

圖9 梢渦的泄渦頻率隨空泡數(shù)的變化情況

6 結(jié) 語

本文通過CFD方法探究了PPTC螺旋槳的敞水水動力性能和空泡特性，主要得到以下結(jié)論：

1）與試驗結(jié)果的對比分析結(jié)果表明，在均流條件下，無論有無空泡產(chǎn)生，螺旋槳推力和扭矩的數(shù)值計算結(jié)果的精度都相對較高，空泡在槳葉表面的分布與試驗結(jié)果吻合較好，定性和定量給出了螺旋槳推力和扭矩系數(shù)隨空泡數(shù)變化的趨勢。

2）利用Q準則對螺旋槳的泄渦結(jié)構(gòu)進行了研究，梢渦、葉根渦和轂渦是三維泄渦的主要組成部分。梢渦的發(fā)放頻率具有隨空泡數(shù)減少而先略有上升再逐漸下降的特征。