一种风机叶片涡流发生器优化方法及系统技术方案

本发明专利技术公开了一种风机叶片涡流发生器优化方法，以风机叶片模型周围区域作为控制区域并网格化；构建网格的控制方程得到风机叶片在来流风速下的第一发电功率；设置涡流发生器的形状及位置，在控制方程中引入涡流发生器源项得到修正后的控制方程，得到涡流发生器处于该位置时风机叶片修正的发电功率；以此得到不同形状的涡流发生器处于不同位置时的多个第二修正发电功率；将每一个第二修正发电功率分别与第一发电功率进行比较，得到不同形状的涡流发生器处于不同位置时的多个无量纲特征功率，多个无量纲特征功率中最大者对应的涡流发生器的位置和形状作为最佳方案。本发明专利技术提供了对涡流发生器进行优化的方法。涡流发生器进行优化的方法。涡流发生器进行优化的方法。

【技术实现步骤摘要】
一种风机叶片涡流发生器优化方法及系统


[0001]本专利技术涉及风力发电
，尤其涉及一种风机叶片涡流发生器优化方法及系统。

技术介绍

[0002]近年来，风力发电行业正在快速发展。由于风向与风速的随机性，风机叶片易在短时间内运行在大攻角工况，导致叶片表面出现气流分离现象，同时叶根多采用厚翼型来满足强度要求，在大攻角下其叶根也容易发生气流分离，导致叶片性能下降。对于风机叶片而言，新技术的研发用以提升叶片的气动性能并实现机组发电量的提升，是必需的。在风机叶片的众多提质增效技术中，涡流发生器(Vortex Generator，简称VG)是相对成熟且广泛应用的技术之一。
[0003]安装涡流发生器是控制流动分离的一种低成本且有效的方法，最早是针对机翼在大攻角条件控制失速而提高机动性的一个经典设计。在大攻角情况 下，涡流发生器能够大幅度增加翼段的升力系数，减小阻力系数。在风力发电机组上合理安装涡流发生器，能够有效地提高风轮的输出功率和轴向推力，从而能够降低风力发电机组叶片的制造成本，使风能较其他能源产业更具竞争力。涡流发生器自身型阻可能导致翼型总体阻力的增加，影响风电机组的发电量。
[0004]公开号为CN105556114A的专利给出了完整的涡流发生器参数与安装方案，但是该专利并没有提出针对不通叶片的涡流发生器的具体设计和优化方法，随着叶片设计的不断更新发展，该专利提出的方案已经难以再适应新型叶片上涡流发生器的设计。公开号为CN108280259A的专利，提出了一种基于计算流体力学CFD方法对涡流发生器进行优化设计的方案，该方案需要对涡流发生器的局部结构进行CAD建模，使计算资源的消耗很大，对设计进行更新优化的难度较大。公开号CN111859651A的专利，通过使用叶素动量法则对叶片的气动性能进行快速计算和优化设计，但是该方法简化过多，无法准确体现出涡流发生器对叶片性能的具体影响，故优化的效果也不明确。公开号为CN112943565的专利，提出了针对一种具体的涡流发生器形式，即波浪形涡流发生器的优化设计方法，无法针对其他的涡流发生器提出设计优化方案。

技术实现思路

[0005]鉴于此，本专利技术提供一种风机叶片涡流发生器优化方法及系统，使用计算流体力学方法对涡流发生器的具体参数进行分析，从而提供对任意设计参数下的涡流发生器性能评估标准。
[0006]为实现上述目的，本专利技术提供一种风机叶片涡流发生器优化方法，所述方法包括步骤：S1、构建风机叶片模型，以环绕该风机叶片模型周围的三维空间区域作为控制区域，将所述控制区域进行网格化；
S2、以控制区域边界处的来流风速和风机叶片表面空气速度作为控制方程的边界条件，构建每一个网格的控制方程，计算得到对应网格的风速和压强，以风机叶片模型表面网格的风速和压强计算得到风机叶片表面的压强和切应力，并得到风机叶片在该来流风速下的第一发电功率；S3、设置涡流发生器的形状，并将该涡流发生器设置在风机叶片表面的任一位置，获取该涡流发生器所占据的网格空间位置，在每一个网格的控制方程中引入涡流发生器源项，以得到每一个网格修正后的控制方程；S4、基于每一个网格修正后的控制方程，计算得到对应网格的修正后的风速和压强，以风机叶片模型表面网格的修正后的风速和压强计算得到风机叶片表面修正后的压强和切应力，并得到当所述涡流发生器处于所述位置时风机叶片在所述来流风速下修正的发电功率；S5、依次设置涡流发生器在风机叶片表面的各个不同位置，以及设置涡流发生器的不同形状，重复执行步骤S3
‑
S4，得到不同形状的涡流发生器处于风机叶片的不同位置时所对应的风机叶片的多个第二修正发电功率；S6、将每一个第二修正发电功率分别与第一发电功率进行比较，得到不同形状的涡流发生器处于风机叶片的不同位置时所对应的多个无量纲特征功率，取多个无量纲特征功率中最大者所对应的涡流发生器的位置和形状作为涡流发生器的最佳设计方案。
[0007]优选的，所述步骤S1包括：构建风机叶片模型的三维空间坐标体系，以风机轮毂的原点中心为坐标体系的原点，x轴方向垂直于风机叶轮指向风机的下游方向，以重力的方向作为z轴方向，y轴方向与x轴方向和z轴方向垂直。
[0008]优选的，所述步骤S1包括：所述控制区域是以风机叶片的质量中心为原点，分别向x、y、z轴三个轴方向上扩充6倍叶片弦长距离的三维空间区域。
[0009]优选的，所述步骤S2包括：以控制区域边界处的来流风速作为控制方程的边界条件： （1）；风机叶片表面空气速度作为控制方程的边界条件，风机叶片表面空气速度为： （2）；其中，为边界点到轮毂中心的矢量位移，η为来流风速下叶片的叶尖速比，R为叶轮半径；对于控制区域中网格的控制方程为：
ꢀ
（3）；其中， （4）； （5）； （6）；其中，t为时间，（u
x
, u
y
, u
z
）为网格风速在x,y,z三个轴向上的速度，p为网格的压强，ρ为空气的密度，ν
a
为空气的有效动力学粘度，ν
t
为空气的湍流动力学粘度，g
z
为z轴方向的重力加速度，U为网格风速的空间导数矩阵。
[0010]优选的，所述步骤S2包括：根据风速随位置变化的关系，计算速度梯度，构建风机叶片表面切应力τ(U)为：
ꢀꢀ
（7）。
[0011]优选的，所述步骤S2包括：将切应力τ(U)和压强p沿着叶片表面和叶片表面单位法向量相点乘，基于面积微元dA进行面积分，求得叶片所受合力矩T
r
(U)为： （8）；其中，ρ为空气的密度，为边界点到风机轮毂中心的矢量位移；
计算得到风机叶片在该来流风速下的第一发电功率P1为： （9）。
[0012]优选的，所述步骤S3包括：涡流发生器源项表示为： （10）； （11）；（12）；其中，为模型常数，通常取值为10，为涡流发生器表面的单位法向向量，为涡流发生器表面的单位切向向量，为和垂直的单位法向向量，（w
x
, w
y
, w
z
）分别为引入涡流发生器后网格风速在x、y、z三个轴向上的速度，V
cell
为当前涡流发生器的单个网格体积，S
VG
为涡流发生器的表面积，L
VG
和H
VG
分别为涡流发生器的长与高，V
s
为引入了涡流发生器修正模型的网格的总体积。
[0013]优选的，所述步骤S3包括：在网格的控制方程中引入涡流发生器源项，以控制区域边界处的来流风速和风机叶片表面空气速度作为控制方程的边界条件，以得到修改正后的控制方程为： （13）；其中， （14）；
ꢀ
（15）。
[0014]将解得的W代入公式（8），得到有涡流发生器时的叶片所受合力矩T
r
(W)，基于公式（9）计算得到第二修正发电功率为P2： （16）。
[0015]优选的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种风机叶片涡流发生器优化方法，其特征在于，所述方法包括步骤：S1、构建风机叶片模型，以环绕该风机叶片模型周围的三维空间区域作为控制区域，将所述控制区域进行网格化；S2、以控制区域边界处的来流风速和风机叶片表面空气速度作为控制方程的边界条件，构建每一个网格的控制方程，计算得到对应网格的风速和压强，以风机叶片模型表面网格的风速和压强计算得到风机叶片表面的压强和切应力，并得到风机叶片在该来流风速下的第一发电功率；S3、设置涡流发生器的形状，并将该涡流发生器设置在风机叶片表面的任一位置，获取该涡流发生器所占据的网格空间位置，在每一个网格的控制方程中引入涡流发生器源项，以得到每一个网格修正后的控制方程；S4、基于每一个网格修正后的控制方程，计算得到对应网格的修正后的风速和压强，以风机叶片模型表面网格的修正后的风速和压强计算得到风机叶片表面修正后的压强和切应力，并得到当所述涡流发生器处于所述位置时风机叶片在所述来流风速下修正的发电功率；S5、依次设置涡流发生器在风机叶片表面的各个不同位置，以及设置涡流发生器的不同形状，重复执行步骤S3
‑
S4，得到不同形状的涡流发生器处于风机叶片的不同位置时所对应的风机叶片的多个第二修正发电功率；S6、将每一个第二修正发电功率分别与第一发电功率进行比较，得到不同形状的涡流发生器处于风机叶片的不同位置时所对应的多个无量纲特征功率，取多个无量纲特征功率中最大者所对应的涡流发生器的位置和形状作为涡流发生器的最佳设计方案。2.如权利要求1所述的风机叶片涡流发生器优化方法，其特征在于，所述步骤S1包括：风机叶片模型具有三维空间坐标体系，以风机轮毂的原点中心为坐标体系的原点，x轴方向垂直于风机叶轮指向风机的下游方向，以重力的方向作为z轴方向，y轴方向与x轴方向和z轴方向垂直。3.如权利要求2所述的风机叶片涡流发生器优化方法，其特征在于，所述步骤S1还包括：所述控制区域是以风机叶片的质量中心为原点，分别向x、y、z轴三个轴方向上扩充6倍叶片弦长距离的三维空间区域。4.如权利要求3所述的风机叶片涡流发生器优化方法，其特征在于，所述步骤S2包括：以控制区域边界处的来流风速作为控制方程的边界条件： （1）；风机叶片表面空气速度作为控制方程的边界条件，风机叶片表面空气速度为：
ꢀꢀ
（2）；其中，为边界点到轮毂中心的矢量位移，η为来流风速下叶片的叶尖速比，R为叶轮半径；对于控制区域中网格的控制方程为：
ꢀ
（3）；其中，
ꢀꢀ
（4）； （5）； （6）；其中，t为时间，（u
x
, u
y
, u
z
）为网格风速在x,y,z三个轴向上的速度，p为网格的压强，ρ为空气的密度，ν
a
为空气的有效动力学粘度，ν
t
为空气的湍流动力学粘度，g
z
为z轴方向的重力加速度，U为网格风速的空间导数矩阵。5.如权利要求4所述的风机叶片涡流发生器优化方法，其特征在于，所述步骤S2包括：根据风速随位置变化的关系，计算速度梯度，构建风机叶片表面切应力τ(U)为： （7）。6.如权利要求5所述的风机叶片涡流发生器优化方法，其特征在于，所述步骤S2包括：将切应力τ(U)和压强p沿着叶片表面和叶片表面单位法向量相点乘，基于面积微元dA进行面积...

【专利技术属性】
技术研发人员：水沛，尹旭烨，
申请(专利权)人：浙江中自庆安新能源技术有限公司，
类型：发明
国别省市：