Savonius风力机横纵向重叠比对功率因数影响.docVIP

Savonius风力机横纵向重叠比对功率因数影响 　　摘要：为研究叶片重叠比对Savonius风力机功率性能的影响，提出横向重叠比与纵向重叠比2个结构参数。建立二维有限元分析模型，应用FLUENT进行非定常数值模拟，基于RNG k-ε湍流模型，采用滑移网格模拟风力机转动，得到不同横向重叠比和纵向重叠比时的平均功率因数。结果表明：随着横向重叠比和纵向重叠比的增加，风力机的平均功率因数均先增大后减小；横向重叠比的最优值为0.10，纵向重叠比的最优值为-0.05，对应的平均功率因数为0.251 5。 　　关键词：Savonius风力机；功率因数；横向重?B比；纵向重叠比；滑移网格；数值模拟 　　中图分类号：TK83 文献标志码：B 　　0引言 　　Savonius风力机属于垂直轴阻力型风力机，该风力机主要由2～3个半圆弧形叶片组成，其工作原理是依靠叶片凹面和凸面的气流阻力差驱动风力机旋转。Savonius风力机具有很多优点：叶片设计简单；不受风向限制；启动力矩大；噪音小；发电机可置于风轮下或地面，因而安装和维护的费用都较低。为提高风力机的风能转化效率，有学者从实验和数值计算2方面研究叶片展弦比、叶片重叠比、叶片数以及风力机分级数等对风力机功率性能的影响。SHELDAHL等指出增加叶片展弦比可提高风力机功率；SHELDAHL等和AKWA等均指出2叶片的最佳重叠比为0.10～0.15；SAHA等指出2叶片风力机功率最高，两级风力机的功率比单级和三级的高。部分学者采用改变风力机结构的方式提高风力机的功率。赵振宙等、KAMOJI等和GUPTA等对螺旋形的Savonius风力机进行实验和数值研究。GOLECHA等和ALTAN等提出利用前置导流板方案，该方案通过导流板的遮挡效应减小叶片回程阻力，同时增加叶片进程方向的流速，进而提高风力机平均功率。 　　本文重点研究叶片重叠比对风力机功率的影响。叶片重叠比包括横向重叠比和纵向重叠比，以往对于叶片重叠比的研究都是针对叶片横向重叠比，对纵向重叠比的研究甚少。本文综合研究叶片横向重叠比和纵向重叠比，以完善该参数对风力机性能影响的规律。 　　2数值模拟方法 　　2.1几何模型 　　以往对于叶片重叠比的研究都只针对叶片横向重叠比，没有对叶片纵向重叠比进行研究分析，因此，本文在传统重叠比定义的基础上，提出横向重叠比和纵向重叠比2种定义，横向重叠比为X/d，纵向重叠比为Y/d，其中：X为横向间隙，Y为纵向间隙，d为叶片直径。 　　Savonius风力机的二维结构示意见图1，其中：θ为叶片方位角（图示位置为θ=90°），ω为风力机旋转角速度，d为叶片直径，D为风力机叶轮直径。本文采用的原模型为美国Sandia实验室采用的风力机模型（D=1 m，d=0.5 m，X/d=0，Y/d=0）。 　　2.2计算域和边界条件 　　为使数值模拟结果不受计算域划分的影响，选择计算域的尺寸为18D×12D。为实现风力机的转动，采用滑移网格，并将所设置的计算域分为静止域与旋转域。静止域为外流场空气域，旋转域为模拟风力机转动的子域，静止域与旋转域之间的数值模拟数据通过交界面传递。 　　设置流体材料为常温空气。计算域边界条件设置为：入口边界为恒定风速入口，大小为美国Sandia实验室采用的风速7 m/s；出口边界条件为压力出口；上下边为静止壁面；叶片为旋转壁面；静止域与转动域的交界面为滑移边界。计算域和边界条件见图2。 　　2.3 　　网格划分 　　选用ANSYS Workbench 15.0自带的网格划分工具进行计算域网格划分。由于在二维计算中四边形网格具有占用内存少的优点，且对边界层计算非常有利，所以整个计算域采用四边形网格划分，见图3a）。为更加准确地模拟Savorous风力机周围空气流动情况，风力机叶片壁面边界划分为边界层网格，并对边界层和旋转域网格进行加密，见图3b）。 　　2.4计算设置 　　风力机的非定常数值模拟计算基于压力求解器，采用SIMPLEC算法，选择RNG k-ε湍流模型。为使数值模拟具有较好的收敛性和计算精度，动量方程、紊动能后方程和耗散率ε方程均采用二阶迎风格式。 　　根据风力机的旋转速度设定不同的计算时间步长，保证每个计算时间步风力机转动1°，每个计算时问步的最大迭代步数为100；流场的连续性、速度分量、湍动能k和耗散率ε的残差收敛标准均设为1×10-5，以保证计算的准确性。 　　3仿真结果验证 　　为验证数值模拟方法的准确性，计算不同尖速比下的平均力矩因数，并与美国Sandia实验室试验得到的平均力矩因数进行对比，结果见图4。由此可以看到，数值模拟计算结果与美国Sandia实验室结果基本吻合，从而验证数值模拟方法的准确性。 　　4仿真结果分析 　　4.1横向