尾缘襟翼偏转角对不同翼型的垂直轴风力机气动影响研究

doi:10.16183/j.cnki.jsjtu.2022.110

尾缘襟翼偏转角对不同翼型的垂直轴风力机气动影响研究

戴孟祎¹, 张志豪¹, 涂佳黄², 韩兆龙^,¹, 周岱¹, 朱宏博¹

1.上海交通大学船舶海洋与建筑工程学院, 上海 200240

2.湘潭大学土木工程学院, 湖南湘潭 411105

Aerodynamic Effect of Deflection Angle of Trailing Edge Flap on Vertical Axis Wind Turbine with Different Airfoils

DAI Mengyi¹, ZHANG Zhihao¹, TU Jiahuang², HAN Zhaolong^,¹, ZHOU Dai¹, ZHU Hongbo¹

1. School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

2. College of Civil Engineering, Xiangtan University, Xiangtan 411105, Hunan, China

通讯作者: 韩兆龙,男,副教授,博士生导师;E-mail:han.arkey@sjtu.edu.cn.

责任编辑: 孙启艳

收稿日期: 2022-04-15

基金资助:

国家自然科学基金项目(42076210)
国家自然科学基金项目(51879160)
国家自然科学基金项目(52088102)
上海市科委重大研究计划(2019-01-07-00-02-E00066)
湖南省自然科学基金项目(2021JJ50027)
湖南省教育厅科学研究项目(21A0103)
政府间国际科技创新合作重点专项(2018YFE0125100)

Received: 2022-04-15

作者简介 About authors

戴孟祎(1997-),女,湖南省长沙市人,硕士生,从事垂直轴风力机相关研究.

摘要

风能转化率偏低是阻碍垂直轴风力机市场化发展的重要原因.尾缘襟翼的设计能够改变叶片表面的流场结构,从而提高垂直轴风力机的气动性能.目前关于不同翼型垂直轴风力机的气动性能随尾缘襟翼的变化规律尚不明确.基于计算流体动力学方法,采用转捩剪切应力输运湍流模型,对3种不同分离式尾缘襟翼的翼型(NACA0018、NACA0021和NACA0024)叶片的H型垂直轴风力机气动性能进行数值研究.验证算例与已有的实验结果对比,结果吻合较好,证实本方法的可靠性.进一步考虑3种基础翼型与5组襟翼偏转角(-16°、-8°、0°、8°、16°)参数,探究垂直轴风力机的气动性能差异,分析其内在机理.研究结果表明:逆风区正向襟翼偏转角可以有效提高叶片的弯矩系数,顺风区负向襟翼偏转角对叶片的弯矩系数产生有利影响.在负向襟翼偏转角下,风能利用率受偏转影响的程度与翼型厚度呈正相关;在正向襟翼偏转角下,风能利用率受偏转影响的程度与翼型厚度呈负相关.研究成果可以为垂直轴风力机尾缘襟翼的应用提供有效参考.

关键词： 垂直轴风力机; 计算流体动力学; 尾缘襟翼; 翼型; 气动性能

Abstract

Low power efficiency is a critical factor that restricts marketization development of the vertical axis wind turbine (VAWT). The proposal of the trailing edge flap can change flow structure on blade surface, so as to improve the aerodynamic performance of VAWT. At present, the variation law of aerodynamic performance of different airfoil VAWT with trailing edge flaps is not clear. Based on the computational fluid dynamics (CFD) method and the shear stress transport (SST) model, a numerical simulation of 3 H-type VAWTs with different airfoils (NACA0018, NACA0021, and NACA0024) with separated trailing edge flap is conducted. It is found that the results of the validation case are in good agreement with experimental results, which verifies the reliability of this method. Afterwards, 3 basic airfoils and 5 groups of flap deflection angle (-16°, -8°, 0°, 8°, and 16°) parameters are selected to explore the difference in the aerodynamic performance of VAWTs. The results indicate that the positive flap deflection angle in the upwind region can effectively improve blade moment coefficient, and the negative flap deflection angle in the downwind region has a beneficial effect. For the negative flap, the degree of wind energy utilization affected by deflection is positively correlated with airfoil thickness, while for the positive flap, the opposite is true. The research results of this paper can provide an effective reference for application of trailing edge flaps of vertical axis wind turbines.

Keywords： vertical axis wind turbines (VAWT); computational fluid dynamics (CFD); trailing edge flap; airfoil; aerodynamic performance

PDF (5904KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

戴孟祎, 张志豪, 涂佳黄, 韩兆龙, 周岱, 朱宏博. 尾缘襟翼偏转角对不同翼型的垂直轴风力机气动影响研究[J]. 上海交通大学学报, 2022, 56(12): 1619-1629 doi:10.16183/j.cnki.jsjtu.2022.110

DAI Mengyi, ZHANG Zhihao, TU Jiahuang, HAN Zhaolong, ZHOU Dai, ZHU Hongbo. Aerodynamic Effect of Deflection Angle of Trailing Edge Flap on Vertical Axis Wind Turbine with Different Airfoils[J]. Journal of Shanghai Jiaotong University, 2022, 56(12): 1619-1629 doi:10.16183/j.cnki.jsjtu.2022.110

近年来,风能以其清洁的特性、丰富的储量以及日趋成熟的风力发电技术,成为未来最具前景的新能源形式之一^[1].风力机作为捕获风能的重要装置,其构型对于提高风力机的风能转化率、降低生产和运行维护成本具有重要意义^[2].与传统水平轴风力机相比,垂直轴风力机具有设计简单、对风向不敏感等优点^[3],但垂直轴风力机较低的风能利用率成为了其大规模商业化应用的瓶颈.

在旋转过程中,垂直轴风力机叶片的攻角会出现周期性变化,当叶片的攻角超过失速攻角时,叶片表面气体就会流动分离,从而使得叶片的气动性能急剧下降^[4].因此,研究人员提出利用不同的流动控制手段改善叶片的气动性能.流动控制手段主要可分为被动控制^[5⇓-7]和主动控制^[8⇓-10]两大类.其中,可变尾缘襟翼技术作为一种常见的主动控制手段,可通过改变翼型改变叶片的气动性能.

Chen等^[8]通过对比固定浆距和变桨距垂直轴风力机的性能,发现叶片在旋转过程中做适当幅度的俯仰运动,能够有效抑制功率输出、转速和转矩输出的波动,并提高风力机的功率系数.Li等^[9]为探索最优桨距控制空间解,搭建由遗传算法和计算流体力学(CFD)仿真模块组成的变桨距自动优化平台,得到能在较广的叶尖速比范围内使风力机平均功率系数均得以提升的优化叶片浆距.向斌等^[10]为减小动态失速对垂直轴风力机气动性能的影响,提出在翼型尾缘处布置动态格尼襟翼的主动控制方法.研究结果表明,采用主动式格尼襟翼提高了垂直轴风力机的运行稳定性,并增大风力机在低风速下的启动力矩.缪维跑等^[11]对NACA0012两段式襟翼翼型进行数值模拟,利用弦线变换得到襟翼摆角与攻角的关系,发现襟翼翼型在静态条件下由摆角引起的攻角迁移现象.研究结果表明,随襟翼摆角增大,有效攻角范围减小,翼型攻角产生迁移,受力亦发生变化.祖红亚等^[12]以NACA0018为基准翼型,采用数值模拟方法研究尾缘襟翼相对长度对翼型气动性能的影响.研究结果表明,襟翼对翼型周围主涡发展和变化具有影响,不仅改善了翼型的失速特性,同时也提高了翼型的气动性能,襟翼翼型的失速攻角在研究范围内均大于基准翼型.

上述研究均是针对某一特定翼型下襟翼对叶片或风力机气动性能的影响所展开.然而,对于不同的翼型,襟翼对叶片和风力机气动性能的影响存在一定差异.因此,针对不同翼型下带有尾缘襟翼的风力机进行气动性能分析和优化设计具有重要工程意义.

本文采用计算流体动力学数值模拟方法,建立添加分离式尾缘襟翼的H型垂直轴风力机二维模型,基于NACA0018、NACA0021和NACA0024共3种基础翼型与-16°、-8°、0°、8°、16°共5组襟翼偏转角参数,从风能利用率、叶片气动力系数、风力机荷载和流场分布等方面对比分析了尾缘襟翼对不同翼型的H型垂直轴风力机的气动性能影响,为进一步研究和提高H型垂直轴风力机气动性能提供参考.

1 计算模型和方法

1.1 几何模型

采用文献[13]研究的三叶片H型垂直轴风力机作为原始模型,具体模型如图1所示.该垂直轴风力机使用NACA0021对称翼型,弦长c=85.8 mm,风轮旋转直径D=1 030 mm,转轴直径为c/6,自由来流速度U_∞=9 m/s,叶片安装角为0°.3个叶片在同一个平面内,叶片长度l=1.456 4 m, 互为120°夹角.风力机模型的展弦比(叶片长度与弦长之比)接近17,相对较大.根据Rezaeiha等^[14]的结论,该风力机叶片的三维效应对叶片中部截面影响不显著,适用于二维CFD模拟.因此,为减小计算量,本研究使用二维CFD模拟.为提高计算效率,忽略支撑杆对风力机整体性能的影响.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 原始H型垂直轴风力机3D模型

Fig.1 3D model of original H-type VAWT

图2为风力机的旋转示意图.其中, α为叶片相对局部速度与弦线之间的夹角,即攻角;R为风轮旋转半径; $ω^{'}$ 为角速度; $F_{T}$ 为叶片所受切向力; $F_{N}$ 为法向力; $F_{L}$ 为升力; $F_{D}$ 为阻力.假定顺风区的自由来流速度与逆风区的自由来流速度相同^[15],即不考虑流动诱导的影响,相对局部速度可表示为

(1)

\begin{array}{l} W = \sqrt{(U_{\infty} + λ U_{\infty} {c o s θ)}^{2} + (λ U_{\infty} {s i n θ)}^{2}} = \\ U_{\infty} \sqrt{1 + 2 λ c o s θ + λ^{2}} \end{array}

式中:λ为叶尖速比;θ为叶片所在方位角.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 风力机旋转示意图

Fig.2 Schematic diagram of wind turbine rotation

攻角可表示为

(2)

\begin{matrix} α = a r c t a n (\frac{s i n θ}{λ + c o s θ}) \end{matrix}

根据叶片所在方位角θ的不同,将风轮旋转区域划分为逆风区(45°≤θ<135°)、背风区(135°≤θ<225°)、顺风区(225°≤θ<315°)和向风区(315°≤θ<45°)4个区域.功率系数C_p可由单叶片的弯矩系数C_m决定.

(3)

\begin{matrix} C_{p} = C_{m} N λ \end{matrix}

(4)

\begin{array}{l} C_{m} = \frac{M}{0.5 ρ U_{\infty}^{2} A_{C} R} \end{array}

式中: $N 为叶片个数; M 为风力机弯矩; ρ 为空气密度; A_{C}$ 为参考面积.

此外,选取切向力系数 $C_{T} 、法向力系数 C_{N} 、推力系数 C_{t h} 和侧向力系数 C_{l a}$ 共4个关键气动力参数,对作用在叶片上的气动力进行分析.

(5)

\begin{matrix} C_{T} = \frac{F_{T}}{0.5 ρ c H W^{2}} \end{matrix}

(6)

C_{N} = \frac{F_{N}}{0.5 ρ c H W^{2}}

(7)

C_{t h} = \frac{F_{t h}}{0.5 ρ U_{\infty}^{2} A}

(8)

C_{l a} = \frac{F_{l a}}{0.5 ρ U_{\infty}^{2} A}

式中: $H$ 为叶片沿展向长度;F_th为叶片所受推力;F_la为侧向力;A为风力机的扫风面积.在二维计算中, $F_{T} 、 F_{N} 、 F_{t h} 和 F_{l a}$ 可根据计算软件监测到的叶片x和y方向上的力 $F_{x} 和 F_{y}$ 与方位角 $θ$ 的关系计算得到.

在原始模型的基础上,建立添加分离式尾缘襟翼的垂直轴风力机几何模型,叶片结构如图3所示.叶片弦长与原始模型的弦长相同,为c=85.8 mm.将叶片前缘至襟翼旋转中心的距离定义为主翼长度c₁,襟翼旋转中心至叶片尾缘的距离定义为襟翼长度c₂,主翼与襟翼的间隙定义为翼缝宽度d.参照祖红亚等^[16-17]的工作,设置襟翼长度为0.3倍弦长,即c₂=0.3c,设置翼缝宽度d=1.7%c.为了研究不同襟翼偏转角对垂直轴风力机气动性能的影响,襟翼偏转角β可绕其旋转中心在 ±16° 范围内偏转,定义偏转角为正表示襟翼向上偏转,偏转角为负表示襟翼向下偏转.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 分离式尾缘襟翼叶片示意图

Fig.3 Schematic diagram of blade with separated trailing edge flap

1.2 数值方法

运用以有限体积法为基础的商业软件Fluent对垂直轴风力机气动性能和流场特性进行分析.考虑到风力机的性能在很大程度上取决于叶片边界层的发展,因此对转捩点的准确预测至关重要.采用4方程转捩剪切应力输运(SST)湍流模型求解纳维-斯托克斯(N-S)方程.转捩SST湍流模型除了使用k-ω SST湍流模型中采用的湍流动能和耗散率方程外,还结合另外两个考虑间歇因子 $γ 和当地转捩雷诺数 R e_{θ t}$ 的输运方程^[18-19].其优点在于,一方面避免了一般情况下对平均场进行积分的过程,计算周期短且计算要求较低;另一方面对自由来流湍流度、分离和压力梯度等影响转捩的因素敏感,可以更准确地捕捉层流到湍流之间的过渡^[20-21].

间歇因子表示在转捩区某一固定点流动处于有脉动和无脉动状态的时间比例,其输运方程为

(9)

\begin{array}{l} \frac{\partial (ρ γ)}{\partial t} + \frac{\partial (ρ U_{j} γ)}{\partial x_{j}} = \\ \frac{\partial}{\partial x_{j}} [(μ + \frac{μ_{t}}{σ_{γ}}) \frac{\partial γ}{\partial x_{j}}] + P_{γ} - E_{γ} \end{array}

(10)

P_{γ} = F_{l e} c_{a 1} ρ S {(γ F_{o n})}^{0.5} (1 - γ)

(11)

E_{γ} = c_{a 2} ρ Ω γ F_{t u} (c_{e 2} γ - 1)

(12)

t = \frac{500 μ}{ρ U_{\infty}^{2}}

式中: t为时间尺度;μ为动力黏度系数; $U_{j}$ 为y轴方向的速度分量; $x_{j}$ 为y轴方向的位移分量; $μ_{t}$ 为转捩常数;S为应变率; $F_{l e}$ 为转捩长度函数; $F_{o n} 和 F_{t u}$ 为转捩控制函数;Ω为涡量; $c_{a 1} 、 c_{a 2} 、 c_{e 2} 和 σ_{γ}$ 为转捩常数.

动量厚度雷诺数 $\tilde{R} e_{θ t}$ 的输运方程为

(13)

\begin{array}{l} \frac{\partial (ρ \tilde{R} e_{θ t})}{\partial t} + \frac{\partial (ρ U_{j} \tilde{R} e_{θ t})}{\partial x_{j}} = \\ \frac{\partial}{\partial x_{j}} [σ_{θ t} (μ + μ_{t}) \frac{\partial \tilde{R} e_{θ t}}{\partial x_{j}}] + P_{θ t} \end{array}

(14)

P_{θ t} = c_{θ t} ρ (R e_{θ t} - \tilde{R} e_{θ t}) (1 - F_{θ t}) / t

式中:混合函数 $F_{θ t}$ 为边界层中关闭源项; $c_{θ t} 和 σ_{θ t}$ 为转捩常数.

在垂直轴风力机非定常气动特性的计算过程中,计算时间步长设置为T/1 440,其中,T为一个旋转周期的时长.确保数值模拟过程中风力机转子每转动0.25°进行一个时间步计算.单次时间步长包含20次迭代,保证模拟的湍动能残差ξ始终控制在ξ=10^-6附近波动.基于Rezaeiha等^[14]的研究结果,湍动能残差控制为ξ=10^-6可以兼顾模拟的时间复杂度和结果精确度.整个模拟过程设置为20个周期,总时长T_max=2.72 s,前10个周期的初始化及流场发展可以保证模拟结果的充分收敛,第11~20个周期的模拟用于分析风力机气动性能的相关参数.

2 计算域及有效性验证

2.1 计算域和网格分布

计算流体动力学分析要求高精度流场,本文建立的二维模拟区域保证流场内流体能够充分发展.为减小进出口边界和侧边界对计算精度的影响^[22],将风力机布置在靠近上游的位置,入口到旋转中心的距离为10D,旋转域直径为1.5D,计算域沿入流方向的长度为30D,宽度为20D.模型的阻塞比(在二维计算中定义为风轮直径与计算域宽度之比)为5%,由对称边界条件导致的流动加速非常微弱,对风力机计算空气动力学性能的影响可以忽略^[23].在计算中,忽略连接杆件的影响,将整个流场分成4个区域:外流域F、转轴控制域Z1、环形旋转域Z2和叶片控制域Z3,各区域之间设置交界面,通过滑移网格技术实现风轮旋转,具体计算域划分如图4所示.考虑叶片表面的粗糙度和扰流效应,叶片布置为无滑移壁面.边界设置为前端速度入口 U_∞=9 m/s.为了统一模型参数的计算与分析,湍流强度设置为3%.空气密度ρ=1.225 kg/m³,动力黏度 μ=1.789 4×10^-5 Pa·s,压力出口为0 Pa,基于弦长的雷诺数Re=5.28×10⁴.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 计算流场区域划分

Fig.4 Division of computational flow field region

为保证旋转域与外流场的良好过渡和下游流场的充分发展与稳定,对外流域的网格以旋转轴为中心进行十字形加密处理,加密宽度为3D,促进实现流场的进一步精确模拟.分别采用二维结构网格与非结构网格进行拓扑,转轴控制域Z1和外流场域F采用结构网格,环形旋转域Z2和叶片控制域Z3采用非结构网格进行网格加密,贴近壁面处使用边界层网格,首层网格高度设置为 0.01 mm,网格增长率为1.2%,使叶片壁面处y⁺≈1,以满足准确刻画边界层流动的要求.详细网格布置如图5所示,单个叶片控制域网格数约为17万,总网格数约为68万.所有仿真计算都在配置为Intel(R) Xeon(R) CPU E5-2676 V3@2.4 Hz的小型服务器上进行,一个计算周期大约需要消耗1 h,约20 h完成一个算例.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 计算域的网格拓扑

Fig.5 Mesh topology of computational domain

2.2 网格独立性测试

网格质量通常对模拟结果的准确性有显著影响.因此,在正式计算前,进行网格独立性测试.改变叶片的最小表面尺寸,共设置3类网格方案,分别为粗网格、中网格和细网格.对叶尖速比λ=2.64时不同网格方案计算的功率系数进行比较,具体网格设置和结果如表1所示.

表1 基于功率系数误差的网格独立性测试

Tab.1 Mesh independence test based on errors of power coefficient

网格方案	总网格数/个	功率系数	相对误差/%
粗网格	450 679	0.312	11.86
中网格	644 226	0.350	1.13
细网格	940 388	0.354

注:“空白”表示无此项.

新窗口打开| 下载CSV

结果表明,细网格和中网格方案的功率系数差别非常小,分别为0.354和0.350,但两者的总网格数量相差较大,约19万网格.另一方面,粗网格方案虽然网格数更少,但计算得到的功率系数的相对误差较大.需要注意的是,网格数量越多并不意味着结果越准确,因为还会受网格形状和分布的影响.因此,综合考虑计算时长与计算精度的要求,本研究基于中网格方案来执行其他仿真.

2.3 模型有效性验证

为评估CFD计算结果的准确性,选用文献[13]的实验结果与当前数值模拟结果进行对比验证.该实验数据已被广泛用于研究垂直轴风力机的气动性能^[24-25],计算结果与实验数据对比如图6所示.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 现有结果与实验和数值数据对比

Fig.6 Comparison of current results with experimental and numerical data

可知,在叶尖速比λ达到额定叶尖速比2.64前,功率系数 $C_{p}$ 随λ增大而增大;而λ达到2.64后,功率系数 $C_{p}$ 随λ增加呈减小态势.在低叶尖速比下,数值模拟的结果较好地表现了功率系数 $C_{p}$ 随叶尖速比λ变化的趋势,而在高叶尖速比时,数值模拟结果对风力机的功率系数存在一定的高估,这与文献[25]中的数值计算结果类似.在高叶尖速比下,C_p被高估主要有两个原因:一是在二维CFD模拟中忽略了叶片高度,无法描述叶片两端的局部流动特征,不存在叶片尖端损失;二是由于CFD模拟过程中的几何简化与建模,例如叶片辐条和连接至风力机塔架的支撑杆会导致较大的阻力,并降低风洞测量时的风力机性能,而这些在本研究的数值模拟中被忽略.相关研究表明,由于支撑臂摩擦而导致的风力机性能损失可以达到20%左右,在高叶尖速比下,这一损失将更明显^[13,15].对于低叶尖速比(存在动态失速),由于几何简化,同样存在类似高估的问题.而在中等叶尖速比下,可能由于二维转捩SST湍流模型无法准确模拟风力机气动效应的固有三维流动复杂性,使这种预期的高估被抵消一部分.

总之,本文所采用的数值模型能够较好地反映垂直轴风力机在风洞中的气动特性,准确地再现实验结果.因此,目前的数值模型可以作为后续研究的一种可靠模拟方法.

3 结果与分析

NACA对称翼型族具有较高的升阻比和良好的失速特性,在H型垂直轴风力机的叶片选型上得到了广泛应用.本文选用3种不同厚度的基础翼型NACA0018、NACA0021和NACA0024,作为不同H型垂直轴风力机的叶片主体,叶片弦长与原始风力机模型的弦长一致.为分析尾缘襟翼对不同垂直轴风力机气动性能的影响规律,针对添加分离式尾缘襟翼的3种基础翼型NACA0018、NACA0021和NACA0024,分别选取襟翼偏转角为 -16°、-8°、0°、8°和16°的5种情况,取原风力机模型的额定叶尖速比2.64进行计算模拟,讨论风力机叶片的气动性能与叶片附近区域的流场特性.

3.1 气动性能与气动荷载

叶片弯矩系数 $C_{m}$ 是决定风能利用系数 $C_{p}$ 的关键参数.图7给出了翼型为NACA0018、NACA0021和NACA0024共3种垂直轴风力机在不同襟翼偏转角下单个叶片的弯矩系数随方位角的变化.可知,在逆风区,正向襟翼偏转角能够有效提高叶片的弯矩系数,且β=16°时,叶片的瞬时弯矩系数达到最优,此时NACA0018风力机的最大弯矩系数由0.179提高至0.246,增长了37.4%,NACA0021和NACA0024风力机的最大弯矩系数也分别增长了25.7%和30.4%.原因是在此方位角区间内,正向偏转的襟翼增大了叶片的有效攻角,使叶片升阻比得到提升.同样,在顺风区,负向襟翼偏转角对叶片的瞬时弯矩系数产生有利影响,且在β=-16° 时弯矩系数达到最优.考虑到中等叶尖速比下的垂直轴风力机本身气动性能较优,并不存在较严重的动态失速问题^[26],因此可知,在其他区域,襟翼偏转角的提升效果有限,甚至低于无偏转时的风力机模型,对叶片弯矩系数造成不利影响.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 不同襟翼偏转角下风力机的单叶片弯矩系数

Fig.7 Moment coefficient of blade at different flap angles

单个叶片的切向力系数和法向力系数对比分别如图8和图9所示.切向力方向与叶片绕旋转轴旋转的切线方向相同为正;法向力垂直于叶片绕转轴旋转的切线,指向风轮旋转中心为正.与单叶片弯矩系数随方位角的变化相似,逆风区的正向襟翼偏转角和顺风区的负向偏襟翼转角可以使叶片的切向力系数得到提高,且β=16°时,叶片切向力系数在θ=90°处达到最大值.此时NACA0018风力机的最大切向力系数由1.757提高至2.532,增长了44.1%, NACA0021和NACA0024风力机的最大切向力系数也分别增长了30.6%和37.9%.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 不同襟翼偏转角下风力机的单叶片切向力系数

Fig.8 Tangential force coefficient of blade with different flap angle

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 不同襟翼偏转角下风力机的单叶片法向力系数

Fig.9 Normal force coefficient of blade at different flap angles

从叶片的法向力系数对比可知,在整个旋转周期范围内,负向襟翼偏转角能够有效降低法向力系数.在β=-16° 时,叶片法向力系数在θ=90°处达到最大值.此时NACA0018风力机的最大法向力系数由10.310下降到7.202,降低了30.1%,NACA0021和NACA0024风力机的最大法向力系数也分别降低了24.6%和27.2%,表示风力机在旋转过程中支撑杆受到的荷载可以得到大幅降低.

不同垂直轴风力机在一个稳定的旋转周期内平均功率系数和平均弯矩系数的对比如图10(a)和图10(b)所示.显然,分离式尾缘襟翼的偏转角对风力机的动力性能有重要影响.对3种垂直轴风力机而言,在某些方位角区间内尾缘襟翼的偏转能够使风力机叶片的瞬时弯矩系数得到提高,但风力机整机的平均功率系数和平均弯矩系数均随襟翼偏转角的绝对值的增大而减小,说明简单地设置固定偏转襟翼并不能使整机动力性能得到明显改善,这与图7分析结果一致.

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 不同垂直轴风力机的对比

Fig.10 Comparison of different VAWTs

对于大型垂直轴风机的设计与运用,叶片周期性的转动将导致整机受到周期性的疲劳荷载,这些疲劳荷载是需要重点考量的因素.风力机的气动载荷是衡量其气动性能的主要参数,图10(c)和10(d)分别为不同垂直轴风力机在一个稳定的旋转周期内平均推力系数和平均侧向力系数的对比,推力和侧向力分别指风力机整体受到的沿来流方向的力和垂直于来流方向的力.可知,负向襟翼偏转角会在一定程度上降低风力机的推力系数;而偏转角为正时,推力系数随偏转角度的增大而增大.对于侧向力系数,整体呈现为随偏转角绝对值的增大而增大的趋势.说明单一地改变襟翼偏转角对改善风机整体荷载的作用效果不突出,与图8和图9的分析结果相符.

为进一步分析尾缘襟翼对不同厚度翼型的垂直轴风力机动力性能的影响,将各风力机在不同襟翼偏转角下的功率系数与 $β = 0 °$ 时的功率系数作差得到 $Δ C_{p}$ ,并求功率系数的变化率,即 $\frac{Δ C_{p}}{C_{p} (β = 0 °)}$ ,可直观得到不同风力机 $Δ C_{p}$ 随尾缘襟翼的变化分布,如图11所示.对任意风力机而言,襟翼偏转角的绝对值越大,功率系数变化越明显,且变化速率逐渐增大.对比3种垂直轴风力机,当襟翼偏转角为负时,NACA0024风力机的功率系数随偏转角变化最明显,结合图10(a)可知,偏转角从0°减小至 -16°,其功率系数由0.301降低至0.129,下降了57.1%,NACA0021和NACA0018风力机的功率系数也分别下降了50.3%和39.2%;当襟翼偏转角为正时,NACA0018风力机的功率系数随偏转角变化最明显,偏转角从0° 增大至16°,其功率系数由0.352降低至0.054,下降了84.7%,NACA0021和NACA0018风力机的功率系数也分别下降了65.3%和55.3%.由此可见,在负向襟翼偏转角下,风力机功率系数受偏转角影响的程度与翼型厚度呈正相关;而在正向襟翼偏转角下,功率系数受偏转角影响的程度与翼型厚度呈负相关.

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 不同垂直轴风力机ΔC_p随尾缘襟翼变化分布

Fig.11 Distribution of ΔC_p with trailing edge flap deflection angle for different airfoils

3.2 流场分布

动态失速与叶片表面发生的流动分离直接相关^[27].图12给出了NACA0018和NACA0024风力机在不同襟翼偏转角下翼型周围的涡量分布对比,以观察流动分离的形成和发展,并进一步解释风力机气动性能特性的内在机理.

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12 不同风力机翼型周围的涡量分布对比

Fig.12 Comparison of vorticity distributions around airfoil for different VAWTs

在方位角θ=60° 时,对于β=0° 的模型,NACA0018风力机的气流基本附着在叶片表面,未发生流动分离现象,而NACA0024风力机尾缘下游区域存在显著的周期性旋涡结构;对于β=16° 的模型,两种不同翼型工况下,襟翼的偏转导致襟翼处产生明显的流动分离,尾缘处形成周期性的卡门涡街.当叶片旋转到θ=90° 位置时,对于β=0° 的模型,由于攻角的增大,尾缘产生轻微流动分离,形成周期性尾涡.同时由于翼缝的存在,襟翼靠近翼缝的位置出现小尺度的涡;对于β=16° 的模型,襟翼的偏转延缓了尾缘处的流动分离,并使翼缝处的涡得到小范围扩散,解释了图7所示的叶片产生更大弯矩系数的结果.当叶片旋转到θ=120° 位置时,由于攻角增大,叶片表面的气流高度分离.后缘涡和前缘涡的强度大大增加,表明在这一阶段发生了严重的动态失速,对应图7中弯矩系数的急剧下降.在θ=150° 位置,由于攻角逐渐减小,从前缘产生的气流分离逐渐受到抑制,后缘产生的旋涡结构向下游处进一步发展.对于β=16° 的工况,叶片襟翼的偏转加剧了后缘附近处的流动分离.当叶片旋转到θ=180° 位置时,从叶片前缘和后缘脱落的涡逐渐消散,攻角在此位置再次归零,气流逐渐重新附着到叶片表面,流动分离强度减弱,因而导致图7中弯矩系数激增.对比不同方位角下,NACA0018和NACA0024风力机有无襟翼偏转工况下的涡量图可知,β=16° 的偏转襟翼对NACA0018风力机流场分布的影响比对NACA0024风力机的影响更显著,尤其在θ=120° 和θ=150° 的位置,这与图11得到的分析结果一致.

4 结论

针对添加尾缘襟翼的垂直轴风力机的空气动力学问题,采用计算流体动力学转捩SST湍流模型,开展3种不同分离式尾缘襟翼的翼型(NACA0018、NACA0021和NACA0024)叶片的H型垂直轴风力机气动性能的数值研究.对比不同襟翼偏转角 -16°、-8°、0°、8°、16° 下3种风力机模型的风能利用率、叶片气动力系数、风力机荷载和流场分布等结果,主要结论如下:

(1) 襟翼偏转角度变化能够导致风力机气动性能变化.在本文的研究模型中,与尾缘襟翼不发生偏转的模型相比,风力机的平均功率系数随襟翼偏转角的绝对值增大而减小.

(2) 在逆风区(45°≤θ<135°),正向襟翼偏转角可以有效提高叶片的弯矩系数,当β=16° 时,叶片的瞬时弯矩系数达到最优,NACA0018、NACA0021和NACA0024风力机的弯矩系数分别增长了37.4%、25.7%和30.4%;在顺风区(225°≤θ< 315°),负向襟翼偏转角对叶片的瞬时弯矩系数产生有利影响,在β=-16° 时弯矩系数达到最优.

(3) 对不同翼型的H型垂直轴风力机,当襟翼偏转角为负时,风能利用率受偏转角影响的程度与翼型厚度呈正相关,即NACA0024风力机受影响程度最大,在β=-16° 时风能利用率下降了57.1%;当襟翼偏转角为正时,风能利用率受偏转角影响的程度与翼型厚度呈负相关,即NACA0018风力机受影响程度最大,在β=16° 时风能利用率下降了84.7%.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[21]

LAM

H F

, PENG

H Y

Study of wake characteristics of a vertical axis wind turbine by two-and three-dimensional computational fluid dynamics simulations

[J]. Renewable Energy, 2016, 90: 386-398.

DOI:10.1016/j.renene.2016.01.011 URL [本文引用: 1]

[22]

ZHANG

L X

, LIANG

Y B

, LIU

X H

, et al.

Aerodynamic performance prediction of straight-bladed vertical axis wind turbine based on CFD

[J]. Advances in Mechanical Engineering, 2013, 5: 905379.

DOI:10.1155/2013/905379 URL [本文引用: 1]

[23]

REZAEIHA

, KALKMAN

, BLOCKEN

CFD simulation of a vertical axis wind turbine operating at a moderate tip speed ratio: Guidelines for minimum domain size and azimuthal increment

[J]. Renewable Energy, 2017, 107: 373-385.

DOI:10.1016/j.renene.2017.02.006 URL [本文引用: 1]

[24]

SOBHANI

, GHAFFARI

, MAGHREBI

M J

Numerical investigation of dimple effects on darrieus vertical axis wind turbine

[J]. Energy, 2017, 133: 231-241.

DOI:10.1016/j.energy.2017.05.105 URL [本文引用: 1]

[25]

REZAEIHA

, KALKMAN

, BLOCKEN

Effect of pitch angle on power performance and aerodynamics of a vertical axis wind turbine

[J]. Applied Energy, 2017, 197: 132-150.

DOI:10.1016/j.apenergy.2017.03.128 URL [本文引用: 1]

[26]

SAGHARICHI

, MAGHREBI

M J

, ARABGOLARCHEH

Variable pitch blades: An approach for improving performance of Darrieus wind turbine

[J]. Journal of Renewable and Sustainable Energy, 2016, 8(5): 053305.

DOI:10.1063/1.4964310 URL [本文引用: 1]

[1]

MITTAL

, MITRA

Determining layout of a wind farm with optimal number of turbines: A decomposition based approach

[J]. Journal of Cleaner Production, 2018, 202: 342-359.

DOI:10.1016/j.jclepro.2018.08.093 URL [本文引用: 1]

[2]

HAND

, KELLY

, CASHMAN

Aerodynamic design and performance parameters of a lift-type vertical axis wind turbine: A comprehensive review

[J]. Renewable and Sustainable Energy Reviews, 2021, 139: 110699.

DOI:10.1016/j.rser.2020.110699 URL [本文引用: 1]

[3]

HAND

, CASHMAN

A review on the historical development of the lift-type vertical axis wind turbine: From onshore to offshore floating application

[J]. Sustainable Energy Technologies and Assessments, 2020, 38: 100646.

DOI:10.1016/j.seta.2020.100646 URL [本文引用: 1]

[4]

Q A

, MAEDA

, KAMADA

, et al.

Study on stall behavior of a straight-bladed vertical axis wind turbine with numerical and experimental investigations

[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 164: 1-12.

DOI:10.1016/j.jweia.2017.02.005 URL [本文引用: 1]

[5]

ZHU

H T

, HAO

W X

, LI

, et al.

Effect of geometric parameters of Gurney flap on performance enhancement of straight-bladed vertical axis wind turbine

[J]. Renewable Energy, 2021, 165: 464-480.

DOI:10.1016/j.renene.2020.11.027 URL [本文引用: 1]

[6]

WANG

H P

, ZHANG

, QIU

Q G

, et al.

Flow control on the NREL S809 wind turbine airfoil using vortex generators

[J]. Energy, 2017, 118: 1210-1221.

DOI:10.1016/j.energy.2016.11.003 URL [本文引用: 1]

[7]

DAM

C P V

The aerodynamic design of multi-element high-lift systems for transport airplanes

[J]. Progress in Aerospace Sciences, 2002, 38(2): 101-144.

DOI:10.1016/S0376-0421(02)00002-7 URL [本文引用: 1]

[8]

CHEN

, SU

S S

, VIOLA

I M

, et al.

Numerical investigation of vertical-axis tidal turbines with sinusoidal pitching blades

[J]. Ocean Engineering, 2018, 155: 75-87.

DOI:10.1016/j.oceaneng.2018.02.038 URL [本文引用: 2]

[9]

, XIAO

Y Q

, XU

Y L

, et al.

Optimization of blade pitch in H-rotor vertical axis wind turbines through computational fluid dynamics simulations

[J]. Applied Energy, 2018, 212: 1107-1125.

DOI:10.1016/j.apenergy.2017.12.035 URL [本文引用: 2]

[10]

向斌, 缪维跑, 李春, 等.

垂直轴风力机叶片尾缘主动式格尼襟翼气动效率研究分析

[J]. 热能动力工程, 2020, 35(4): 242-250.

[本文引用: 2]

XIANG

Bin

, MIAO

Weipao

, LI

Chun

, et al.

Research of aerodynamic efficiency of active gurney flaps on the trailing edge of vertical axis wind turbine blades

[J]. Journal of Engineering for Thermal Energy and Power, 2020, 35(4): 242-250.

[本文引用: 2]

[11]

缪维跑, 李春, 聂佳斌, 等.

襟翼翼型位置对气动性能的影响研究

[J]. 能源研究与信息, 2015, 31(4): 242-246.

[本文引用: 1]

MIAO

Weipao

, LI

Chun

, NIE

Jiabin

, et al.

Influence of the flap airfoils with different positions on the aerodynamic performance

[J]. Energy Research and Information, 2015, 31(4): 242-246.

[本文引用: 1]

[12]

祖红亚, 李春, 李润杰, 等.

襟翼相对长度对翼型气动性能的影响

[J]. 动力工程学报, 2015, 35(8): 666-673.

[本文引用: 1]

Hongya

, LI

Chun

, LI

Runjie

, et al.

Effect of relative flap length on aerodynamic performance of the airfoil

[J]. Journal of Chinese Society of Power Engineering, 2015, 35(8): 666-673.

[本文引用: 1]

[13]

RACITI

CASTELLI M

, ARDIZZON

, BATTISTI

, et al.

Modeling strategy and numerical validation for a darrieus vertical axis micro-wind turbine

[C]//Proceedings of the ASME 2010 International Mechanical Engineering Congress and Exposition. Vancouver, British Columbia, Canada: ASME, 2010: 409-418.

[本文引用: 3]

[14]

REZAEIHA

, MONTAZERI

, BLOCKEN

Characterization of aerodynamic performance of vertical axis wind turbines: Impact of operational parameters

[J]. Energy Conversion and Management, 2018, 169: 45-77.

DOI:10.1016/j.enconman.2018.05.042 URL [本文引用: 2]

[15]

PARASCHIVOIU

Wind turbine design with emphasis on Darrieus concept

[M]//Ion paraschivoiu. Canada: Presses inter Polytechnique, 2002.

[本文引用: 2]

[16]

祖红亚, 李春, 陆云凤, 等.

襟翼翼缝相对宽度对翼型气动性能影响研究

[J]. 能源工程, 2015(3): 12-19.

[本文引用: 1]

Hongya

, LI

Chun

, LU

Yunfeng

, et al.

Study on effect of relative width of flap slot on airfoil aerodynamic performance

[J]. Energy Engineering, 2015(3): 12-19.

[本文引用: 1]

[17]

李润杰, 祖红亚, 李春, 等.

襟翼翼缝相对宽度对翼型动态气动性能的影响

[J]. 热能动力工程, 2016, 31(4): 38-44.

[本文引用: 1]

Runjie

, ZU

Hongya

, LI

Chun

, et al.

Effect of the relative width of wing flap slit on the aerodynamic performance of airfoil

[J]. Journal of Engineering for Thermal Energy and Power, 2016, 31(4): 38-44.

[本文引用: 1]

[18]

GHASEMIAN

, ASHRAFI

Z N

, SEDAGHAT

A review on computational fluid dynamic simulation techniques for Darrieus vertical axis wind turbines

[J]. Energy Conversion and Management, 2017, 149: 87-100.

DOI:10.1016/j.enconman.2017.07.016 URL [本文引用: 1]

[19]

MENTER

F R

, LANGTRY

R B

, LIKKI

S R

, et al.

A correlation-based transition model using local variables: Part I. Model formulation

[J]. Journal of Turbomachinery, 2006, 128(3): 413.

DOI:10.1115/1.2184352 URL [本文引用: 1]

[27]

POST

M L

, CORKE

T C

Separation control using plasma actuators: Dynamic stall vortex control on oscillating airfoil

[J]. AIAA Journal, 2006, 44(12): 3125-3135.

DOI:10.2514/1.22716 URL [本文引用: 1]

[20]

ARAB

, JAVADI

, ANBARSOOZ

, et al.

A numerical study on the aerodynamic performance and the self-starting characteristics of a Darrieus wind turbine considering its moment of inertia

[J]. Renewable Energy, 2017, 107: 298-311.

DOI:10.1016/j.renene.2017.02.013 URL [本文引用: 1]

Study of wake characteristics of a vertical axis wind turbine by two-and three-dimensional computational fluid dynamics simulations

2016

... 运用以有限体积法为基础的商业软件Fluent对垂直轴风力机气动性能和流场特性进行分析.考虑到风力机的性能在很大程度上取决于叶片边界层的发展,因此对转捩点的准确预测至关重要.采用4方程转捩剪切应力输运(SST)湍流模型求解纳维-斯托克斯(N-S)方程.转捩SST湍流模型除了使用k-ω SST湍流模型中采用的湍流动能和耗散率方程外,还结合另外两个考虑间歇因子

γ 和 当 地 转 捩 雷 诺 数 R e_{θ t}

的输运方程^[18-19].其优点在于,一方面避免了一般情况下对平均场进行积分的过程,计算周期短且计算要求较低;另一方面对自由来流湍流度、分离和压力梯度等影响转捩的因素敏感,可以更准确地捕捉层流到湍流之间的过渡^[20-21]. ...

Aerodynamic performance prediction of straight-bladed vertical axis wind turbine based on CFD

2013

... 计算流体动力学分析要求高精度流场,本文建立的二维模拟区域保证流场内流体能够充分发展.为减小进出口边界和侧边界对计算精度的影响^[22],将风力机布置在靠近上游的位置,入口到旋转中心的距离为10D,旋转域直径为1.5D,计算域沿入流方向的长度为30D,宽度为20D.模型的阻塞比(在二维计算中定义为风轮直径与计算域宽度之比)为5%,由对称边界条件导致的流动加速非常微弱,对风力机计算空气动力学性能的影响可以忽略^[23].在计算中,忽略连接杆件的影响,将整个流场分成4个区域:外流域F、转轴控制域Z1、环形旋转域Z2和叶片控制域Z3,各区域之间设置交界面,通过滑移网格技术实现风轮旋转,具体计算域划分如图4所示.考虑叶片表面的粗糙度和扰流效应,叶片布置为无滑移壁面.边界设置为前端速度入口 U_∞=9 m/s.为了统一模型参数的计算与分析,湍流强度设置为3%.空气密度ρ=1.225 kg/m³,动力黏度 μ=1.789 4×10^-5 Pa·s,压力出口为0 Pa,基于弦长的雷诺数Re=5.28×10⁴. ...

CFD simulation of a vertical axis wind turbine operating at a moderate tip speed ratio: Guidelines for minimum domain size and azimuthal increment

2017

Numerical investigation of dimple effects on darrieus vertical axis wind turbine

2017

... 为评估CFD计算结果的准确性,选用文献[13]的实验结果与当前数值模拟结果进行对比验证.该实验数据已被广泛用于研究垂直轴风力机的气动性能^[24-25],计算结果与实验数据对比如图6所示. ...

Effect of pitch angle on power performance and aerodynamics of a vertical axis wind turbine

2017

Variable pitch blades: An approach for improving performance of Darrieus wind turbine

2016

... 叶片弯矩系数

C_{m}

是决定风能利用系数

C_{p}

的关键参数.图7给出了翼型为NACA0018、NACA0021和NACA0024共3种垂直轴风力机在不同襟翼偏转角下单个叶片的弯矩系数随方位角的变化.可知,在逆风区,正向襟翼偏转角能够有效提高叶片的弯矩系数,且β=16°时,叶片的瞬时弯矩系数达到最优,此时NACA0018风力机的最大弯矩系数由0.179提高至0.246,增长了37.4%,NACA0021和NACA0024风力机的最大弯矩系数也分别增长了25.7%和30.4%.原因是在此方位角区间内,正向偏转的襟翼增大了叶片的有效攻角,使叶片升阻比得到提升.同样,在顺风区,负向襟翼偏转角对叶片的瞬时弯矩系数产生有利影响,且在β=-16° 时弯矩系数达到最优.考虑到中等叶尖速比下的垂直轴风力机本身气动性能较优,并不存在较严重的动态失速问题^[26],因此可知,在其他区域,襟翼偏转角的提升效果有限,甚至低于无偏转时的风力机模型,对叶片弯矩系数造成不利影响. ...

Determining layout of a wind farm with optimal number of turbines: A decomposition based approach

2018

... 近年来,风能以其清洁的特性、丰富的储量以及日趋成熟的风力发电技术,成为未来最具前景的新能源形式之一^[1].风力机作为捕获风能的重要装置,其构型对于提高风力机的风能转化率、降低生产和运行维护成本具有重要意义^[2].与传统水平轴风力机相比,垂直轴风力机具有设计简单、对风向不敏感等优点^[3],但垂直轴风力机较低的风能利用率成为了其大规模商业化应用的瓶颈. ...

Aerodynamic design and performance parameters of a lift-type vertical axis wind turbine: A comprehensive review

2021

A review on the historical development of the lift-type vertical axis wind turbine: From onshore to offshore floating application

2020

Study on stall behavior of a straight-bladed vertical axis wind turbine with numerical and experimental investigations

2017

... 在旋转过程中,垂直轴风力机叶片的攻角会出现周期性变化,当叶片的攻角超过失速攻角时,叶片表面气体就会流动分离,从而使得叶片的气动性能急剧下降^[4].因此,研究人员提出利用不同的流动控制手段改善叶片的气动性能.流动控制手段主要可分为被动控制^[5⇓-7]和主动控制^[8⇓-10]两大类.其中,可变尾缘襟翼技术作为一种常见的主动控制手段,可通过改变翼型改变叶片的气动性能. ...

Effect of geometric parameters of Gurney flap on performance enhancement of straight-bladed vertical axis wind turbine

2021

Flow control on the NREL S809 wind turbine airfoil using vortex generators

2017

The aerodynamic design of multi-element high-lift systems for transport airplanes

2002

Numerical investigation of vertical-axis tidal turbines with sinusoidal pitching blades

2018

... Chen等^[8]通过对比固定浆距和变桨距垂直轴风力机的性能,发现叶片在旋转过程中做适当幅度的俯仰运动,能够有效抑制功率输出、转速和转矩输出的波动,并提高风力机的功率系数.Li等^[9]为探索最优桨距控制空间解,搭建由遗传算法和计算流体力学(CFD)仿真模块组成的变桨距自动优化平台,得到能在较广的叶尖速比范围内使风力机平均功率系数均得以提升的优化叶片浆距.向斌等^[10]为减小动态失速对垂直轴风力机气动性能的影响,提出在翼型尾缘处布置动态格尼襟翼的主动控制方法.研究结果表明,采用主动式格尼襟翼提高了垂直轴风力机的运行稳定性,并增大风力机在低风速下的启动力矩.缪维跑等^[11]对NACA0012两段式襟翼翼型进行数值模拟,利用弦线变换得到襟翼摆角与攻角的关系,发现襟翼翼型在静态条件下由摆角引起的攻角迁移现象.研究结果表明,随襟翼摆角增大,有效攻角范围减小,翼型攻角产生迁移,受力亦发生变化.祖红亚等^[12]以NACA0018为基准翼型,采用数值模拟方法研究尾缘襟翼相对长度对翼型气动性能的影响.研究结果表明,襟翼对翼型周围主涡发展和变化具有影响,不仅改善了翼型的失速特性,同时也提高了翼型的气动性能,襟翼翼型的失速攻角在研究范围内均大于基准翼型. ...

Optimization of blade pitch in H-rotor vertical axis wind turbines through computational fluid dynamics simulations

2018

垂直轴风力机叶片尾缘主动式格尼襟翼气动效率研究分析

2020

垂直轴风力机叶片尾缘主动式格尼襟翼气动效率研究分析

2020

襟翼翼型位置对气动性能的影响研究

2015

襟翼翼型位置对气动性能的影响研究

2015

襟翼相对长度对翼型气动性能的影响

2015

襟翼相对长度对翼型气动性能的影响

2015

Modeling strategy and numerical validation for a darrieus vertical axis micro-wind turbine

2010

... 采用文献[13]研究的三叶片H型垂直轴风力机作为原始模型,具体模型如图1所示.该垂直轴风力机使用NACA0021对称翼型,弦长c=85.8 mm,风轮旋转直径D=1 030 mm,转轴直径为c/6,自由来流速度U_∞=9 m/s,叶片安装角为0°.3个叶片在同一个平面内,叶片长度l=1.456 4 m, 互为120°夹角.风力机模型的展弦比(叶片长度与弦长之比)接近17,相对较大.根据Rezaeiha等^[14]的结论,该风力机叶片的三维效应对叶片中部截面影响不显著,适用于二维CFD模拟.因此,为减小计算量,本研究使用二维CFD模拟.为提高计算效率,忽略支撑杆对风力机整体性能的影响. ...

... 可知,在叶尖速比λ达到额定叶尖速比2.64前,功率系数

C_{p}

随λ增大而增大;而λ达到2.64后,功率系数

C_{p}

随λ增加呈减小态势.在低叶尖速比下,数值模拟的结果较好地表现了功率系数

C_{p}

随叶尖速比λ变化的趋势,而在高叶尖速比时,数值模拟结果对风力机的功率系数存在一定的高估,这与文献[25]中的数值计算结果类似.在高叶尖速比下,C_p被高估主要有两个原因:一是在二维CFD模拟中忽略了叶片高度,无法描述叶片两端的局部流动特征,不存在叶片尖端损失;二是由于CFD模拟过程中的几何简化与建模,例如叶片辐条和连接至风力机塔架的支撑杆会导致较大的阻力,并降低风洞测量时的风力机性能,而这些在本研究的数值模拟中被忽略.相关研究表明,由于支撑臂摩擦而导致的风力机性能损失可以达到20%左右,在高叶尖速比下,这一损失将更明显^[13,15].对于低叶尖速比(存在动态失速),由于几何简化,同样存在类似高估的问题.而在中等叶尖速比下,可能由于二维转捩SST湍流模型无法准确模拟风力机气动效应的固有三维流动复杂性,使这种预期的高估被抵消一部分. ...

Characterization of aerodynamic performance of vertical axis wind turbines: Impact of operational parameters

2018

... 在垂直轴风力机非定常气动特性的计算过程中,计算时间步长设置为T/1 440,其中,T为一个旋转周期的时长.确保数值模拟过程中风力机转子每转动0.25°进行一个时间步计算.单次时间步长包含20次迭代,保证模拟的湍动能残差ξ始终控制在ξ=10^-6附近波动.基于Rezaeiha等^[14]的研究结果,湍动能残差控制为ξ=10^-6可以兼顾模拟的时间复杂度和结果精确度.整个模拟过程设置为20个周期,总时长T_max=2.72 s,前10个周期的初始化及流场发展可以保证模拟结果的充分收敛,第11~20个周期的模拟用于分析风力机气动性能的相关参数. ...

Wind turbine design with emphasis on Darrieus concept

2002

... 图2为风力机的旋转示意图.其中, α为叶片相对局部速度与弦线之间的夹角,即攻角;R为风轮旋转半径;

ω^{'}

为角速度;

F_{T}

为叶片所受切向力;

F_{N}

为法向力;

F_{L}

为升力;

F_{D}

为阻力.假定顺风区的自由来流速度与逆风区的自由来流速度相同^[15],即不考虑流动诱导的影响,相对局部速度可表示为 ...

... 可知,在叶尖速比λ达到额定叶尖速比2.64前,功率系数

C_{p}

随λ增大而增大;而λ达到2.64后,功率系数

C_{p}

随λ增加呈减小态势.在低叶尖速比下,数值模拟的结果较好地表现了功率系数

C_{p}

襟翼翼缝相对宽度对翼型气动性能影响研究

2015

... 在原始模型的基础上,建立添加分离式尾缘襟翼的垂直轴风力机几何模型,叶片结构如图3所示.叶片弦长与原始模型的弦长相同,为c=85.8 mm.将叶片前缘至襟翼旋转中心的距离定义为主翼长度c₁,襟翼旋转中心至叶片尾缘的距离定义为襟翼长度c₂,主翼与襟翼的间隙定义为翼缝宽度d.参照祖红亚等^[16-17]的工作,设置襟翼长度为0.3倍弦长,即c₂=0.3c,设置翼缝宽度d=1.7%c.为了研究不同襟翼偏转角对垂直轴风力机气动性能的影响,襟翼偏转角β可绕其旋转中心在 ±16° 范围内偏转,定义偏转角为正表示襟翼向上偏转,偏转角为负表示襟翼向下偏转. ...

襟翼翼缝相对宽度对翼型气动性能影响研究

2015

襟翼翼缝相对宽度对翼型动态气动性能的影响

2016

襟翼翼缝相对宽度对翼型动态气动性能的影响

2016

A review on computational fluid dynamic simulation techniques for Darrieus vertical axis wind turbines

2017

γ 和 当 地 转 捩 雷 诺 数 R e_{θ t}

A correlation-based transition model using local variables: Part I. Model formulation

2006

γ 和 当 地 转 捩 雷 诺 数 R e_{θ t}

Separation control using plasma actuators: Dynamic stall vortex control on oscillating airfoil

2006

... 动态失速与叶片表面发生的流动分离直接相关^[27].图12给出了NACA0018和NACA0024风力机在不同襟翼偏转角下翼型周围的涡量分布对比,以观察流动分离的形成和发展,并进一步解释风力机气动性能特性的内在机理. ...

A numerical study on the aerodynamic performance and the self-starting characteristics of a Darrieus wind turbine considering its moment of inertia

2017

γ 和 当 地 转 捩 雷 诺 数 R e_{θ t}

〈

〉