風(fēng)帆陣列氣動干擾特性數(shù)值分析

胡以懷， 李妍嫣， 唐娟娟

(1.上海海事大學(xué) 商船學(xué)院， 上海 201306； 2.上海船舶研究設(shè)計院， 上海 201203)

風(fēng)帆陣列氣動干擾特性數(shù)值分析

胡以懷1， 李妍嫣1， 唐娟娟2

(1.上海海事大學(xué) 商船學(xué)院， 上海 201306； 2.上海船舶研究設(shè)計院， 上海 201203)

針對某遠(yuǎn)洋散貨船設(shè)計一套風(fēng)帆陣列。為了研究風(fēng)帆之間的相互干擾現(xiàn)象，利用CFD軟件對單帆和風(fēng)帆陣列的空氣動力學(xué)特性進(jìn)行數(shù)值模擬，分析風(fēng)帆之間的干擾特征。結(jié)果表明，盡管變化趨勢一致，但由于風(fēng)帆陣列存在帆間的流場干擾，每只風(fēng)帆的升/阻力系數(shù)都較單個風(fēng)帆要小，特別是在攻角大于25°時，最大值相差30%左右，且阻力系數(shù)下降較快，這為風(fēng)帆的實船布置和應(yīng)用提供了參考。

風(fēng)帆助航；風(fēng)帆陣列；氣動干擾；CFD計算

0 引言

為獲得較大裝帆面積，充分利用船舶的甲板面積，風(fēng)帆在實際安裝使用時大都采用組合式風(fēng)帆模式。由于風(fēng)帆間會存在相互干擾的現(xiàn)象，不能簡單地將各帆的推力之和作為風(fēng)帆陣列的總推力，因此，有必要研究風(fēng)帆陣列的氣動干擾特性，以進(jìn)一步優(yōu)化風(fēng)帆的布置。

本文以我校遠(yuǎn)洋實習(xí)船為例，在船上設(shè)計安裝了8只風(fēng)帆，并建立了“育明”輪風(fēng)帆陣列模型，以研究風(fēng)帆陣列的空氣動力性能及其帆間的干擾特性。

1 風(fēng)帆陣列模型

“育明”輪是一艘載重48 000 t大型遠(yuǎn)洋散貨船，船體寬值較大。為了充分利用風(fēng)能，設(shè)計了雙排四列布置的風(fēng)帆陣列，舷側(cè)過道預(yù)留1 m的距離，風(fēng)帆中心離船體中心線13.13 m，風(fēng)帆的寬度設(shè)為20 m，風(fēng)帆陣列布置如圖1所示。

考慮到流場的對稱性，取一半的流場進(jìn)行計算，風(fēng)帆陣列的矩形計算域如圖2所示。在風(fēng)帆表面采用結(jié)構(gòu)化矩形網(wǎng)格，然后由非結(jié)構(gòu)化網(wǎng)格過渡到邊界，在靠近風(fēng)帆表面處網(wǎng)格劃分細(xì)密以獲得較為精確的流場模擬。為了易于區(qū)分，如圖3所示給風(fēng)帆進(jìn)行了編號。

圖1 “育明”輪風(fēng)帆陣列布置

圖2 風(fēng)帆陣列3D網(wǎng)格計算域(一半)

圖3 風(fēng)帆陣列附近網(wǎng)格劃分

2 風(fēng)帆陣列氣動特性數(shù)值模擬

類似單帆，在小攻角范圍內(nèi)，風(fēng)帆陣列周圍的流場是定常的，升/阻力系數(shù)最終會回歸到一個穩(wěn)定值。在大攻角狀態(tài)下，流場由于尾部漩渦而呈現(xiàn)周期性，升/阻力系數(shù)近似取一個波動周期內(nèi)的平均值。如圖4和圖5所示。

圖4 15°攻角下的升力系數(shù)和阻力系數(shù)計算迭代值

圖5 90°攻角下的升力系數(shù)和阻力系數(shù)計算迭代值

不同攻角下，風(fēng)帆陣列的速度分布云圖如圖6～圖8所示。

分析上圖可知，小攻角時帆翼上表面流體速度大于帆翼下表面的流體速度，符合儒可夫斯基升力定理中關(guān)于風(fēng)帆空氣動力特性的分析。另外，隨著攻角的逐漸變大，風(fēng)帆周圍的流場越來越復(fù)雜，帆翼前緣流體速度與其他位置速度的差距越來越明顯。當(dāng)攻角大于20°后，帆翼上表面存在流體速度幾乎為0的薄層，且帆1、帆3、帆5、帆7這種情況尤其明顯。帆2、帆4、帆6、帆8由于受到帆1、帆3、帆5、帆7尾流的影響，順著來流方向存在一段流體低速區(qū)，并逐漸變大。

圖9和圖10分別為風(fēng)帆陣列中8個帆各自的升/阻力系數(shù)曲線與單帆升/阻力系數(shù)對比曲線。由圖9可知，單帆和風(fēng)帆陣列的升力系數(shù)均近似呈兩頭小、中間大的特征，升力系數(shù)最大值均出現(xiàn)在攻角40°左右。但由于風(fēng)帆陣列存在帆間干擾，8個風(fēng)帆各自的升力系數(shù)都比單帆來得小。在攻角α=0°～25°范圍內(nèi)，由于帆2、帆4、帆6、帆8基本上處于帆1、帆3、帆5、帆7的尾流影響區(qū)中，受到的干擾較大，其升力系數(shù)比單帆要小很多。帆8升力系數(shù)最大，但也僅為單帆的45%左右。帆1、帆3、帆5、帆7的升力系數(shù)在攻角α=0°～20°范圍內(nèi)基本與單帆相同，帆間的干擾影響較小。在攻角α=25°～90°范圍內(nèi)，8個風(fēng)帆的升力系數(shù)接近，但總體呈下降趨勢，與單帆相比差距較大，帆間相互干擾作用更加顯著。

圖6 0°攻角下風(fēng)帆陣列速度云圖

圖7 20°攻角下風(fēng)帆陣列速度云圖

圖8 40°攻角下風(fēng)帆陣列速度云圖

圖9 風(fēng)帆陣列升力系數(shù)曲線

圖10 風(fēng)帆陣列阻力系數(shù)曲線

由圖10可知，單帆、風(fēng)帆陣列的阻力系數(shù)曲線也呈兩頭小、中間大的特征，但阻力系數(shù)最大值均在攻角60°左右，且8個風(fēng)帆的阻力系數(shù)均比單帆小。在攻角α=0°～25°范圍內(nèi)，帆1、帆3、帆5、帆7的阻力系數(shù)與單帆時的阻力系數(shù)較接近，帆間干擾影響較小，而帆2、帆4、帆6、帆8與單帆的相差稍大。隨著攻角的繼續(xù)增大，尤其在攻角α>50°后，帆間干擾作用越來越明顯，8個帆與單帆的差值愈來愈大。

計算時假設(shè)風(fēng)帆轉(zhuǎn)動時每個風(fēng)帆轉(zhuǎn)動的角度是相同的，因此風(fēng)帆陣列產(chǎn)生的總升力即為每個風(fēng)帆產(chǎn)生的升力之和，風(fēng)帆列陣的理論合成結(jié)果與風(fēng)帆陣列數(shù)值模擬結(jié)果比較如圖11所示。

觀察圖11可以發(fā)現(xiàn)，8個風(fēng)帆產(chǎn)生的理論總升力/阻力并非是單面帆產(chǎn)生的升力/阻力相加。理論計算結(jié)果(單帆×8)遠(yuǎn)大于風(fēng)帆列陣的數(shù)值計算結(jié)果，最大值相差約30%，但整體變化趨勢仍保持一致。理論升力最大值依然出現(xiàn)在攻角40°左右，阻力最大值依然出現(xiàn)在攻角60°左右。當(dāng)攻角α>50°后，轉(zhuǎn)帆角比較大時，帆之間的重疊面積增大，帆間干擾明顯增大，風(fēng)帆陣列的阻力下降比單帆快；在α=90°時，風(fēng)帆陣列的總阻力僅為單帆阻力總和的25%左右。

圖11 風(fēng)帆列陣?yán)碚摵铣膳c計算總升力系數(shù)、阻力系數(shù)對比

3 結(jié)論

雖然實船應(yīng)用中每艘船可以安裝多個風(fēng)帆，但各風(fēng)帆之間存在空氣動力學(xué)的相互干擾，實際的總升力和阻力要小于理論總和。數(shù)值分析表明，風(fēng)帆陣列中每只風(fēng)帆的升/阻力系數(shù)都較單個風(fēng)帆的小，特別是在攻角大于25°時，最大值相差30%左右，且阻力系數(shù)下降較快。在實際應(yīng)用時，建議各帆的轉(zhuǎn)角可以不盡相同，應(yīng)采用一種更優(yōu)化的風(fēng)帆間隔或轉(zhuǎn)動方式，以避免風(fēng)帆陣列之間相互干擾。

[ 1 ] 王立明.風(fēng)帆輔助推進(jìn)裝置及控制系統(tǒng)的建模與仿真研究[D].武漢:武漢理工大學(xué)，2010.

[ 2 ] 王福軍.計算流體動力學(xué)分析[M].北京:清華大學(xué)出版社，2004.

[ 3 ] 吳光中,宋婷婷,張毅.Fluent基礎(chǔ)入門與案例精通[M].北京:電子工業(yè)出版社，2012.

[ 4 ] 朱紅鈞,林元華,謝龍漢.Fluent流體分析及仿真實用教程[M].北京：人民郵電出版社，2010.

[ 5 ] 楊龍霞.風(fēng)帆助航遠(yuǎn)洋船的翼帆性能及其機(jī)槳配合研究[D].上海:上海交通大學(xué),2013.

[ 6 ] 張軍,趙峰.大攻角翼型繞流的流場測試及數(shù)值計算[J].中國造船，2000,41(2):22-27.

[ 7 ] 王元靖,范召林,呂全洲,等.大攻角流動的非定常特性分析[J].空氣動力學(xué)學(xué)報，2007,25(2):145-149.

[ 8 ] 袁新,江學(xué)忠.翼型大攻角低速分離流動的數(shù)值模擬[J].工程熱物理學(xué)報，1999,20(2):161-165.

[ 9 ] 李艷凱,范秦寅,董非,等.風(fēng)機(jī)葉片翼型大攻角下氣動性能的數(shù)值分析[J].太陽能學(xué)報，2010,31(12):1639-1643.

[10] 王獻(xiàn)孚.船用翼理論[M].北京：國防工業(yè)出版社，1998.

[11] 劉文霞.帆板搖帆運(yùn)動帆氣動力特性的數(shù)值模擬研究[D].青島：中國海洋大學(xué)，2007.

[12] 張福成.基于Fluent的風(fēng)力發(fā)電機(jī)流場仿真研究[D].大連：大連理工大學(xué)，2012.

Numerical Analsis of Aerodynamic Characteristics between Aligned Sails

HU Yihuai1, LI Yanyan1, TANG Juanjuan2

(1.Merchant Marine College, Shanghai Maritime University, Shanghai 201306, China；2.Shanghai Merchant Ship Design and Research Institute, Shanghai 201203, China)

A set of aligned sails is designed for an ocean-going bulk carrier on the main deck. To analysis the disturbance between these aligned sails the dynamic characteristics of the sails were simulated with CFD software. It is proved that the lift and drag coefficients have the same tend against attack angle with that of a single sail, but the disturbance between the aligned sails does exist. This disturbance could induce less lift and drag coefficients of the sail by 30 percent compared to a single sail, especially under larger attack angles beyond 25°. These results could provide reference for the practical application of sails on board ship.

sail assisted ship； aligned sail set； aerodynamic disturbance； CFD calculation

上海市科委重點(diǎn)科研項目(項目編號：08210511800)。

胡以懷(1964-)，男，教授，研究方向為船舶新能源利用和輪機(jī)仿真與診斷。

1000-3878(2017)01-0006-06

U662

A