风力发电机械中机组叶片小翼三维场数值模拟研究

2023-10-24李擎雯田克琴唐钰铭

机械与电子 2023年10期

程乐,李擎雯,田克琴,唐钰铭

(国家电投集团广西兴安风电有限公司,广西桂林 541000)

0 引言

近年来,全球资源愈发紧张,工业化发展加速了资源的干涸,造成了大部分国家的资源短缺,也破坏了原有的生态平衡。由于某些能源不可再生,使资源再利用成为了现今最需要重视的问题之一,因此,人们将目光放在了可循环利用的风能上,风能作为一种可再生能源,不仅可以就地利用,还不受地域环境影响,可以代替大多数能源。研究发现,风能不仅能够代替机械能,还能将动能转换成电能,满足人类对电力的需求,而风力发电不仅取之不竭,其发电后的产物还不对环境造成危害,风能的持续发展可能会开启一个新的时代。

风力发电虽然环保卫生,但受到多方面的限制,其取决于天气与光照的强度,如果用之不慎,还会造成一定的事故。由于机械组叶片旋转会受到一定的压力,产生的能量与损耗不同,为使所有参数的变化在可控范围内,达到最佳效果,需要对风力发电机械中机组叶片小翼三维场数值模拟,以进一步保证风力发电机组的稳定性,解决能源短缺的问题。

1 风机效率上限计算

以风能作为发电的源头是环保且合理的,也可能是未来主要的发电模式,根据气候的改变推测出风力的大小和风的方向,原理是利用季节的不同来确定风的走向。风力发电机械中机组叶片小翼结构如图1所示。

图1 机组叶片小翼结构

在合适的位置将机械的叶片调整至阻力较小的角度来进行旋转,使风能转化为电能,传输到线路中,整个转变过程和力学与能量守恒规律相似,因此以能量守恒作为标准并模拟其能量转化[1]。假设在三维场景中的叶片体积不变,当光照正常,风力较大时,空中的气压与阻力就会占据一定空间,使场景变得拥挤,而这些阻力并不是一成不变的,而是随着叶片的旋转实现气体的交换,将气流维持在一个范围内,并计算体积大小,将高速旋转的叶片看作是一个圆形气流交换通道[2],那么则有以下规则:

a.气体交换产生的阻力不具备干扰性,进出的体积不变,但存在质量的改变,属于稳定的流体。

b.叶片旋转异常快速,不存在摩擦力,产生的空气阻力也可以忽略不计[3]。

c.叶片的数量越多,交换气流的速度越快,两者呈正比。

d.空气流形成的圆形通道每个角度受到的力都是相同的,可以相互代替。

e.控制速度的主要是叶片的上部,尾部几乎没有作用。

f.场景中的出入口面积都是相同的,其中的气压也是随着入口的变化而变化的。

按照上述规则来预测整个系统中动能与风能的能量转换,假设在三维场景中气体交换时不存在空气阻力,忽略自身能量的影响,保证质量守恒,也就是说空间中气流是不断进出的,叶片的旋转速度相同,当出现压力时,反方向会同时出现一个阻力,与压力相抗衡,力的分析如图2所示。

图2 贝兹理论分析示意图

根据图2中的提示,气流量交换[4]的过程用公式表示为

A=mV=(ρSV)V

(1)

A为气流量计算值;m为气体的质量;V为交换的速度;S为气体占据空间的体积;ρ为气体本身的密度,则阻力与压力的守恒公式为

P=m1V1-m2V2=(ρS1V1)V1-(ρS2V2)V2

(2)

P为压力。在能量保持平衡时,一定时间内气体的交互率与速率是相同的[5],用公式表示为

(3)

当出现在机械叶片上时,气体流速就自然小于叶片旋转的速率,不仅能够保持能量守恒[6],还不会有能量损失,则测量出的圆形通道面积可表示为

(4)

假设叶片前后部分受到的压力与空气阻力相等,即P1=P2,那么旋转时速度也相等,即V1=V2,则产生的摩擦力可以表示为

F=S(P1-P2)

(5)

比较2个公式中产生的能量[7],则有

(6)

V4为叶片上下部分的平均速率。

基于上述过程对动能与风能的能量转换,计算出风机效率上限,为叶片小翼三维场数值模拟提供基础。

2 叶片三维建模

将计算所得到的各参数平移至三维场景中,将发电机与其他设备的位置用坐标来描述[8],并记为[x′,y′,z′,1]=[x,y,z,1]·T,T为周期,那么在1个周期内空间中所有坐标的变换矩阵H为

(7)

矩阵中的所有元素都代表叶片的坐标点,按照其排列顺序确定变换前后的位置[9],假设以叶片前段作为原点,坐标为(x1,y2),转换后的坐标就为(x2,y1),原点坐标为(X,Y),那么转换过程的表达式为

(x2,y2)=(x1,y1)×N

(8)

N为倍数。叶片旋转后导致空间中所有点变换了各自的位置[10],利用不变的参数来描述出实际的坐标,其公式为

(9)

α为叶片之间的空间角度。

在上述建模的基础上,机械按照风力的等级将风能分别使用在叶片不同的发电机上,风力机源源不断地吸收风力,一方面转换成动能促使叶片持续旋转,另一方面用来抵抗空气阻力[11],一般利用特征系数来描述风力机吸收风能的多少,用公式表示为

(10)

B为机械吸收的全部风能;σ为特征系数。为了更加精确地表示风力等级对风轮组的影响,保证在转换过程中能量守恒,将叶片与风速之间的比例用公式表示为

(11)

λ为速度比;n为风轮的数量;R为风轮旋转时的半径;r为单位时间内风轮的转数。

3.绿色低碳成为新的经济发展方向。深圳的空气质量已经是国内1000万级人口城市的最优水平，是国家发展循环经济的示范城市，在经济不断快速发展的同时，深圳不断增加环境保护方面的投资，从1996年的不到20亿元增加到如今的298.65亿元，在国民生产总值中所占比重保持在2.8%左右，位于全国前列，体现了深圳对环境保护和绿色低碳产业的重视。深圳正处在改善环境质量的重要时期，绿色低碳成为新的经济发展方向。

当打开发电机械时会立即产生扭矩[12],这是由于叶片在发动过程中加速了气流的交换,导致出现摩擦阻力,这种情况直接影响了特征系数与转速比,因此,在研究发电设备时要着重考虑扭矩的大小,保证在阻力最小的前提下,得到最大的电功率[13],那么系数的公式就变为

(12)

L为扭矩的长度。

经过上述计算能够实现风能转换电能的最大化,减少耗能。

3 压力修正

利用风能发电已经逐渐成为我国建立发电厂的备选技术之一,通常要保证机械中叶片保持匀速运动,还要考虑到光照的影响。使用SIMPLE算法控制压力的变化范围,赋予一个可以迭代的数值,其原理是将压力分解[14],通过方向不同的力建立离散方程,确定各方向受到的压力与风速,由于压力看不见且不能测量,那么得到的风速就不一定遵循方程规律,所以对给予的压力逐一修正,修正的目的是对流场控制,具体如图3所示。