风力机叶片三维数值计算方法确认研究 the validate and research on 3d numerical simulation of the aerodynamic perormance of wind turbine blades.pdfVIP

第3l卷第3期 太阳能学报 V01．31．No．3 2010年3月 ACrAENERGIAESOLARISSINICA Mar．．20lO 风力机叶片三维数值计算方法确认研究 范忠瑶1，康 顺1，王建录2 (1．华北电力大学电站设备状态监测与控制教育部重点实验窀，北京102206； 2．东方电气集团东方汽轮机有限公司，德阳618000) 摘要：采用CFD软件包FINETM／TURBO，以两叶片NREI。PhaseVI风力机的风轮为对象，进行了风力机风轮叶片 三维绕流的定常粘性数值模拟研究。通过详细对比计算结果与实验数据(包括功率、叶片展向5个截面压力系数 分布及沿叶展方向载荷系数分布)，确认在大部分风速条件下数值模拟可以很好的预计风力机气动性能。然后分 析r计算域尺度、边界条件和湍流模型等对数值模拟结果的影响，为采用CFD技术对实际风力机叶片三维气动性 能进行精确数值模拟提供参考。 PhaseVI 关键词：风力机；叶片；数值模拟；NUMECA；NREL 中图分类号：TK83 文献标识码：A O引言 度的非定常空气动力学(UAE)实验【10|，因实验数据 众所周知，旋转的风力机叶片周围存在复杂的 全面、可靠而被国内外学者广泛应用于数值模拟技 三维粘性流场。基于动量．叶素理论及其修正的方 术的验证和发展。其风轮由两叶片组成，直径为 法在附着流动时可以得到较好的气动性能计算结 10．058m。采用变截面扭叶片，失速控制方式，额定 果。然而随着风速的提高，流动分离增大，风轮旋转 功率19．8kW，转速72r／min。实验风力机做了改装， 的三维流动效应更加明湿。在离心力和科氏力的作 进行了上、下风向，轮毂刚性、翘板型连接等35种风 用下…，径向的二次流动增强。这时，基于二维翼型 力机结构在不同风速条件下的1700多个实验。本 气动特性的计算模型需要结合流动机理进一步改 文选择上风向、轮毂刚性连接结构、拆除了五孔探针 进[2’3J。这方面的实验也由于受到各种条件的限制 以减小对流场扰动的UAE—S系列实验模型进行数值 而不易展开。因此，目前CFD方法已成为分析叶片 计算对比分析。 气动性能的主要工具。研究人员从不同角度展开了 广泛研究，并获得了许多有益结果【4圳。然而，要将 2数值方法 CFD直接应用于工程实际并获得可信赖的数值计算 结果，还存在合理的计算域尺度与边界条件的选取、 能够满足计算精度要求的网格质量与网格数目选择 以及合适的湍流模型选择等问题。 程，中心节点的有限体积法离散，显式Runge—Kutta Phase 法求解，全多重网格初场处理和多重网格迭代加速 本文以NREL VI风力机为例，在确认数值 计算可靠性的基础I二，分析了改变计算域尺度、计算 以及低速流动的预处理技术等。选用的湍流模型为 域中转动部分的尺度、边界条件和湍流模型等对数 Stress 值计算结果的影响。 k．(1，(ShearTransport)(以下简称SsT)模型。 2．1网格自动生成 1研究模型