自由面條件下水平軸潮流能葉輪水動力研究

王樹齊，張亮，徐剛，朱仁慶

(1. 江蘇科技大學(xué) 船舶與海洋工程學(xué)院，江蘇 鎮(zhèn)江 212003； 2. 哈爾濱工程大學(xué) 船舶工程學(xué)院，黑龍江 哈爾濱 150001)

?

自由面條件下水平軸潮流能葉輪水動力研究

王樹齊1，張亮2，徐剛1，朱仁慶1

(1. 江蘇科技大學(xué) 船舶與海洋工程學(xué)院，江蘇 鎮(zhèn)江 212003； 2. 哈爾濱工程大學(xué) 船舶工程學(xué)院，黑龍江 哈爾濱 150001)

為了利用近水面較高的潮流流速，潮流能葉輪一般安裝在接近自由面的位置，導(dǎo)致葉輪的水動力特性受到自由面的影響。本文基于CFD方法，通過改變?nèi)~尖浸沒水深，對水平軸葉輪在自由面條件下的水動力特性展開研究。研究結(jié)果表明：葉輪軸向載荷和能量利用率平均值隨葉尖浸沒水深的增加而增加，并逐漸接近于無自由面時的計算結(jié)果；葉輪軸向載荷和能量利用率瞬時值產(chǎn)生波動，波動幅值隨葉尖浸沒水深的增加而減小；當(dāng)葉尖浸沒水深為0.79D時，軸向載荷系數(shù)及能量利用率已基本不隨葉片方位角的變化而改變，即自由面已幾乎無影響。研究結(jié)果可為潮流電站設(shè)計及工程應(yīng)用提供依據(jù)和參考。

葉尖浸沒水深；潮流能；水平軸葉輪；自由面；水動力載荷

在實際水平軸潮流電站中，為了充分利用表面較高的流速，水平軸潮流能水輪機(jī)一般由漂浮式載體或固定式載體支撐其遠(yuǎn)離海底而靠近自由面，水輪機(jī)葉輪面臨水中空化、自由面以及重力水波引起的速度梯度等不同程度的影響，這種影響一方面導(dǎo)致葉輪功率特性變差，另一方面則可能導(dǎo)致干擾載荷累積、使水輪機(jī)或支撐結(jié)構(gòu)疲勞甚至斷裂破壞，釀成重大事故。

水平軸潮流能葉輪[1-5]近水面運行時，葉輪的運動會造成自由面的變形，自由面的變形反過來會影響葉輪的性能。葉輪和自由面相互影響的過程會涉及到旋渦、湍流與自由面的非定常、非線性相互作用，使水平軸潮流能葉輪周圍流場變得非常復(fù)雜。因此，準(zhǔn)確預(yù)報水平軸葉輪在自由面條件下的水動力性能顯得尤為重要。

Bahaj等[6-9]在不同翼型、槳距角以及兩種不同葉尖浸沒水深條件下，采用模型試驗方法，研究了水輪機(jī)性能、葉片空化以及雙轉(zhuǎn)子干擾問題，得出了水輪機(jī)在不同葉尖浸沒水深時的功率及載荷規(guī)律。Coiro等[10]在拖曳水池中對S805翼型水翼進(jìn)行了測試。水輪機(jī)直徑0.8 m,水輪機(jī)浸沒水深在自由面下1.3 m。拖車速度范圍1.3～3.5 m/s，槳距角范圍-4～4°，試驗過程中沒有研究空化問題。Maganga等[11]針對來流湍流度8%～25%，偏航角-10～20°，自由面以下深度0.94D、1.37D和2.04D(D為葉輪直徑)等條件下在一個封閉的循環(huán)水槽中對3葉片水輪機(jī)進(jìn)行實驗研究。實驗揭示了來流有較高湍流度時，效率與推力系數(shù)下降約9%，在較大偏航角時，水輪機(jī)性能明顯下降，并給出了不同葉尖浸沒水深下的效率和推力系數(shù)。國內(nèi)針對不同葉尖浸沒水深下水平軸潮流能葉輪水動力性能鮮有研究。

從以上研究可以看出，針對自由面條件下水平軸水輪水動力問題的研究主要采用模型試驗的方法，該方法主要不足之處在于：模型試驗一般測試的是水輪機(jī)旋轉(zhuǎn)過程中的平均效率及載荷，無法體現(xiàn)自由面對水輪機(jī)瞬時性能的影響；模型試驗無法給出復(fù)雜的流場細(xì)節(jié)。因此，本文基于ANSYS CFX方法，研究葉尖浸沒水深對水平軸葉輪水動力性能及載荷的影響規(guī)律，為潮流電站設(shè)計提供設(shè)計依據(jù)，以及為工程應(yīng)用提供參考和借鑒。

1 CFD數(shù)值模擬

1.1 計算模型建立

葉輪計算模型采用哈爾濱工程大學(xué)深海中心可再生能源的實驗?zāi)Ｐ?。該實驗?zāi)Ｐ偷娜~輪直徑D為0.7 m，葉片數(shù)為2。葉片采用變弦長及變槳距角形式，除葉片根部外，展向各位置處的截面翼型均為S809。葉輪模型如圖1所示。計算模型建好之后，為了提交計算精度，需建立一個足夠大的計算域來模擬流場，計算域的出口與入口都是正方形，邊長是10倍葉輪直徑；入口到葉輪旋轉(zhuǎn)平面的距離是5倍葉輪直徑；出口到葉輪旋轉(zhuǎn)平面的距離為15倍葉輪直徑。

1.2 網(wǎng)格劃分

網(wǎng)格劃分是CFD模擬過程中較為關(guān)鍵的環(huán)節(jié)，也是比較耗時的環(huán)節(jié)，網(wǎng)格質(zhì)量的好壞直接影響數(shù)值模擬精度和計算效率。本文采用結(jié)構(gòu)化網(wǎng)格和非結(jié)構(gòu)化網(wǎng)格相結(jié)合的網(wǎng)格形式對流場進(jìn)行離散，旋轉(zhuǎn)域由于葉片外形復(fù)雜，采用非結(jié)構(gòu)化網(wǎng)格形式，靜止域采用結(jié)構(gòu)化網(wǎng)格形式。文獻(xiàn)[12]對網(wǎng)格數(shù)量及收斂性進(jìn)行了研究，結(jié)果表明：當(dāng)網(wǎng)格數(shù)量達(dá)到231萬，模型表面網(wǎng)格第一層高度0.000 5 m時，繼續(xù)加密網(wǎng)格，計算結(jié)果已基本上沒有變化，因此該文計算中采用此網(wǎng)格數(shù)類型。如圖1所示。

1.3 邊界條件設(shè)置

模擬自由面的難點是需要模擬兩相流，即水和空氣，水汽交界面采用體積分?jǐn)?shù)法(VOF)處理。在CFD數(shù)值模擬中，多相流模擬的穩(wěn)定性及收斂性是比較差的，對網(wǎng)格要求比較高。具體設(shè)置如下：設(shè)置大氣壓為參考壓力，給定重力加速度的方向；將入口邊界設(shè)置為速度入口條件，給定均勻來流速度 、湍流參數(shù)及水和空氣的體積分?jǐn)?shù)，體積分?jǐn)?shù)根據(jù)葉尖浸沒水深H(葉尖浸沒水深是指葉尖距水面最短距離)得到；出口邊界設(shè)置為開放的壓力邊界，給定相對壓力及水和空氣的體積分?jǐn)?shù)，相對壓力設(shè)為水的靜壓，水和空氣的體積分?jǐn)?shù)與入口一致；流體計算域的左右兩側(cè)和底面設(shè)置為自由滑動壁面；流體計算域的頂部設(shè)置為開放的壓力邊界，相對壓力設(shè)為0，水的體積分?jǐn)?shù)為0，空氣的體積分?jǐn)?shù)為1；葉片及輪轂表面設(shè)置為不可滑移壁面；給定旋轉(zhuǎn)域旋轉(zhuǎn)角速度；靜止域和旋轉(zhuǎn)域之間通過交界面連接。計算中湍流模型采用SST模型，求解器為瞬態(tài)求解器，時間步長為葉輪旋轉(zhuǎn)2°所用的時間。

多相流模擬時經(jīng)常會出現(xiàn)發(fā)散的現(xiàn)象，特別在剛開始計算的時候。因為在剛開始計算時，由于流體及物體的突然運動，導(dǎo)致相與相之間的交界面變化顯著，很容易導(dǎo)致計算發(fā)散。為了避免出現(xiàn)類似錯誤，先計算得到葉輪不旋轉(zhuǎn)時穩(wěn)定的自由表面條件，讓自由面適應(yīng)網(wǎng)格模型，以此為初始條件，計算葉輪旋轉(zhuǎn)時的水動力性能。

圖1 模型及網(wǎng)格示意圖Fig.1 Schematic and mesh model

2 計算結(jié)果分析

2.1 計算值與試驗值對比

采用上述網(wǎng)格模型及邊界條件設(shè)置對哈爾濱工程大學(xué)深海中心可再生能源所試驗狀態(tài)(流速U=1.5 m/s，H=0.79D)進(jìn)行模擬，得到葉輪在自由面條件下的能量利用率及軸向載荷系數(shù)。水輪機(jī)的能量利用率CP是表示葉輪吸收潮流動能能力的參數(shù)，葉輪軸向載荷系數(shù)CZ是葉輪沿旋轉(zhuǎn)軸方向載荷的無量綱化，葉尖速比是葉尖旋轉(zhuǎn)線速度與來流速度的比值，表達(dá)式為

(1)

(2)

(3)

式中：P為葉輪功率，W；ρ為流體密度，kg/m3；U為流體流速，m/s；FZ為葉輪旋轉(zhuǎn)軸方向載荷,N；ω為葉輪旋轉(zhuǎn)角速度，rad/s；R為葉輪半徑，m；n為葉輪轉(zhuǎn)速，r/min。

計算結(jié)果的平均值與試驗值對比如圖2所示，圖中實驗值是在哈爾濱工程大學(xué)船舶拖曳水池中進(jìn)行測試的，平均值是指葉輪旋轉(zhuǎn)一圈時計算結(jié)果的算術(shù)平均。從圖2中可以看出：在高速比時計算值與試驗值吻合比較好，基本一致，但在低速比時計算值高于實驗值，速比越低，差異越明顯。這是因為能量利用率最大值所對應(yīng)的速比為最優(yōu)速比，而在葉片設(shè)計時，最優(yōu)速比一般與翼型的最優(yōu)攻角對應(yīng)。當(dāng)速比大于最優(yōu)速比時，葉輪旋轉(zhuǎn)時葉片各截面翼型的實際攻角小于翼型的最優(yōu)攻角，流動不發(fā)生分離，計算值與試驗值吻合比較好；但當(dāng)速比小于最優(yōu)速比時，葉片各截面翼型的實際攻角大于翼型的最優(yōu)攻角，流動產(chǎn)生分離，甚至出現(xiàn)失速現(xiàn)象，計算值與試驗值存在一定偏差。整體上來說，計算值與試驗值吻合較好，驗證了該數(shù)值方法研究水平軸葉輪在自由面條件下的水動力的可行性與準(zhǔn)確性。

圖2 計算值與試驗值對比Fig.2 Comparison between calculated and experimental values

2.2 自由面影響分析

在上述數(shù)值計算的基礎(chǔ)上，分別對另外兩個不同葉尖浸沒水深(H=0，H=0.21D)和無自由面進(jìn)行模擬，分別計算三種不同轉(zhuǎn)速(150、230、310 r/min)下的能量利用率及軸向載荷系數(shù)。不同轉(zhuǎn)速下計算結(jié)果的平均值如圖3所示，圖例中∽表示無自由面，流速1.5 m/s。從圖3可以看出：隨著葉尖浸沒水深的增加，能量利用率和軸向載荷系數(shù)逐漸增加，且轉(zhuǎn)速越大，增加越明顯；葉尖浸沒水深為0.79D時的計算結(jié)果與無自由面時的計算結(jié)果非常接近，說明葉尖浸沒水深達(dá)到0.79D時，自由面對葉輪的水動力已幾乎無影響。這是因為葉輪旋轉(zhuǎn)使自由面產(chǎn)生興波，且葉尖浸沒水深越小，興波越明顯，從而使葉輪水動力性能及載荷隨葉尖浸沒水深的減小而減小。

不同葉尖浸沒水深下軸向載荷系數(shù)及能量利用率瞬時值隨葉片方位角的變化規(guī)律如圖4所示，方位角為0°時，葉片軸線與水平面平行。從圖4中可以得知：自由面使葉輪軸向載荷系數(shù)及能量利用率的瞬時值產(chǎn)生波動，波動幅值隨葉尖浸沒水深的增加而減??；當(dāng)葉尖浸沒水深達(dá)到0.79D時，葉輪軸向載荷系數(shù)及能量利用率的瞬時值的波動已經(jīng)很小了，說明自由面對葉輪的水動力已幾乎沒影響。這是因為葉尖浸沒水深越小，葉輪旋轉(zhuǎn)過程中水面興波越大，能量耗散越多，從而導(dǎo)致葉輪水動力性能及載荷降低。載荷的波動影響葉輪結(jié)構(gòu)的疲勞壽命，能量利用率的波動不利于電能的穩(wěn)定輸出，在結(jié)構(gòu)和電控系統(tǒng)設(shè)計中需要特別考慮。

圖3 不同葉尖浸沒水深下的平均計算值Fig.3 The average calculated values of different tip-immersion

圖4 不同葉尖浸沒水深的瞬時計算值Fig.4 The instantaneous calculated values of different tip-immersion

3 流場特性分析

圖5給出了計算模型在轉(zhuǎn)速為230 r/min，葉尖浸沒水深為0時，葉片在不同方位角下表面壓力分布。從圖5中可以看出：葉片迎流面的壓力明顯高于背流面的壓力，從而形成來流方向的壓力差，即軸向載荷；當(dāng)葉片軸線與水平面平行時，兩個葉片的壓力分布相同；但當(dāng)葉片軸線與水平面垂直時，兩個葉片的壓力分布明顯不同，即靠近自由液面的葉片的表面壓力比較低。由于葉片在不同位置時表面壓力不同，從而導(dǎo)致葉片載荷在旋轉(zhuǎn)過程中發(fā)生波動。

圖6給出了計算模型在轉(zhuǎn)速230 r/min時，不同葉尖浸沒水深下自由面示意圖，由圖6可以看出：當(dāng)葉尖浸沒水深為0時，自由液面出現(xiàn)非常明顯的興波，在向下游傳播過程中，興波逐漸減小并向兩側(cè)擴(kuò)散；隨著葉尖浸沒水深的增加，自由液面興波逐漸減小，能量耗散降低，從而導(dǎo)致水輪機(jī)軸向載荷及能量利用率增加。

圖5 葉片表面壓力分布Fig.5 Pressure distribution of blade

圖6 自由面示意圖Fig.6 The diagram of free surface

4 結(jié)論

采用ANSYS CFX對三維水平軸葉輪在自由面條件下的水動力特性進(jìn)行了計算分析，研究結(jié)果表明：

1)自由面使葉輪軸向載荷系數(shù)及能量利用率的瞬時值產(chǎn)生波動，波動幅值隨葉尖浸沒水深的增加而減小，載荷的波動影響葉輪結(jié)構(gòu)的疲勞壽命，能量利用率的波動不利于電能的穩(wěn)定輸出，在結(jié)構(gòu)和電控系統(tǒng)設(shè)計中需要特別考慮；

2)軸向載荷系數(shù)及能量利用率的平均值隨葉尖浸沒水深的增加而增加，并逐漸接近于無自由面時的計算結(jié)果；

3)當(dāng)葉尖浸沒水深為0.79D時，軸向載荷系數(shù)及能量利用率已基本不隨葉片方位角的變化而改變，即自由面已幾乎無影響。

[1]沈云, 李龍, 朱多彪. 水平軸潮流水輪機(jī)轉(zhuǎn)輪尾流特性數(shù)值分析[J]. 水電能源科學(xué), 2013, 31(10): 149-151.

SHEN Yun, LI Long, ZHU Duobiao. Numerical analysis of wake performance of horizontal-axis tidal current turbine[J]. Water resources and power, 2013, 31(10): 149-151.

[2]MILNE I A, DAY A H, SHARMA R N, et al. Blade loading on tidal turbines for uniform unsteady flow[J]. Renewable energy, 2015, 77: 338-350.

[3]MILNE I A, DAY A H, SHARMA R N, et al. Blade loads on tidal turbines in planar oscillatory flow[J]. Ocean engineering, 2013, 60: 163-174.

[4]BAHAJ A S, MOLLAND A F, CHAPLIN J R, et al. Power and thrust measurements of marine current turbines under various hydrodynamic flow conditions in a cavitation tunnel and a towing tank[J]. Renewable energy, 2007, 32(3): 407-426.

[5]LEE J H, PARK S, KIM D H, et al. Computational methods for performance analysis of horizontal axis tidal stream turbines[J]. Applied energy, 2012, 98: 512-523.

[6]BAHAJ A S, BATTEN W M J, MCCANN G. Experimental verifications of numerical predictions for the hydrodynamic performance of horizontal axis marine current turbines[J]. Renewable energy, 2007, 32(15): 2479-2490.

[7]BATTEN W M J, BAHAJ A S, MOLLAND A F, et al. Experimentally validated numerical method for the hydrodynamic design of horizontal axis tidal turbines[J]. Ocean engineering, 2007, 34(7): 1013-1020.

[8]BATTEN W M J, BAHAJ A S, MOLLAND A F, et al. The prediction of the hydrodynamic performance of marine current turbines[J]. Renewable energy, 2008, 33(5): 1085-1096.

[9]BATTEN W M J, BAHAJ A S, MOLLAND A F, et al. Hydrodynamics of marine current turbines[J]. Renewable energy, 2006, 31(2): 249-256.

[10]COIRO D P, MAISTO U, SCHERILLO F, et al. Horizontal axis tidal current turbine: numerical and experimental investigations[C]//Offshore Wind and Other Marine Renewable Energies in Mediterranean and European Seas. Rome, Italy: European Seminar, 2006.

[11]MAGANGA F, GERMAIN G, KING J, et al. Experimental characterisation of flow effects on marine current turbine behaviour and on its wake properties[J]. IET renewable power generation, 2010, 4(6): 498-509.

[12]王樹齊, 肖鋼, 張亮, 等. 潮流能水平軸水輪機(jī)支撐立柱干擾研究[J]. 華中科技大學(xué)學(xué)報: 自然科學(xué)版, 2014, 42(4): 81-85.

WANG Shuqi, XIAO Gang, ZHANG Liang, et al. Supporting column influence analysis of Horizontal axis tidal current turbine[J]. Journal of huazhong university of science and technology: natural science edition, 2014, 42(4): 81-85.

Hydrodynamic analysis of a tidal current impeller in a horizontal axis under the condition of a free surface

WANG Shuqi1, ZHANG Liang2, XU Gang1, ZHU Renqing1

(1. School of Naval Architecture and Ocean Engineering，Jiangsu University of Science and Technology, Zhenjiang 212003, China；2. College of Shipbuilding Engineering，Harbin Engineering University, Harbin 150001, China)

To take advantage of the high tidal current velocity near the surface of water, a tidal current impeller is generally installed near the free surface, which influences the hydrodynamic characteristics of the impeller. Based on the CFD numerical method, hydrodynamic characteristics of horizontal axis tidal current impellers were studied in this paper, in consideration of free surface and by changing the tip-immersion. The results show that the average axial load and energy utilization ratio of the impeller increase with the increase of the tip-immersion, and are gradually close to the calculated results without regard to the free surface. The instantaneous axial load and energy utilization ratio of the impeller have a fluctuation. The more the tip-immersion, the smaller the wave amplitude of instantaneous value will be. Basically, the axial load and energy utilization ratio will not change with the azimuthal change of the impeller when the tip-immersion is 0.79D. Hence, the free surface has almost no effect on hydrodynamic characteristics of impeller. The results of the study can provide a basis and reference for the design and engineering application of a tidal power station.

tip-immersion; tidal current energy; horizontal axis impeller; free surface; hydrodynamic characteristics

2015-09-02.

日期：2016-08-29.

國家自然科學(xué)基金項目(51179077,51579120, 51309125).

王樹齊(1986-) ，男，講師，博士；

徐剛(1981-) ，男，副教授.

徐剛，E-mail：me_xug@qq.com

10.11990/jheu.201509007

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20160829.1421.050.html

TK730，O352

A

1006-7043(2016)10-1330-05

王樹齊，張亮，徐剛，等. 自由面條件下水平軸潮流能葉輪水動力研究[J]. 哈爾濱工程大學(xué)學(xué)報， 2016, 37(10): 1330-1334.

WANG Shuqi, ZHANG Liang, XU Gang， et al. Hydrodynamic analysis of a tidal current impeller in a horizontal axis under the condition of a free surface[J]. Journal of Harbin Engineering University, 2016, 37(10): 1330-1334.