失速限制型风轮叶片流场的数值模拟 Numerical simulation of stall regulated wind turbine.pdfVIP

第27卷第4期 空气动力学学报 V01．27。No．4 2009年08月 ACTAAERoDYNAMICASINICA Aug．．2009 文章编号t0258—182512009)04加405-06 失速限制型风轮叶片流场的数值模拟 肖中云，周 铸，陈作斌，刘 剐 (中国空气动力研究与发展中心．四川绵阳 621000) 摘要：采用双时间步方法在旋转坐标系下求解非定常方程，构造了零偏航条件下风轮叶片流场的数值模拟方法， 采用该方法对失速限制型风轮叶片NRELPhase-Ⅵ进行了模拟．其中网格采用重叠网格进行模拟，控制方程为可 压缩N．S方程．同时加入了低速预处理技术．模拟的风速范围从13m／s到25m／s．这时侯风轮叶片发生失速以限制 扭矩输出．本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析，从风轮扭矩、载荷分布等 方面对数值计算结果进行了对比验证。 关键诩：风轮叶片，零偏航；失速，旋转坐标系#低速预处理 中图分类号：TK8 文献标识码：A 机叶片的流场速度通常很低，就意味着非定常对应的 O 引 言 频率较高，因此只适合于用非定常方程进行模拟。 风能作为一种可再生的清洁能源，具有煤电等一 国内刘思永[2]等对风力涡轮的流场算法进行过 些传统能源开发不可比拟的环境优势，同时也具有大 研究，流动控制方程采用柱坐标系下的可压缩N—S 规模开发和商业化发展的前景。传统的风轮机设计 方程，考虑的来流风速为85m／s。实际上自然风要远 都采用基于叶素理论、动量理论等一些工程方法，相 低于这个速度，流场中存在着大片的不可压区域，在 比之下CFD方法具有更高的精度，也是适合于设计 采用可压缩方程求解时会碰到很强的刚性问题。这 未来大型风轮设备的主流方法。随着计算技术的发 一问题可以通过低速预处理方法加以解决，具体来说 展，CFD技术逐渐开始在风轮叶片的优化设计、性能 是向控制方程中加入预处理矩阵，使方程满足条件数 预估分析甚至整机载荷评估方面起到重要作用；在相 为1的量级，从而消除系统刚性，使可压缩方法在计 关的流动机理研究方面，用CFD方法研究与非定常 算低速流场时既保持精度，同时又具有很好的收敛性 流动有关的动态失速，动态载荷也得到了越来越多的 质。低速预处理矩阵有多种不同形式[3]，大致可以分 关注。 为Choi—Merkle、Turkel和VanLeer型预处理矩阵。 NREL Phase．VI风机是美国再生能源实验室的 第六阶段的风机试验模型’，采用失速限制的方法对能 阵演化而来，和Turkel的预处理矩阵又具有类似的 量输出进行控制，其原理是当风速超过风机的额定风 形式，矩阵的构造相对简单，该方法在工程数值模拟 速以后，从叶片根部开始会出现失速，以此限制扭矩 中得到了较广泛的应用。 输出和保持发电系统的稳定。风轮叶片类似于直升 零偏航条件下的风轮叶片流场具有旋转对称性， 机旋翼，其不同之处在于翼型厚度和弯度较大，容易 通过在非惯性旋转坐标系下求解方程可以使计算得 形成气流分离、再附等复杂的流动现象。气流分离使 到简化。本文的计算条件取自NRELPhase-Ⅵ实 流场的非定常特征十分明显，可以用如下减缩频率来 验，来流风速范围为13m／s～25m／s，对应了风轮叶 定量表示[1]，七=眦／2y，其中∞表示非定常频率，c代片失速以后的流动状态，本文通过数值模拟对风轮叶 表特征长度(用叶片的弦长表示)，y表示平均来流速 片的绕流特点、扭矩和载荷分布等方面进行了分析， 度。文献[1]