论文摘要
随着能源问题和环境问题的日益突出,作为一种清洁的可再生能源,风能的开发利用已得到世界各国的高度重视,我国风力发电事业近十年也得到了快速发展。在风力机研制和生产过程中还存在一些问题,其中一个关键问题是适应于不同风场的风力机的气动特性。目前广泛应用的气动性能研究方法是基于传统的叶素理论以及二维翼型风洞实验数据,而实际的风力涡轮是三维及旋转的,采用该方法低估了实际风轮的动力产生。随着现代计算机技术的发展和三维湍流模拟技术的提高,CFD方法在风力机气动性能研究中的作用越来越明显。本文选取在传统水平轴风力机中使用较多的NACA63-215翼型建立了绕翼型流动的二维和三维物理、数学模型,采用NUMECA公司的商业软件FINE/TURBO软件包进行计算,湍流粘度分别采用基于RANS的Baldwin-Lomax,Spalart-Allmaras和k-ε三种湍流模型来处理,采用时间相关法求解,Jameson的中心差分格式,四阶Runge-Kutta法进行时间推进。为了加速收敛,使用了多重网格技术和隐式残差平均技术,结果以等值线、云图、速度矢量图、流线、二维曲线等方式较为形象直观地描述了流场的运动规律。数值模拟结果表明,基于RANS的Spalart-Allmaras湍流模型计算出的二维翼型气动性能与实验测试值更接近。而后基于叶片的旋转坐标系建立了叶片三维物理和数学模型,湍流模型选取在翼型气动性能计算中所得结果与实验值较为接近的Spalart-Allmaras湍流模型,结果显示结合定量的压力分布曲线和定性的壁面极限流线清晰地反映不同工况的流场。通过对比二维非旋转与三维旋转模型的计算结果发现,三维旋转风轮截面的静压分布与二维非旋转条件下翼型表面静压分布有较大差别,从而验证了静态失速问题的存在并分析出现此现象的机理。本论文研究表明,在旋转的作用下,风力机叶片的气动特性大大改变,所以研究三维旋转效应,进一步探索认清其规律有一定的工程价值;另一方面,CFD方法应用于求解风力机流场特性是可行的,可作为实验研究的一种补充和预测手段,随着计算机的发展以及湍流模型的不断完善,计算流体力学的应用必将成为风力机气动性能计算的发展趋势。