网格单元类型对立轴风力机数值模拟影响的研究.pdfVIP

网格单元类型对立轴风力机数值模拟影响的研究 梅毅，曲建俊，许明伟 (哈尔滨工业大学机电工程学院，哈尔滨150001；) 恭蔓：本文基于非定常Navier-Stokes方程和滑移网格技术，分别采用三角形和四边形两种网格单元类型，模拟了10m／s 风速下垂直轴风力机处于最优尖速比时的风轮流场，通过比较计算最大功率系数值和风洞实验数据，分析风轮旋转过 程中两种模拟方法得到压力系数分布，研究了网格单元类型对垂直轴风力机气动性能数值计算精度的影响。结果表明， 两种模拟计算值与实验值之间都存在差别，四边形网格计算最大功率系数的精度比三角形网格高8．26％。。与四边形网 格相比，采用三角形网格模拟时叶片气流分离的位置更靠近前缘，负压面压力峰值更小。在垂直轴风力机流场数值模 拟中，采用四边形网格可以得到更为精确的结果。 关■阅：垂直轴风力机：数值模拟：滑移网格：网格单元：精度 O引言 风能是一种具有较好开发利用前景的可再生资源。水平轴风力机是目前技术比较成熟，市场占有 率较高的风能利用设备，但水平轴风力机并不适应所有的场合。直叶片垂直轴风力机由于具有无需对 风装置，噪声小以及加工制造相对简单等优势，更适合作为城市高楼楼顶、边远山村、边防哨所等地 域的风能利用设备，引起了研究人员对垂直轴风力机相关技术问题的广泛关注【1．3】。 随着计算机技术和数值算法迅速发展，基于N—S方程的CFD方法可以清晰、准确地揭示出非定 常流场的流动机理，并获得相应的气动性能数据，逐渐成为研究垂直轴风力机的重要工具。杨从新[4】 则采用四边形网格应用于CFD计算。关于哪种网格单元类型应用于垂直轴风力机时计算精度更高， 现有文献中尚未有相关研究的报道。 垂直轴风力机流场复杂，攻角变化范围大，叶片动态失速和流动分离频繁，对数值模拟提出了较 高要求。网格是CFD模型的几何表达形式，网格单元类型对计算精度有重要影响，需要进行有针对 性的研究。本文通过与风洞实验数据对比研究了网格单元类型对CFD计算精度的影响，研究结果可 为垂直轴风力机数值模拟技术在工程中的应用提供参考。 1研究模型 本文以文献[8】中的风力机为研究模型，将其风洞实验数据与CFD计算结果进行对比，以验证何 种网格单元类型更适合垂直轴风力机模拟。风力机回转直径2．5m，风轮转轴直径O．1m，叶片高3m， 对比研究。 2数值方法 2．1控制方程与湍流模型 直叶片垂直轴风力机叶片截面相同，可将三维流场简化为二维流场，以减小计算量。在以恒角速 度绕风轮主轴旋转的相对直角坐标系下，时间平均控制方程组为 望+旦邋：o at 眈 (1) 良，I—I ， J 旦逊at+塑a≯x=一鱼Ox,+毒卜(考0x+丝Ox,]]+s 式中，D是流体密度，P是压强，lae是有效粘性系数，Si是广义源项，xi和yi是各坐标分量，ui和 uj是平均相对速度分量。 k．u模型在壁面处采用 经过时均后，方程中因出现了雷诺应力等脉动关联项未封闭。两方程SST 标准k．u模型，在边界层边缘和自由剪切层采用标准k．￡模型。由于综合了k-w模型在近壁区和k-E k-w模型对逆压梯度较大的流场和分离流的预测有更高的精度，因此在 模型在远场计算的优点，SST 本文的研究中选用SSTk-w湍流模型来封闭时均方程组。 2．2几何建模与网格划分 采用Gambit软件进行几何建模和网格划分。建立图1所示的风轮叶片和转轴二维模型以及计算 为150，风轮旋转区域直径为2①。分别用四边形网格和三角形网格划分计算域(如图2)，网格总数 网格处理旋转区域和非旋转区域之间的耦合问题。 图1计算域示意图 (a)四边形网格 (b)三角形网格 图2计算域网格划分示意图 2．3边界条件与求解设定 利用商用软件fluent求解时均方程组，进口边界AC设置为速度进口边界条件，来流速度为10m／