螺旋槳四象限水動力性能數值計算

張文璨，董國祥，陳偉民，杜云龍，任海奎

(上海船舶運輸科學研究所 航運技術與安全國家重點實驗室，上海 200135)

?

螺旋槳四象限水動力性能數值計算

張文璨，董國祥，陳偉民，杜云龍，任海奎

(上海船舶運輸科學研究所 航運技術與安全國家重點實驗室，上海 200135)

采用求解RANS方程的方法并運用計算流體力學(Computational Fluid Dynamics, CFD)軟件Fluent，對荷蘭MARIN船模水池B系列螺旋槳的敞水四象限水動力性能進行CFD數值模擬；同時，將數值模擬結果與B系列螺旋槳敞水四象限試驗圖譜進行比較。結果表明，計算誤差較小，驗證了計算策略的正確性，也為后續進行螺旋槳四象限水動力性能研究打下了基礎。

B系列螺旋槳；四象限；水動力性能；數值模擬

0　引　言

敞水槳四象限(前進中正車、前進中倒車、后退中正車、后退中倒車)水動力性能數值模擬是當前螺旋槳計算流體力學(Computational Fluid Dynamics, CFD)研究熱點，但目前相關研究較少。由于螺旋槳四象限動態模型領域在表達上有困難，因此迄今為止船舶運動仿真一般僅限于前向，即第一象限內的航行。螺旋槳在反轉和反向制動時將進入推力和扭矩的第二象限，而在船舶倒航時將進入第三象限，這超出了一般螺旋槳圖譜的表達范圍，為螺旋槳逆轉工況下的運動仿真帶來了一定的困難。進行敞水槳四象限水動力性能研究不僅具有諸多理論意義，而且是解決諸多艦船操縱性和速航性課題(離靠泊、艦船避碰、各種特殊機動)的關鍵所在[1-4]。在該背景下，選取MARIN B系列螺旋槳中有代表性的螺距比P/D=1.0、盤面比AE/AO=0.7的四葉槳B4-7010，對其敞水四象限水動力性能進行數值模擬，并將分析后的結果與B系列螺旋槳敞水四象限試驗圖譜進行對比。結果表明：計算值與圖譜值有較高的貼合性，采取的計算策略是成功的。

1　CFD計算模型與計算策略

1.1計算模型

以選取MARIN B系列螺旋槳中有代表性的螺距比P/D=1.0、盤面比AE/AO=0.7的四葉槳B4-7010為研究對象，計算采用模型尺度，槳模直徑D=0.2 m。B4-7010模型的幾何要素見表1。

表1　B4-7010模型幾何要素

圖1　B4-7010螺旋槳幾何模型

利用Fluent前處理軟件Gambit建立B4-7010螺旋槳的幾何模型(見圖1)，具體步驟為：

(1) 螺旋槳二維幾何數值進行坐標轉化，得到螺旋槳三維空間坐標點，并將三維空間坐標點文件導入到Gambit中；

(2) 遵循“由點生線、由線生面、由面生體”的原則生成螺旋槳實體。

1.2計算策略

1.2.1計算域及網格劃分

計算域為與螺旋槳同軸的圓柱體(圓柱體的高度和直徑分別)為螺旋槳直徑的30倍及10倍)，其中：速度入口距離螺旋槳約7倍螺旋槳直徑；壓力出口距離螺旋槳約20倍螺旋槳直徑。

網格劃分是指計算區域的離散化，即將空間上連續的計算區域劃分為多個子區域，并確定每個區域中的節點。網格劃分的本質是用有限個離散的點代替原來的連續空間。數學上，生成網格后(離散化后)連續的控制方程將被離散化，描寫流動和傳熱的偏微分方程將轉化為各個節點上的代數方程組。網格劃分是CFD數值計算過程中最重要的環節，直接影響模擬的精度和效率。若網格過疏，則會導致模擬精度不高；若網格過密，則又會導致計算量增大及計算難以收斂[5-6]。

為便于進行網格劃分，需要將無限流域流場空間劃分為多個小空間。計算空間域被分為外部計算域、過渡區和螺旋槳旋轉域等3個部分。

(1) 最外層為外部計算域，該區域中壓力和速度分布變化并不劇烈，主要反映導管螺旋槳的來流和尾流情況，因此可劃分較粗的結構化網格；

(2) 在外部計算域和螺旋槳旋轉域的過渡區中，速度壓力分布變化較為劇烈，需要加入較為細化的非結構網格；

(3) 在螺旋槳旋轉域中，速度壓力分布變化非常劇烈，應采用比過渡區更加細致的非結構化網格進行計算。

在這3個區域中，外部計算域和過渡區均為靜態域，而螺旋槳旋轉域采用多重旋轉坐標系(Multiple Rotational Frame, MRF)模型，設置為繞x軸正方向旋轉，角速度為20 r/s。

在進口邊界處設置為Velocity-inlet條件；出口邊界定義為Pressure-outlet邊界；外部計算域圓柱體表面設為Symmetry邊界；所有壁面均設為無滑移固壁條件。

由于第一、二象限為前進，第三、四象限為后退，因此第一、二象限與第三、四象限的來流方向相反，同時螺旋槳旋轉域及外部計算域的網格劃分也會有些許差異。第一、二象限計算域網格劃分總數為1 817 278，其中內部計算域網格劃分數為909 515，偏斜率>0.8的網格數為6，網格質量達到要求；第三、四象限計算域網格劃分總數為1 594 806，其中內部計算域網格劃分數為810 301，偏斜率>0.8的網格數為0，網格質量同樣達到要求。

第一、二象限和第三、四象限計算螺旋槳旋轉域及全部計算空間域網格劃分情況見圖2和圖3。圖2、圖3中：Gx，Gy，Gz為O-xyz坐標系中Ox，Oy，Oz坐標軸的方向。

a) 螺旋槳旋轉域網格劃分情況

b) 全部計算空間域網格劃分情況

1.2.2計算流體力學基本方程

連續性方程及動量方程分別為式(1)和式(2)。

(1)

(2)

1.2.3湍流模型

鑒于所模擬流場的特性，非定常湍流計算采用RNGk-ε雙方程模型。其形式如下[7]

(3)

(4)

1.2.4計算參數設置

計算時運用SIMPLE半隱式連接壓力方程；梯度插值采用Green-Gauss Cell Based方案；壓力插值采用PRESTO!格式；動量、湍流動能及湍流耗散度的離散均選用二階迎風差分格式[8]。

a) 螺旋槳旋轉域網格劃分情況

b) 全部計算空間域網格劃分情況

2　計算及數據分析

2.1計算策略正確性驗證

無因次推力系數[9-10]為

(5)

無因次扭矩系數為

(6)

這些系數是水動力螺距角β的函數。

(7)

式(5)～式(7)中：T為螺旋槳推力；Q為螺旋槳扭矩；ρ為水的密度，取998 kg/m3；VA為進速，單位為m/s；n為螺旋槳轉速，取20 r/s；D為螺旋槳模型直徑，取0.2 m。

依據計算結果及式(5)～式(7)，根據不同象限計算出每個進速VA所對應的水動力螺距角β、推力系數CT和扭矩系數CQ，見表2～表5。

表2　一象限計算數據分析結果

表3　二象限計算數據分析結果

表4　三象限計算數據分析結果

表5　四象限計算數據分析結果

根據以上數據，繪制出B系列螺旋槳四象限水動力性能圖譜值和計算值曲線(見圖4)。

由圖4可知，圖譜值與計算值的貼合性較高，計算結果較為滿意，證明了所采取的計算策略的正確性。

2.2螺旋槳四象限水動力性能的影響分析

以VA=1.4 m/s時的工況為例，對4個象限的螺旋槳葉背、葉面壓力分布情況進行分析，研究不同象限螺旋槳推力的變化情況。4個象限的螺旋槳葉背、葉面壓力云圖見圖5和圖6。

圖4　四象限水動力性能曲線圖

a) 一象限

b) 二象限

c) 三象限

d) 四象限

由圖5可知，螺旋槳葉背對一、四象限而言為吸力面；而對二、三象限而言為壓力面。吸力面壓力由葉根到葉梢及由導邊到隨邊均逐漸降低；而壓力面壓力由葉根到葉梢逐漸升高，二象限壓力由導邊到隨邊逐漸升高，三象限壓力則由導邊到隨邊先降低再升高，在葉面中部有一段徑向分布的低壓區。

由圖6可知，螺旋槳葉面對一、四象限而言為壓力面；而對二、三象限而言為吸力面。壓力面壓力由葉根到葉梢及由隨邊到導邊逐漸升高；而吸力面壓力由葉根到葉梢及由隨邊到導邊均逐漸降低。

3　結　語

利用RANS雷諾時均方程組對B4-7010螺旋槳的敞水四象限水動力性能進行了數值模擬。首先，利用GAMBIT對螺旋槳進行了三維建模；其次，對螺旋槳的計算域進行了設計和網格劃分，并設置了合適的邊界條件；最后將其導入FLUENT進行計算。由計算結果與圖譜值的對比可知，計算結果與圖譜有著較高的貼合性。這既證明了計算結果的精確性，又證明了所采用的計算策略的正確性。同時，為后續導管螺旋槳、可調螺距螺旋槳四象限水動力性能的研究打下了基礎。

a) 一象限

b) 二象限

c) 三象限

d) 四象限

[1]李理，劉可，李超，等. 螺旋槳四象限水動力性能數值模擬及應用[J]. 艦船科學技術, 2012, 34(7): 8-14.

[2]孫洪波，洪碧光，谷元全. 螺旋槳的四象限推力系數和轉矩系數的MATLAB實現方法[J]. 大連海事大學學報，2006，32(2): 42-44.

[3]孫永華. 荷蘭船模水池B系列螺旋槳新舊圖譜的比較[J]. 船舶, 1997(3): 30-36.

[4]孫洪波. 螺旋槳逆轉工況下船舶運動建模與仿真[D]. 大連:大連海事大學，2007.

[5]盛振邦，劉應中.船舶原理[M].上海：上海交通大學出版社,2004.

[6]王獻孚, 周樹信, 陳澤梁, 等. 計算船舶流體力學[M]. 上海:上海交通大學出版社，1992.

[7]黃勝, 王超, 王詩洋. 不同湍流模型在螺旋槳水動力性能計算中的應用與比較[J]. 哈爾濱工程大學學報，2009，30(5): 481-485.

[8]王詩洋, 解學參, 王超, 等. 參數選擇對螺旋槳水動力性能的影響[J]. 船舶工程，2010，32(6): 4-8.

[9]付頤鑫.船舶螺旋槳敞水性能CFD模擬[D].大連:大連海事大學,2012.

[10]王亮.船舶螺旋槳變參數水動力性能研究[D].大連:大連理工大學,2012.

Numerical Simulation of Four-Quadrant Hydrodynamic Performance of Propeller

ZHANG Wencan,DONG Guoxiang,CHEN Weimin,DU Yunlong,REN Haikui

(StateKeyLaboratoryofNavigationandSafetyTechnology，ShanghaiShip&ShippingResearchInstitute，Shanghai200135，China)

Four-quadrant hydrodynamic performance numerical simulation of open-water propellers is a hotspot in the field of propeller numerical simulation. However, there has not been much research work on the four-quadrant hydrodynamic performance of propellers published. Through solving RANS equations and using CFD software Fluent, CFD numerical simulation of MARIN series B propeller is carried out to investigate its open-water four-quadrant hydrodynamic performance. The simulation results are compared with the test-graph results to prove the validity of calculation strategy. This research may inspire further research on the four-quadrant hydrodynamic performance of propellers.

MARIN series B propeller; four-quadrant; hydrodynamic performance; numerical simulation

2015-07-31

張文璨(1990—)，男，遼寧沈陽人，碩士生，主要從事船舶推進器水動力性能研究。

1674-5949(2016)01-001-07

U661.31+3；U664.33

A