一种风力机涡尾迹修正方法技术

本发明专利技术涉及风力机测试技术领域，公开了一种风力机涡尾迹修正方法，通过获取风力机叶片数量、风轮半径、风轮转速和叶片弦长，根据风力机运行当地的主风向的风速、空气密度及空气粘度，确定风力机叶片的来流风速和迎角，计算涡尾迹形状，根据尾迹角区域使用不同涡核模型更加准确合理地描述尾迹流场。在初始尾迹角附近时，采用Scully涡模型；向前、向后的尾迹角范围，采用Lamb‑Oseen涡模型；根据流场尾迹角区域选用不同的有效粘性扩散因子。本发明专利技术根据尾迹流场时间历程和空间分布的差异，对一般涡尾迹计算的尾流模型进行修正，使涡尾迹计算的结果更合理、准确。

A Method of Wind Turbine Vortex Wake Correction

The invention relates to the technical field of wind turbine testing, and discloses a method for correcting wind turbine vortex wake. By obtaining the number of wind turbine blades, the radius of wind turbine, the speed of wind turbine and the chord length of the blade, according to the wind speed, air density and air viscosity of the local main wind direction in which the wind turbine operates, the wind speed and angle of attack of the wind turbine blades are determined, and the shape of the vortex wake is calculated, according to the wake angle region. Different vortex core models are used to describe the wake flow more accurately and reasonably. Scully vorticity model is used near the initial wake angle, Lamb_Oseen vorticity model is used in the forward and backward wake angle range, and different effective viscous diffusion factors are selected according to the wake angle region of the flow field. According to the difference of time history and spatial distribution of wake flow field, the wake model of general vorticity wake calculation is modified to make the result of vorticity wake calculation more reasonable and accurate.

【技术实现步骤摘要】
一种风力机涡尾迹修正方法
本专利技术涉及风力机测试
，具体涉及一种风力机涡尾迹修正方法。
技术介绍
风力机气动特性计算方法主要有三种：叶素动量理论、涡尾迹方法、计算流体力学方法。综合考虑方法的准确性和计算成本，涡尾迹方法更具优势，而且其本质上具有旋涡特性，能更准确地计及尾流场尾涡间相互诱导作用，是模拟风力机气动特性较为灵活的数值工具。自由涡尾迹方法不需要流场的先验数据，直接通过解涡线控制方程得到尾迹结构，有较强的理论根据，是涡尾迹方法中应用最广泛的。自由涡尾迹方法中，涡线从叶片尾缘拖出流向下游，涡线的控制节点随当地流速自由移动，通过求解涡线控制方程得到最终的尾迹几何形状。目前，一般采用自由涡尾迹建模计算风力机气动特性，尾迹形状的描述均是采用上述从叶片尾缘拖出自由移动至远场的方式。这些模型中涡尾迹的计算采用单一涡核模型，并根据某一固定的有效粘性扩散因子作为代表，风力机涡尾迹方法过于简化。因此亟需一种根据尾迹流场时间历程和空间分布的差异，对一般涡尾迹计算的尾流模型进行修正，使涡尾迹计算的结果更合理、准确。
技术实现思路
基于以上问题，本专利技术提供一种风力机涡尾迹修正方法，根据尾迹流场时间历程和空间分布的差异，对一般涡尾迹计算的尾流模型进行修正，使涡尾迹计算的结果更合理、准确。为解决以上技术问题，本专利技术提供了一种风力机涡尾迹修正方法，包括如下步骤：S1：获取风力机叶片数量、风轮半径、风轮转速和叶片弦长；S2：根据风力机运行当地的主风向的风速、空气密度及空气粘度，确定风力机叶片的来流风速和迎角；S3：计算涡尾迹形状，在初始尾迹角附近，尾迹角的范围为30°～320°采用Scully涡模型；尾迹角30°以前、320°以后尾迹角范围，采用Lamb-Oseen涡模型；涡核模型计算诱导速度公式为：其中r为风轮半径，Γ为涡强，rc为对应于周向诱导速度最大位置的涡核半径；n的取值为1或2；S4：根据流场尾迹角区域选用不同的有效粘性扩散因子；计算涡周围的切向旋转速度，计算公式为：rc为粘性核半径，为用rc无量纲化的半径；rc随时间t变化关系为：rc随尾迹角ζ变化关系为：其中，γ为运动粘性系数，γ＝1.4607×10-5m2/s；α＝1.25643；δ为有效粘性扩散因子，其取值为从叶尖涡“卷起”阶段结束位置(尾迹角为240°)开始到其后120°范围，取值为0～2，尾迹角超过360°开始，增长到10以上；Ω为风轮旋转速度，r0为叶尖涡“卷起”阶段结束位置的初始涡核半径，ζ0为叶尖涡“卷起”阶段结束位置的尾迹角。进一步地，步骤S4中的有效粘性扩散因子基于试验数据，具体来说是在流场中播撒示踪粒子，采用粒子图像测速仪(也即PIV)测量两幅间隔非常短的(通常为微秒量级)粒子图像，通过对图像进行互相关计算得到流场中的速度矢量分布，对其进行后处理得到有效粘性扩散因子的实际测量值或平均值。与现有技术相比，本专利技术的有益效果是：(1)根据尾迹角区域使用不同涡核模型更加准确合理地描述尾迹流场。在初始尾迹角附近时，采用Scully涡模型；向前、向后的一定尾迹角范围，采用Lamb-Oseen涡模型。(2)根据流场尾迹角区域选用不同的有效粘性扩散因子，更准确地描述涡核内部涡量向外的输运。有效粘性扩散因子取值为从叶尖涡“卷起”阶段结束位置(尾迹角为240°)开始到其后120°范围，取值为0～2，尾迹角超过360°开始，增长到10以上。(3)基于试验数据或根据试验规律确定的平均值，有效粘性扩散因子的量值具有充分依据。附图说明图1为实施例中风力机涡尾迹修正方法的流程图；具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本专利技术作进一步的详细说明，本专利技术的示意性实施方式及其说明仅用于解释本专利技术，并不作为对本专利技术的限定。实施例：参见图1，一种风力机涡尾迹修正方法，包括如下步骤：S1：获取风力机叶片数量、风轮半径、风轮转速和叶片弦长；S2：根据风力机运行当地的主风向的风速、空气密度及空气粘度，确定风力机叶片的来流风速和迎角；S3：计算涡尾迹形状，在初始尾迹角附近，尾迹角的范围为30°～320°采用Scully涡模型；尾迹角30°以前、320°以后尾迹角范围，采用Lamb-Oseen涡模型；涡核模型计算诱导速度公式为：其中r为风轮半径，Γ为涡强，rc为对应于周向诱导速度最大位置的涡核半径；n的取值为1或2，n＝1对应的是Scully涡核模型，n＝2时对应的是Lamb-Oseen涡核模型。S4：根据流场尾迹角区域选用不同的有效粘性扩散因子；计算涡周围的切向旋转速度，计算公式为：为用rc无量纲化的半径；rc为粘性核半径，rc随时间t变化关系为：rc随尾迹角ζ变化关系为：其中，γ为运动粘性系数，γ＝1.4607×10-5m2/s；α＝1.25643；δ为有效粘性扩散因子，其取值为从叶尖涡“卷起”阶段结束位置(尾迹角为240°)开始到其后120°范围，取值为0～2，尾迹角超过360°开始，增长到10以上；Ω为风轮旋转速度，r0为叶尖涡“卷起”阶段结束位置的初始涡核半径，ζ0为叶尖涡“卷起”阶段结束位置的尾迹角。步骤S4中的有效粘性扩散因子基于试验数据，具体来说是在流场中播撒示踪粒子，采用粒子图像测速仪(也即PIV)测量两幅间隔非常短的(通常为微秒量级)粒子图像，通过对图像进行互相关计算得到流场中的速度矢量分布，对其进行后处理得到有效粘性扩散因子的实际测量值或平均值。如上即为本专利技术的实施例，上述实施例以及实施例中的具体参数仅是为了清楚表述专利技术验证过程，并非用以限制本专利技术的专利保护范围，本专利技术的专利保护范围仍然以其权利要求书为准，凡是运用本专利技术的说明书及附图内容所作的等同结构变化，同理均应包含在本专利技术的保护范围内。本文档来自技高网...

【技术保护点】
1.一种风力机涡尾迹修正方法，其特征在于，包括如下步骤：S1：获取风力机叶片数量、风轮半径、风轮转速和叶片弦长；S2：根据风力机运行当地的主风向的风速、空气密度及空气粘度，确定风力机叶片的来流风速和迎角；S3：计算涡尾迹形状，在初始尾迹角附近，尾迹角的范围为30°～320°采用Scully涡模型；尾迹角30°以前、320°以后尾迹角范围，采用Lamb‑Oseen涡模型；涡核模型计算诱导速度公式为：

【技术特征摘要】
1.一种风力机涡尾迹修正方法，其特征在于，包括如下步骤：S1：获取风力机叶片数量、风轮半径、风轮转速和叶片弦长；S2：根据风力机运行当地的主风向的风速、空气密度及空气粘度，确定风力机叶片的来流风速和迎角；S3：计算涡尾迹形状，在初始尾迹角附近，尾迹角的范围为30°～320°采用Scully涡模型；尾迹角30°以前、320°以后尾迹角范围，采用Lamb-Oseen涡模型；涡核模型计算诱导速度公式为：其中r为风轮半径，Γ为涡强，rc为对应于周向诱导速度最大位置的涡核半径；n的取值为1或2；S4：根据流场尾迹角区域选用不同的有效粘性扩散因子；计算涡周围的切向旋转速度，计算公式为：rc为粘性核半径，为用rc无量纲化的半径；rc随时间t变化关系为：rc随...

【专利技术属性】
技术研发人员：武杰，陈立，肖京平，王强，史雨，李国强，金华，王桥，赵献礼，
申请(专利权)人：中国空气动力研究与发展中心低速空气动力研究所，
类型：发明
国别省市：四川,51