风力机翼型动态失速数值模拟 investigation of dynamic stall of wind turbine airfoil.pdfVIP

第33卷第10期 VoL 太阳能学报 33．No．10 20]2年10月 ACTAENERGIAESOLARISSINICA Oct．。2012 文章编号：0254-0096(2012)10—1641--07 风力机翼型动态失速数值模拟 满洪海，杜 刚，陈 江 (北京航空航天大学航空发动机数值仿真中心，北京100191) 度、湍流模型对俯仰振荡翼型数值模拟的影响，并与实验结果进行对比验证。数值模拟得出的气动力迟滞曲线与 实验结果变化趋势符合较好，表明了数值方法的有效性。采用此方法研究了北航自主开发的风力机专用翼型的动 态失速特性。通过分析表明该翼型在动态失速过程中升力系数下降平缓，俯仰力矩振荡变化不显著，适合应用于 兆瓦级大型风力机叶片。 关键词：动态失速；风力机；翼型；数值模拟 中图分类号：TK89 文献标识码：A k-oJ模型)，空间、时间离散格式均采用二阶迎风格 0 引 言 式。为提高收敛速度和计算精度，选用基于密度的 动态失速是指翼型在做非定常运动时，其扰流 耦合隐式求解器。 流场的分离和失速发生超过其临界迎角，形成动态 非结构网格与结构网格相比，舍弃了网格连接 延迟这一复杂的非定常气动现象。根据粘性作用的 性和正交性的限制，且具有易于控制网格单元大小、 大小，动态失速可分为轻失速和深失速。轻失速不 形状及网格点位置等诸多优点【2]。因此非结构网格 仅表现出一般静态失速的特征，还表现有非定常分 具有更大的几何灵活性，可有效防止网格运动过程 离的强粘性或无粘相互作用的性质，此时边界层厚 中的扭曲和变形，对于复杂外形和运动外形具备更 度为翼型厚度的量级；深失速则表现为全粘性现象， 好的普遍适应性。基于这些特点，本文数值计算中 存在高度非线性的压力脉动，在翼型表面上有大涡 使用的网格是采用ANSYS 运动，此时边界层厚度可达到翼型弦长的量级¨。。 网格，外边界距离壁面大于等于30倍弦长。 目前，国外对于翼型的动态失速性能无论是数 网格的运动采用弹簧近似法和局部网格重构法 值计算还是实验均有广泛研究，但国内针对风力机 相结合的办法，使用弹簧近似法时，网格运动过程 翼型的动态失速研究还较少。本文旨在通过数值方 中，其拓扑始终不变，无需插值，保证了计算精度。 法模拟低雷诺数下，NACA0012翼型和风力机翼型而局部网格重构法主要对运动翼型边界附近区域的 s809的动态失速特性。研究不同网格密度、湍流模网格起作用，它对复杂的外形具有很好的适应性。 型对数值模拟的影响，并与实验结果进行对比，验证 翼型绕1／4弦线轴做俯仰振荡，迎角变化规律 数值方法的有效性。采用此方法研究北航自主开发 为： 的风力机专用翼型的动态失速特性。 理=ao+△d×sin(2Kt’) 式中，d。——平衡迎角；△d——振幅角；卜减缩 1 研究方法 频率；t’——无量纲化时间，t+=tU。／c。此外本文 采用商用CFD软件进行二维翼型的数值模拟。 中雷诺数均基于无穷远处绝对速度u。和翼型弦长 控制方程组为Navier．Stokes方程，采用3种湍流模c，定义为Re=U。c／v。 型(一方程s—A模型、两方程标准k-e模型和标准 计算域如图1所示，由翼型表面和远场边界组 收稿日期：2010-09—10 万方数据 太 阳 能 学 报