基于ANSYS的风力机双向流固耦合分析.docVIP

小型风力机风轮流线分析研究 杨柳1，汪建文1,2，马剑龙1，魏海姣，张彦奇1，李佩林1 内蒙古工业大学 能源与动力工程学院，内蒙古 呼和浩特 010051；2.风能太阳能利用技术省部共建教育部重点实验室，内蒙古 呼和浩特010051） 摘要：基于非稳态算法针对某一小型风力机风轮进行风轮流线分析，研究发现随着距离风轮中心30%R、60%R、90%R截面变化翼型前缘处速度越来越大；风轮叶尖处流线密度最大，并且风轮流线有明显内卷的现象。距离风轮中心截面30%R翼型尾部出现了明显的失速涡，可得风轮根部会造成大量的风量损失，主要做功截面分布于叶片中部以及叶梢。随着距离风轮中心30%R、60%R、90%R0 引言 风力机运行于开放的大气环境中， 图1风轮模型 2计算域的建立与设置 利用流场计算软件CFX进行非稳态流场模拟，由于风力机表面是无滑移，无渗透，不可压缩风场是有旋流场运动，加之FSI面速度压力分布均相等[7]，故将计算域分为两个域：静止域与旋转域，静止域的主要作用是用于进出口风速和风压的设置，其网格类型为普通的滑移网格；由于风力机旋转过程中旋转域的网格会发生变形，故旋转域的网格类型设定为变形网格，用于与风力机叶片与外界环境的互相传递。如图所示2计算域。 对于ANSYS的CFX中非稳态的设置主要包括出口压力的设置、入口风速设置、静止域和旋转域设置、交界面设置。静止域的进口风速设置为：8m/s，出口相对压力为：0；旋转域的旋转速度的为：57.14rad/s，同时FSI面必须在旋转域中设置，并且FSI面为风力机的所有叶片表面。计算域设置如图2计算域设置所示。 图2 计算域 3结果分析与处理 3.1风轮流线分布 经过非稳态算法计算模拟可得风轮流线分布情况如下。 （a）风轮侧面流线图 （b）风轮垂直迎风面流线图 图3 风轮的流线图 从3(a)风轮侧面流线图中可以看到在风力机叶尖处流线密度最大，因此可得叶尖周围的流体流动速度最大，反之可知，风力机叶尖处的运动速度最大，从3(b)风轮垂直迎风面流线图，可以观察到有明显流线内卷的现象，可能由于在风力机叶片根部的流动现象比较复杂，同时出现分流的现象，这种流线内卷现象会造成风能损失。 3.2风轮叶片三截面速度、流线云图分析 为了比较分析距离风轮不同截面的速度、流线的分布，在风力机同一叶片下、相同攻角处，分别选取叶片展向方向上距离风轮中心三个截面分别是30%R、60%R、90%R，如图8所示为三截面处速度、流线分布云图。 如图所示4速度分布云图可以观察到在三个截面中随着沿叶展方向距离风轮中心30%R、60%R、90%R的变化前缘处速度越来越大。 30%R 60%R 90%R 图4 截面速度分布云图 30%R 60%R 90%R 图5截面流线分布云图 在图5流线分布云图中从距离风轮中心30%R的截面中可以看到翼型尾部出现了明显的失速涡，因此风轮主要做功截面分布于叶片中部以及叶梢，根部会造成大量的风量损失，从而可以得到在此处发生了较大的分离现象。在距离风轮中心60%R、90%R两截面中，可以看到靠近翼形表面的流线曲率半径逐渐减小，根据公式，其中为沿流线法线方向的积分常数，为流线曲率半径，为流体速度，所以沿叶片展向，流体速度越来越大，从而又可以得出叶梢处的速度达到最大，从而看出与图4的结果相符合。同时也可以观察到随着距离风轮中心30%R、60%R、90%R三个截面的变化，尾缘分离的现象越来越弱，流线逐渐趋于紧密，根据流体力学的分析，发生尾缘分离现象主要是由于逆压梯度的产生。 4结论 利用非稳态算法针对直径为1.4m的小型水平轴风力进行风轮尾迹流场分析（速度、流线）。以来流风速为8m/s ，风轮转速为57.14rad/s为例，分析风轮速度、流线变化规律： 三个截面中随着叶展方向距离风轮中心30%R、60%R、90%R截面变化翼型前缘处速度越来越大。 通过风轮流线分布图可知风轮叶尖处的流线密度最大，并且风轮流线有明显内卷的现象，造成能量损失。 从距离风轮中心30%R截面中可以看到翼型尾部出现了明显的失速涡，因此可得风轮主要做功截面分布于叶片中部以及叶梢，根部会造成大量的风量损失。 随着距离风轮中心30%R、60%R、90%R三个截面的变化，尾缘的分离现象越来越弱，并且该流线逐渐趋于紧密。 参考文献： [1] Arakawa C, Fleig O, Iida M, et al. N