风力机叶片翼型的气动性能在很大程度上影响着风力机的发电效率[1]。翼型在不同攻角下的气动性能有较大差异,尤其是当翼型处于失速状态时,风力机会产生较大的振动,导致其安全性能降低。因此,准确地预测翼型的气动性能是十分有必要的[2-3]。
由于不同湍流模型对翼型气动性能参数的计算结果会有较大不同,因此选择一种合适的模型来进行数值计算尤为关键[4],为此,许多学者都进行了相关的研究。丁勤卫等[5] 采用S-A、T-SST这2 种湍流模型对翼型的气动性能进行了数值计算,结果表明,当翼型出现失速时,T-SST 模型可以更好地预测其气动性能参数。彭续云等[6] 采用3 种湍流模型对NACA63-421 翼型进行了气动性能数值计算,结果发现,k-KL-ω 模型对阻力系数的预测效果最好;在大攻角范围内,SST k-ω 模型计算得到的升力系数值与实验值最为吻合。
目前对于湍流模型影响翼型气动性能参数计算精度的研究已较为成熟,而将S-A、SST k-ω和RNG k-ε 这3 种湍流模型进行综合比较分析的文献较少,甚至大多数文献都未关注到离散格式对翼型绕流数值计算精度的影响。因此,本文选取S-A、SST k-ω 和RNG k-ε 这3 种湍流模型对风力机叶片NACA0018 翼型进行气动性能参数计算,得到一定攻角范围内的气动性能参数情况,然后将其与文献[7] 给出的风洞实验值进行对比,并讨论了一阶迎风、二阶迎风、QUICK 这3 种离散格式对计算精度的影响。
1 翼型的气动性能参数当来流流过翼型表面时,由于上、下翼面的来流速度不同,导致上翼面的速度大、压力小,下翼面的速度小、压力大。上、下翼面的压差即为翼型所受到的升力,升力沿圆周切线方向的分力推动风轮转动。升力系数Cl、阻力系数Cd、升阻比Cl/Cd 均是表示翼型气动性能优劣的参数,其中,升阻比越高,说明翼型的综合气动性能越好。
升力系数Cl 的表达式为:
式中,Fl 为翼型所受升力;ρ 为空气密度;u为来流风速;c为翼型弦长。阻力系数Cd 的表达式为: