风力机翼型等速上仰动态失速数值模拟.docx

风力机翼型等速上仰动态失速数值模拟作者简介：周正（1989-）男（汉）硕士研究生摘要：采用k-ωSST模型利用CFD软件模拟了NRELS809翼型正弦振荡动态失速并将结果和俄亥俄州立大学（OSU）风洞试验值对比显示出较好的一致性验证了所用方法的有效性.在此基础上对该翼型在雷诺数Re=1.0×106时以攻角变化率α・=34.54（°）・s-1等速上仰动态失速过程进行了数值模拟详细描述了等速上仰动态失速过程涡的发展以及翼型周围流场的分布.结果表明动态失速现象是由前缘主涡和尾缘逆向涡交替作用引起；其气动特性曲线的分析结果表明其失速前气动性能较静态时有较大提升.关键词：S809翼型；等速上仰；动态失速；数值模拟中图分类号：TK83文献标志码：A对于水平轴风力机风剪切、垂直风、大气湍流、塔影效应以及偏航运行等因素使得风轮叶片的来流攻角快速变化表现出与风洞实验中完全不同的气动特性升力系数和俯仰力矩系数均出现明显差别[1].翼型作为风力机叶片的基本要素显著影响后者气动性能.为充分利用充满能量的流动非常有必要理解和计算各种条件下的动态失速过程[2].许多学者进行了相关研究[3-6].准确计算翼型在动态失速时的气动特性系数对于提高风力机性能及确保风力机安全设计都具有十分重要的意义.为此本文选用NRELS809翼型对其进行等速上仰动态失速数值模拟.1数值计算方法1.1控制方程1.2数值模拟方法本文数值模拟计算域如图1所示长度为弦长的34倍宽度为弦长的28倍充分考虑了尾迹效应的影响.在Gambit2.2中对翼型计算域分块生成非结构网格内部圆形旋转部分较密外部静止部分相对较疏.图2为翼型边界层网格前缘部分视图其周围共布置500个节点边界层首层厚度设为0.0001m增长因子为1.08无量纲壁面高度y+值小于1以捕捉边界层附近流动现象.2计算结果及分析2.1气动力系数曲线2.2流场及动态失速涡分析图6反映了等速上仰动态失速时不同攻角下翼型周围流线图和涡量图.攻角为6°和10.01°时翼型周围为附着流动状态升力系数随着攻角几乎线性增大.攻角增大到13.67°时升力系数达到第一个峰值翼型尾缘处开始产生逆向小涡.攻角继续增大到20.93°过程中尾缘小涡的发展使得升力系数逐渐减小并在20.93°处前缘一次涡开始诱导出二次涡此后直至21.83°二次涡快速发展（如攻角21.43°所示）与一次涡合并成一主涡同时尾缘逆向涡逐渐脱落两者共同作用使得升力系数达到最大值.攻角仰至22.28°时主涡逐渐远离翼型表面流动发生大范围分离并在尾缘处诱导出二次逆向小涡翼型气动特性严重下降阻力系数和力矩系数快速增长.攻角继续上仰到22.93°前缘二次涡诱导的三次涡快速发展尾缘逆向涡渐渐脱落使升力系数增加.并在攻角仰至23.54°时二次涡和三次涡合并成一充分发展的主涡在尾缘处再次诱导出逆向小涡升力系数又一次到达峰值.此后的上仰过程基本与上述过程一致不同之处在于前缘未能诱导出三次涡.由此可知：等速上仰翼型动态失速现象为前缘主涡和尾缘逆向涡交替产生、发展、脱落的过程其中伴随着气动特性系数的快速变化主导着失速后翼型的气动性能.3结论（1）采用k-ωSST模型计算所得结果与OSU风洞试验值吻合较好验证其可有效模拟动态失速过程.（2）翼型等速上仰动态失速扩展了升力系数线性增长段延迟了失速攻角失速前气动性能较静态时有较大提升.（3）等速上仰动态失速为前缘主涡和尾缘逆向涡交替作用的过程涡的发展密切影响气动特性曲线的变化.参考文献：[1]李春叶舟高伟等.现代陆海风力机计算与仿真[M].上海：上海科学技术出版社2012：179-290.[2]LEISHMANJG.DynamicstallexperimentsontheNACA23012aerofoil[J].ExperimentsinFluids19909（1/2）：49-58.[3]SHIHCLOURENCOLMVANDOMMELENLetal.Unsteadyflowpastanairfoilpitchingataconstantrate[J].AIAAJournal199230（5）：1153-1161.[4]LEISHMANJGBEDDOESTS.Asemiempiricalmodelfordynamicstall[J].JournaloftheAmericanHelicopterSoc