布置粗糙带的垂直轴风力机输出特性研究

《布置粗糙带的垂直轴风力机输出特性研究》论文发表期刊：《可再生能源》；发表周期：2020年12期

《布置粗糙带的垂直轴风力机输出特性研究》论文作者信息：通讯作者： 陈永艳（1976-），女，博士，副教授，研究方向为风能利用技术。

　　摘要：为提高垂直轴风力机的风能利用率，改善其输出特性，文章首先在NACA0012翼型压力面后部布置粗糙带，以粗糙元高度、间距高度比以及叶片安装角作为3个影响因素进行正交设计;然后通过风洞实验和数值模拟相结合的方法进行对比验证。验证结果表明：粗糙元高度对风力机功率系数的影响比间距高度比大;粗糙带可显著提升风力机的扭矩系数，最高可提升33%;粗糙带可有效提高风力机的起动性能，当风速为6m/s时，与光滑翼型风力机相比，平均静态扭矩系数可提高49%。

　　关键词：垂直轴风力机;粗糙带;输出特性;扭矩系数

　　0引言

　　由于化石能源的不可再生性，以及世界范围内环境问题的日益严重，大力发展包括风能在内的可再生能源已经成为人们的共识。垂直轴风力机以其结构简单、无需对风等独特优势受到业界普遍关注，国内外学者在优化翼型改善气动性能方面进行了诸多研究，并取得了较大进展1。

　　Nagao SP和Boese M"将沟槽分别布置于叶片压力面和吸力面进行研究，结果表明，吸力面布置沟槽可以减小阻力，压力面布置沟槽则会增加阻力，严重影响风力机的气动性能。汪泉问分别用对称翼型和非对称翼型建立仿真模型进行二维模拟计算，结果表明：在低叶尖速比下，非对称翼型的气动特性较好;在高叶尖速比下，对称翼型的气动特性较好。胡丹梅 对融合小翼风力机在额定工况下的气动性能进行了数值模拟，得到了不同倾角的小翼会改变叶尖流场分布，从而改变叶尖部分和小翼的输出功率的结论。郭少真n将不同的小翼布置于风力机的后缘，结果表明，风力机可通过小翼消除叶尖涡，进而提升扭矩系数。董世充对翼型表面粗糙度变化时，在两个特定攻角下翼型气动性能的变化情况进行了研究，结果表明，在翼型吸力面、前缘和尾缘布置粗糙带时，翼型的气动性能变化较大。李德顺网运用数值模拟的方法对DU96-W-180翼型前缘不同磨损形貌进行探究，结果显示：形貌为小坑时，其对升力系数和阻力系数的影响均较小;形貌为脱皮时，叶片会提前发生流动分离，与光滑翼型相比，升阻比最高可减少60%左右。张旭[0针对考虑粗糙度敏感位置的风力机翼型钝尾缘改型前后的气动性能进行了研究，揭示了钝尾缘改型对表面粗糙翼型增升效果的影响规律，钝尾缘改型使升力系数和最大升阻比均明显升高，显著改善了表面粗糙翼型的气动性能。上述研究表明，表面粗糙度对翼型气动性能有较大影响。目前，研究的翼型集中在水平轴风力机翼型上，对垂直轴风力机翼型的研究鲜有涉及。

　　基于此，本文将粗糙带应用于垂直轴风力机，将其布置于翼型压力面后部0.1C-0.6C处(C为弦长)，通过风洞实验和数值模拟相结合的方法进行对比验证，探究粗糙带对垂直轴风力机输出特性的影响，为进一步改善垂直轴风力机的气动性能提供参考。

　　1正交试验设计

　　粗糙带的主要参数为粗糙元高度h和间距高度比d，又因为叶片安装角对垂直轴风力机的气动性能影响较大，本次正交试验选取了3个因素，分别为安装角a，h和d，其中d为两个粗糙元之间的距离与h的比值。a和d均选取了4个水平，根据以往试验数据，a选择了效果较好的4个角度。试验表明，当风速为15m/s，d为7时，增阻效果最好，结合设定风速，本次正交试验选择了该4个水平，详细数值如表1所示。

　　考虑到粗糙元太高可能对风力机原有气动性能产生反作用，故h选取了3个水平，因为没有现成的混合水平正交表，故须要用到“拟水平法"，具体为根据实际经验，在h现有的3个水平中选一个较好的水平作为h的第四水平，这样便可以采用Lo(4)正交表来安排试验。根据以上各因素和各水平的要求，采用La(4)正交表，该表为3因素4水平，试验次数为16次。

　　2风洞实验

　　本次风洞实验采用NACA0012翼型叶片，风力机布置如图1所示，叶片参数如表2所示。进行风洞实验的目的是确定适合布置粗糙带的最佳a和粗糙带参数，这不仅可以为后续数值模拟提供指导和依据，还可以与数值模拟结果进行对比验证。

　　2.1风洞实验方案

　　实验所用粗糙元截面为长方形，考虑到可操作性，粗糙元截面长度为5mm。为了减缓叶片压力面边界层的气流速度，增大内、外表面压差，进而提高升力系数，故将其布置于叶片压力面后部0.1C-0.6C处。实验在内蒙古工业大学风能基地低速风洞进行，通过在风洞口手持热线式风速仪标定风速，输出风速为9m/s，通过调节负载箱改变转速进而控制叶尖速比A，通过观察Fluke高精度功率分析仪来实时记录风力机瞬时功率。

　　2.2风洞实验结果

　　将测得的实验数据进行处理后，得到无粗糙带叶片风力机与布置不同粗糙带风力机在不同转速下的功率值，发现风力机在不同工况下测得的功率随入的变化趋势基本一致，均在A为1.6时取得最大值。图2为无粗糙带风力机和布置1-3型粗糙带风力机的功率随入的变化曲线，1-3型粗糙带具有普适性，一般性。

　　由图2可知，当入小于1.6时，粗糙带对风力机功率的提升作用较为明显。

　　以A为1.6时进行分析，根据经验a设为20，得到此时的h与d的交互作用表(表3)。表中数值为该工况下的功率值与光滑翼型风力机功率值的差值。

　　依据表3的分析结果，当d设为3时，得到h与a的交互作用表(表4)。

　　由表3，4可知：d无明显规律可循，总体来看d较大时功率提升效果较好;h对功率的影响为1 mm>2 mm>3 mm，当h为1 mm时，效果更好;除了安装角之外，h对风力机功率值的影响比d大。

　　3计算模型和数值模拟方法

　　3.1几何模型

　　数值模拟可以灵活地设定参数，为确定更加合理的粗糙带尺寸，在风洞实验的基础上进行数值模拟，图3为翼型NACA0012表面粗糙带分布图。

　　将粗糙元为方形的粗糙带布置于翼型压力面后面0.1C-0.6C处，基于风洞实验结论，h分别为0.8，12 mm和1.4 mm，粗糙元d分别为4，4.5，5和5.5，h与d依次组合匹配，研究h和d对风力机输出特性的影响。为确保充分流动，计算域长为20D，宽为10D，风轮距出口的距离为L=15D。图4为计算区域示意图。

　　由于网格数量与质量会直接影响模拟的效率与模拟结果的准确性，针对翼型布置粗糙带的网格区域，采用三角形非结构化网格对其进行加密，第一层网格高度Ah=0.002 mm，Y=1.5，时间步长为0.004 7 s，为风轮旋转100所需时间。并设计多种网格划分方案，进行网格无关化验证，验证结果如图5所示。布置不同高度粗糙带的风力机最佳网格数量：光滑表面模型为38.9万;h为0.8mm时，网格数量为41.3万;h为1.0mm时，网格数量为40.5万;h为1.2 mm时，网格数量为42.1万;h为1.4mm时，网格数量为43.6万。

　　3.2边界条件

　　壁面附近采用增强型壁面函数。差分格式选取二阶迎风模式，速度与压力的迭代求解采用SIMPLEC算法，收敛条件设置为残差水平下降到10，本流场选取速度入口边界条件，速度方向与入口边界垂直，出口为压力出口，上、下边界和翼型表面均采用无滑移壁面边界条件，静止域与旋转域交界面设置为滑移交界面。

　　4结果与分析

　　4.1粗糙带对风力机转矩系数影响数值模拟转矩系数指风产生的旋转扭矩有多少被风力机利用。

　　Ca-M

　　2PARV21)

　　式中：M为风力机转矩，N-m;V为风速，m/s;R为风轮半径，m;A为风轮扫略面积，m'sp为空气密度，kg/m'。

　　图6所示为不同粗糙带的风力机的C随A的变化曲线。

　　由图6可知：因为经过了正交试验对粗糙带参数的筛选，故模拟采用的不同粗糙带对风力机扭矩系数均产生了积极影响;当A小于0.5时，布置有粗糙带的风力机的Ca较光滑翼型风力机略有提高;当入为0.5-2.5时，不同的粗糙带均对风力机产生了不同程度的改善，并且当A为1.5左右时，风力机的C.取得了最大值，其中粗糙带布置为1.2-4.5的风力机表现最为突出，其最大Ca较光滑翼型风力机提高了33%;当入超过1.5后，由于动态失速严重，C迅速下降。

　　4.2粗糙带对风力机起动性性的影响在此采用光滑翼型和布置1.2-4.5粗糙带的NACA0012翼型风力机为计算模型，在风速为6，8，10 m/s时，分别计算风力机在静止状态下一个转动周期内Ca随旋转角的变化情况，计算角度间隔为10"，计算结果如图7所示。

　　由图7可知：在一个旋转周期内，无论布置粗糙带与否，风力机Ca的总体变化基本趋于一致，除风速为10m/s时布置1.2-4.5型粗糙带的风力机外，其余工况均显示，当方位角为90时，风力机C.最低;当风速为6m/s时，与光滑翼型风力机相比，布置12-45型粗糙带的风力机的C.略有提高，并且这一现象在方位角超过900之后更加明显;当风速分别为8m/s和10m/s时，在方位角为0-600时，布置粗糙带的风力机的Ca略有提升，在方位角为60-1200时明显增大。图8所示为不同风速下有无粗糙带的风力机的平均Ca的柱状图。

　　由图8可知：当风速为6m/s时，无粗糙带凤力机的平均Ga为0.039，布置粗糙带的风力机的平均Ca为0.058，提升了49%;当风速为8m/s时，光滑翼型风力机的平均Ca为0.061，布置粗糙带的风力机的Ca为0.078，提升了28%;当风速为10 m/s时，光滑翼型风力机的平均Ca为

　　0.082，布置粗糙带的风力机的平均C为0.119提升了45%。因此，叶片外表面后部布置粗糙带对风力机起动性影响明显，当风速为6m/s时，粗糙带对风力机Ca的提升效果更加明显。

　　5结论

　　本文通过风洞实验和数值模拟对比验证的方法，研究了不同尺寸的粗糙带对风力机的C、功率以及风力机起动性的影响，得出如下结论。

　　0d无明显规律可循，但可以看出d为4.5时，功率提升效果较好，h对风力机C值的影响较d更大。

　　②不同粗糙带对风力机C均可产生积极影响，其中布置1.2-4.5型粗糙带的风力机，其最大C较光滑翼型风力机提高了33%

　　③粗糙带可有效提高风力机的起动性能，当

　　风速为6m/s时，与光滑翼型风力机相比，布置粗

　　糙带的风力机的平均静态Ca可提高49%。

　　参考文献：

　　[1]Mohamed M H，Ali A M，Hafiz A A.CFD analysis forH-rotor Darrieus turbine as a low speed wind energy converter[DJ]Engineering Science and Technology，2015，18(1)：1-13.

　　[2]舒超，蔡新，潘盼，等.垂直轴风力机改进翼型气动及

　　功率特性研究J.水电能源科学，2015，33(1);162-165.

　　[3]Nagao S，Breugelmans F A E.TS-25 investigation of riblets in a CDB，DCA and 65-S compressor cascade(Session B-1：Cascade Design 1)[A].TransporlationFrontiers for the Next Millennium：Mecting of the Institute of Transportation Engineers[C].Tokyo:Gas Turbine Society of Japan，1999.

　　[4]Boese M，Folttner L.Effects of riblets on the loss behavior of a highly loaded compressorcascade[A].Proceedings of ASME TURBO EXPO 2002[C].Amsterdam：Power for Land，Sea，and Air，2002.

　　[5]汪泉，陈晓田，黄攀，等.垂直轴风力机两种翼型气动性能比较研究[J.可再生能源，2020，38(2)：119-204.

　　[6]胡丹梅，吴志祥，张开华，等，融合小翼对风力机气动性能的影响[U].可再生能源，2019，37(9)：136-

　　1361.

　　[7]郭少真，陈永艳，田瑞，等，加装后缘小翼的垂直轴风力机输出特性研究1J可再生能源，2020，38(2)：193-198.

　　[8]董世充，王晓东，康顺表面粗糙度对翼型气动性能影响的数值模拟[J]工程热物理学报，2015，36(6)：

　　1238-1241.

　　9]李德顺，王成泽，李银然，等.叶片前缘磨损形貌特征对风力机翼型气动性能的影响[J农业工程学报，2017，33(22)：269-275.

　　[10]张旭，刘海龙，王格格，等，考虑粗糙度敏感位置的钝尾缘翼型气动性能研究[J.农业工程学报，2017，33(8)：82-89.