孫 科，張 亮，何環(huán)宇

(哈爾濱工程大學 船舶工程學院，150001 哈爾濱)

?

三葉片垂直軸水輪機自啟動性能數值計算

孫 科，張 亮，何環(huán)宇

(哈爾濱工程大學 船舶工程學院，150001 哈爾濱)

為尋求固定偏角垂直軸水輪機的啟動性能規(guī)律，應用CFX流體-剛體耦合運動求解方法實現了垂直軸潮流能水輪機自啟動過程的瞬態(tài)數值模擬，結合靜態(tài)啟動力矩系數的計算結果，分析在不同初始方位角下水輪機的啟動時間以及角速度的穩(wěn)定情況.研究結果表明:對于三葉片固定偏角垂直軸水輪機而言，方位角在60°～130°范圍內啟動性能較好，但在90°位置靜態(tài)力矩會出現突降，不利于自啟動；方位角在140°～170°范圍內啟動性能較差，但系統(tǒng)阻尼為0時，在所有位置均能夠實現自啟動，且啟動過程穩(wěn)定后水輪機轉速均達到飛逸轉速.

潮流能；垂直軸水輪機；啟動位置；自啟動性能；方位角；系統(tǒng)阻尼

垂直軸潮流能水輪機是獲取潮流能的重要裝置[1-2],其自啟動性能是關系到水輪機年發(fā)電量的關鍵參數，也是評價發(fā)電機組性能優(yōu)劣的重要指標之一[3-4].對于垂直軸葉輪的自啟動性能研究最早始于風力發(fā)電機，1985年Rajaoferson等[5]通過實驗證明Darrieus型風力機具有一定的自啟動性能.Kirke[6-7]等詳細闡述了垂直軸風力機的自啟動性能，指出固定偏角垂直軸風力機一個最大的局限性就是自啟動性能差，并提出了幾種改進方式.隨著計算機的快速發(fā)展以及CFD軟件的日漸成熟，2010年開始,研究者們紛紛采用數值算法來研究垂直軸葉輪的啟動性能.例如Beri等[8]利用CFD軟件計算得到NACA2415翼型變弧固定葉片垂直軸風力機具有自啟動性能的結論；李鳳來等[9]通過加裝葉柵附體使潮流能垂直軸水輪機的自啟動性能得到進一步加強;Batista等[10-11]通過研究EN0005新型葉片和NACA系列不對稱翼型的流體性能使垂直軸風力機具有一定的自啟動性能；張新建等[12]研究了固定偏角垂直軸風力機的安裝半徑、葉片弦長和來流速度對自啟動性能的影響.以上數值研究都是基于穩(wěn)態(tài)研究，即在給定均勻流的狀態(tài)下，研究風(水)輪靜止不動的狀態(tài)下，葉片對主軸產生的靜力矩;2011年，張學偉等[13]利用Fluent軟件實現了潮流能水輪機自啟動過程的瞬態(tài)數值模擬；Untaroiu等[14]應用CFX CEL Expression語言實現了垂直軸風力機二維和三維的瞬態(tài)自啟動過程模擬，并對比了湍流度、時間步等參數設置對計算精度的影響，但上述研究者們都沒有分析葉輪的啟動位置(初始方位角)對于自啟動過程的影響.

本文著重研究垂直軸水輪機的初始方位角對于啟動性能的影響.首先分析三葉片固定偏角垂直軸水輪機處于不同方位角時葉片對主軸的穩(wěn)態(tài)啟動力矩；然后采用CFX軟件的rigid-body模塊實現了水輪機自啟動過程的剛體-流體耦合瞬態(tài)數值模擬；最后對應穩(wěn)態(tài)啟動力矩曲線，探討啟動位置對于瞬態(tài)啟動時間和轉速穩(wěn)定性的影響.

1 數值模型

1.1 模型參數及坐標系

本文采用李志川[15]研發(fā)的垂直軸水輪機試驗模型作為研究對象，模型參數見表1.水輪機坐標系定義如圖1所示，圖1中笛卡爾坐標系原點位于主軸中心，x軸方向平行于來流方向，y軸方向與x軸正交；θ為葉片的方位角，取x軸為起點逆時針旋轉為正；Q為葉片對主軸的合力矩.

表1 李志川[15]垂直軸水輪機模型參數

圖1 水輪機坐標系

1.2 控制方程

CFX軟件通過求解NS方程解決黏性流體問題，不可壓縮黏性流體的連續(xù)性方程和RANS方程組為：

水輪機葉片剛體運動控制方程組為：

(1)

(2)

(3)

(4)

(5)

(6)

(7)

1.3 網格、計算域和模型設置

采用CFX13.0全隱式耦合算法求解黏性流體控制方程，流體和剛體運動控制方程的耦合采用Rigid body模塊完成.計算域和網格如圖2所示，葉輪所在區(qū)域為旋轉域，外部區(qū)域為靜止域，采用結構化網格，內外域交界面采用滑移網格技術.流場入口設定為均勻流速1 m/s，5%湍流強度，上下邊界和出口設定為開放出口，絕對壓力為101 325 Pa，葉片表面設定為無滑移壁面.湍流模型為SST模型.穩(wěn)態(tài)的最大迭代步為200，物理時間為0.3 s；瞬態(tài)計算的時間步長為0.01 s.

教學模式指在一定的教學思想指導下，反映特定教學理論邏輯輪廓的，為保持某種教學任務的相對穩(wěn)定而具體的教學活動結構．目前，在全球范圍內的高等教育過程中采用的教學模式主要有2類4種[16-19]：一是以“教師”為中心的教學模式，包括LBL（Lecture-based learning）和CBL（Case-based learning）；二是以“學生”為中心的教學模式，包括PBL（Problem-based learning）和RBL（Resource-based learning）．

圖2 模型計算域及網格

2 模型驗證

由于李志川[15]水輪機模型沒有進行啟動性能的測試，為了驗證計算模型的精度，本文采用文獻[14]中英國Durham大學的垂直軸風力機的二維計算模型作為驗證算例(模型參數見表2).在保證計算域和邊界條件相同的基礎上，計算了啟動過程中轉速隨時間的變化曲線，并與文獻[14]中的二維結果進行了對比.

表2 Durham垂直軸風力機模型參數

圖3給出啟動過程的轉速變化對比曲線，可以看到，在初始階段自啟動角速度數值接近，在10 s之后數值出現偏差，兩條曲線的最大誤差為7.3%.產生誤差的原因主要有：

1)Untaroiu的網格總數為13.6萬，Y+=1.7 ；本文受限于計算條件，網格總數降到9.8萬，Y+=1.9.

2)文獻[14]采用的是CEL Expression語言來實現流體-剛體運動耦合，是一種顯式算法；而本文采用的是rigid body模型來實現流體-剛體運動耦合，是一種全隱式算法.

3)顯式算法離散格式的穩(wěn)定性取決于時間步和空間步的選取，文獻[14]對比了時間步對計算結果的影響，選取dt=0.001 s；由于全隱式算法無條件穩(wěn)定，為了節(jié)約計算成本，本文在保證收斂性的基礎上選取dt=0.01 s.

圖3 啟動過程中葉輪瞬時轉速隨時間的變化

整體來看，兩條曲線的變化規(guī)律一致，都是隨著時間的增加轉速逐漸爬升，0～10 s之內幾乎呈線性變化，10～15 s之間轉速爬升加速，超過15 s之后轉速爬升回緩，接近20 s時轉速接近最大值.這說明無論應用CFX軟件的rigid body模塊還是應用CEL Expression語言，在定性上討論瞬態(tài)自啟動過程的變化規(guī)律都是可行的.

3 結果和討論

1)利用NACA0018翼型升阻力系數的靜態(tài)試驗值，推導出葉片在不同方位角下的切向力，然后疊加3個葉片切向力總和從而計算出葉片對主軸的合力矩.

2) 采用穩(wěn)態(tài)模型，應用CFX軟件模擬水輪機葉片在不同初始位置時的流場，待計算收斂后統(tǒng)計3個葉片對主軸的合力矩.如圖1所示，由于3個葉片沿圓周對稱分布，本文取葉片1在θ= 60°～170°之間的位置進行分析，每間隔10°為一個工況.為了數據分析的方便，定義水輪機的速比和力矩系數如下：

式中：ρ為介質密度；b為葉片展長(二維計算展長取單位為1)；R為水輪機半徑.

圖4給出兩種計算方法得到的水輪機主軸的靜態(tài)力矩系數對比.可以看到水輪機力矩推導值相較于CFX的穩(wěn)態(tài)計算值偏大，這是因為用升阻力系數推導葉片受力時，假定葉片之間無干擾；但用CFX建立數值模型時，可以充分考慮上游葉片對下游葉片流場的影響，這種影響在大部分的位置角時會導致下游葉片的入流速度損失，因此CFX的計算值更符合真實的流場特性.整理來看，固定偏角垂直軸水輪機的靜態(tài)力矩值都很小，但在阻尼為零的假定下，只要力矩不出現零值或負值，根據式(1)～(7)，理論上就可以實現葉輪啟動過程的瞬態(tài)模擬.

圖4 不同方位角下水輪機的靜態(tài)力矩系數

圖5為水輪機在不同初始方位角時啟動過程中角速度(ω)隨時間(t)的變化曲線.圖5中曲線證實了公式分析的結論，在系統(tǒng)零阻尼的假定條件下，水輪機在各個初始方位角時都具有自啟動性能，且旋轉角速度最終都穩(wěn)定在6.5 rad/s左右.無論初始方位角在什么位置，曲線的變化趨勢相同，葉輪的角速度從線性爬升階段過渡到爬升速度加劇階段，經過速度回落階段最終達到穩(wěn)定值.由于初始方位角的不同，角速度的爬升速度和達到的最大值也有所不同.方位角在θ=70°時水輪機自啟動角速度達到最大值用時最短.在θ=60°、70°、80°、100°、110°、120°以及130°的位置，自啟動角速度的變化曲線相近，到達各自最大值的用時相對其他5種工況少，屬于優(yōu)選自啟動方位角范圍.在θ=90°、140°、150°、160°和170°的位置，自啟動角速度到達最大值用時相對較長，其中θ=140°的自啟動時間最長，這5個位置屬于次選自啟動方位角.

在t=13～24 s階段各角速度從最大值回落到穩(wěn)定值.θ=140°用時最長且角速度最大值時間節(jié)點排在最后，θ=70°用時最短且角速度最大值時間節(jié)點排在最前面.優(yōu)選自啟動方位角回穩(wěn)的用時相對次選自啟動方位角回穩(wěn)用時要小，會更快達到穩(wěn)定值.對應圖4的CFX靜態(tài)力矩曲線看出，140°是靜態(tài)力矩最小的位置，所以啟動速度最慢；150°、160°和170°位置的靜態(tài)啟動力矩均低于優(yōu)選啟動方位角的靜態(tài)力矩.因此可以判斷，水輪機的靜態(tài)啟動力矩是影響水輪機啟動性能的重要因素.在優(yōu)選自啟動方位角的范圍內，即θ=60°～130°之間，θ=90°位置的靜態(tài)力矩出現了一個突降的階躍.如圖6所示，此時葉片位置的分布比較特殊，葉片3完全處于葉片2的尾流區(qū)域內，葉片2脫落的尾渦對于葉片3近前方的合速度產生很大的影響，導致葉片3對啟動力矩做出零貢獻甚至是負貢獻，這是在初始啟動位置中應避免的一個特殊方位角.

圖5 啟動過程中角速度隨時間的變化曲線(θ=60°～170°)

從圖5中可以看出，任何初始方位角下啟動，葉輪的角速度最終都穩(wěn)定在6.5 rad/s左右.在來流速度不變的情況下，水輪機以角速度ω=6.5 rad/s做均速旋轉，計算結果顯示平均力矩系數Cq為0，證明了ω=6.5 rad/s為該流速下葉輪的飛逸轉速.圖7給出水輪機以恒定轉速ω=6.5 rad/s做均速旋轉時的力矩系數與自啟動狀態(tài)下轉速達到6.5 rad/s時力矩系數的對比，兩條曲線基本吻合，恒定轉速計算的曲線波動比自啟動瞬態(tài)計算的波動略高，均值都為0.由式(6)可以推斷，當無阻尼狀態(tài)下葉輪自由旋轉時，轉速達到穩(wěn)定即為葉輪的角加速度為0，此時葉片對主軸的合力矩均值必為0.也就是說，無阻尼狀態(tài)下水輪機自啟動過程達到的角速度穩(wěn)定值必為水輪機的飛逸轉速值.以往在進行轉速恒定假設數值模擬的時候，需要不斷改變速比點來尋求飛逸轉速的位置，應用本文的自啟動數值模擬方法可以估算出在不同流速下水輪機模型的飛逸轉速.

圖6 穩(wěn)態(tài)流場速度分布圖 (θ=90°)

圖7 力矩系數隨方位角變化曲線

4 結 論

1)應用CFX的rigid body模塊，可實現固定偏角垂直軸水輪機的剛體-流體耦合瞬態(tài)自啟動過程的數值模擬，應用該方法可以估算出在不同流速下水輪機的飛逸轉速.

2)三葉片固定偏角垂直軸水輪機的靜態(tài)啟動力矩在所有方位角下均大于零，因此在系統(tǒng)零阻尼的條件下，水輪機在所有的初始方位角下都具有自啟動性能.

3)方位角θ=60°～130°屬于優(yōu)選自啟動范圍，其中θ=90°由于流場的特殊性，自啟動性能較差；方位角θ=140°～170°屬于次選自啟動范圍角，其中θ=140°為最差自啟動方位角，靜態(tài)啟動力矩最小，自啟動用時最長，應盡量避免在啟動時刻葉片處于該位置.

[1] 李志川, 張亮, 孫科,等.垂直軸潮流水輪機數值模擬研究[J].太陽能學報,2011,32(9):1321-1326.

[2] AJEDEGBA J O, NATERER G F, ROSEN M A, et al．Effects of blade configurations on flow distribution and power output of a zephyr vertical axis wind turbine[C] //Proceedings of the 3rdIASME/WSEAS International Conference on Energy & Environment. Stevens Point, Wisconsin: World Scientific and Engineering Academy and Society, 2008: 480-486．

[3] 王世澎.波浪對潮流影響的數值模擬研究[D].大連:大連理工大學,2007：20-22.

[4] HE Huanyu, SUN Ke, MA Yong,et al. Self-starting performance numerical analysis of fixed-pitch vertical axis hydro-turbine[C]// Proceedings of the 2ndInternational Conference on Energy Engineering and Environment Engineering. Pfaffikon: Trans Tech Publications, 2014:102-105.

[5] RAJAOFERSON S, COMOLET R. Self-starting process of the Darrius wind-turbine[C] //European Wind Energy Conference. [S.l.]:[s.n.], 1985:523-528.

[6] KIRKE B K.Evaluation of self-starting vertical axis wind turbines for seand alone applications[D].Australia: School of Engneering,Griffith University,1998:41-50.

[7] KIRKE B K,LAZAUSKAS L. Limitations of fixed pitch Darrieus hydrokinetic turbines and the challenge of variable pith[J]. Renewable Energy,2011,36(3):893-897.

[8] BERI H, YAO Yingxue. Effect of camber airfoil on self-starting of vertical axis wind turbine[J]. Journal of Enviromental Science and Technology, 2011,4(3): 302-312.

[9] 李鳳來,姜德志，孟巧.豎軸固定直葉片潮流水輪機自啟動性能研究[J].工程與實驗, 2011, 51(3): 23-26.

[10]BATISTA N C , MELICIO R. MATICS J C O, et al.New blade profile for darrieus wind turbines capable to self-start[C]//IET Conference on Renewable Power Generation, RPG 2011. Edinburgh: Institution of Engineering and Technology, 2011:1-5.

[11]BATISTA N C, MELICIO R, MATICS J C O, et al. Self-start evaluation in lift-type vertical axis wind turbines: methodology and Computational tool applied to asymmetrical airfoils[C]// International Conference on Power Engineering, Energy and Electrical Drives (POWERENG). Malaga: IEEE Computer Society, 2011: 1-6.

[12]張新建,蔡新,潘盼.H型垂直軸風力機啟動性能分析[J]. 水電能源科學, 2013, 31(5): 243-246.

[13]張學偉,張亮, 李志川, 等. 垂直軸自由變偏角水輪機的多體耦合數值求解方法[J]. 計算力學學報, 2013, 30(3): 370-375.

[14]UNTAROIU A, WOOD H G, ALLAIRE P E, et al. Investigation of self-starting capability of vertical axis wind turbines using a computational fluid dynamics approach[J]. Journal of Solar Energy Engineering, 2011, 133(4): 112-120.

[15]李志川.垂直軸潮流能水輪機水動力特性數值模擬與實驗研究[D].哈爾濱：哈爾濱工程大學，2011: 35-36.

(編輯 張 紅)

Numerical computation on self-starting position of a three-blade vertical axis hydro-turbine

SUN Ke，ZHANG Liang，HE Huanyu

(College of Shipbuilding Engineering, Harbin Engineering University,150001 Harbin, China)

For studying on start-up performance rules of fixed pitch tidal current turbine, the transient self-starting process of vertical axis hydro-turbine is numerical simulated by using CFX fluid-rigid body coupled motion calculating method. Combining static start-up torque coefficient results, the self-starting time and the angular speed stability of the turbine are analyzed in different initial azimuth angles. It is shown that, as for the three-blade fixed-pitch vertical axis hydro-turbine, the turbine has a better self-starting performance at azimuth angles of 60°～130°, but the angle of 90°is not a good choice because the static start-up torque will drop suddenly at this point. The turbine has a worse self-starting performance at azimuth angles of 140°～170°. However, the three-blade vertical axis hydro-turbine has self-starting performance at any initial position when the system damping is zero and the angular speed of the main shaft can arrive at runaway speed finally.

tidal current energy; vertical axis hydro-turbine; start-up position; self-starting performance; azimuth angle; system damping

10.11918/j.issn.0367-6234.2016.04.024

2014-10-21.

國家自然科學基金(51209060); 高等學校博士學科點專項科研基金(20122304120035).

孫 科(1978—)，女，講師，博士;

張 亮(1959—)，男，教授，博士生導師.

張 亮, heusunke@126.com.

O352

A

0367-6234(2016)04-0144-05