«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2019, Vol. 40 Issue (8): 1387-1392  DOI: 10.11990/jheu.201805056
0

引用本文  

马勇, 胡超, 郑雄波, 等. Π型竖轴双转子水轮机纵荡水动力特性研究[J]. 哈尔滨工程大学学报, 2019, 40(8), 1387-1392. DOI: 10.11990/jheu.201805056.
MA Yong, HU Chao, ZHENG Xiongbo, et al. Hydrodynamic characteristics of Π-type vertical-axis parallel dual-rotor turbine in surge motion[J]. Journal of Harbin Engineering University, 2019, 40(8), 1387-1392. DOI: 10.11990/jheu.201805056.

基金项目

国家自然科学基金项目(51779062,51679053);中央高校基本科研业务费项目(HEUCFP201714);黑龙江省自然科学基金项目(A2016001)

通信作者

郑雄波, E-mail:zhengxiongbo@hrbeu.edu.cn

作者简介

马勇, 男, 副教授, 博士生导师;
胡超, 男, 博士研究生;
郑雄波, 男, 讲师

文章历史

收稿日期:2018-05-14
网络出版日期:2019-04-15
Π型竖轴双转子水轮机纵荡水动力特性研究
马勇 1,2, 胡超 1, 郑雄波 1, 邓锐 1,2     
1. 哈尔滨工程大学 深海工程技术研究中心, 黑龙江 哈尔滨 150001;
2. 中山大学 海洋工程与技术学院, 广东 珠海 519082
摘要:为提高潮流能装置发电功率及可靠性,加快潮流能开发利用,本文提出了一种Π型竖轴双转子水轮机的设计,并对Π型竖轴双转子水轮机的纵荡运动特性进行了研究。基于CFD方法,通过将水轮机旋转运动与振荡运动相结合的方法模拟水轮机的纵荡运动,对竖轴双转子水轮机在不同振荡频率、振荡幅值及叶尖速比下的纵荡运动进行模拟计算,同时与均匀流无振荡工况下的水轮机作对比,分别对振荡运动下双转子水轮机受力、获能效率及流场的变化特征进行分析。研究结果表明:双转子并列布置可提高水轮机获能效率、纵荡幅值及频率对平均获能效率影响较小。纵荡下双转子水轮机所受推力及侧向力的波动较大,且推力及侧向力波动幅值与振幅及振荡频率成正比。振荡频率影响泄涡频率,振荡幅值影响尾涡区长度。
关键词潮流能    竖轴双转子水轮机    水动力特性    强迫振荡    CFD方法    
Hydrodynamic characteristics of Π-type vertical-axis parallel dual-rotor turbine in surge motion
MA Yong 1,2, HU Chao 1, ZHENG Xiongbo 1, DENG Rui 1,2     
1. Deepwater Engineering and Technology Research Center, Harbin Engineering University, Harbin 150001, China;
2. School of Marine Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, China
Abstract: To improve the power output and reliability of the tidal current device and accelerate the use of tidal current energy, in this paper, we present our design of a Π kind of vertical-axis dual rotor, and we investigate the hydrodynamics of the surge motion. Using computational fluid dynamics, we combine the rotary and oscillation motions of the turbine to simulate the turbine oscillation movement. We simulate and perform calculations of the vertical-axis dual-rotor turbine operating at different oscillation frequencies, oscillation amplitudes, and tip speed ratios, and compare our results with a turbine working at uniform flow without oscillation. We then analyze the variation characteristics of the force, energy efficiency, and flow field of the dual-rotor turbine under oscillating motion. The results shows that a parallel dual-rotor arrangement can improve the power output efficiency, but the oscillation frequency and amplitude have less influence on the average power output. The fluctuation amplitudes of the thrust and lateral force of the dual-rotor turbine under oscillation movement are greater and proportional to the oscillation amplitude and frequency. The oscillation frequency affects the vortex discharge frequency and the oscillation amplitude affects the length of the tail vortex region.
Keywords: tidal power    vertical-axis dual rotors turbine    hydrodynamic characteristics    forced oscillation    computational fluid dynamics(CFD)    

潮流能作为一种清洁可再生能源,因其具有储量大、分布广、可预测性强的优点,近年来成为国内外研究的热点[1-3]。竖轴水轮机作为潮流能水轮机中一种重要的水轮机形式,具有设计简单、无需换向机构、密封要求低、噪音小等特点[4-5]。为了提高潮流能装置发电功率,提高经济效益比,采用单装置双机组布置可行性较强。漂浮式载体形式是竖轴水轮机重要的载体形式之一,具有甲板面积大,方便布置运维等优点[6]。本文提出了一种漂浮式载体搭载Π型竖轴双转子水轮机的型式。该型式水轮机结构简单,生产制造成本较低,可全水深工作。连接叶片和主轴的轮辐位于水面以上,降低了轮辐能量耗散,进一步提高了水轮机的获能效率。同时叶片的水下部分没有轮辐和其他连接构件,解决了渔网和海草缠绕问题,便于叶片安装、维护、检修及海生物清除等海上作业。除叶片外,所有构件位于水面以上,降低了水密难度和水轮机故障率,便于安装维护和运行管理。

为了更加准确地预报双转子水轮机水动力性能,近年来研究者针对均匀流下双转子水轮机的相位角[7]、相对旋向[8]、轴间距[9]、排布方式[10]等进行了研究,研究成果表明这些参数的改变,会对水轮机水动力性能产生有利影响或负面干扰。实际上,潮流电站在海洋环境中工作时,会由于波浪的作用而产生周期性运动,这对双转子水轮机水动力特性会产生较大影响[11-13]。当竖轴水轮机叶片展弦比大于7时,可不考虑水轮机三维效应的影响。本文基于粘性CFD方法,通过将转子的旋转运动与转子的振荡运动相结合来模拟双转子水轮机在海洋环境下的工作情况。本文不考虑轮辐干扰及三维效应的影响,将转子纵荡的三维问题简化为二维问题,建立了Π型竖轴双转子强迫振荡下的水动力性能模拟方法,研究了不同叶尖速比、振荡幅值、振荡频率对双转子水动力性能的影响。

1 双转子无振荡数值模拟

竖轴双转子水轮机结构形式如图 1所示。

Download:
图 1 Π型双转子结构 Fig. 1 Structure of Π type dual rotor

在展开数值模拟时,定义$\overline{\mu}=\left(\int_{0}^{T} \mu \mathrm{d} t\right) / T $代表叶轮载荷均值,$ \hat{\mu}=\left(\mu_{\max }-\mu_{\min }\right)$代表载荷波动幅值。

水轮机在实际工作过程当中,影响其性能的因素很多,本文着重考虑主要因素,忽略次要因素,有利于理清研究方向,更进一步优化水轮机性能。本文规定以下无因次量:

$ \lambda=\omega R / V_{A} $ (1)
$ C_{\overline{F}_{X}}=\overline{F}_{X} /\left(0.5 \rho V_{A}^{2} D b\right) $ (2)
$ C_{\overline{F}_{Y}}=\overline{F}_{Y} /\left(0.5 \rho V_{A}^{2} D b\right) $ (3)
$ C_{P}=Q_{\omega} /\left(0.5 \rho V_{A}^{3} D b\right) $ (4)

2个水轮机旋转中心的间距为12 m(1.5DD为转子直径),并采用内对转方式以抵消侧向力的影响,相比其他旋转方式,内对转的形式也具有最高的获能效率。模拟时为减少受力波动,便于分析,将2个水轮机的初相位均设为0°。为提高网格质量,计算域包括旋转域和静止域,并规定坐标系X方向与来流一致,如图 2所示。

Download:
图 2 双转子无振荡模拟 Fig. 2 Simulation of dual rotor without oscillation

给定来流Ve=3 m/s,设定2个水轮机的角速度分别为ω和-ω,计算时叶尖速比λ分别为1.0、1.5、2.0、2.5、3.0、3.5、4.0。定义CP为双转子的能量利用率,计算方法为2个单转子获能效率的算术平均值。

提取不同速比下的CP,并与单转子进行对比,如图 3所示。可以看出,当两轮心间距为12 m(1.5D),低速比时, 双转子布置对效率提升效果不明显,随着速比逐渐趋近最优速比,双转子获能效率较单转子有着明显提升,当λ=2.0时,此时的双转子比单转子在平均能量利用率上提升了15.3%,因此合理的双转子布置能有效提高叶轮获能效率。

Download:
图 3 单/双转子获能效率对比 Fig. 3 Comparison of power output of single and dual rotor

λ=2.0时,考察水轮机在实现最大获能效率时的受力情况,得到无因次化的推力、转矩系数,如图 4所示。由图可以看出,双转子水轮机受到的XY方向的合力均呈周期性变化。X方向推力合力的峰值波动幅值很小,峰值包络线可近似为一条直线。Y方向侧向力合力历时曲线峰值波动很大,但侧向力均值$\overline{C}_{F r Y} $ =0.006 8,双转子内对转布置可以起到抵消侧向力的作用。经过计算可以发现,侧向力CFrY峰值包络线变化周期与转子周期保持一致。

Download:
图 4 双转子无振荡时推力系数和侧向力系数变化曲线 Fig. 4 Curves of thrust and lateral force coefficient without oscillation
2 双转子振荡数值模拟

双转子水轮机的振荡运动CFD模拟所确定计算域如图 5所示,规定X正方向为来流方向,包括静止域、旋转域及振荡域,三者之间通过交界面连接。其中,旋转域含有2个尺度相同且旋向相反的水轮机;振荡域为沿X方向的正弦运动,表达式为ξ=ξ0sin(ωt),ζ0为振荡幅值,ω为振荡圆频率。

Download:
图 5 双转子振荡模拟 Fig. 5 Simulation of dual rotor with oscillation

CFD模拟时设定双转子水轮机的内对转,还要实现双转子水轮机的振荡运动。采用CEL语言控制水轮机的旋向,以及动网格和滑移网格相结合的网格运动形式,网格如图 6所示。振荡域和旋转域同时受到X方向的简谐运动ξ=ξ0sin(ωt),采用动网格形式。旋转域在振荡域的内部,作为子域进行旋转运动,旋转域与振荡域通过滑移网格进行连接,以此来同步实现水轮机叶轮振荡与水轮机叶轮旋转,同时减小网格变形,保证模拟精度。

Download:
图 6 双转子振荡模拟全流场网格 Fig. 6 Simulation mesh of dual rotor with oscillation
2.1 叶尖速比的影响

首先研究不同叶尖速比对于叶轮获能效率的影响。分别计算当叶尖速比λ为1.0、1.5、2.0、2.5、3.0、3.5等7个工况下,2个转子的平均获能效率,计算结果如图 7所示。从图中可以看出,获能效率随叶尖速比的变化趋势同无振荡运动工况下相比并未发生改变,最佳速比处于λ=2.0左右。

Download:
图 7 CP随叶尖速比的变化趋势 Fig. 7 Variation trend of CP vs. tip speed ratio

然后研究振荡运动下,不同叶尖速比下的受力情况,计算结果如图 8所示。可以看出,X方向合力系数的CFrX$\hat{C}_{F r X} $与叶尖速比λ成正相关关系,均随着叶尖速比的提高而提高。Y方向合力系数$\overline{C}_{F r Y} $并未随着叶尖速比λ的变化而发生改变,且不同速比下λ$\overline{C}_{F r Y} $值均围绕在0进行微弱波动,说明振荡运动下,速比并没有对Y方向合力均值发生影响,还是处于一个相互抵消的状态;而Y方向合力系数$\hat{C}_{F r Y} $与叶尖速比λ的变化成正相关,随着速比的升高而缓慢升高。

Download:
图 8 推力及侧向力系数变化幅值及均值随速比变化曲线 Fig. 8 Curves of fluctuation amplitude and mean value of thrust and lateral force coefficient along with λ
2.2 振荡幅值的影响

为分析不同振荡幅值对双转子水轮机的受力、获能效率及流场的影响,选取不同的振荡幅值ξ0为0.4、0.8、1.0、1.2 m。在Ve=3 m/s,λ=2,振荡圆频率ω=1.0 rad/s,振荡位移ξ=ξ0sin(ωt)下进行CFX模拟。

图 9可知,仅改变振荡幅值时,水轮机振荡运动和旋转运动之间相位关系保持不变,仅推力及侧向力幅值发生改变,且推力及侧向力波动幅值与振荡幅值成正比,但推力均值并未发生改变。相比于非振荡条件,振荡条件下的水轮机推力及侧向力波动更大,双转子布置方式对主轴疲劳特性提出了更高的要求。

Download:
图 9 不同振荡幅值时推力及侧向力系数历时曲线 Fig. 9 Curves of thrust and lateral force coefficient at different amplitude along with time

图 10为获能效率随振荡幅值变化曲线。可以看出,随着纵荡幅值的增加,叶轮获能效率并未发生较大改变,纵荡运动对水轮机获能效率并未产生较大影响。这是因为纵荡时叶轮盘面相对速度按照振荡频率周期性变化,因此仅造成获能效率周期性波动,但对均值无影响。

Download:
图 10 不同纵荡幅值时CP变化曲线 Fig. 10 Curve of CP at different surge amplitude

不同振荡幅值下的双转子水轮机的流场变化如图 11所示。中间射流的存在会把中间流道内原本连续的涡给冲散,因此随着振荡幅值的提升,中间的涡列比外侧2排涡列连续性更差且都是间断的。更多的是集中涡而不是连续涡,此时对叶片受力的脉动程度更强,因此对其疲劳程度要求也越高。

Download:
图 11 尾流场随振荡幅值的变化趋势 Fig. 11 Variation trend between wake flow field and oscillation amplitude
2.3 振荡频率的影响

为分析不同振荡频率对双转子水轮机的受力、获能效率及流场的影响,选取不同的振荡频率ωo为0.4、0.8、1.0、1.2 rad/s。在Ve=3 m/s,λ=2,振荡幅值ξ=1.0 m,振荡规律为ξ=ξsin(ωot)下进行CFX模拟。

图 12为在3个不同振荡频率下的推力及侧向力波动历时曲线。当振荡频率发生改变时,旋转和振荡相位之间相位关系应发生改变。但在0.4、0.8、1.2 rad/s时,推力及侧向力历时曲线波动趋势和相位特性几乎一致。

Download:
图 12 不同频率时侧向力及推力系数历时曲线 Fig. 12 Curves of thrust and lateral force coefficient along with time

双转子水轮机的受力均值和波动值随频率发生改变趋势如图 13所示。在均值方面,不同振荡频率并没有对$\overline{C}_{F r X} $$\overline{C}_{F r Y} $产生明显的影响,$\overline{C}_{F r X} $$\overline{C}_{F r Y} $均在围绕在某个值微幅波动。在波动幅值方面,$\hat{C}_{F r X} $$\hat{C}_{F r Y} $均随着振荡频率的提高而升高。

Download:
图 13 推力和侧向力系数波动幅值及均值随频率变化关系 Fig. 13 Curves of fluctuation amplitude and mean value of thrust and lateral force coefficient along with frequency

图 14给出了不同纵荡频率下叶轮获能效率曲线,0 rad/s表示无振荡。从图中可以看出纵荡频率对叶轮效率影响较小。同理,不论是是纵荡幅值还是纵荡频率发生改变,叶轮盘面相对速度的改变可以看成是叶轮盘面均匀来流速度Ve和振荡速度u=ξωcos(ωt)的叠加,对叶轮平均效率没有较大影响。

Download:
图 14 获能效率随频率变化关系 Fig. 14 Curve of CP along with frequency

最后研究不同振荡频率下的双转子水轮机的流场变化特点,将本节所涉及到的3个工况下的尾涡变化,按由小到大列举如图 15所示。随着振荡频率的增加,影响了流场的泄涡频率,泄涡频率也会增加,涡的分布更加分散,但是振荡频率的变化并没有明显影响尾流区长度。可见振荡频率更多地影响泄涡频率,而振荡幅值更多地影响尾流区长度和速度补偿的快慢。

Download:
图 15 尾流场随振荡频率的变化趋势 Fig. 15 Variation trend between wake flow field and oscillation frequency
3 结论

1) 双转子水轮机布置形式可有效提高水轮机获能效率,相比于单转子水轮机获能效率提升12.6%左右。

2) 纵荡运动工况下,双转子水轮机获能效率变化趋势与无振荡情况下变化趋势相同。纵荡幅值及纵荡频率对水轮机平均获能效率影响较小,相比于无振荡情况,双转子推力及侧向力波动更为剧烈,这对叶轮结构强度及叶片疲劳寿命会产生较大影响。推力及侧向力波动幅值同纵荡频率及幅值成正比。

3) 纵荡情形下,双转子尾流区更多呈现的是集中涡,叶片及主轴受力脉动性更强,对结构疲劳要求更高。

4) 纵荡幅值主要是影响尾流区长度及尾流区速度补偿快慢,纵荡频率主要是影响泄涡频率,涡的分布更加分散,对尾流区长度影响较小。

参考文献
[1]
张亮, 李新仲, 耿敬, 等. 潮流能研究现状2013[J]. 新能源进展, 2013, 1(1): 53-68.
ZHANG Liang, LI Xinzhong, GENG Jing, et al. Tidal current energy update 2013[J]. Advances in new and renewable energy, 2013, 1(1): 53-68. (0)
[2]
LIU Hongwei, MA Shun, LI Wei, et al. A review on the development of tidal current energy in China[J]. Renewable and sustainable energy reviews, 2011, 15(2): 1141-1146. DOI:10.1016/j.rser.2010.11.042 (0)
[3]
MA Yong, LI Tengfei, ZHANG Liang, et al. Experimental study on hydrodynamic characteristics of vertical-axis floating tidal current energy power generation device[J]. China ocean engineering, 2016, 30(5): 749-762. DOI:10.1007/s13344-016-1001-y (0)
[4]
汪鲁兵, 张亮, 曾念东. 一种竖轴潮流发电水轮机性能优化方法的初步研究[J]. 哈尔滨工程大学学报, 2004, 25(4): 417-422.
WANG Lubing, ZHANG Liang, ZENG Niandong. An optimization method for improving hydrodynamic performance of the vertical-axis turbine for tidal streams energy conversion[J]. Journal of Harbin Engineering University, 2004, 25(4): 417-422. DOI:10.3969/j.issn.1006-7043.2004.04.004 (0)
[5]
李志川, 张亮, 孙科, 等. 垂直轴潮流水轮机数值模拟研究[J]. 太阳能学报, 2011, 32(9): 1321-1326.
LI Zhichuan, ZHANG Liang, SUN Ke, et al. Numerical simulation of vertical axis tidal turbine[J]. Acta energiae solaris sinica, 2011, 32(9): 1321-1326. (0)
[6]
SHENG Qihu, JING Fengmei, ZHANG Liang, et al. Study of the hydrodynamic derivatives of vertical-axis tidal current turbines in surge motion[J]. Renewable energy, 2016, 96: 366-376. DOI:10.1016/j.renene.2016.04.074 (0)
[7]
王凯, 孙科, 张亮, 等. 相位角对双机组立轴水轮机的水动力性能影响[J]. 哈尔滨工程大学学报, 2016, 37(1): 104-109.
WANG Kai, SUN Ke, ZHANG Liang, et al. Impact of phase angle on the hydrodynamic performance of a bi-unit vertical axis tidal current energy water turbine[J]. Journal of Harbin Engineering University, 2016, 37(1): 104-109. (0)
[8]
LI Ye, CALISAL S M. Modeling of twin-turbine systems with vertical axis tidal current turbine:Part Ⅱ-torque fluctuation[J]. Ocean engineering, 2011, 38(4): 550-558. DOI:10.1016/j.oceaneng.2010.11.025 (0)
[9]
孙科, 李岩, 王凯, 等. 串列竖轴水轮机尾流场影响CFD模拟分析[J]. 哈尔滨工业大学学报, 2018, 50(5): 185-191.
SUN Ke, LI Yan, WANG Kai, et al. CFD simulation analysis on the wake effect of tandem vertical axis tidal turbines[J]. Journal of Harbin Institute of Technology, 2018, 50(5): 185-191. (0)
[10]
BAI Guanghui, LI Wei, CHANG Hao, et al. The effect of tidal current directions on the optimal design and hydrodynamic performance of a three-turbine system[J]. Renewable Energy, 2016, 94: 48-54. DOI:10.1016/j.renene.2016.03.009 (0)
[11]
张亮, 王树齐, 马勇, 等. 潮流能水平轴叶轮纵摇运动水动力分析[J]. 哈尔滨工程大学学报, 2015, 36(3): 307-311.
ZHANG Liang, WANG Shuqi, MA Yong, et al. The pitch hydrodynamic analysis of tidal current energy horizontal axis impeller[J]. Journal of Harbin Engineering University, 2015, 36(3): 307-311. (0)
[12]
SUFIAN S F, LI Ming, O'CONNOR B A. 3D modelling of impacts from waves on tidal turbine wake characteristics and energy output[J]. Renewable energy, 2017, 114: 308-322. DOI:10.1016/j.renene.2017.04.030 (0)
[13]
胡健, 李聪慧, 张维鹏, 刘立超.基于CFD的桨舵水动力干扰研究[J].应用科技, 2017, 44(3): 5-11, 18.
HU Jian, LI Conghui, ZHANG Weipeng, et al. Investigation of hydrodynamic interaction between propeller and rudder based on computational fluid dynamics[J]. Applied science and technology, 2017, 44(3): 5-11, 18. (0)