基于风洞试验和仿真的垂直轴风力发电机的湍流模型分析.docx

1、基于风洞试验和仿真的垂直轴风力发电机的湍流模型分析基于风洞试验和仿真的垂直轴风力发电机的湍流模型分析 杨燕昭 1，2，郭志平 1，张艳锋 1，2，李庆安 3 （1.内蒙古工业大学 机械学院，内蒙古 呼和浩特 010051；2.日本国立三重大学 共生环境工学研究科，三重县津市 5148507；3.日本国立三重大学 机械工学研究科，三重县津市 5148507） 摘 要：对于湍流的计算一直是流体力学领域的难点，而垂直轴风力发电机周围流场的湍流非常复杂，对其CFD仿真模拟选取湍流模型一直各抒己见。选取目前使用最为广泛的SST k-、Realize k-、RNG k-和标准k-湍流模型进行CFD仿真，模

2、拟双叶片直线翼垂直轴风力发电机的叶片表面压力分布和单叶片转矩系数，并与采用高速多点压力测量仪和力矩传感器测量的实验数据进行对比分析。其结论为，在叶片失速区，采用标准k-和Realize k-湍流模型CFD仿真模拟的结果与实验值拟合度很低；当叶片处于深度失速区时，采用RNG k-湍流模型CFD仿真模拟的结果与实验结果拟合度最高，而在轻失速区时，采用SST k-湍流模型模拟得到的值与实验测量值拟合度最好。从单叶片转矩系数角度分析可得，采用SST k-湍流模型模拟得到的值与实验测量值拟合度最好。 关键词：垂直轴风力发电机；SST k-；Realize k-；RNG k-；标准 k-湍流模型 1 引言

3、 与技术非常成熟的水平轴风力发电机相比，垂直轴风力发电机的发明时间最早，却研发较晚1-2。其原因是垂直轴风力发电机的叶片的相对风速和功角随回转角不断变化，导致其流场极为复杂3-4。如果采用人工计算其流场，一次至少需要几年时间，在没有高性能计算机和实验条件较为落后的时代，垂直轴风力发电机的研究由此被束缚了前进的脚步5。 近年来，随着高性能计算机的快速发展和实验条件的不断改善，且垂直轴风力发电机具有无需对风装置和低风速启动的优良特性，重新获得了科研人员的青睐6。为了节约科研经费和减少研发周期，科研人员大都在实验之前利用CFD仿真技术进行前期预研7。但是，CFD仿真模拟所采用的湍流模型却是各异。目前

4、采用较多的是 SSTk-、Realizek-、RNGk- 和标准 k- 湍流模型，对这四个模型应用在CFD仿真模拟垂直轴风力发电机流场的优缺点各抒己见。文献8采用RNGk-湍流模型进行CFD仿真模拟了小型垂直轴风力发电机的空气动力学特性，认为RNGk-湍流模型计算大窝量分离的流场会更加准确；文献9则采用SSTk-湍流模型模拟小型风力发电机翼型修改的影响，认为此模型结合了k-湍流模型和标准SSTk-湍流模型的优势，对于可变加权函数可以提供更大的权重；文献10利用Realizek-和SSTk-湍流模型进行CFD仿真计算H型小型垂直轴风力发电机的能量转，认为SSTk-湍流模型对于入口自由流过于敏感，

5、且花费仿真时间较多，不如Realizek-湍流模型性能优良；文献11利用标准k-湍流模型模拟了垂直轴风力发电机的三维流场特性。 到底哪个模型计算垂直轴风力发电的气动特性更加准确。对双叶片直线翼垂直轴风力发电机分别采用RNG k-湍流模型和SST k-湍流模型进行CFD仿真模拟，并利用风洞试验进行验证，对比分析得到最适合CFD仿真模拟垂直轴风力发电机的湍流模型。 2 物理模型 垂直轴风力电机的叶轮简图，如图1所示。此风轮采用双叶片直线翼模型，叶轮的回转直径D=2.0m，叶轮高H=1.2m，叶轮的叶片采用NACA0021对称翼型，其弦长c=0.265m。 图1 风机叶轮简图Fig.1Fan Imp

6、eller Diagram 叶片受力和速度矢量图，如图2所示。角度、和分别为叶片的回转角、功角、安装角和总角，叶片的总角为叶片功角和安装角之和。叶片顺时针旋转。沿回转圆半径方向的力FN为法向力，沿回转圆切线方向的力FT为切向力。风速U、V和W分别为叶片的局部风速、切向风速和相对风速，由局部风速U和切向风速V合成得到相对风速W。 图2 叶片受力与速度适量图Fig.2The Force and Velocity of Blade in the Blade CFD仿真的数值域，如图3所示。此域的长度为20D，宽度为10D，高度为D。由于此模型采用滑移网格技术，此数值域分为旋转数值域（Move）和静止

7、数值域（Static）。CFD仿真的旋转数值域，如图4所示。其旋转数值域的外圆半径为1.2m，内圆的半径为0.8m，旋转数值域与静止数值域的交界面采用interface边界条件，入口采用inlet边界条件，出口采用oulet边界条件，分别采用SST k-、Realize k-、RNG k-和标准k-湍流模型。叶片表面附近对网格进行了加密处理，如图5所示。尤其叶片前缘的网格在纵向与横向均得到加密，其原因是叶片前缘的压力梯度大于叶片尾缘所致。 图3 CFD仿真数值域Fig.3CFD Simulation Number Range 图4 CFD仿真旋转数值域Fig.4CFD Simulation R

8、otation Counting Range 图5 叶片附近网格分布Fig.5Grid Distribution Near Blade 3 实验设备 图6 垂直轴风力发电机风洞实验原理图Fig.6Experimental Schematic of Wind Tunnel in Vertical Axis Wind Turbine 实验任务是测量叶片回转一周过程中叶片表面的压力分布。实验采用开口半径为1.5m的大型开口回流型风洞，如图6所示。风力发电机采用1：1比例，测压设备采用多点测量仪，其压力传感器的型号为ZOC22B，测量功率系数采用的力矩传感器的型号为TS-2700。压力补正原理图，如图

9、7所示。叶片表面的压力分布与风洞吹风实验的压力分布相同，叶片表面设置32个压力点，其压力点的布置采用从尾缘至前缘密度逐渐增大，叶片表面压力孔的直径为0.4mm。 图7 压力测量补正系统原理简图Fig.7Schematic Diagram of Pressure Measurement Complement System 4 CFD仿真结果与实验结果对比分析 4.1 叶片表面压力分布对比分析 几乎所有由CFD仿真模拟所得叶片表面压力差都大于实验测量所得，如图8所示。 图8 在各湍流模型下仿真计算和实验测量得到的压力分布对比Fig.8Shows the Comparison of Pressure

10、 Distribution Under Each Turbulent Model and the Experimental Measurement 但总体拟合度较好，上流域的拟合度高于下流域，在下流域CFD仿真模拟误差较大。SSTk-、Realizek-、RNGk-和标准k-湍流模型在叶片尾缘的拟合度非常高。但是，在叶片前缘的拟合度较低，其计算的偏差也主要来源于叶片的前缘。此现象的原因是叶片前缘压力梯度大，相对于叶片尾缘对湍流模型的选择较为敏感。当回转角在 0、30、180、210、240、370、300和 330时，上述四个湍流模型模拟所得压力分布值拟合度很高，在这些回转角度下，叶片表面流场

11、基本处于层流状态，所以，其拟合度很高。但是，当回转角在60、90、120、150时，采用上述四个湍流模型模拟所得叶片表面的压力分布拟合度较低，尤其是标准k-和Realizek-湍流模型模拟的结果差异最大，而此时的叶片正处于失速状态，叶片表面的气流出现分离，叶片附近出现湍流，正是对湍流模型选择优良的验证。有文献12可知，当回转角在90附近时，叶片处于深度失速状态，此时，标准k-湍流模型模拟结果的精度却很低，说明采用标准k-湍流模型模拟深度失速的叶片流场误差很大，而在回转角下，标准k-湍流模型模拟的结果与实验测量的结果拟合度最高，说明，RNGk-湍流模型对于深度失速的模拟更加合适。但是，当回转角在

12、60、120的轻度失速区，SSTk-湍流模型模拟的结果拟合度最好。 4.2 单叶片转矩系数分析 分别在 SST k-、Realize k-、RNG k- 和标准 k- 湍流模型下，如图9所示。对垂直风力发电机进行CFD仿真模拟和实验测量得到的单叶片转矩系数进行对比，五条单叶片转矩系数曲线均在回转角为100附近出现波峰，出现这种现象的原因是被测量叶片正处于叶轮的上流域，由文献7可知由于上流域对下流域流场的影响，使得下流域的出现风能亏损和复杂流场。五条曲线均在回转角为30和180附近出现小波谷，在30附近出现小波谷的原因是此时叶片的相对风速达到最小值，在180附近出现小波谷的原因是此时叶片的切向速

13、度与风向相同。由图还可得到，当回转角在（070）区间，实验测量值大于所以CFD仿真模拟所得值，SST k-和Realize k-湍流模型CFD仿真模拟所得结果更加贴近实验值，当回转角在（70180）区间，SST k- 湍流模型预实验值拟合度最好，Realize k-和标准k-湍流模型的拟合度最差。当回转角在（180360）区间，SST k-、Realize k- 和RNG k-湍流模型的模拟结果很相近，其值小于实验测量值，这是主要因为上流域对下流域的影响，使得下流域气流非常复杂，并且还有相对风速的影响，这些原因都造成了回转角在（180360）时，CFD仿真模拟值小于实验测量值。综上分析，从单叶

14、片转矩系数角度分析可得，采用SST k-湍流模型模拟得到的值与实验测量值拟合度最好。 图9 在各湍流模型下CFD仿真计算和实验计算得到的回转一周单叶片转矩系数对比Fig.9Shows the Comparison of Single-Blade Torque Coefficients in the Rotation of CFD Simulation and Experimental Calculation Under Each Turbulent Model 5 结论 通过上述分析可得以下结论： （1）上述四个湍流模型模拟所得的叶片表面压力差的差异主要发生在叶片前缘；当叶片处于层流区时，上述

15、四个湍流模型模拟所得的叶片表面压力差值拟合度非常好。 （2）在叶片失速区，采用标准k-和Realize k-湍流模型CFD仿真模拟的结果与实验值拟合度很低；当叶片处于深度失速区时，采用RNG k-湍流模型CFD仿真模拟的结果预实验结果拟合度最高，而在轻失速区时，采用SST k-湍流模型模拟得到的值与实验测量值拟合度最好。 （3）从单叶片转矩系数角度分析可得，采用SST k-湍流模型模拟得到的值与实验测量值拟合度最好。 参考文献 1崔民选.中国能源发展报告R.北京：社会科学文献出版社，2010（5）：152-170.（Cui Min-xuan.Chinas energy development

16、reportR.Beijing：Social Science Literature Publishing House，2010（5）：152-170.） 2IEC C D V.61400-12-1：Wind turbines-Part 12-1：Power performance measurements of electricity producing wind turbinesJ.Geneva，Switzerland：IEC，2005. 3Li Q，Maeda T，Kamada Y，Murata J，Furukawa K and Yamamoto M.Measurement of the

17、Flow Field around Straight-bladed Vertical Axis Wind Turbine.Journal of Wind Engineering&Industrial Aerodynamics，2016（151）：70-78. 4Li Q，Maeda T，KamadaY，Murata J，Kawabata T，Furukawa K.Analysis of Aerodynamic Load on Straight-bladed Vertical Axis Wind TurbineJ.Journal of Thermal Science，2014，23（4）：314

18、-324. 5GWEC.Global Wind Outlook 2013，http：/ 9-April-2014，Access 2014.06.28. 6何志宁，刘座铭，张宗珍.垂直轴风力机的发展历史及现状J.机械工程师，2014（2）：39-40.（He Zhi-ning，Liu Zuo-ming，Zhang Zhong-zhen.Development history and current situation of vertical axis wind turbineJ.Mechanical engineer，2014（2）：39-40.） 7Li Q，Maeda T，Kamada Y.

19、Wind tunnel and numerical study of a straight-bladed vertical axis wind turbine in three-dimensional analysis（Part I：For predicting aerodynamic loads and performance）J.Energy，2016（106）：443-452. 8Howell R，Qin N，Edwards J.Wind tunnel and numerical study of a small vertical axis wind turbineJ.Renewable

20、 Energy，2010，35（2）：412-422. 9Ismail M F，Vijayaraghavan K.The effects of aerofoil profile modification on a vertical axis wind turbine performanceJ.Energy，2015（80）：20-31. 10Mohamed M H，Ali A M，Hafiz A A.CFD analysis for H-rotor Darrieus turbine as a low speed wind energy converterJ.Engineering Scienc

21、e and Technology，an International Journal，2015，18（1）：1-13. 11王巧红，王洛明.小型垂直轴风力发电机风轮流场的三维数值模J.机械设计与制造，2010（4）：10-13.（Wang Qiao-hong，Wang Luo-ming.Small vertical axis wind turbine wind turbine wind turbine wind turbine wind power rotation field 3d numerical modelJ.Mechanical Design and Manufacture，2010（

22、4）：10-13.） 12Dunne R，McKeon B J.Dynamic stall on a pitching and surging airfoilJ.Experiments in Fluids，2015，56（8）：1-15. Analysis of Turbulence Model of Vertical Axis Wind Turbine with Wind Tunnel Experiments and Simulations YANG Yan-zhao1，2，GUO Zhi-ping1，ZHANG Yan-feng1，2，LI Qing-an3（1.College of Me

23、chanical Engineering，Inner Mongolia University of Technology，Inner Mongolia Hohhot 010051，China；2.Graduate School of Bioresources，Mie University，Tsu city 5148507，Japan；3.Division of System Engineering，Mie University，Tsu city 5148507，Japan） Abstract：The calculation of the turbulence intensive has alw

24、ays been differences in the field of the fluid mechanics，and the turbulence model of the flow field around the Vertical Axis Wind Turbine is very complicated.It is always different idea that the turbulence model for the simulation of the vertical axis wind turbine is selected.In this paper，the CFD s

25、imulations is carried out by the most widely used turbulence models of SST k-，Realize k-，RNG k- and standard k- to get the blade surface pressure distribution Single-blade torque coefficient.And the simulated results and the experimental data that get by the torque meter and the multi-port pressure

26、devices were compared and analyzed.It is concluded that quality of the CFD simulation results with the standard k- and Realize k- turbulence model are very low in the stall region.When the blades are in the deep stall，the results of CFD simulation with RNG k- turbulence model were in good agreement

27、with the experimental results.The results of CFD simulation with SST k- turbulence model were in good agreement with the experimental results in the light stall region.Through the perspective of single-blade torque coefficient，the simulated value obtained by SST k- turbulence model is the best fit with the experimental value. Key Words：Vertical-Axis Wind Turbine；SST k-；Realize k-；RNG k-；Standard k- Turbulence Model 中图分类号：TH16；TK83 文献标识码：A 文章编号：1001-3997（2017）10-0133-03 来稿日期：2017-04-08 作者简介：郭志平，（1957-），男，内蒙古丰镇人，硕士研究生，博士生导师，教授，主要研究方向：动力机械； 李庆安，（1985-），男，山东廊坊人，博士研究生，研究员，主要研究方向：小型风车的标准化高效化研究