导流型垂直轴风力机气动特性及其计算方法的数值研究.doc

中国工程热物理学会 流体机械 学术会议论文 编号： 中国工程热物理学会 流体机械 学术会议论文 编号：097004 导流型垂直轴风力机气动特性及其 计算方法的数值研究 王鑫 童正明 王企鲲 (上海理工大学 动力工程学院，上海 200093) 联系电话 E-mail:wangxin_usst@139.com 摘要：本文分析了传统垂直轴风力机效率低的原因，在此基础上提出了一种带有导叶的导流型垂直轴风力机，并数值研究了该种风力机的气动性能，对导流型风力机内部旋转区域的计算方法进行了分析。导叶不仅可以有效地降低因来流对动叶轮吸力面的直接冲击而造成的阻力扭矩，而且还有助于改善来流对动叶轮压力面的有效冲击，这些均使该风力机动叶轮的旋转扭矩得到显著增加。通过比较内部区域的计算结果，给出不同计算方法在工程应用中的所适用的范围与场合。 关键词：垂直轴风力机、导流叶片、数值研究、多参考坐标系、滑移网格  0 引言 风能作为一种清洁可持续利用能源，在现今社会得到了广泛的重视，风力发电作为一种储量丰富、无污染的清洁能源倍受关注，已成为构成电力的一个重要部分[1]，而风力机本身的性能对风能利用起到了至关重要的作用。 （a）导流型垂直轴风轮 图1. 垂直轴风轮示意图 （b）无导叶垂直轴风轮 风力机按结构形式大体上分为水平轴和垂直轴两种。水平轴风力机发展较完善，由于其具有较高的效率而作为当今主流的风力机类型得到了广泛的应用。然而由于水平轴风力机安装难度高、投资成本大且不便移动等缺点，它的使用目前往往主要限于并网发电中，而对于一些小型离网型发电系统中，目前使用最为广泛的仍是垂直轴风力机。 与水平轴风力机相对，垂直轴风力机具有投资成本低、安装方便、便于移动、运行噪音小等特点，其缺点是效率较低。这些在阻力型垂直轴风力机上表现地尤为突出[2]。 由于阻力型风力机迎风面上一侧风阻的存在，大大降低了垂直轴叶轮的有效受风性能。因此，只要降低该侧上风阻，即有望能提高垂直轴风力机的效率。 本课题组提出一种导流型垂直轴风力机，图1（a）是其结构的平面示意图。它主要是在传统的垂直轴风轮如（图1（b）所示）的外围添加一个由数个导流叶片沿圆周方向均匀排列而构成的导风轮。导风轮是静止部件，主要起导流作用。在其导流作用下，使来流完全导向左侧风叶的压力面，同时对右侧风叶而言，导流叶片能起到挡流作用，使来流不直接冲击风叶的吸力面，这样就有望通过简单的二元叶轮而实现对整个垂直轴风力机的作功效率的提高[3,4]。 采用CFD（计算流体动力学，为Computational Fluid Dynamics的缩写）技术对导流型垂直轴风力机的内部流场进行数值计算是获得其气动性能的一种快捷而经济的方法。尽管CFD计算结果尚须利用实验进行校核，但在设计的初期预先采用CFD技术进行计算显然可以使随后的实验更具有针对性，这已在众多的叶轮机械的设计过程中已被证实[5]。 目前，对于旋转式叶轮机械内部流场的CFD计算，要处理动、静两部分区域交界面上计算信息的传递，这主要有两种计算方法：多参考坐标系法（Multiple Reference Frame, 以下简称“MRF法”)和滑移网格法（Sliding Mesh, 以下简称“SM法”）。其中，MRF法作为一种准定常的近似求解方法，以简便适用的特点在工程中被广泛应用。相比较而言,滑移网格法是一种非定常方法，通过旋转区域与静止区域的相对滑动，在每个时间步内不断产生新的网格，从而给出随时间变化的数值解[6,7]。 由于整个导流型垂直轴风力机绕流流场的数值模拟是典型的对多尺度、大分离、复杂流场的数值模拟。对于此种流场用何种计算方法较为适宜目前尚未可知。基于此，本文对同一种导流型垂直轴风力机的绕流场，分别运用MRF和SM两种方法对其进行模拟计算，从数值计算角度对两种方法的实用性和准确性进行对比分析。 1 计算模型与方法 为降低计算时间，本文选取的导流型垂直轴风力机的物理模型是一种二维叶轮，其具体结构如图1（a）所示。计算区域是以导流型垂直轴风力机为中心，上延其导流叶轮外径的4倍、下延12倍，左右方向延伸5倍的二维矩形范围。 计算网格采用结构与非结构