水平轴风力机叶片附近区域流场的数值分析-numerical analysis of flow field near horizontal axis wind turbine blades.docx

摘要风轮附近区域流场的速度分布和压力分布一直是国内外学者关注的热点。了解该区域流动的分布规律对掌握水平轴风力机功率的获取和提高，有着极为重要的意义。本文利用CFD软件Fluent，按照由浅入深、循序渐进的方式对风轮附近流场进行了数值研究，以等值线图、云图、流动矢量、流线、迹线、二维曲线图等方式较为形象直观地描述了流场运动规律。模拟过程中，将主要控制参数——升阻力系数以及部分流场视图，分别与相关文献和PIV实验进行了对照，误差在工程允许范围内。文章主要分两部分进行。第一部分，奠定研究基础。首先简要介绍第二代100W水平轴风力机的优化设计，设计基于叶素动量理论，利用Matlab编程的方式完成。其次，对设计风力机选用翼型NACA4415的绕流问题进行了数值研究和湍流模型选用的对比研究。第二部分，主体研究部分。先分别研究了二维、三维单叶片准动态绕流时计算区域的流场特性，尤其是尾流区的流场分布情况。最后又进一步探索性研究了三维旋转流场的特性，分别给出了风洞模型比和叶尖小翼对风力机尾流区的影响。数值模拟过程中，数学模型基于那维斯托克斯方程，采用基于压力求解器和离散求解器，稳态隐式方程解法，维数与湍流模型的选取由实际计算过程决定，压力速度耦合采用SIMPLE（求解压力耦合方程组的半隐式方法）算法，离散格式为二阶迎风格式；物理模型造型以及网格离散由前处理软件Gambit完成。为了尽可能减少总体网格数以及改善网格质量，采用了“自下而上”的方式，即“点－线－面－体”的步骤进行网格划分。论文研究证明：一方面，数值研究应用于求解风力机流场特性可行，可作为实验研究的一种补充和预测手段。另一方面，湍流模型、翼型攻角、风洞风轮模型比以及小翼对风力机叶片附近的流场都有很大的影响，进一步探索认清其规律有一定的意义。关键词：水平轴风力机；数值模拟；尾流；升力系数；阻力系数AbstractThevelocitydistributionandpressuredistributionintheflowfieldlocatedatthenearbyregionoftherotorhasbeenafocusofresearchwhichattractsbothdomesticexpertsandforeignscholars.Abetterunderstandingoftheflowinthisareahelpsmastertheacquirementrulesandimprovementmethodsofpoweroutput.ThenumericalinvestigationontheflowfieldaroundtherotoriscarriedonstepbystepbyvirtueofFluent—CommercialCFDsoftwareinthispaper.Theflowisintuitionsticlydescribedbythenumericalresults,whicharepresentedintheformofcontour/colormaps/vector/line/rake/XYcurves.Intheprocessofsimulation,themainparameters-lift/dragcoefficientsandpartialflowvisualizationarecomparedwiththerelatedliteratureandPIVresults.Finally,wefindtheerrorisacceptableintheengineeringaspect.Ithastwomajorpartsinthisarticle.Inthefirstpart,theresearchbasisisset.Atthebeginning,theoptimaldesignofthesecondarygeneration100Watthorizontalaxiswindturbineisintroduced.Thedesignbasesonbladeelementmomentum(BEM)method.Itiscarriedoutbyself-madeprogramwhichbasesonMatlab.Then,theseparatedflowofchosenairfoilNACA4415isresearched.What’more,fivekindsofturbulencemodelsarecompared.Inthesecondpart,themainbodyresearchofthearticleispresented.Westudytheflowcharacteristicoftheseparatedflowaround2Dand3