2017, Vol. 33
可再生能源的开发和利用日益受到人们越来越多的重视, 其中海洋能因经济性、储量多以及环保性等优点成为国内外研究的热点.根据轴线与来流方向可将水轮机分为水平轴式叶轮和垂直轴式叶轮.水平轴式叶轮发电效率较高, 垂直轴式叶轮可以吸收任何方向的水流且结构简单.目前比较多的研究机构主要研究轮机水动力学性能, Guerri等[1]仿真了二维竖直轴风力机的叶片受力情况, 采用的湍流模型为k-ω模型.研究发现, 在较低速比时, 叶片附近会产生很强的涡流, 模拟叶片的受力将有助于研究流体与风机结构间的相互影响.Myers等[2]针对水平轴潮流能水轮机进行了多次实验分析和理论研究, 分析了叶片安装角对水轮机水动力学性能的影响; Saha等[3]通过水洞实验对Savonius叶轮的结构参数进行了优化设计.Mohamed等[4]研究了导流罩对Savonius叶轮发电性能的影响.黄鹏等[5]对比了相同叶尖速比、不同叶片半径的风力机气动特性.本文设计了一种卧式水轮机, 它既有水平轴式水轮机发电效率高的优点, 又有垂直轴式水轮机可以捕获多方向水流的优点.本文重点研究叶片展长对水轮机水动力学性能的影响, 分析不同叶片展长的水轮机叶片对其转矩系数和获能系数的影响, 为实际的设计工作提供指导.
1 卧式水轮机机械结构和工作原理 1.1 卧式水轮机机械结构图 1为卧式水轮机整体结构.图 1(a)为卧式水轮机系统, 该系统由安全防护沉浮系统保护并漂浮于海面上, 通过锚链固定.由波浪能集成能量推动卧式水轮机转动, 完成由海洋能向旋转的机械能的转化, 并通过增速传动机构传递给旁侧发电机, 实现电能转化.于此同时, 安装于装置顶部并露出水面的远程监控系统将对装置发电数据进行实时保存和传输, 岸上用户可直观监控.卧式水轮机的关键是卧式水轮机的设计, 水轮机样机如图 1(b)所示.参考风力机设计理论[6-7], 采用5片叶片NACA63-018对称翼型结构, 结合浪流在变速运动下最优叶片设计理论, 对叶片弦长、展长等几何尺寸进行初步设计, 具体参数如表 1所示, 叶片界面为对称截面.本文主要研究叶片展长b对水轮机水动力学性能的影响.
|
图 1 卧式水轮机整体结构 Fig.1 Structure of the horizontal hydro turbine |
|
表 1 水槽实验样机模型参数 Table 1 Parameters of prototyping model for the flume experiment |
图 2为水轮机叶片受力分析.由图可知:水流以速度v向右流动, 水轮机叶片则以角速度ωo绕坐标轴原点旋转, t为时间, 定义图示位置为零时刻; 水流对叶片有升力Fl和阻力Fd作用, 定义叶片所受的阻力Fd沿x轴正方向, 升力Fl沿y轴正方向, 将它们分别沿着叶片的切线方向和法线方向进行分解, 即可得到单个叶片的切向力和法向力; 切向力促使水轮机进行旋转运动, 法向力则组成水轮机的径向载荷.
|
图 2 水轮机叶片受力分析 Fig.2 Force analysis of the turbine blade |
为了计算方便, 定义如下无量纲量.
尖速比
|
(1) |
式中:R为水轮机的半径; v为来流速度.
叶片升力系数
|
(2) |
式中:ρ为水流密度; A为参考面积.
叶片阻力系数
|
(3) |
水轮机所有叶片总的瞬时扭矩系数
|
(4) |
|
(5) |
式中:d为水轮机转子直径; M为水轮机所有叶片受到的总扭矩; Ft为推力; Ct为推力系数; h为叶片的弦长.
因为M随叶片的相位角α改变, 故Ctm也随α改变.水轮机在达到稳定状态时, 叶片受力呈周期性变化, 周期为2π.计算在一个周期内所有叶片总的平均扭矩系数Cm, 即
|
(6) |
同样, 定义水轮机瞬时的能量利用率
|
(7) |
式中, m为水轮机的质量.
定义水轮机的能量利用率
|
(8) |
实验在国家海洋局东海标准计量中心的水槽中进行, 水槽尺寸为200 m×4 m×3 m, 水深为2 m, 可产生的速度范围为0.5~3 m·s-1.水槽安装图如图 3所示.实验水轮机布置于槽体下方, 并保证水轮机轴平行于槽底, 水轮机主轴连接JN-338转矩转速测量仪和WT-230功率仪等测量仪器.当行车以不同的速度运动时, 带动水轮机主轴转动, 从而可以得到不同转速下的功率.根据式 (7) 可以得出能量利用率Cp.实验中通过造波器改变来流速度, 制作不同展长的样机模型.浪流工况实验方案如表 2所示.
|
图 3 水槽安装图 Fig.3 Installation of the flume |
|
|
表 2 浪流工况实验方案 Table 2 Experimental scheme for the working conditions of wave and current |
水轮机转矩为其水动力学性能至关重要的一个参数, 其测量值的精确性显著影响卧式轮机获能系数.参考文献[8-9], 测量时选择基于精确扭矩测量路线.该路线主要基于电阻应变片在外界力的作用下产生机械变形从而引起电阻变化的原理.所需实验设备有电阻应变片、YD-28A型动态电阻应变仪、信号采集仪以及计算机.实验过程中, 电阻应变片被牢固地粘贴在卧式轮机主轴顶端测量点上.由于水轮机主轴转动过程中产生机械扭转变形引起电阻的变化, 通过信号采集仪采集相应的电压变化, 并传递给计算机.经电阻与扭矩Q的前期标定公式得出扭矩.JN-338转矩测量图如图 4所示.
|
图 4 JN-338转矩测量图 Fig.4 Torque measurement by JN-338 |
实验中发电机功率的测量设备包括WT 230功率仪、整流器、电子负载以及与其配套的软件.电子负载及其配套监控软件界面如图 5所示.发电机与功率仪采用Y-Y三相三线连接, 通过整流器后接入由美尔诺电子公司设计的直流电子负载系统.该系统实现两个功能, 其一是通过WT 230功率仪测量并保存发电机功率; 其二是通过电子负载调节系统负载电阻, 以改变卧式轮机的转速, 保证其在所需标准叶尖速比工况条件下运行.
|
图 5 电子负载及其配套监控软件界面 Fig.5 Electronic load and the interface of its monitoring software |
本文采用Fluent软件进行水动力学数值计算.首先对该软件分析的理论基础进行概述.
3.1 理论基础 3.1.1 基本守恒定律流体流动需遵守的基本守恒定律包括质量守恒定律、动量守恒定律、能量守恒定律.若流动为湍流, 则必须同时遵守湍流输运方程.
质量守恒方程
|
(9) |
式中:u为速度矢量; p为压力.
动量守恒方程
|
(10) |
式中:τ为微元表面黏性应力; f为作用在微团上单位质量流体的质量力.
能量守恒方程
|
(11) |
式中:e为内能; q为交换热量; φ为广义变量.
3.1.2 SST k-ω湍流模型SST k-ω湍流模型综合了边界层外部模型的独立性和近壁面模型的稳定性的优点.该模型可以适用于速度压力梯度变化的各种物理现象, 模拟计算性能优良.因此, SST k-ω湍流模型对处理本文中雷诺数不同的问题时具有良好的通用性.
湍流强度方程
|
(12) |
式中:Gk为由平均速度梯度引起的湍动能产生项; Yk为湍流引起的k的耗散; Γk为k的有效扩散项.
湍流频率方程
|
(13) |
式中:Yω为由湍流引起的ω的耗散; Dω为交叉扩散项; Gω为ω产生项; Γω为ω的有效扩散项.
3.1.3 滑移网格模型运用滑移网格模型计算旋转区域可以达到很高的精度.通过求解两个或多个计算区域, 每个计算区域之间沿区域分界面移动, 网格分界面可以进行两个或多个区域的数据传递.在计算过程中, 每个计算区域沿网格分界面相互平移, 但是不会像动网格那样发生网格的变形.为了计算交界面的数据交换, 在新的时间步长时要确定交界面的重合面.理论上, 网格交界面的通量是根据交界面的重合面来计算[7].
3.2 水轮机数值计算为了进一步研究卧式水轮机叶轮的工作机理, 重点研究叶片展长对水轮机能量利用率的影响.本文采用滑移网格技术进行非定常数值计算.
3.2.1 物理模型简化选择叶片展长分别为100、200、400、600、800 mm时进行三维建模.在数值计算中, 为了在降低计算量的同时又不影响计算精度, 对叶轮结构进行了相应的简化, 仅保留5片叶片.
3.2.2 计算域选取及边界条件设置为了使来流充分发展, 避免由于计算区域过小使计算结果产生很大误差, 将计算区域设定为由前端为半径5 m、高为10 m的半圆柱, 以及10 m×10 m的长方体组合而成, 叶轮位于计算域中间且距上游3 m.计算域和边界条件如图 6所示.计算域由1个外部的静止计算域、1个内部的旋转计算域组成.5个叶片位于内部旋转计算域中.计算域的边界条件:静止计算域和滑动计算域的交界面设定为滑移边界; 出口边界设为自由流出口; 叶片表面设定为无滑移壁面; 考虑实际工作情况下的水流速度, 将入口边界设为恒定速度入口.
|
图 6 计算域和边界条件 Fig.6 Computational region and boundary conditions |
利用Gambit软件进行网格划分, 网格可分为结构化网格和非结构化网格.结构化网格具有生成速度快、数据结构简单、网格质量好、网格划分光滑和更能反映实际模型的特征等优点.但是, 结构化网格的使用范围相对较窄, 复杂模型则比较难生成结构化网格.非结构化网格是一种无规则随机的网格结构, 可以很好地弥补结构化网格的缺陷, 但是计算存储量非常大.因此, 本文采用结构化网格和非结构化网格复合的方法进行网格划分, 并在壁面上设置边界层进行局部加密, 最终获得的网格数大于1×106, 壁面网格为0.1 mm.水轮机网格如图 7所示.
|
图 7 水轮机网格 Fig.7 Turbine mesh |
计算过程基于SST k-ω模型.SST k-ω模型适用于具有涡流、快速应变的多重剪切流动, 尤其适用于计算回旋机械的数值计算.采用非平衡壁面函数的近壁面处理和二维稳态分解法的隐式解法.在控制方程的离散格式中, 压力插值采用Standard方式, 压力速度耦合采用SIMPLE算法, 动量方程、紊动能k方程和耗散率ε方程均采用二阶迎风格式[10].
3.2.5 计算结果验证SST k-ω湍流模型经常被用来模拟旋转叶轮.为了确保SST k-ω湍流模型和滑移网格技术对流场模拟的准确性, 使用之前将在某研究所进行的实验结果和数值模拟数据进行对比, 结果如图 8所示.实验结果与计算结果的变化趋势一致, 而且实验值整体小于模拟值, 最大误差在15%左右.这是由于实验中模型尺寸较小, 质量小, 会有能量损失.因此, 两者误差在合理范围, 证明SST k-ω模型和滑移网格技术的模拟数据可以反映卧式水轮机水动力学特性的真实信息.
|
图 8 实验值和模拟值对比 Fig.8 Comparison between the test and numerical simulation results |
叶片展长为800 mm的水轮机叶片压力云图和速度矢量图如图 9所示.由图 9(a)可知, 叶片中心位置所受的压力最大, 然后向两端递减, 尾端压力最小, 很容易产生三维效应, 影响水轮机能量捕获效率.图 9(b)为叶片速度矢量图, 可见, 最接近来流的叶片所受到的压差最大, 动能也最大, 且随着展长的增加压差慢慢减小.这也解释了叶片展长越大其启动速度越高的原因.
|
图 9 叶片压力云图和速度矢量图 Fig.9 Pressure contours and velocity vectors of the blade |
由式 (8) 可得叶轮能量利用率随着尖速比的变化关系, 结果如图 10所示.随着叶轮尖速比的增加, 叶轮能量利用率先增加后减小.当λ=3时, Cp出现最大值, 为0.449.能量利用率的这种变化趋势可以解释为:该水轮机叶片为升力型叶片, 随着尖速比λ的增加, 叶轮上平均扭矩系数呈近似直线递减, 而能量利用率是平均扭矩系数和尖速比之积, 因此平均功率系数与尖速比之间是二次函数的关系, 所以存在极大值Cp, max.叶轮的驱动力为水轮机叶轮的升力, 所以叶轮转速可以大于来流速度.
|
图 10 λ-Cp图 Fig.10 Relationship between λ and Cp |
(1) 基于SST k-ω模型建立的计算流体力学模型, 在网格数大于1×106、壁面网格为0.1 mm时的计算结果可以很好地预测叶片水动力的变化趋势.
(2) 实验结果和数值模拟结果对比证明:结合滑移网格技术的数学模型可以反映卧式轮机的流场特征.
(3) 叶片展长对叶轮转矩和能量利用率影响较大, 主要是由于叶片的切向力起主要作用.叶片展长越大, 叶片受到的切向力越大, 因此, 水轮机对能量的吸收能力越强.
(4) 随着叶片展长的增加, 叶轮能量利用率随着尖速比的增加先增加后减小.在低尖速比情况下, 叶片展长的影响不是很大, 但在尖速比大于2.5时, 叶片展长对能量利用率的影响显著.
| [1] | GUERRI O, SAKOUT A, BOUHADEF K. Simulations of the fluid flow around a rotating vertical axis wind turbine[J]. Wind Turbine , 2007, 31 (3) :149–163. |
| [2] | MYERS L E, BAHAJ A S. An experimental investigation simulating flow effects in first generation marine current energy converter arrays[J]. Renewable Energy , 2012, 37 (1) :28–26. DOI:10.1016/j.renene.2011.03.043 |
| [3] | SAHA U K, THOTLA S, MAITY D. Optimum design configuration of Savonius rotor through wind tunnel experiments[J]. Journal of Wind Engineering and Industrial Aerodynamics , 2008, 96 (8/9) :1359–1375. |
| [4] | MOHAMED M H, JANIGA G, PAP E, et al. Optimization of Savonius turbines using an obstacle shielding the returning blade[J]. Renewable Energy , 2010, 35 (11) :2618–2626. DOI:10.1016/j.renene.2010.04.007 |
| [5] | 黄鹏, 王宏光. 相同叶尖速比不同转速的垂直轴风力机气动性能分析[J]. 能源研究与信息 , 2014, 30 (3) :150–156. |
| [6] | ESTEBAN M, LEARY D. Current developments and future prospects of offshore wind and ocean energy[J]. Applied Energy , 2012, 90 (1) :128–136. DOI:10.1016/j.apenergy.2011.06.011 |
| [7] | 马舜, 李伟, 刘宏伟, 等. 25 kW独立运行式水平轴潮流能发电系统[J]. 电力系统自动化 , 2010, 34 (14) :18–22. |
| [8] | 王树杰, 鹿兰帅, 李冬, 等. 潮流发电帆翼式柔性叶片水轮机实验研究[J]. 海洋工程 , 2009, 27 (1) :83–89. |
| [9] | 王树杰, 李冬, 赵龙武, 等. 柔性叶片水流发电模型实验研究[J]. 太阳能学报 , 2010, 31 (3) :362–366. |
| [10] | 韩占忠, 王敬, 兰小平. FLUENT:流体工程仿真计算实例与应用[M]. 北京: 北京理工大学出版社, 2004 . |