基于cfd的定桨距垂直轴风力机变桨距控制研究

基于cfd的定桨距垂直轴风力机变桨距控制研究

新能源的开发和利用受到国内外科学家的高度重视。风能作为清洁可再生能源,已成为最具潜力的替代能源之一。随着研究的深入以及风力发电技术的不断发展,垂直轴风力机凭借噪音小、安装维护方便以及翼型制造容易等优势,成为关注的焦点,并具有了与水平轴风力机相媲美的潜力。然而,垂直轴风力机启动性能差甚至无法自启动,以及风能利用率相比水平轴风力机较低等问题,极大限制了垂直轴风力机的发展和推广。研究表明,变桨距技术是解决垂直轴风力机存在问题的有效方法。Pryor发明了一种应用W型机构的垂直轴风力机变桨距方法,Pawsey设计了一款基于气动力与惯性力原理的离心锤式变桨距垂直轴风力机。Dress研制了“Cycloturbine”风力机,采用凸轮机构实现了变桨距控制,通过理论分析发现该类风力机的最大风能利用率可达0.45。张亮等研究了桨距角对水轮机性能的影响。张松等则设计了一款同步变桨距垂直轴风力机,由1台伺服电机实现多组叶片的变桨距控制,样机在尖速比0.6附近获得最大风能利用率。Schönborn等研究了一种基于正弦规律的变桨距水轮机,该水轮机在尖速比3附近获得最大发电效率,约为38%。Hwang等提出了具有摆线规律的可独立变桨距控制的水轮机,变桨距控制下的性能相比定桨距水轮机有70%的提高。虽然变桨距垂直轴风力机已经得到初步发展,但变桨距控制规律却一直处于实验室阶段。本文作者针对风力机变桨距控制规律展开深入研究,以最大瞬时功率为目标函数设计桨距角优化程序,并在CFD中对设计的变桨距规律进行验证,分析变桨距控制下风力机的气动性能。1风力机性能分析为了获得变桨距控制规律,需要以叶片的气动载荷和风力机的气动性能为基础建立风力机的数学模型,分析桨距角与风能利用率的关系。由于流管理论可以有效地预测风力机的整体性能,是一种很好的辅助风力机设计的模型,故采用双盘面多流管理论研究风力机的性能。1.1垂直轴风力机个体来流速度模型双盘面多流管模型由Paraschivoiu提出,其在流管理论的基础上又将每个流管分割为上、下2个盘面,可以实现上、下2个盘面的不同性能分析,如图1所示。依据流管理论,风力机局部风速为:式中:V∞为无穷远处来流风速;Vu和Vd分别为上、下盘面局部风速;κu和κd分别为上、下盘面速度因子。现定义变量κθ为垂直轴风力机0°~360°范围内的速度因子,式(1)可以简化为:其中:当叶片位于上盘面时κθ=κu;当叶片位于下盘面时κθ=κd(2κu-1)。垂直轴风力机的相对速度W与局部风速Vθ的关系如图2所示。式中:U为风力机切向速度,且U=ω×R。风力机局部尖速比λθ可以表示为:由图2可知风力机的来流角γ为根据动量理论,单流管内的轴向力为其中:ρ为空气密度;AS为流管面积。基于叶素理论,单流管内作用在叶片上的气动力为则联立式(7)和(9),采用数值解法即可得到风力机的速度因子。1.2叶片桨距角度为分析叶片桨距角与风能利用率的关系,建立了叶片的气弹动力学模型,如图3所示。作用在叶片上的气动载荷简化为力矢量FA和相对于简化中心A的力偶MA,点P和点G分别为叶片的转动中心和重心。如图3所示,升力Fl与阻力Fd的关系为并定义垂直轴风力机的切向力系数Cτ为式中:Cl为升力系数;Cd为阻力系数。值得注意的是:简化中心A随方位角θ变化而变化,并记作rA。由简化中心位置变化导致的额外气动力矩会影响叶片桨距角,进而影响风力机的气动性能,严重时会导致叶片颤振。因此,有必要分析作用在叶片转轴处的俯仰力矩,即由式(12)可知:叶片转轴处的俯仰力矩会对叶片桨距角产生影响,这样就会引入额外的桨距角βex,其中:Kα为叶片转动中心的刚度。则叶片实际的桨距角βa为单叶片产生的瞬时力矩Mθ为单叶片的力矩系数Cm为:则变桨距垂直轴风力机的瞬时功率Pθ为式中:N为叶片数量;Mθi为第i个叶片产生的瞬时力矩,且i=1,2,…,N。变桨距垂直轴风力机的风能利用率Cp可以表示为由式(15)和式(17)可知:垂直轴风力机的瞬时功率是关于升力系数,阻力系数以及局部尖速比λθ的函数,而叶片桨距角直接决定了叶片的气动性能,故可通过改变桨距角方法以实现最大功率输出,提高垂直轴风力机的风能利用率。2垂直轴风力机cfd仿真研究为了获得不同尖速比下的桨距角变化规律,先对定桨距垂直轴风力机进行CFD仿真研究,分析不同桨距角(对于定桨距垂直轴风力机即为安装角)下的气动性能以及桨距角与风能利用率的关系,进而得到桨距角变化趋势。2.1平均数据处理对于直叶片垂直轴风力机,在不考虑叶片端部湍流或者叶片无限长的情况下,翼型周围的流场很相似,可以将三维流场简化为二维流场。根据质量守恒,动量守恒以及连续性理论建立流场模型,采用二维不可压缩N-S方程作为控制方程,并对流体方程进行了RANS雷诺平均处理,即瞬时速度ui和压力p分别为各自时均值和脉动值的和。在笛卡儿坐标系下流体方程表示为式中:项为雷诺应力,并将速度脉动的二阶关联量表示成平均速度梯度与湍流黏性系数的乘积,即式中:μt和k分别为动态粘涡系数和湍流动能。由Howell等的研究发现,重整化群的k-ε模型修正了耗散率方程,保证风力机的数值模拟有较高的精度,故采用k-εRNG模型,即由于k-εRNG模型可使得流体方程封闭,有效解决了高雷诺数问题,而低雷诺数问题则采用壁面函数法处理。同时采用压力-速度耦合以及标准壁面方程的方法,研究大攻角下流场的分离情况。本文建立了4叶片垂直轴风力机的CFD模型,风轮半径为2m,叶片展长为2m,采用NACA0015对称翼型,弦长为400mm,来流风速为10m/s,并采用速度进口、压力出口和无滑移壁面等边界条件。2.2有效攻角的影响由式(17)可知:翼型的气动参数决定了风力机的最大瞬时功率,获得翼型的气动参数很重要。对NACA0015翼型进行CFD仿真研究,获得了翼型的静态升阻力系数,如图4所示。由图4可知:随着弦长的增加,升力系数Cl增加,阻力系数Cd下降,翼型的升阻比变大,有助于风能利用率的提高;随着弦长C的增加,阻力系数曲线右移,表明翼型的静态失速角增加,削弱了翼型边界层的分离;当翼型失速后,压差阻力相比摩擦阻力起主导作用,阻力系数急剧上升;对于升力系数,失速前升力系数随有效攻角变大而变大,失速后升力系数快速下降;最大升力系数随弦长增大而增大。经CFD研究还得到了不同攻角下翼型的压力分布,如图5所示。由图5可知:翼型前缘附近的负压梯度对升力的形成贡献很大;在低压区形成吸力使翼型绕叶片转轴从高压区转向低压区,即负压形成了负升力。这表明与负升力对应的有效攻角为负值。2.3垂直轴风速对风力的影响在获得翼型的气动性能之后,还需研究叶片桨距角与风能利用率的关系,以便规划变桨距控制规律。由于尖速比决定了风力机的运转特性,直接影响风能的利用,所以,研究风速与尖速比的对应关系也很重要。经CFD仿真研究,得到不同风速下尖速比与风能利用率的关系,如图6所示。由图6可知:随着风速的增大,风能利用率的峰也增大,且风能利用率峰值所对应的尖速比也增加,说明高尖速比下的风能利用率较高,且风力机尖速比在2.5~3.5范围内具有较高的风能利用率。此外,在尖速比小于2的情况下,风力机启动后的风能利用率差别不大,垂直轴风力机在尖速比λ<2.5的区间内的气动性能有待改善。固定垂直轴风力机叶片桨距角进行CFD仿真,通过改变桨距角得到了桨距角与风力机性能的关系曲线,如图7所示。由图7可以看出:桨距角为-4°时的风能利用率最高,并且具有较宽的有效工作范围。经分析可知桨距角-4°是定桨距风力机性能的转折点,风能利用率呈现先增加后下降的趋势;当桨距角由2°向-4°变化时,风能利用率峰值朝低尖速比方向移动,且对应的最大风能利用率提高。由此可知:外偏置桨距角可以增大风能利用率的峰值和有效作用范围。现对桨距角进行如下解释:垂直轴风力机的运行轨迹为圆,若翼型前缘位于圆内则定义为内偏置桨距角,且为正值;若翼型前缘恰在圆上,则桨距角为0°;否则为外偏置桨距角,且为负值。图8所示为定桨距垂直轴风力机在不同桨距角下的力矩系数。由图8可知:上盘面相比下盘面的力矩系数较大。这是因为在上盘面大部分风能转化为机械能,下盘面风速几乎是上盘面风速的一半,加之风能与风速的三次方呈正比,所以,下盘面的力矩系数很小。从图8还可以看出:当桨距角为正值时,风力机在较小的方位角就进入有效力矩作用范围,也最先退出有效力矩作用区;当桨距角为负值时,相比正值桨距角较晚进入有效力矩作用范围,且随着桨距角幅值的增大,相应的方位角也增大。由图8可知:桨距角为-4°时力矩系数较高且有效力矩作用区较宽,具有较好的气动性能;当桨距角为正值时,风力机进入有效力矩作用区时的力矩系数基本一致;当桨距角为负值时,风力机退出有效力矩作用区时的力矩系数差别也较小。由此可得出:设置内偏置桨距角或者小幅值的外偏置桨距角,使风力机在小方位角进入有效力矩作用区;且设置外偏置桨距角使风力机延缓退出有效力矩作用区,通过该方法可以实现风力机气动性能的提高。3修改桨距离检测规律和模拟分析3.1不同尖速比下的桨距角变化规律图9所示为确定垂直轴风力机桨距角变化规律的程序框图。其原理为:根据结构参数、流场和运行条件以及气动参数求解速度因子;在考虑局部风速和局部尖速比的情况下,根据式(17)建立风力机瞬时功率;对初始气动参数进行修正,求解瞬时功率最大值及对应的桨距角,最终得到风力机0°~360°范围内的桨距角控制规律。应用MATLAB软件编写优化程序,获得不同尖速比下桨距角变化规律,如图10所示。图10所示为不同尖速比下的桨距角控制规律。图10中曲线关于180°中心对称,与风力机流场特性基本吻合;随着尖速比的增加,桨距角变化量减小,表明高尖速比下是否采用变桨距控制还需全局衡量。3.2垂直轴风力机的变桨距控制由式(11)可知:切向力系数在一定程度上可以反映风力机性能,经CFD仿真得到了定桨距和变桨距风力机切向力系数,如图11所示。由图11可以看出:无论是定桨距还是变桨距风力机,都在上盘面取得切向力系数的峰值,而下盘面的切向力系数很小。这是因为在上盘面流体不会或很少发生分离,湍流损失较少,大部分气动力转化为升力;而下盘面流体速度较低,或因忽略下盘面流线方向变化以及尾涡影响,致使下盘面的切向力系数较小。此外,与定桨距垂直轴风力机相比,在变桨距控制下风力机在方位角90°~180°范围内的性能得到了改善,风力机切向力有效作用范围约占风力机整周的1/2,而定桨距垂直轴风力机却只占有1/4圆周。结果表明该主动式变桨距控制规律显著提高了垂直轴风力机的整体性能。另外,经CFD仿真还得了作用在叶片上的俯仰力矩,在变桨距控制下叶片的俯仰力矩与定桨距风力机的对比结果,如图12所示。由图12可知:定桨距垂直轴风力机上盘面俯仰力矩变化剧烈,而在变桨距控制下作用在叶片上的俯仰力矩变化比较平缓。这是因为变桨距控制提高了叶片的动态失速角,有效地抑制了动态失速现象的发生,而且叶片桨距角的动态调节降低了压力梯度。在变桨距控制下,作用在叶片上的俯仰力矩大部分是负值,且与叶片运动方向一致,对垂直轴风力机的变桨距控制有利,可以降低由伺服电动机在主动变桨距控制下功率消耗;而作用在定桨距垂直轴风力机上的俯仰力矩是必须要克服的阻力矩。图13所示为作用在风力机上的合力矩系数。从图13可见:与定桨距垂直轴风力机相比,变桨距控制下的合力矩系数较大,并且力矩系数波动较小,表明变桨距控制可以提供更加稳定的功率输出;叶片在相同方位角下具有一致的运动规律,且呈周期性变化,说明合力矩系数只与方位角有关,与具体哪个叶片起作用无关。图14所示为变桨距控制下垂直轴风力机的风能利用率。由图14可知:随着尖速比的增大风能利用率提高,且在尖速比为2.6附近取得峰值。主动式变桨距控制极大地提高了尖速比λ<2.6区间下的风能利用率。在尖速比λ=2,风速为10m/s的情况下,在变桨距规律控制下垂直轴风力机的风能利用率达到39.12%,而定桨距风力机(β=0°)的风能利用率只有19.80%,风力机利用率提高近2倍。对于定桨距垂直轴风力机,桨距角为-4°时具有较高的风能利用率,但也