Factors of Simulation Accuracy for Wind Turbine Airfoil Aerodynamic Characteristics
-
摘要: 本文以风力机翼型的气动力特性为研究对象,分别建立了多面体和六面体两种体网格模型,并分别在S-A和SST k-ω湍流模型下进行CFD模拟,还与风洞试验结果进行对比。计算结果表明SST k-ω湍流模型的计算结果相对于S-A模型与风洞试验值更为接近,六面体网格模型的计算结果相对于多面体网格模型与风洞试验值更为接近,并且所划分的网格质量越好CFD模拟得越准确。将等截面模型的风洞实验结果或CFD模拟结果应用于真实叶片的气动力特性计算是比较合理的。Abstract: In this paper, we respectively established the polyhedral and hexahedral volume mesh models, and used CFD method to calculate the aerodynamic characteristics of a wind turbine airfoil under the S-A and SST k-ω turbulence model.The calculated result was compared with wind tunnel test data.It was discovered that the calculation results under SST k-ω turbulence model were more close to the wind tunnel test than S-A turbulence model.And the results showed that choosing better quality hexahedral mesh could get more accurate simulated result than polyhedral mesh.Wind tunnel test or CFD method could reasonably be applied to the aerodynamic characteristics calculation for a real wind turbine airfoil.
-
Keywords:
- CFD Simulation /
- Airfoil /
- Aerodynamic Force /
- Turbulence Model /
- Mesh
-
随着世界能源危机的日益严重,风能作为一种清洁的可再生能源越来越受到各国政府的重视[1]。风电技术的快速发展也使得未来风力发电机风轮直径有着变大的发展趋势,这就要求能更加准确地预测风力机翼型的气动力系数。Christopher等人[2]认为在用CFD模拟翼型气动性能时,二维与三维得到的结果不相同,三维模拟结果比二维的精度要高。潘小涛等人[3]认为风力机三维等截面模型比二维模型的CFD计算结果略精确,特别是失速后的升力系数,且三维模型模拟得到的最大升力系数比二维翼型的要大。但是不管是三维还是二维模型,其计算结果均与风洞试验值存在一定的偏差。因此本文以某风力机专用翼型为研究对象, 采用不同的计算网格与湍流模型,研究不同CFD模型的计算差异,以确定能最准确获得翼型气动力特性的CFD模拟方法。
1. 计算模型
1.1 网格划分
为了研究不同体网格和湍流模型对风力机翼型三维气动力特性的影响,本文针对某风力机专用翼型划分了非结构化和结构化两种体网格,并用其进行运算和后处理分析。首先对展长为10m的等截面翼型三维模型添加外流域,外流域采用半径为30m、高为10.4m的圆柱形空间几何体来表示。需要说明的是翼型位于外流域的正中位置(如图 1.1所示),上下端各有0.2m余量以方便体网格的划分,圆柱形的流场可以在不改变体网格的前提下就不同攻角进行计算分析。然后对三维模型进行体网格的划分工作,并且分别建立了非结构化多面体网格和结构化六面体网格两种模型用以比较(如图 1.2所示),划分网格的具体参数设定已列于表 1.1中。最终多面体总网格数约273万,六面体总网格数约257万。模型采用的翼型弦长设定为1,且弦长方向平行于0度攻角方向。
表 1.1 网格划分的参数设置基本尺寸 边界层 叶片 外流域 目标尺寸 最小尺寸 目标尺寸 最小尺寸 多面体 总厚0.0065m 0.015m 0.005m 0.5m 1.25m 0.5m 共6层 六面体 递增比例1:1.6 0.015m 0.005m 0.25m 1.0m 1.2 湍流模型
目前,工程计算中常用的基于RANS的湍流模型有基于Boussinesq涡粘性假设的零方程模式、一方程模式、两方程模式。由于一方程模型中的S-A模型相对于其它湍流模型而言,不仅计算量小、稳定性好,而且不要求设置壁面函数,所以其计算误差对于网格的划分质量不太敏感,对于翼型问题的计算精度较高[4],参照文献[3]中采用的S-A模型进行分析。但是据资料[5]显示,现在多数航空航天与机械汽车产业中也常用k-ω湍流模型对飞机机翼或是赛车零部件进行CFD分析。于是在保证其他条件相同的情况下,本文也同时采用两方程的SST k-ω湍流模型与S-A模型进行对比。
1.3 边界条件
本文选用的翼型为某公司风力机专用翼型,并已完成3106雷诺数条件下的全攻角风洞试验,对应马赫数为0.126 1。风速由马赫数确定,攻角定义如图 1.3所示。本文CFD模拟选用理想可压缩气体,叶片边界条件为wall,流场周边边界条件为free stream,外流域的上下圆形表面边界条件为symmetry plane,攻角范围为0°至30°,间隔2.5°,且均采用分离求解器求解。
2. 结果分析
2.1 湍流模型不同的影响
与风洞试验采用的三维等截面模型相似,为了方便得到升阻力系数且与风洞试验结果进行对比,本节针对1.1节中所提到的以非结构化多面体网格来划分的三维模型,分别就S-A湍流模型和SST k-ω湍流模型做CFD模拟。收敛标准为所有的残差曲线达到水平,由于S-A模型计算结果中的残差值在大于15°的攻角下无法平稳,故只对0°至20°攻角下的CFD计算结果进行比较。对于本文采用的等截面三维CFD模型,升阻力系数可通过在reports中设置参考面积为10m2并定义升、阻力系数Cl、Cd相关参数后直接获得,结果如图 2.1所示。可以看出在0°至10°攻角下两种湍流模型的升力系数模拟结果非常接近,且与风洞试验结果相差不多。但是在大于10°的攻角下S-A模型的模拟结果明显偏离了试验值,相反的k-ω模型的模拟结果则相对较好。因此可以认为针对此类风力机翼型,更适合采用SST k-ω湍流模型对其进行CFD模拟。
2.2 体网格划分形式的影响
体网格划分得合理与否往往决定了CFD模拟结果的准确性,本节将着重考虑体网格的划分方式对风力机翼型气动力特性CFD模拟的影响。
2.2.1 多面体网格与六面体网格的比较
为了研究不同的体网格形式对CFD模拟结果的影响,本节分别针对1.1节中所提及的用非结构化多面体网格和结构化六面体网格划分出来的两种体网格模型进行模拟分析,湍流模型的选取上采用2.1节中效果较好的k-ω模型。计算结果如图 2.2所示,可以看出在0°至12.5°攻角下,两种划分网格方式的计算结果十分接近且与试验结果比较吻合。但是在15°至30°攻角下,以多面体网格划分的模型计算结果(升力系数Cl)呈现出先平稳后急剧下降的曲线,不论是数值还是曲线变化趋势都与实验结果差异较大;而以六面体网格划分的模型计算结果来看,则与试验相似的有升力系数在超过某临界攻角后随其增加而递减的趋势。不过在15°攻角下升力系数突然下降,即明显的失速状态,这是在风洞试验结果中没有体现出来的。为此对12.5°至17.5°之间的攻角做加密处理,以更好的反映出翼型气动攻角增加到一定程度时,翼型所产生的升力突然减小的失速状态,图 2.2中的局部放大部分就体现了这种状态。而升力系数骤降后能较快回升,并且随攻角增加只有小幅度下降的趋势,这说明此翼型具有良好的失速性能[6]。
2.2.2 网格质量好坏的影响
参照文献[3]对同一种风力机翼型做CFD模拟,且物理模型的参数设定均相同。但是发现由于所生成体网格质量的好坏存在差异,同样采用六面体网格对等截面三维叶片模型进行划分的情况下,模拟结果却有着较大的区别(如图 2.3所示)。文献[3]中的图 2显示划分的六面体网格在叶尖等位置较狭长,而本文通过合理设置边界层网格和局部位置加密的方式将体网格划分得比较均匀(接近正方体)。因此在网格总数偏少的情况下,本文的模拟结果相较于文献[3]仍更为准确。具体来看在0°至15°攻角下,本文的升力系数模拟结果更接近试验值,且在攻角大于12.5°之后与风洞试验结果一致有明显的下降趋势。而就阻力系数而言,本文的模拟结果已经与风洞试验结果基本吻合。这就说明划分的六面体网格若可以尽可能避免狭长的长方体形式,而多采用正方体或者近正方体的网格形式,就能在一定程度上提高CFD模拟计算结果的准确性。
2.3 叶片整体模型与等截面模型的比较
随着计算机技术的发展,CFD数值风洞技术在翼型空气动力特性预测中起着越来越重要的作用。但是目前风力机气动特性的预测,多还是通过对真实叶片各特征段的等截面模型进行风洞试验,进而叠加合成来完成。不过真实风力机的叶片不仅仅是变截面的,还存在扭角。风洞试验又或者是等截面模型CFD模拟得到的气动力参数能否应用于真实三维叶片,以及对实际工程的影响还有待进一步研究。本节对某公司提供的风力机专用真实整体翼型做CFD模拟,使用六面体网格进行划分,物理模型设定与1.1节一致。方便起见,就真实叶片根部位置按图 1.3规定将攻角定义为0°,取与1.1节中等截面翼型一致的截断面(距叶片根部30.5m处)进行对比。此时需要将等截面模型风向角调整为9.347°,比较两者的风压分布(如图 2.4所示)可以看出叶片整体模型与等截面模型的CFD模拟结果十分相近,所以将等截面模型的风洞实验结果应用于真实叶片的空气动力特性计算是比较合理的。
3. 结论
由于风力机翼型的几何外形复杂,所以目前工程实践中一般通过等截面三维模型进行风洞试验,或者通过CFD模型计算模拟来得到其气动力系数。然而风洞试验结果与CFD模拟结果往往存在着一定的差异,并且以上述方法获得的气动力系数是否能适用于真实叶片也有待研究。本文分别建立了非结构化多面体和结构化六面体两种体网格模型,并分别在S-A和SST k-ω湍流模型下对风力机翼型的气动力特性进行CFD模拟,还与风洞试验结果进行对比,得到的主要结论如下:
(1) SST k-ω湍流模型的计算结果相对于S-A模型与风洞试验值更为接近。虽然在0°至10°攻角下两种湍流模型的计算结果基本一致,但是在失速后S-A模型的模拟结果明显偏离了风洞实验值;
(2) 结构化六面体网格模型的计算结果相对于非结构化多面体网格模型与风洞试验值更为接近。同样的在0°至12.5°攻角下两者的计算结果也基本一致,但是采用结构化六面体网格能更好地模拟失速之后翼型升力系数的变化趋势,且数值也与风洞试验值更为相近;
(3) 所划分的结构化六面体网格若可以尽可能避免狭长的长方体形式,而多采用正方体形式,就能在一定程度上提高CFD模拟计算结果的准确性;
(4) 对比相同翼型截面的压力分布图,可以看出等截面三维模型与真实整体叶片模型的CFD模拟结果十分接近。所以综合上述结论,将等截面模型的风洞实验结果或等截面三维模型的CFD模拟结果应用于真实叶片的空气动力特性计算是比较合理的。
-
表 1.1 网格划分的参数设置
基本尺寸 边界层 叶片 外流域 目标尺寸 最小尺寸 目标尺寸 最小尺寸 多面体 总厚0.0065m 0.015m 0.005m 0.5m 1.25m 0.5m 共6层 六面体 递增比例1:1.6 0.015m 0.005m 0.25m 1.0m -
[1] 潘艺, 周鹏展, 王进.风力发电机叶片技术发展概述[J].湖南工业大学学报, 2007, 21(3):48-51. http://www.doc88.com/p-036414426255.html [2] Stone C P, Tebo S M, Duque E P N. Computational fluid dynamics of flat-back airfoils for wind turbine applications[C]//44th AIAA Aerospace Sciences Meeting and Exhibit, (Reno, NV), AIAA-2006-0194. 2006.
[3] 潘小涛, 黄铭枫, 楼文娟. 风力机翼型三维气动力特性CFD模拟[OL]. [2013-09-11]. 中国科技论文在线. [4] 韩占忠, 王敬, 韩小平.FLUENT流体工程仿真设计实例与应用[M].北京:北京理工大学出版社, 2004:98-121. [5] 尹耀安. 大型风力机叶片设计与CFD分析[D]. 华中科技大学, 2008. https://wenku.baidu.com/view/bdacb74ffe4733687e21aab6.html [6] 叶枝全, 黄继雄, 陈严, 等.风力机新系列翼型气动性能研究[J].太阳能学报, 2002, 23(2):211-216. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_tynxb200202018