1. 引言
近年来,随着国民经济的快速发展,整个社会对能源的需求量正在急剧上升,然而对于我国来说,传统化石能源不论是从可持续性以及环境角度,均已经难以满足社会发展的要求,发展新能源已逐步成为能源发展的重要目标。在风电行业,我国的大型风电场大多位于“三北”地区,这给基础建设、电力输送、监控管理以及运行维护都带来相当的难度,对于中东部以及南部等人口密集地区,经济基础好,能源及电力需求量大,风能密度普遍不高,年平均风速大多数处于低风速环境,如何高效利用低品位风能进行发电,得到人们越来越广泛的重视。
在低品位风能的高效利用方面,张臣等人通过对渐缩型聚风装置的研究,探索了渐缩型风筒对于微风汇聚的积极作用 [1];田德等人对渐缩–渐扩型聚风筒进行了相关实验研究,提出一种聚风型风筒的方法,对风能进行浓缩,达到了聚风提速的效果 [2] [3] [4] [5];大屋裕二与张兴联合研制了紧凑型聚风发电装置,利用出口处的负压抽吸作用对聚风筒风能进行提速,表明渐扩型聚风筒有着优秀的聚风能力 [6] [7] 。在对低品位风能聚风型风筒的研究过程中,本研究团队注意到风筒的尺寸对于风机的效率有着较大影响,且不同的模型尺寸对模型聚风能力的影响程度也各不相同。因此,本文基于CFD数值模拟仿真方法 [8] ,采用参数化设计与实验研究相结合的方式,针对聚风筒尺寸效应对于低品位聚风效果的影响规律以及变化趋势进行了重点研究。
2. 物理模型
本文研究的对象为聚风筒模型,如图1(a)和图1(b)所示,该模型主要由三段构成,前端为收缩管型进口,后端为扩张管型出口,中间为直筒状喉口。前端收缩管长度为lin,进口张角为θin;后端扩张管长度为lout,出口张角为θout;中间喉口段长度lmid,喉口段直径为dm。
采用控制变量的方法,主要研究尺寸变化对聚风筒聚风效率的影响规律。文中选择喉口直径为500 mm的模型作为参照模型,并设置喉口直径分别为800 mm,1200 mm,1600 mm模型作为放大实验对照组,设置喉口直径为300 mm,200 mm,100 mm,50 mm的模型作为缩小实验对照组,各组对照模型以其喉口直径大小为依据,等比例进行缩放。
(a)
(b)
Figure 1. Three-dimensional model and structural parameters. (a) Three-dimensional model; (b) Structural parameters
图1. 三维模型及结构参数。(a) 三维模型;(b) 结构参数
3. 聚风筒数值仿真模型及实验验证
3.1. 数值仿真模型
本文采用全流场数值模拟的方法,对聚风筒的内流场运动规律及外流场运动规律进行研究。为提高计算精度同时降低网格总数,文中采用远场模型与近场模型相配合形成混合结构的数值计算模型,如图2(a)所示。远场模型为大长方体流体域,并且宽度和高度在尺寸上均高于风筒最大直径的10倍,而长度方面则高于风筒总长的20倍 [7] 。
(a)
(b)
Figure 2. Full flow field model and mesh generation of the air duct. (a) Numerical simulation model; (b) Mesh generation
图2. 聚风筒全流场模型及模型网格剖面图。(a) 数值计算模型;(b) 网格划分
在网格划分方法上,近场模型采用非结构化网格的方式进行划分,并在近壁面处生成3层结构化三棱柱网格 [7] ,如图2(b)所示。模型的远场采用结构化网格进行划分,控制各个位置的网格疏密度,使得内部网格密度较大,越靠近边缘网格密度越小。采用混合网格的方式将结构化远场模型与非结构化近场模型进行连接,获得整体网格数663万。
3.2. 实验测试系统及数值结果验证
将上述对照实验中的参照模型(喉口直径500 mm)按照1:1大小做出实物模型,并依据实验要求放置于实验风场中,具体实物模型与测试方法如下图3所示。对风洞出口进行风场区域量化,选择该风洞出口处面平均风速最接近3 m/s的位置,进行实验测试工作,经测算实验风场的平均风速为2.45 m/s。
![](//html.hanspub.org/file/4-2810105x14_hanspub.png)
Figure 3. Test system of the air duct
图3. 聚风筒实验测试系统
按照上述数值模拟方法进行仿真计算,为表征流体沿流线的速度分布规律,数值仿真的采点方法为沿着垂直于轴线方向等间距采点,如图4(a)所示。针对参照模型(喉口直径500 mm)各测点位置按照来流风速为2.45 m/s进行测试,分别进行实验测试和数值仿真,得到喉口处的风速值,如图4(b)所示,结果显示,各个测点处数值仿真结果与实验测试结果有较好的吻合度,表明本文构建的聚风筒数值模型的有效性。另外,对仿真模型喉口处的风速值取加权平均,可知在来流风速为2.45 m/s时,模型的喉口平均风速达到3.26 m/s,其喉口处平均风速相对于来流风速值提升了33%,这说明本文提出的聚风筒模型具有明显的低风速风能汇聚效果。
(a)
(b)
Figure 4. Location of different measuring points and wind speed value curves. (a) Location of measuring points; (b) Wind speed value curves
图4. 不同测点位置及风速分布曲线。(a) 测点位置;(b) 风速值曲线
4. 聚风筒尺寸效应计算结果及分析
图5和图6分别为不同尺寸下的聚风筒速度云图和压力云图。分析可得,在模型喉口直径不低于500 mm时缩小或放大参照模型的尺寸,喉口速度均保持在较高水平,而且随着尺寸由大到小,喉口风速也是呈现出正相关线性变化趋势,但是当喉口直径低于500 mm时,喉口风速会急剧下降,喉口风速大小依然呈现与模型尺寸大小成正相关的变化趋势。分析原因为:空气具有一定的粘性,根据粘性力方程τ = μ∙du/dz,在模型尺寸较小时(即z值较小时),模型内流场的速度梯度du/dz较大,粘性应力以及摩擦损失较高,进而阻止了风筒内气体流速的提高,而当模型尺寸较大时,这种阻碍作用相对减弱。而且模型较小时,边界层厚度相对于模型尺寸来说就相对较大,边界层效应带来的流动损失也会更大,因而在模型喉口直径低于500 mm后,喉口风速的下降速度也会加快。当模型尺寸放大到参照组模型的2.4倍(喉口直径1200 mm)及以上时,喉口的风速并不能带来较为明显的提升,这也与边界层效应和空气粘性效应相对于模型整体尺寸来说达到饱和有关。
![](//html.hanspub.org/file/4-2810105x17_hanspub.png)
Figure 5. Velocity contours of air duct at different size
图5. 不同尺度下的风筒速度云图
![](//html.hanspub.org/file/4-2810105x18_hanspub.png)
Figure 6. Pressure clouds of air duct at different size
图6. 不同尺度下的风筒压力云图
![](//html.hanspub.org/file/4-2810105x19_hanspub.png)
Figure 7. Relationship between the wind speed ratios and air duct size
图7. 聚风筒喉口风速比随风筒尺寸的变化曲线
根据边界层理论可知,流体贴近固体壁面的速度为零,并沿着径向逐渐提升至和喉口最大风速值相接近,由于不同尺寸的模型喉口风速值近似为3 m/s,变化幅度小,而模型的喉口直径却有较大变化,根据粘性应力理论,流体在模型喉口处速度沿着径向变化将有很大差距,粘性应力的大小也相差甚远。本文为表征粘性力针对不同规格聚风筒的尺寸效应,将模型喉口平均风速与喉口直径的比值作为速度梯度du/dz的近似代替,绘制成如下曲线,考察该特征曲线同时对比图7可以看出,随着模型尺寸的减小,喉口风速急剧下降的同时,喉口处的速度梯度也急剧上升,由此侧面说明了上述尺寸效应分析的正确性。
针对不同尺寸的数值模型,对其喉口截面风速值取加权平均,绘制成如图8所示特征曲线,当喉口直径低于500 mm时,模型尺寸对模型气动特性的影响明显,且尺寸越小,喉口风速的下降趋势越大,而高于此值时,逐渐增大模型尺寸,模型喉口风速的增加趋势将逐渐趋缓,此结果与上述云图分析结果吻合。
![](//html.hanspub.org/file/4-2810105x20_hanspub.png)
Figure 8. Curve of average wind speeds at different sizes of the air duct
图8. 不同尺度下的喉口平均风速变化曲线
5. 结论
本文研究了聚风筒在低品位风能汇聚下的尺寸效应,考察了不同尺寸下聚风结构内外的流动分布和特征,同时也量化分析了尺寸效应对于汇聚特性的影响规律,得到结论如下:
1) 针对本文研究的聚风筒结构,喉口处的聚风能力与风筒尺寸成正相关关系,变化曲线符合指数型变化规律,当喉口直径低于500 mm时,风筒的聚风能力会随着尺寸的减小急剧下降,而风筒直径高于500 mm时,风筒聚风能力的上升幅度将减弱,并随尺寸的增大而逐渐趋于平稳。
2) 聚风筒的尺寸效应与聚风筒内粘性力摩擦损失的大小密切相关,风筒尺寸较小时,筒内流体速度梯度大,摩擦损失高,聚风效率低;风筒尺寸较大时,筒内流体速度梯度小,摩擦损失低,聚风效率高。
基金项目
湖北省技术创新专项重大项目(编号:2017AAA035);低品位能源利用技术及系统教育部重点实验室开放基金(No. LLEUTS-201905)。
NOTES
*通讯作者。