Savonius型垂直轴风力机叶型优化设计.pdf

第 卷第期湖北工业大学学报 年 月V o l N o J o u r n a l o fH u b e iU n i v e r s i t yo fT e c h n o l o g yA u g 收稿日期 基金项目国家自然科学基金()第一作者汪泉(),男,湖北武汉人,湖北工业大学副教授,研究方向为机械系统优化设计,风力发电机叶片关键技术,机器人技术.通信作者王环均(),男,湖北荆州人,湖北工业大学硕士研究生,研究方向为机械工程.文章编号 ()S a v o n i u s型垂直轴风力机叶型优化设计汪泉,王环均,杨书益,王冯云(湖北工业大学机械工程学院,湖北 武汉 )摘要为探索叶型形状对S a v o n i u s型垂直轴风力机性能的影响,得到较优叶型以提高其风能利用率,以回旋式S a v o n i u s型垂直轴风力机的风轮主体为研究对象,提出了基于三角函数的垂直轴风力机叶型廓线表达方法,并根据该方法设计出两种新型垂直轴风力机叶型.采用C F D数值计算方法来计算优化叶型的功率系数,从功率系数、力矩系数、压力分布、速度分布和涡量分布这个方面讨论优化叶型较初始叶型的优越性.结果表明:在来流风速为 m/s的环境中,当叶尖速比为时,风轮最大功率系数值为 ,提升了约,并且在叶尖速比高达时,仍然保持较高的功率系数.本研究对垂直轴风力机叶片选型提供参考依据,对于提高垂直轴风力机的功率特性具有一定的理论意义与工程应用价值.关键词S a v o n i u s型垂直轴风力机;优化叶型;功率系数;三角函数 中图分类号T K 文献标识码A S a v o n i u s型垂直轴风力机具有低风速下自启能力强、无需考虑对风损失等优点,它的缺点也十分明显,如工作速比范围很小、风能利用率较低.为提高风能利用率,更好地利用S a v o n i u s型垂直轴风力机自身的优势,国内外学者从各种几何参数、转子的设计、导流装置等对风力机性能的影响进行了研究.李岩等以重叠比为研究对象,利用P I V测试系统进行试验,发现当S a v o n i u s风力机叶片之间存在适当的重叠比时,可以提高S a v o n i u s型风力机的静态起动性能.王伟等提出一种双侧外形不同的叶轮方案,并对双侧外形参数进行优化,发现相比于常规的S a v o n i u s叶 轮,优 化 风 轮 的 发 电 效 率 提 高 了.I r a b u等发现通过应用导向叶片,S a v o n i u s型垂直轴风力机的性能可以提高约.上述研究很少分析叶型不同对风能利用率的影响,且并没有提出一种通用的函数表达回旋式S a v o n i u s垂直轴风力机叶型,通过函数系数的优化得到性能俱佳的垂直轴风力机.因此,本文以回旋式S a v o n i u s型垂直轴风力机的风轮主体为研究对象,提出基于正弦三角函数的垂直轴风力机叶型廓线表达方法,并根据该方法设计出两种新型垂直轴风力机叶型,最后通过功率系数和力矩系数展现优化叶型对性能的改善程度,并通过压力分布、速度分布和涡量分布这个方面的流场分析探讨其优化原理.S a v o n i u s型垂直轴风力机建模C F D二维模型建立如图所示,参照S h e l d a h l等进行的风洞试验实物模型,建立回旋式S a v o n i u s型垂直轴风力机二维模型.模型的详细参数如表所示.表回旋式S a v o n i u s型垂直轴风力机模型参数模型参数数值叶片数N叶片弦长L/m叶片厚度t/mm风轮直径D/m 风轮旋转半径R/m 风轮高度H/m重叠距离O/mm 叶片展弦比S重叠比e 由于垂直转轴对风力机功率的影响微乎其微,故最终风力机计算域简化模型如图所示,风力机计算域大小为 R R(R为风力机旋转半径).将计算域分为中心旋转区域和远场静止区域,A表示静止区域,B表示旋转区域.旋转区域与静止区域之间的边界为交界面,能实现数值模拟过程中旋转区域与静止区域之间数据传递.设定图中风轮方位角为,叶片按空间位置分布,即风轮所处方位角的变化分为两个部分:前进叶片()和回转叶片().图回旋式S a v o n i u s型垂直轴风力机模型图计算域网格划分由于叶片形状是两个对称的半圆形,叶片周围的网格划分很难采取结构化网格,因此在静止区域采用结构化网格,针对旋转区域则采用非结构化网格.整个计算域网格如图 a所示,旋转区域网格划分如图 b所示,叶片边界层附近网格划分如图 c所示,静止区域与旋转区域交界面处的网格划分如图 d所示.由于计算的是瞬态运动问题,静止区域和旋转区域之间的相对运动必须考虑,非定常的相互作用不可忽略,可以采用滑移网格或者动网格实现静止区域和旋转区域的流场耦合求解.图计算域网格划分湍流模型与求解器设置本文选择剪切应力传输(S S T)k模型来模拟湍流,对近壁面区域的边界层进行了增强处理;选择压力基求解器,使得求解更稳定,收敛更快速.为了提高计算效率计算精度,采用工程常用的S I M P L E算法,离散格式为二阶迎风(S e c o n d o r d e r u p w i n d).计算模型边界条件设置如图所示,根据文献 的计算模型边界条件设置.考虑到垂直轴风力机适合在较高速状态下运行,具有较为稳定良好的性能,在 m/s的风速时,垂直轴风力机的风能利用率较高,将计算域左边设为速 度 入 口(v e l o c i t y i n t l e t),来 流 风 速 大 小 为 m/s,方向从左到右,湍流强度为.右边设为压力出口(P r e s s u r e o u t l e t),因两个叶片相同,载荷和边界条件对称,为了减少计算量,上下边设为对称边界(S y mm e t r y).旋转区域与静止区域的交界处设置为交界面(I n t e r f a c e),便于数值模拟时旋转区域与静止区域之间的数据传递.叶片边界设置为移动壁面(M o v i n g w a l l,n o s l i p),叶片近壁面处边界层网格节点间隔宽度设置为mm,初始高度设为 mm,增长率为,共 层.网格独立性验证本文建立种不同网格密度的计算域网格,网格数量分别为 W、W和 W.在叶尖速比为下,种不同网格数量的模型经数值计算后,得到力矩系数曲线如图所示.经计算得出 W、W和 W网格数量对应的平均力矩系数分别为 、和 .从图中可 以 看 出 W网格的力矩系数曲线与其他两套网格相差较大,而 W网格和 W网格的力矩系数曲线几乎重合,其平均力矩系数差值小于.此外,因 W网格的计算时长要远多于 W网格,考虑到计算时长的因素,选择 W网格作为后续风力机功率的计算.060120 180 240 300 3600 6.0 4.0 2.?Cm66W84W102W图不同网格数量下单个叶片力矩系数曲线 叶型优化设计叶型优化普通型S a v o n i u s型垂直轴风力机的S I N 型和双圆弧形叶型均在不同方面性能优于传统半圆形叶片.为了结合这两种叶型的优势,采用复合正弦三湖北工业大学学报 年第期角函数表达类似双圆弧结构的叶型(简称S I N 型)本文进一步提出一种采用正弦三角函数表达回旋式S a v o n i u s型垂直轴风力机叶片叶型廓线的参数化方法,因为研究对象是回旋式S a v o n i u s型垂直轴风力机叶型,重叠比e大于出现重叠部分,此时的S I N 型和S I N 型的数学模型表现为分段函数,以原点呈对称分布.优化叶型表达式如下:yais i n(bixc)x ,ais i n(bixc)x,i,n参照S a n d i风洞试验实物模拟的主要几何参数,构建符合条件的回旋式S a v o n i u s型垂直轴风力机S I N 型和S I N 型的风轮转子简化模型.在构建S I N 型 的 二 维 转 子 模 型 时,需 要 使 用 软 件MAT L A B按照优化叶型表达式计算参数ai、bi和ci,并通过细微调整优化,发现当i为时,能够拟合得到一条形状类似双圆弧结构的复合正弦三角函数模型.参照S a n d i风洞试验的数据,得到的优化参数ai、bi和ci比较准确.因为分段函数的两部分(即前进叶片和回转叶片)以原点对称,只需要计算出一部分的参数即可.以前进叶片为例,当i时,得到符合条件的优化叶型S I N 型的简化模型(图 a);当i时,得到优化叶型S I N 型(图 b).计算所得到的S I N 型和S I N 型的具体参数见表和表.表S I N 型参数参数数值参数数值a c /b/表S I N 型参数参数数值参数数值a b a c a c b c b 图三角函数表达叶型的简化模型叶型优化结果对S I N 型和S I N 型进行数值计算并分析其结果,再通过与半圆形叶型对比来验证优化叶型的优越性.由于叶片形状的不同,不同叶型经网格划分得到的网格数量并不会相同.为忽略网格数量对计算 结 果 造 成 的 影 响,选 择 网 格 数 量 在 范 围 为 W均可.将符合条件的m e s h文件导入软件F L U E NT中进行求解.为避免时间步长对计算结果的影响,选择旋转 为一个时间步长.计算收敛后在第周期达到稳定状态,故计算数据均选自第周期.S a v o n i u s型垂直轴风力机的平均功率系数可以反映风力机风轮在多个旋转周期内的风能转换率,是衡量垂直轴风力机发电性能的关键指标.接下来针对不同叶尖速比下改进叶型的S a v o n i u s型垂直轴风力机平均功率系数的改善进行讨论,两种新型叶型和传统半圆形叶型的风力机风轮的平均功率系数曲线如图所示,倒三角表示的是S h e l d a h l进行风洞试验所得出的功率系数曲线,正方形、圆形和正三角分别表示传统半圆形、S I N 型和S I N 型经过C F D数值计算得出的风力机平均功率系数曲线.由于C F D数值计算半圆形叶型的平均功率系数时做了简化处理,计算的风能利用率比试验值大,从图中可以看出,相同叶型的半圆形与试验值相比,其风轮平均功率系数曲线变化趋势相同而数值整体偏大.不同叶型风轮对应的平均功率系数随叶尖速比的增大都是呈现先增大后减少的趋势,但是曲线数值变化上的差异性较大.从整体上看,S I N 型的性能略优于半圆形叶型,S I N 型风轮的平均功率系数远大于S I N 型和半圆形叶型.当叶尖速比较小时(),优化叶型S I N 型和S I N 型风轮的平均功率系数均小于半圆形叶型.图不同叶尖速比下的的平均功率系数可以明显看出,S I N 型风轮的最大平均功率系数最大,出现在时值为 ,增长值为 ,增长率为,较传统半圆形风轮的最大平均功率系数 ,增长了 ,大约提升了.当时,S I N 型风轮的平均功率系数略低于传统半圆形;当时,其风轮平均功率系数高于传统半圆形,并且增长值随着叶尖速比的增大而提高;当时,增长率高达;在较高速运转()时,风力机仍然保持较高的平均功率系数,甚至是在时,其平均功率系数值都保持在较高水平上.第 卷第期汪泉,等S a v o n i u s型垂直轴风力机叶型优化设计力矩系数直接决定了功率系数的大小,平均力矩系数表示风力机风轮在多个旋转周期内的风力机输出力矩能力的系数.图展示的是不同叶尖速比下叶片上的平均力矩系数曲线,倒三角表示的是根据参考风洞试验功率系数计算出来的力矩系数曲线,正方形、圆形和正三角分别表示传统半圆形、S I N 型和S I N 型经过C F D仿真计算出来的叶片平均力矩系数曲线.从图中可以看出,各叶型叶片平均力矩系数的变化大致相同,都呈现出随叶尖速比增大而降低的趋势.相对其他叶型,S I N 型叶片的平均力矩系数曲线显得更加平缓.图不同叶尖速比下的叶片平均力矩系数流场特性分析压力分布图展示了传统半圆形、S I N 型和S I N 型在叶尖速比为时风轮附近的压力分布云图.从整体趋势上分析,当方位角 时,最大压力出现在前进叶片的叶尖位置附近;随着风轮的旋转,最大压力出现的位置转移到回转叶片凸侧附近;最后当方位角为 时,回转叶片变为前进叶片,最大压力又回到前进叶片的叶尖处.前进叶片的凸侧和回转叶片的凹侧均受到较大的负压,而这些负压在数值上均大于最大压力,两个叶片凹凸侧压力差的大小不同,使得转子上产生力矩为风轮的旋转提供了转矩.优化叶型S I N 型的负压范围更大,且其回转叶片的凸侧最大压力范围也更大,造成的压力差就更大,也使得S I N 型转子上的转矩更大.当方位角 时,回转叶片的凹侧和前进叶片凸侧靠近叶尖处的负压较大,回转叶片凹侧的负压只有分布范围的区别,而数值大小的范围相同,个叶型前进叶片凸侧靠近叶尖处的负压大于 P a的范围大小排序为:S I N 型S I N 型半圆形.可以明显看出,S I N 型前进叶片叶尖处的负压最大,S I N 型其次,而传统半圆形叶型最小.为了更好地分析叶片上压力分布的影响,建立压力系数C P分布散点图.如图所示,实线、划线和点线分别代表S I N 型、半圆形和S I N 型叶片上的压力系数分布散点曲线.由图可见,无论方位图下的压力分布云图角如何变化,S I N 型前进叶片处压力系数的上部分推力侧曲线略低于半圆形,而下部分吸力侧曲线则远低于半圆形,推力侧与吸力侧的差值大于半圆形.在回转叶片处,随着方位角的增大,S I N 型推力侧和吸力侧曲线均呈现出,从低于半圆形到高于半圆形最后又低于半圆形的变化趋势,即在负压较大的叶片凹侧,S I N 型前进叶片处的负压值和压差大小明显大于半圆形叶型.而在回转叶片处,S I N 型叶片凹侧的负压值不总是大于半圆形叶型,压差大小没有过于明显的差距.优化叶型S I N 型叶片整体上受到的压差更大,风轮产生的旋转力矩也就湖北工业大学学报 年第期更大,相应的功率系数也就更大.图下的压力系数分布散点图速度分布对比分析对风轮及周围的流场速度分布进行分析,有助于进一步了解风力机叶片动态失速的特性.动态失速是一种随时间变化而产生的气流分离和失速现象,主要包括流体在叶片边界层与吸力侧表面的分离,以及随后被卷入叶尖的过程,常在非定常运动中发生,可能会导致运动部件疲劳和结构失效.动态失速现象是垂直轴风力机气流的主要特征,是导致其发电效率低的主要因素之一.应用计算流体力学的方法,计算过程中发现,在叶尖速比较低时,风力机叶片的气动性受到叶片失速现象严重.图 中展示的是叶尖速比为时的速度分布云图.该叶尖速比下的功率系数最大,气动性能最好,叶片产生的动态失速现象相对较弱,但是在叶片附近还是能明显看到涡旋分离.从整体上分析,气流在叶片边界层的分离主要发生在叶根处、前进叶片的叶尖处和回转叶片的凸侧;随叶片旋转,前进叶片叶尖处形成的涡旋逐渐脱落汇入下游尾流,叶根处的涡旋只有少量涡脱落汇入回转叶片叶尖处,而回转叶片凸侧的涡旋在逐渐分离过程中形成一道涡旋轨迹.可以明显看出,相比于半圆形叶型,优化叶型S I N 型的前进叶片叶尖处产生的分离现象较小,叶根处也没有明显的涡脱落现象,回转叶片叶尖处的涡旋轨迹更为连贯,产生的涡脱落现象减弱.S I N 型涡旋轨迹更为连贯,脱落的涡量与S I N 型大致相同,但叶根处的涡旋还存在有少量涡脱落现象.S I N 型较为明显地减少了涡旋的脱落,有效改善了动态失速现象.=0?=30?=60?=90?=120?Velocitym/s0 4 7.9 4.14 1.18 8.23 5.28 2.32 9.37 6.42 3.47第 卷第期汪泉,等S a v o n i u s型垂直轴风力机叶型优化设计()a?()b SIN1?=150?=180?()c SIN2?图 下的速度分布云图涡量分布对比分析涡量是描述涡旋运动最重要的物理量之一,即流体速度矢量的旋度,可以用来度量涡旋的强度和方向.图 展示的是种叶型风力机在叶尖速比为时的涡量分布情况.从图中可以看出,涡旋源主要位于前进叶片凸侧叶尖处和回转叶片凹侧的叶根处.当方位角为 时,新涡旋产生,接着随风轮的旋转,涡旋慢慢从叶片脱离;当旋转至方位角为 时,前进叶片处的涡旋从右下方分离;当方位角为 时,上一个周期的回转叶片上形成的涡旋分离.由于旋转区域是逆时针旋转,涡旋在脱离过程中,都会有一部分随风轮的旋转产生尾流,而这部分尾流在逆时针旋转过程中会重新融入回转叶片上方的涡旋中.与图 a中的半圆形叶型相比较,图 b中的S I N 型和图 c中的S I N 型前进叶片的凹侧存在较大范围的涡,回转叶片凹侧的涡旋范围也更大,尤其是S I N 型更加明显.不同于S I N 型回转叶片处的涡旋在凹侧靠近叶根位置,S I N 型回转叶片叶根处的涡旋在凸侧,凹侧中心位置也存在更明显且范围更大的涡旋,而且在前进叶片凸侧的叶根处涡旋也特别明显.从以上的分析来看,在整旋转周期内,优化叶型S I N 型的叶片凹侧始终保留有大范围且涡量较大的涡旋,这意味着优化叶型有效减少了涡分离带来的力矩损失.=0?=30?3268104 140 176 212 248320 356 392 428 464 500-40-4284Velocitys1-()a?=60?=90?=120?=150?()b SIN1?=180?()c SIN2?图 下的涡量分布云图 结论本文以回旋式S a v o n i u s型垂直轴风力机的风轮主体为研究对象,提出了采用三角函数的数学模型来表达叶型,并成功得到了优化叶型.通过对比分析传统半圆形叶型和优化叶型在风力机风轮功率系数、力矩系数以及风轮在计算收敛后的旋转过程中速度分布、压力分布和涡量分布情况,发现优化叶型有效提高了回旋式S a v o n i u s型垂直轴风力机的风能利用率.对回旋式S a v o n i u s型垂直轴风力机的研究中,得到了如下主要结论.)从整体上看,S I N 型的性能要略强于半圆形叶型,而S I N 型的性能远强于S I N 型和半圆形叶型.从S a v o n i u s型垂直轴风力机的平均功率系数、平均力矩系数来分析种叶型风力机的气动性能.通过对比传统半圆形叶型与优化叶型S I N 型和S I N 型,发现S I N 型只是在较大叶尖速比情况下优于半圆形叶型,而S I N 型的气动性能最优,且在叶尖速比为时功率系数最大.最后,在压力分布、速度分布和涡量分布这三个方面,对比分析种湖北工业大学学报 年第期叶型风轮的流场特性,从而印证了优化叶型S I N 型的性能优越性.)在来流风速为 m/s的环境中,当叶尖速比为时,优化叶型S I N 型风轮的功率系数最大,较半圆形叶型风轮的最大功率系数增长了 ,提升了约,并且在叶尖速比高达时,仍然保持优较高的功率系数.)相比于半圆形叶型,优化叶型S I N 型的叶片受到的压差更大,风轮产生的旋转力矩更大;涡旋分离现象减弱,动态失速现象有所改善;涡量流失较小,叶片上的力矩损失有所降低.参考文献李岩,赵守阳,曲春明,等 S a v o n i u s风力机静态流场P I V可视化试验研究J排灌机械工程学报,():王伟,宋保维,毛昭勇,等 S a v o n i u s风机叶轮双侧外形优化设计J哈尔滨工程大学学报,():I R A B U K,R OYJN C h a r a c t e r i s t i c so fw i n dp o w e ro nS a v o n i u s r o t o ru s i n gag u i d e b o x t u n n e lJ E x p e r i m e n t a lT h e r m a la n dF l u i dS c i e n c e,():S HE L D AHLRE,F E L T ZLV,B L A C KWE L LBFW i n d t u n n e l p e r f o r m a n c ed a t a f o r t w o a n d t h r e e b u c k e tS a v o n i u sr o t o r sJ J o u r n a lo fE n e r g y,():陈晓田 H型垂直轴风力机翼型气动外形设计与性能评价研究D武汉:湖北工业大学,汪泉,甘笛,杨书益,等考虑攻角范围的垂直轴风力机叶片翼型优化设计J农业工程学报,():张仲柱,王会社,赵晓璐,等水平轴风力机叶片气动性能研究J工程热物理学报,():雷航垂直轴风力机非定常空气动力特性数值模拟分析D上海:上海交通大学,P A B L O O,THOR S T E NS,L U I SR E f f e c to fb l a d ec a m b e r i n go nd y n a m i cs t a l l i nv i e wo fd e s i g n i n gv e r t i c a la x i st u r b i n e sJ J o u r n a lo fF l u i d sE n g i n e e r i n g,():B U CHN E RA,L OHR Y M,MA R T I N E L L IL,e ta l D y n a m i cs t a l l i nv e r t i c a la x i sw i n dt u r b i n e s:C o m p a r i n ge x p e r i m e n t sa n dc o m p u t a t i o n sJ J o u r n a lo fW i n dE n g i n e e r i n ga n dI n d u s t r i a lA e r o d y n a m i c s,:O p t i m i z a t i o nD e s i g no fS a v o n i u sV e r t i c a lA x i sW i n dT u r b i n eB l a d eWANGQ u a n,WANG H u a n j u n,YANGS h u y i,WANGF e n g y u n(S c h o o l o fM e c h a n i c a lE n g i n e e r i n g,H u b e iU n i vo fT e c h,W u h a n ,C h i n a)A b s t r a c t:I no r d e r t oe x p l o r e t h e i n f l u e n c eo f b l a d e s h a p eo n t h ep e r f o r m a n c eo f S a v o n i u sv e r t i c a l a x i sw i n dt u r b i n e,ab e t t e rb l a d es h a p ew a so b t a i n e dt oi m p r o v e i t sw i n de n e r g ye f f i c i e n c y I nt h i ss t u d y,S a v o n i u st y p er o t a r yv e r t i c a l a x i sw i n d t u r b i n ew a s t a k e na s t h e r e s e a r c ho b j e c t,a n d t h e e x p r e s s i o nm e t h o do f v e r t i c a l a x i sw i n dt u r b i n eb l a d ep r o f i l eb a s e do nt r i g o n o m e t r i c f u n c t i o nw a sp u t f o r w a r d,a n dt w on e wv e r t i c a la x i sw i n d t u r b i n eb l a d ep r o f i l e sw e r ed e s i g n e da c c o r d i n g t o t h i sm e t h o d C F Dn u m e r i c a l c a l c u l a t i o nm e t h o di su s e d t oc a l c u l a t e t h ep o w e r c o e f f i c i e n t o f t h eo p t i m i z e db l a d ep r o f i l e F i n a l l y,t h e a d v a n t a g e so f t h eo p t i m i z e db l a d ep r o f i l eo v e rt h ei n i t i a lo n ea r ed i s c u s s e df r o mf i v ea s p e c t s:p o w e rc o e f f i c i e n t,t o r q u ec o e f f i c i e n t,p r e s s u r ed i s t r i b u t i o n,v e l o c i t yd i s t r i b u t i o na n dv o r t i c i t yd i s t r i b u t i o n T h er e s u l t ss h o wt h a t i nt h ee n v i r o n m e n to f i n c o m i n gw i n ds p e e do f m/s,t h ep o w e r c o e f f i c i e n t r e a c h e s t h em a x i m u mw h e nb l a d e t i ps p e e dr a t i o i s,a n dt h em a x i m u mp o w e rc o e f f i c i e n t i s ,a ni n c r e a s eo fa b o u t Wh e nt h et i ps p e e dr a t i o i sa sh i g ha s,t h ep o w e rc o e f f i c i e n t i ss t i l l h i g h T h i ss t u d yp r o v i d e sar e f e r e n c e f o r t h es e l e c t i o no fv e r t i c a l a x i sw i n d t u r b i n eb l a d e s,a n dh a s a c e r t a i n t h e o r e t i c a l r e s e a r c hs i g n i f i c a n c e a n de n g i n e e r i n ga p p l i c a t i o nv a l u e f o r i m p r o v i n gt h ep o w e rc h a r a c t e r i s t i c so fv e r t i c a l a x i sw i n dt u r b i n e K e y w o r d s:s a v o n i u sv e r t i c a la x i sw i n dt u r b i n e;t oo p t i m i z et h eb l a d e;p o w e rf a c t o r;t r i g o n o m e t r i cf u n c t i o n s 责任编校:张众第 卷第期汪泉,等S a v o n i u s型垂直轴风力机叶型优化设计