许波峰^,*,?,2), 朱紫璇^*,?, 戴成军^*,?, 蔡新^?, 王同光^**, 赵振宙^*,?*河海大学可再生能源发电技术教育部工程研究中心, 南京 210000
^?河海大学江苏省风电机组结构工程研究中心, 南京 210000
^**南京航空航天大学江苏省风力机设计高技术研究重点实验室, 南京 210000

INFLUENCE OF WIND SHEAR ON AERODYNAMIC CHARACTERISTICS AND WAKE SHAPE OF WIND TURBINE BLADES¹⁾

Xu Bofeng^,*,?,2), Zhu Zixuan^*,?, Dai Chengjun^*,?, Cai Xin^?, Wang Tongguang^**, Zhao Zhenzhou^*,?*Research Center for Renewable Energy Generation Engineering of Ministry of Education, Hohai University, Nanjing 210000, China
^?Jiangsu Province Wind Power Structural Engineering Research Center, Hohai University, Nanjing 210000, China
^**Jiangsu Key Laboratory of Hi-Tech Research for Wind Turbine Design, Nanjing University of Aeronautics and Astronautics, Nanjing 210000, China

通讯作者: 2) 许波峰, 副教授, 主要研究方向: 风力机空气动力学. E-mail:bfxu1985@hhu.edu.cn

收稿日期:2020-08-16接受日期:2020-10-8网络出版日期:2021-02-07

基金资助:

1) 江苏风力发电工程技术中心开放基金(ZK19-03-01)和江苏省研究生科研与实践创新计划资助项目.SJCX20_0179

Received:2020-08-16Accepted:2020-10-8Online:2021-02-07
作者简介 About authors


摘要
风力机通常运行在非定常工况中,其气动性能及尾迹会随着工况的变化而变化. 风剪切是风力机长期所处的环境,它会影响到叶片气动载荷、尾迹形状、总体性能等,分析风剪切作用下的叶片气动性能对风力机的设计有重要意义.本文采用一种时间步进自由涡尾迹(free vortex wake, FVW)方法,耦合FVW方法与风剪切模型,计算不同风剪切因子作用下叶片的气动力系数、推力以及风轮后的尾迹形状变化,研究尾迹形状变化对风轮旋转平面诱导速度及风力机叶片气动性能的影响. 结果表明:在风剪切入流条件下, 随着风剪切因子的增大,风力机的气动力系数随时间做周期性波动的幅度加剧, 推力的平均值逐渐减小,尾迹倾斜程度增大, 尾迹在轮毂下方的倾斜程度更明显;尾迹形状的变化使风轮平面轴向诱导速度因子分布不均匀,同时使风力机的总体性能降低且偏离较大;倾斜尾迹相比于对称尾迹对风轮平面处的诱导影响有明显差别, 波动幅值增大,气动力系数在波谷处的偏差比波峰处大. 尾迹越倾斜,风轮旋转平面处的载荷不对称性越明显.
关键词： 风力机;风剪切;自由涡尾迹;气动性能;尾迹

Abstract
A wind turbine usually operates in unsteady conditions, and the aerodynamic performance and wake of the rotor will change with the change of working conditions. Wind shear is the most common and long-term environment of wind turbine. It often affects the aerodynamic load, the shape of the wake and the overall performance of the wind turbine. So it is important to analyze the aerodynamic performance of blade of the wind turbine under the wind shear condition. In this paper, a time-marching free vortex wake method, which coupled with the wind shear model, is used to calculate the aerodynamic coefficient, thrust and wake shape change under different wind shear factors. At the same time, the influence of the wake shape variation on the induced velocity of the rotor rotating plane and the aerodynamic performance of the wind turbine blade are studied. The results show that: under the condition of wind shear inflow, with the increase of the wind shear factor, the fluctuation amplitude of aerodynamic coefficient of the wind turbine, which fluctuates periodically with time is increased, the average thrust decreases gradually, the tilted degree of wake increases and the tilted degree of wake under the center of hub is more obvious. The distortion of wake shape makes the distribution of the axial induced velocity factor in the rotational plane uneven, and makes the overall performance of wind turbine decrease and deviate greatly. There is a distinct difference of the induced effect on the rotational plane between the tilted wake and the symmetrical wake. The fluctuation amplitude of the aerodynamic coefficient induced by the tilted wake is larger than that induced by the symmetrical wake, and the deviation of the wave trough is more obvious than the wave crest. The more tilted the wake is, the more obvious the load asymmetry in the rotating plane is.
Keywords：wind turbine;wind shear;free vortex wake;aerodynamic performance;wake


PDF (12502KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
许波峰, 朱紫璇, 戴成军, 蔡新, 王同光, 赵振宙. 风剪切对风力机叶片气动性能及尾迹形状的影响¹⁾. 力学学报[J], 2021, 53(2): 362-372 DOI:10.6052/0459-1879-20-289
Xu Bofeng, Zhu Zixuan, Dai Chengjun, Cai Xin, Wang Tongguang, Zhao Zhenzhou. INFLUENCE OF WIND SHEAR ON AERODYNAMIC CHARACTERISTICS AND WAKE SHAPE OF WIND TURBINE BLADES¹⁾. Chinese Journal of Theoretical and Applied Mechanics[J], 2021, 53(2): 362-372 DOI:10.6052/0459-1879-20-289

引言

风力机运行在大气边界层中, 其叶片性能会受到风剪切的影响. 在风剪切入流条件下, 风轮旋转平面速度分布不均匀, 叶片载荷呈周期性变化, 风力机的性能也会有所降低^[1]. 此外, 从叶片尖部拖出的尾流的畸变会反过来作用于风轮平面, 进一步影响叶片气动性能. 随着机组的大型化发展, 风剪切对风力机叶片的影响会更加明显.

风力机叶片的气动性能计算方法主要有动量-叶素理论(blade element momentum, BEM)方法、CFD方法、自由涡尾迹方法. BEM方法在科学研究与工业领域得到了广泛应用, 对于非定常入流条件下的仿真需要加入相关修正模型来提高计算精准度^[2-7]. 唐新姿等^[8]基于修正的BEM和非嵌入式概率配置点法, 提出了一种风力机不确定空气动力学法, 研究了不确定性风速对风力机的影响. Ali等^[9]采用一种经验非定常BEM计算高空风力机气动特性, 在峰值工况运行下的性能与CFD方法结果相比, 其偏差在4.1%以内. 但在风轮平面上的入流分布明显不均匀的情况下, 气动载荷的计算会有较大的偏差^[10]. CFD方法的计算精度优于BEM, 不少****采用CFD方法做了非定常工况下风力机s的相关仿真^[11-12]. 但CFD方法的一些参数设置会对计算结果有一定的影响, 并且其流场区域网格的绘制费时费力, 计算成本高. 自由涡尾迹(free vortex wake, FVW)方法的计算精确度高于BEM且计算量低于CFD方法. 许波峰等^[13-14]研究了涡核尺寸对FVW方法准确预估气动性能的影响, 保证了该方法的精确性. FVW方法在非定常工况下的研究较多^[15-17], 仿真结果较好. 因此采用自由涡尾迹方法来计算风剪切下的气动性能具有一定的优势.

考虑风剪切的作用对风力机的设计有着重要作用, Kavari等^[18]基于BEM研究了风剪切对叶片设计的影响, 叶根所受影响不大, 而大部分变化发生在叶片长度的0.2到0.8之间. 对风剪切作用下气动性能研究十分重要, 张旭耀等^[19]研究了风剪切来流下水平轴风力机流场特性与风轮气动载荷的分布规律; Shen等^[20]采用升力面法和FVW方法计算了风剪切入流下风力机的气动载荷. 这些研究都指出叶片载荷变化呈现周期性波动, Chen等^[21]基于自由涡尾迹模型和精确梁理论的气动弹性模型研究了利用后掠叶片进行被动载荷控制的方法来减轻载荷变化. 除了对气动性能进行分析, 王海鹏等^[22]研究了风剪切对风力机近尾迹流动特性的影响; 也有****研究对比了定常和风剪切情况下的尾迹结构形状, 风剪切下尾迹呈现出不对称的现象^[23-24]. FVW方法不仅能够精确快速的计算出风剪切条件下的非定常气动特性, 还能够反映尾迹与风轮平面之间的作用. Shaler等^[25]分析了风力机尾迹对其他风力机作用的影响, 其功率与叶片弯矩的不稳定性增大. Su等^[26]使用混合FVW方法探讨了利用倾斜尾迹的转向来减轻对下风向风力机发电潜力影响的问题. Jeong等^[27]采用BEM方法与FVW方法研究出风力机在风剪切和紊流条件下会受到不利的空气载荷和叶片变形, 同时得出了在较低的风速下尾迹对叶片的气动和结构性能有明显的影响, 因此需要对尾迹动力学进行准确评估的结论.

不同的风剪切因子会对叶片气动性能和尾迹形状的畸变产生不同程度的影响, 同时, 尾迹的不对称除了对其下游的风力机产生影响外, 它对自身风轮叶片也存在诱导作用, 具有一定的研究价值. 风轮旋转平面速度分布是来流风速与诱导速度之和, 来流风速主要受风剪切影响, 诱导速度来源于叶片附着涡与尾迹涡的作用. 为了研究风剪切作用下倾斜尾迹对叶片气动性能的影响大小, 本文采用一种"D3P3"格式的时间步进自由涡尾迹方法^[28], 重点在于研究风剪切作用下尾迹形状的畸变对叶片流场的影响, 因此建立了两种尾迹模型, 即剪切入流$+$叶片附着涡$+$倾斜尾迹与剪切入流$+$叶片附着涡$+$对称尾迹, 其示意图如图1所示, 并计算两种模型作用下的叶片气动性能, 为风力机安全性设计和优化设计提供一定的依据.

图1

新窗口打开|下载原图ZIP|生成PPT
图1两种尾迹模型示意图

Fig.1Schematic diagram of two wake models

1 计算方法

1.1 自由涡尾迹模型

风轮坐标规定如图2所示, $z$轴为风轮旋转轴, 以下风向为正, $y$轴竖直向下, $x$轴与$y$, $z$轴符合右手法则.

图2

新窗口打开|下载原图ZIP|生成PPT
图2自由涡尾迹模型

Fig.2Free vortex wake model



采用Weissinger-L升力面模型建立叶片的气动模型, 如图2所示. 由置于1/4弦线处的附着涡线代替叶片, 并按"arccos余弦法"将其离散成一系列直涡线段, 每段附着涡环量由对应叶素的气动特性及Kutta-Joukowski定理确定. 附着涡在空间和时间上的变化分别采用尾随涡和脱体涡来表示, 从附着涡线节点处拖出. 叶素控制点位于叶素中线的3/4弦点处, 涡线对叶素控制点和涡线节点的诱导速度由Biot-Savart定律^[29]确定, 涡线周向诱导速度采用Lamb-Oseen涡核模型修正^[29]来避免产生奇点. 叶片尾迹被划分为近尾迹和远尾迹, 近尾迹在寿命角为60$^\circ$^[30]处截断, 远尾迹由叶尖涡组成, 长度为风轮直径的两倍. 涡线的控制方程偏微分形式为

(1)$\begin{eqnarray} \label{eq1} \frac{\partial { r}\left( {\psi ,\zeta } \right)}{\partial \psi }+\frac{\partial { r}\left( {\psi ,\zeta } \right)}{\partial \zeta }=\frac{1}{\varOmega }\left[ {{ V}_{\infty } +{ V}_{\rm ind} \left( {{ r}\left( {\psi ,\zeta } \right),t} \right)} \right] \end{eqnarray}$
式中, $ r$是涡线节点的矢量位置, $\psi $是叶片旋转时的方位角, $\zeta $是寿命角, $\varOmega $是风轮转速, $ V_{\infty}$是自由流速度, $ V_{\rm ind}$是流场中所有涡线对该节点的总诱导速度.

自由涡尾迹的尾迹求解方法分为松弛迭代法、时间步进法两种. 时间步进法能够计算出尾迹随时间不断发生形状和位置的变化, 用于风剪切条件下风力机气动性能的计算更为合适. 对于空间步($\zeta$)采用五点中心差分法, 对于时间步($\psi$), 采用三步三阶预估—校正格式——"D3P3"格式^[28].

计算风剪切作用下的气动性能及尾迹形状时需要考虑到失速延迟现象和动态失速现象, 因此采用Du-Selig三维旋转效应模型^[31]和Leishman-Beddoes动态失速模型^[32]进行修正来提高计算的准确性.

1.2 风剪切入流模型

地表上空0$\sim$2000 m存在大气边界层, 风速在这个范围内会随高度$z$的变化而变化, 根据IEC标准, 建立正常风廓线模型, 速度可表示为

(2)$\begin{eqnarray} \label{eq2} { V}\left( z \right)={ V}_{\rm hub} \left( {\frac{z}{z_{\rm hub} }}\right)^{\alpha } \end{eqnarray}$
式中$ V(z)$为高度$z$处的风速, $ V_{\rm hub}$为轮毂高度处风速, $z_{\rm hub}$风力机轮毂高度, $\alpha $为风剪切因子.

计算过程中方位角步长取10$^\circ$^[33-34], 将风轮旋转平面的径向方向和方位角方向被划分为36份, 每个节点的速度可用式(3)来表示

(3)$\left. \begin{array} \\ z(l,i)=z_{\rm hub} -\dfrac{R\cdot l}{36}\cdot \sin \dfrac{2\pi (i-1)}{36} \\[3mm] { V}(j,l,i)={ V}_{\infty } \cdot \lt(\dfrac{z(l,i)}{z_{\rm hub}})^{\alpha } \\ \end{array} \right\}$
式中$(l,i)$表示风轮平面上某一点的坐标, $l$为径向方向, $i$为方位角方向, $j$表示时间步步数, $z(l,i)$表示某一点的高度, $ V(j$, $l$, $i)$为风轮平面内某一点在$x$, $y$, $z$方向的风速.

尾迹位置使用"D3P3"格式迭代更新的过程中, 尾迹节点的速度计算式为

(4)$\begin{eqnarray} \label{eq4} { U}={ U}_{\rm BOD} +{ U}_{\rm TIP} +{ U}_{\rm TRD} +{ U}_{\rm SHD} +{ V}(z) \end{eqnarray}$
式中$ U_{\rm BOD}$为附着涡对该节点的诱导速度, $ U_{\rm TIP}$为叶尖涡对该节点的诱导速度, $ U_{\rm TRD}$为尾随涡对该节点的诱导速度, $ U_{\rm SHD}$为脱体涡对该节点的诱导速度.

1.3 计算流程

根据自由涡尾迹模型及风剪切模型, 该方法的计算流程图如图3所示. 其中几何残差用新旧尾迹各节点坐标的均方差表示

(5)$\begin{eqnarray} \label{eq5} R_{\rm MS} =\sqrt {\frac{\sum\limits_{j=1}^{j_{\max } } {\sum\limits_{k=1}^{k_{\max } } {\lt(r_{j,k}^{\rm new} -r_{j,k}^{\rm old})^{2}} } }{j_{\max } k_{\max } }} \end{eqnarray}$
式中$j_{\max}=N_{\rm T}$, $k_{\max}=N_{\rm T}N_{\rm C}$, $N_{\rm T}$为每一圈的时间步数, $N_{\rm C}$为整个尾迹的圈数.

图3

新窗口打开|下载原图ZIP|生成PPT
图3时间步进自由涡尾迹法流程图

Fig.3Flow chart of time-marching free vortex wake methodology

2 模型验证

采用自由涡尾迹模型计算稳态入流条件下NREL Phase VI风力机的气动性能及尾迹形状, 并与实验数据^[35]进行对比. 主要计算7 m/s和13 m/s风速下叶片沿展向分布的法向力系数$C_{\rm n}$和切向力系数$C_{\rm t}$以及不同风速下的载荷. 因为目前主流风力机都是变桨变速型, 大风速下会通过变桨来控制叶片的迎角, 一般不会处于深失速状态, 而Phase VI为失速型风力机, 大风速下处于深失速状态, 气动特性难以准确预估, 所以模型验证中不考虑大风速下的气动性能.

7 m/s风速下计算结果如图4所示, 该方法计算得到的法向力系数、切向力系数与实验数据误差较小, 计算结果较为精确. 图5为13 m/s风速下计算得到的气动力系数, 可以看到叶根处的计算结果与实验数据的误差较大, 这是由于叶根相比于叶片其他位置更容易进入失速状态, 其周围的流体分离使得流场变得复杂, 气动载荷的误差则会增大. 从叶中至叶尖, 法向力系数、切向力系数的计算结果与实验结果基本吻合.

图4

新窗口打开|下载原图ZIP|生成PPT
图47 m/s风速下的法向力系数与切向力系数

Fig.4Normal force coefficient and tangential force coefficient under 7 m/s wind speed

图5

新窗口打开|下载原图ZIP|生成PPT
图513 m/s风速下的法向力系数与切向力系数

Fig.5Normal force coefficient and tangential force coefficient under 13 m/s wind speed



图6为不同风速下风力机推力和转矩的计算值与实验值的对比. 可以看到随着风速的增大, 风轮总体气动载荷的变化趋势与实验数据保持一致, 10 m/s风速以下计算结果与实验值高度吻合, 10 m/s风速以上由于叶根出现了部分失速现象, 法向力和切向力预测存在偏差,导致总体气动载荷计算存在一定的误差.

图6

新窗口打开|下载原图ZIP|生成PPT
图6不同风速下的推力和转矩

Fig.6Thrust and torque at different wind speeds



在叶尖速比为10, 桨距角为4$^\circ$ 的条件下, 尾迹形状的输出及其与实验结果的对比如图7所示, 可以看到自由涡尾迹方法能够很好地仿真出尾迹结构, 说明它能够研究风剪切下尾迹形状的变化.

图7

新窗口打开|下载原图ZIP|生成PPT
图7尾迹形状验证

Fig.7Verification of wake shape



总的来数, 上述的计算结果大部分在误差范围内, 所以该自由涡尾迹模型能够用于研究风轮的气动载荷及尾迹形状.

3 风剪切下叶片的气动性能与尾迹形状

采用973项目"大型风力机的空气动力学基础研究"中1.5 MW风力机^[36]. 该风力机的叶片数目为3, 轮毂高度65 m, 风轮直径83 m, 额定转速17.2 r/min. 在风力机运行的环境中, 风剪切因子会受地形、粗糙度、温度、季节、时间等各种因素影响, 它的取值可以在0.1$\sim$0.5之间^[37-38]. 本文计算的入流条件为: 轮毂高度处风速为额定风速10.4m/s, 此时风轮转速为额定转速, 风剪切因子分别取 0.1, 0.2, 0.3, 0.4, 0.5, 其入流风速随高度的变化如图8所示. 随着风剪切因子的增大, 风速在风轮平面最高点和最低点的速度差会随之增大, 当$\alpha $为0.1时, 风速差达1.5 m/s; 当$\alpha $为0.5时, 风速差达7 m/s.

图8

新窗口打开|下载原图ZIP|生成PPT
图8不同风剪切因子下的入流情况

Fig.8Inflow conditions under different wind shear coefficients

3.1 风剪切作用下的气动性能

风力机在风剪切作用下的实际气动性能会影响到发电量和风力机的运转, 因此需要对该工况下气动性能进行分析以便为风力机的叶片设计提供相关依据. 图9为风剪切因子从0.1至0.5条件下, 叶片在80%截面处的法向力系数和切向力系数. 可以看出风剪切的作用使气动力系数随时间做周期性波动, 这与入流风速在高度上的不均匀分布及叶片旋转过程中经过风轮平面的最低点与最高点有关, 这些波动会使风力机承受疲劳载荷从而影响其安全性及使用寿命, 风力机在设计过程中考虑风剪切作用显得十分重要. 随着风剪切因子的增大, 气动力系数的波动幅度明显变大, $\alpha$在0.1, 0.3, 0.5时的波动幅度分别为0.1, 0.29, 0.51, 风力机在设计时应充分考虑所选厂址的风剪切因子来进行设计.

图9

新窗口打开|下载原图ZIP|生成PPT
图9不同风剪切因子下的法向力系数和切向力系数

Fig.9Normal force coefficient and tangential force coefficient under different wind shear coefficients



图10为风剪切因子从0.1至0.5条件下风力机所受的推力. 随着风剪切因子的增大, 风力机所受推力平均值逐渐减小, 推力波动更为剧烈, 说明风剪切较大的叶片机组与无风剪切作用下相比, 无论是平均值还是波动幅值都相差较大. 因此, 风力机的设计更应该准确计算风场的风剪切因子并考虑风剪切因子对载荷的影响.

图10

新窗口打开|下载原图ZIP|生成PPT
图10不同风剪切因子下的推力

Fig.10Thrust under different wind shear coefficients

3.2 风剪切作用下尾迹形状的变化

在给定的入流条件下, 计算了不同风剪切因子作用于风力机后的叶尖涡尾迹形状, 如图11所示, 3种不同颜色分别代表了3个叶片的尾迹. 从图示结果可以看出风剪切的存在会使叶尖涡尾迹呈现出不对称、倾斜的形状. 根据本文规定的坐标系, 在风剪切的作用下, $y$轴负半轴所受的来流风速大于正半轴所受的来流风速, 转动到负半轴的尾迹节点的运动速度更快一些, 因此会出现图示所示的倾斜形状. 从计算结果可以看到随着风剪切因子的增大, 尾迹的不对称现象和倾斜程度更加显著, $\alpha $在0.4和0.5时远尾迹会出现向上卷起的现象, 这些都与风轮平面最高与最低风速差的增大有关.

图11

新窗口打开|下载原图ZIP|生成PPT
图11不同风剪切因子的尾迹形状

Fig.11Wake shape under different wind shear coefficients



图12给出了不同风剪切因子作用下, 叶片1尾迹在第4圈的最低点和最高点的轴向相对位置. 由图可得出, 随着风剪切因子的增大, 最低点的轴向相对位置逐渐减小, 而最高点的轴向相对位置逐渐增大, 即尾迹的倾斜程度越来越明显.

图12

新窗口打开|下载原图ZIP|生成PPT
图12叶片1尾迹的第4圈最低点和最高点的轴向相对位置

Fig.12Axial relative position of the lowest point and the highest point of the wake of blade 1 in turn 4



图13给出了不同风剪切因子作用下, 叶片1尾迹在第1圈和第4圈的最低点距离叶尖涡尾迹中心线的相对位移. 由图可得出, 第1圈最低点距中心线的相对位移基本不变, 而随着风剪切因子的增大, 第4圈最低点距中心线的相对位移逐渐减小, 说明尾迹向上卷起的现象越明显.

图13

新窗口打开|下载原图ZIP|生成PPT
图13叶片1尾迹的第1圈和第4圈最低点距尾迹中心线的相对位移

Fig.13Relative displacement of the lowest point in the first lap and the fourth lap of the wake of blade 1 from the centerline

3.3 对称尾迹与倾斜尾迹下的气动性能

风剪切作为来流风速, 对风力机叶片气动性能产生一定的影响, 同时, 其导致的畸变倾斜的尾迹会反过来作用于风轮平面. 为了研究倾斜尾迹对风力机叶片气动性能的影响, 建立了如图1所示的两个尾迹模型进行对比, 尾迹分别为风剪切作用下的真实倾斜尾迹与假设来流风速均等于轮毂处风速且无风剪切的虚拟对称尾迹, 两种模型下的入流风均为剪切风.

首先计算了$\alpha $为0.3时, 倾斜尾迹和对称尾迹分别作用下风轮旋转平面的轴向诱导速度因子分布情况, 如图14和图15所示. 可以看到倾斜尾迹对风轮平面的诱导速度因子存在明显的不对称现象, 而对称尾迹对风轮平面的诱导速度因子的分布较为均匀.

图14

新窗口打开|下载原图ZIP|生成PPT
图14倾斜尾迹对风轮旋转平面的轴向诱导速度因子

Fig.14Effect of tilted wake on axial induced velocity factor of wind turbine plane

图15

新窗口打开|下载原图ZIP|生成PPT
图15对称尾迹对风轮旋转平面的轴向诱导速度因子

Fig.15Effect of symmetrical wake on axial induced velocity factor of wind turbine plane



图16给出了$\alpha$ 为0.3时, 叶片80%截面处在对称与倾斜尾迹诱导作用下的法向力系数和切向力系数. 从图中的计算结果来看, 无论尾迹是对称的还是倾斜的, 在风剪切来流条件下都会呈现出随时间的周期性波动. 对称尾迹下法向力系数与切向力系数的标准差分别为0.098 48和0.028 45, 倾斜尾迹下法向力系数与切向力系数的标准差分别为0.111 13和0.031 73, 可见, 对称尾迹作用下的气动力系数波动幅度比倾斜尾迹的波动幅度小. 尾迹在轮毂下方的倾斜程度更加明显, 所以可以看到图中气动力系数在波谷的偏差比波峰的大. 图17为两种尾迹作用下气动力系数在波谷处的值的大小, 即叶片转至向下时80%截面处气动力系数的对比. 随着$\alpha $的增大, 对称尾迹与倾斜尾迹对风轮最下方叶片尖部的气动载荷影响差别越明显, 说明尾迹越倾斜, 风轮旋转平面处的载荷不对称性越明显.

图16

新窗口打开|下载原图ZIP|生成PPT
图16对称与倾斜尾迹下的法向力系数和切向力系数

Fig.16Normal force coefficient and tangent force coefficient under the influence of symmetrical and tilted wake

图17

新窗口打开|下载原图ZIP|生成PPT
图17两种尾迹的法向力系数和切向力系数在不同风剪切因子下偏差

Fig.17Deviation of normal force coefficient and tangential force coefficient of two kinds of wake under different wind shear coefficients



图18为$\alpha $为0.3时, 对称与倾斜尾迹作用下的推力. 由图中结果可知, 无论是倾斜尾迹还是对称尾迹的诱导, 在风剪切作用下推力都呈现出随时间的波动. 对称尾迹诱导下风力机所受的推力大于倾斜尾迹下风力机所受的推力, 说明尾迹的变形会使风力机的总体性能减小且偏离较大, 所以在风力机设计时需要考虑到尾迹变形对风力机性能的影响.

图18

新窗口打开|下载原图ZIP|生成PPT
图18对称与倾斜尾迹作用下的推力

Fig.18Thrust under symmetrical and tilted wake



因此, 在计算风剪切作用下叶片的气动性能时, 除了要准确计算风剪切对叶片入流的影响, 还要准确预估倾斜的叶片尾迹形状.

4 结论

本文基于时间步进自由涡尾迹方法仿真了风剪切入流工况及不同风剪切因子作用下, 风轮平面气动性能以及尾迹形状的变化, 同时拟定了一个对称尾迹来研究尾迹形状变化对叶片气动性能的影响. 通过对比分析得到如下结论:

(1)风剪切入流条件下, 风力机叶片叶素的气动特性及叶片所受推力都随时间做周期性波动.

(2)轮毂处风速相同的条件下, 风剪切较大的叶片机组与无风剪切作用下相比, 无论是平均值还是波动幅值都相差较大, 随着风剪切因子的增大, 风力机的气动特性波动的幅度越大, 所受推力越小, 同时, 尾迹不对称、倾斜的形状会更加明显.

(3)相同入流条件下, 倾斜尾迹与拟定的对称尾迹相比, 风轮平面的轴向诱导速度分布不均匀, 气动力系数的波动更大, 风力机所受平均推力更小, 且随着风剪切因子的增大, 对风轮叶片尖部的气动载荷影响差别越明显. 尾迹越倾斜, 风轮旋转平面处的载荷不对称性越明显.

(4)风剪切作用下, 要提高风力机动态载荷计算的精度, 正确预估倾斜的叶片尾迹形状与准确的风剪切入流分布同样重要.

参考文献 View Option 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1] 温斌荣, 魏莎, 魏克湘 等. 风切变和塔影效应对风力机输出功率的影响
机械工程学报, 2018,54(10):124-132
[本文引用: 1]

( Wen Binrong, Wei Sha, Wei Kexiang, et al. Influences of wind shear and tower shadow on the power output of wind turbine
Journal of Mechanical Engineering. 2018,54(10):124-132 (in Chinese))
[本文引用: 1]

[2] 吴斌, 董礼, 廖明夫 等. 水平轴风力机气动计算的叶素动量修正法
机械科学与技术, 2011,30(12):2124-2128, 2134
[本文引用: 1]

( Wu Bin, Dong Li, Liao Mingfu, et al. Amended BEM method for aerodynamic analysis of wind turbine
Mechanical Science and Technology for Aerospace Engineering. 2011,30(12):2124-2128, 2134 (in Chinese))
[本文引用: 1]

[3] Holierhoek JG, de Vaal JB, van Zuijlen AH. et al. Comparing different dynamic stall models
Wind Energy, 2013,16(1):139-158


[4] Liu X, Lu C, Liang S. et al. Improved dynamic stall prediction of wind turbine airfoils
Energy Procedia, 2019,158:1021-1026


[5] Kabir IFSA, Ng EYK. Insight into stall delay and computation of 3D sectional aerofoil characteristics of NREL phase VI wind turbine using inverse BEM and improvement in BEM analysis accounting for stall delay effect
Energy, 2017,120:518-536


[6] Rajan A, Ponta FL. Aeroelastic analysis of the 3-dimensional interference patterns of wind-turbine rotors: The 3-D DRD-BEM model
Renewable Energy Focus, 2018,26:22-39


[7] Pinto RLUD, Goncalves BPF. A revised theoretical analysis of aerodynamic optimization of horizontal-axis wind turbines based on BEM theory
Renewable Energy, 2017,105:625-636
[本文引用: 1]

[8] 唐新姿, 王效禹, 袁可人 等. 风速不确定性对风力机气动力影响量化研究
力学学报, 2020,52(1):51-59
[本文引用: 1]

( Tang Xinzi, Wang Xiaoyu, Yuan Keren, et al. Quantitation study of influence of wind speed uncertainty on aerodynamic forces of wind turbine
Chinese Journal of Theoretical and Applied Mechanics. 2020,52(1):51-59 (in Chinese))
[本文引用: 1]

[9] Ali QS, Kim MH. Unsteady aerodynamic performance analysis of an airborne wind turbine under load varying conditions at high altitude
Energy Conversion and Management, 2020,210:112696
[本文引用: 1]

[10] Wei Y, Ferreira CS, van Kuik G, et al. Verifying the blade element momentum method in unsteady, radially varied, axisymmetric loading using a vortex ring model
Wind Energy, 2017,20(2):269-288
[本文引用: 1]

[11] Cai X, Gu RR, Pan P. et al. Unsteady aerodynamics simulation of a full-scale horizontal axis wind turbine using CFD methodology
Energy Conversion and Management, 2016,112:146-156
[本文引用: 1]

[12] Abbaspour M, Radmanesh AR, Soltani MR. Unsteady flow over offshore wind turbine airfoils and aerodynamic loads with computational fluid dynamic simulations
International Journal of Environmental Science and Technology, 2016,13(6):1525-1540
[本文引用: 1]

[13] 许波峰, 刘冰冰, 冯俊恒 等. 自由涡尾迹方法中涡核尺寸对风力机气动计算的影响
力学学报, 2019,51(5):1530-1537
[本文引用: 1]

( Xu Bofeng, Liu Bingbing, Feng Junheng, et al. Influence of vortex core size on aerodynamic calculation of wind turbine in free vortex wake method
Chinese Journal of Theoretical and Applied Mechanics. 2019,51(5):1530-1537 (in Chinese))
[本文引用: 1]

[14] Xu BF, Liu BB, Cai X. et al. Accuracy of the aerodynamic performance of wind turbines using vortex core models in the free vortex wake method
Journal of Renewable and Sustainable Energy, 2019,11(5):053307
[本文引用: 1]

[15] Shen X, Hu P, Chen JG. et al. The unsteady aerodynamics of floating wind turbine under platform pitch motion
Journal of Power and Energy, 2018,232(8):1019-1036
[本文引用: 1]

[16] 林易, 李晔, 段磊. 漂浮式风力机非定常气动特性分析
船舶与海洋工程, 2019,35(6):8-14


( Lin Yi, Li Ye, Duan Lei. Analysis on the unsteady aerodynamics of floating wind turbines
Naval Architecture and Ocean Engineering. 2019,35(6):8-14 (in Chinese))


[17] Lee H, Lee DJ. Wake impact on aerodynamic characteristics of horizontal axis wind turbine under yawed flow conditions
Renewable Energy, 2019,136:383-392
[本文引用: 1]

[18] Kavari G, Tahani M, Mirhosseini M. Wind shear effect on aerodynamic performance and energy production of horizontal axis wind turbines with developing blade element momentum theory
Journal of Cleaner Production, 2019,219:368-376
[本文引用: 1]

[19] 张旭耀, 杨从新, 李寿图 等. 风剪切来流下风力机流场特性与风轮气动载荷研究
太阳能学报, 2019,40(11):3281-3288
[本文引用: 1]

( Zhang Xuyao, Yang Congxin, Li Shoutu, et al. Study on flow field characteristics and aerodynamic loads of horizontal axis wind turbine in shear inflow
Acta Energiae Solaris Sinica. 2019,40(11):3281-3288 (in Chinese))
[本文引用: 1]

[20] Shen X, Zhu XC, Du ZH. Wind turbine aerodynamics and loads control in wind shear flow
Energy, 2011,36(3):1424-1434
[本文引用: 1]

[21] Chen JG, Shen X, Zhu XC. et al. A study on the capability of backward swept blades to mitigate loads of wind turbines in shear flow
Journal of Energy Resources Technology-Transactions of the ASME, 2019,141(8):081201
[本文引用: 1]

[22] 王海鹏, 邱庆刚, 张博. 风剪切来流对风力机近尾迹流动特性影响
热科学与技术, 2017,16(2):120-126
[本文引用: 1]

( Wang Haipeng, Qiu Qinggang, Zhang Bo. Effects of wind shear on near wake characteristics of wind turbine blade
Journal of Thermal Science and Technology. 2017,16(2):120-126 (in Chinese))
[本文引用: 1]

[23] 周文平, 唐胜利, 吕红. 风剪切和动态来流对水平轴风力机尾迹和气动性能的影响
中国电机工程学报, 2012,32(14):122-127
[本文引用: 1]

( Zhou Wenping, Tang Shengli, Lü Hong. Effect of transient wind shear and dynamic inflow on the wake structure and performance of horizontal axis wind turbine
Proceedings of the CSEE. 2012,32(14):122-127 (in Chinese))
[本文引用: 1]

[24] Sezer-Uzol N, Uzol O. Effect of steady and transient wind shear on the wake structure and performance of a horizontal axis wind turbine rotor
Wind Energy, 2013,16(1):1-17
[本文引用: 1]

[25] Shaler K, Kecskemety KM, McNamara JJ. Benchmarking of a free vortex wake model for prediction of wake interactions
Renewable Energy, 2019,136:607-620
[本文引用: 1]

[26] Su KY, Bliss D. A numerical study of tilt-based wake steering using a hybrid free-wake method
Wind Energy, 2020,23(2):258-273
[本文引用: 1]

[27] Jeong MS, Kim SW, Lee I. et al. Wake impacts on aerodynamic and aeroelastic behaviors of a horizontal axis wind turbine blade for sheared and turbulent flow conditions
Journal of Fluids and Structures, 2014,50:66-78
[本文引用: 1]

[28] Xu BF, Wang TG, Yuan Y. et al. Unsteady aerodynamic analysis for offshore floating wind turbines under different wind conditions
Philosophical Transactions of the Royal Society A-Mathematical Physical and Engineering Sciences, 2015,373(2035):20140080
[本文引用: 2]

[29] 左潞, 唐植懿, 许波峰 等. 涡模型对风力机气动特性的影响研究
可再生能源, 2016,34(10):1491-1496
[本文引用: 2]

( Zuo Lu, Tang Zhiyi, Xu Bofeng, et al. Investigation of effects of vortex models on wind turbine aerodynamic characteristics
Renewable Energy Resources. 2016,34(10):1491-1496 (in Chinese))
[本文引用: 2]

[30] Xu BF, Wang TG, Yuan Y. et al. A simplified free vortex wake model of wind turbines for axial steady conditions
Applied Sciences, 2018,8(6):866
[本文引用: 1]

[31] Elgammi M, Sant T. A new stall delay algorithm for predicting the aerodynamics loads on wind turbine blades for axial and yawed conditions
Wind Energy, 2017,20(9):1645-1663
[本文引用: 1]

[32] 李国强, 张卫国, 陈立 等. 风力机叶片翼型动态试验技术研究
力学学报, 2018,50(4):751-765
[本文引用: 1]

( Li Guoqiang, Zhang Weiguo, Chen Li, et al. Research on dynamic test technology for wind turbine blade airfoil
Chinese Journal of Theoretical and Applied Mechanics. 2018,50(4):751-765 (in Chinese))
[本文引用: 1]

[33] Bhagwat M, Leishman JG. Stability, consistency and convergence of time marching free-vortex rotor wake algorithms
Journal of the American Helicopter Society, 2001,46(1):59-71
[本文引用: 1]

[34] Gupta S, Leishman JG. Accuracy of the induced velocity from helicoidal vortices using straight-line segmentation
AIAA Journal, 2005,43(1):29-40
[本文引用: 1]

[35] Hand MM, Simms DA, Fingersh LJ, et al. Unsteady aerodynamics experiment phase VI: Wind tunnel test configurations and available data campaign
Technical Report NREL/TP-500-29955, National Renewable Energy Laboratory, Golden, 2001
[本文引用: 1]

[36] Wang TG, Wang L, Zhong W. et al. Large-scale wind turbine blade design and aerodynamic analysis
Chinese Science Bulletin, 2012,57(1):1-7
[本文引用: 1]

[37] Werapun W, Tirawanichakul Y, Waewsak J. Wind shear coefficients and their effect on energy production
Energy Procedia, 2017,138:1061-1066
[本文引用: 1]

[38] Okorie ME, Inambao F, Chiguvare Z. Evaluation of wind shear coefficients, surface roughness and energy yields over inland locations in Namibia
Procedia Manufacturing, 2016. 7:630-638
[本文引用: 1]