1、I毕业设计论文院系动力工程系专业班级热能与动力工程学生姓名XXX指导教师XXX题目风机翼型边界层分离的二维数值模拟研究二XX年X月热能与动力工程毕业设计I风机翼型边界层分离的二维数值模拟研究摘要当风机工作时，气体流道的几何形状改变会使流体运动速度的大小和方向发生改变，从而产生流动分离。流动分离产生的冲击会造成流动损失。流体运动速度的大小和方向的改变，也会使得气体在进入叶片入口和从叶轮出来进入压出室时，流动角不等于叶片的安装角，从而产生冲击损失，影响风机的效率和性能。由于气体进入叶片入口时存在着冲击速度，使气体在风机叶片的吸力面上形成旋涡，造成边界层分离现象而会导致能量损失。针对G473风机翼型
2、，利用商业软件FLUENT的前期处理工具GAMBIT建立二维不可压缩湍流模型，再利用FLUENT对翼型在36到8的空气来流攻角下的气动特性进行了相应的数值模拟计算，然后，对不同攻角下模拟所得到的速度矢量图进行比较分析，得出风机翼型边界层分离和攻角的关系。关键词风机翼型；边界层；数值模拟；攻角热能与动力工程毕业设计IITHE2DNUMERICALSIMULATIONOFTHEBOUNDARYLAYERSEPARATIONONAWINDTURBINEAIRFOILABSTRACTWHENTHEFANWORKING,THEGASFLOWCHANNELGEOMETRYWILLCHANGE，WHICHM
3、AKESTHEFLUIDVELOCITYMAGNITUDEANDDIRECTIONCHANGERESULTINGFLOWSEPARATIONTHEFLOWSEPARATIONWILLCAUSETHEFLOWLOSSTHECHANGEOFTHEFLUIDVELOCITYMAGNITUDEANDDIRECTIONMAKESTHEFLOWANGLEBENOTEQUALTOTHEINSTALLATIONANGLEWHENTHEGASISGOINGINTOTHEIMPELLERFROMTHEENTRANCEANDOUTFROMTHEIMPELLER,RESULTINGINSHOCKLOSSTHESHOC
4、KLOSSWILLAFFECTTHEEFFICIENCYANDPERFORMANCEOFFANSWHENTHEGASWITHIMPACTSPEEDIMPORTEDINTOTHEENTRANCEOFIMPELLER,ITWILLBRINGABOUTTHEVORTEXONTHESUCTIONSIDETHISISTHEREASONLEADINGTOBOUNDARYLAYERSEPARATIONWITHTHEHELPOFGAMBIT，APROCESSINGTOOLOFFLUENTSOFTWARE，AINCOMPRESSIBLETURBULENCEMODELOFAKINDOFWINDTURBINEAIR
5、FOILWASBUILTOFCOURSE，THESPECIFICAIRFOILOFTHISSTUDYISG473UNDERTHEDIFFERENTANGLEOFATTACK，THEAERODYNAMICPERFORMANCEOF2DAEROFOILOFWINDTURBINEAIRFOILWASSIMULATEDANDANALYZEDBYUSINGTHEFLUENTSOFTWARETHEAOAOFTHISSTUDYWASCHANGEDFROM36TO8THENCOMPARETHESPEEDVECTORDIAGRAMSOBTAINEDVIATHEFLUENTSOFTWAREANDFINDOUTTH
6、ERELATIONBETWEENTHEBOUNDARYLAYERSEPARATIONONTHEWINDTURBINEAIRFOILANDTHEANGLEOFATTACKKEYWORDSWINDTURBINEAIRFOILBOUNDARYLAYERSEPARATIONNUMERICALSIMULATIONANGLEOFATTACK热能与动力工程毕业设计III目录摘要IABSTRACTII1绪论111研究背景及意义112国内外研究现状和趋势213研究方法及主要内容32翼型基本知识421几何参数422气动特性523影响气动特性的主要因素63数值模拟理论831边界条件的确定832K模型84数值模拟结果
7、及分析1041利用GAMBIT建立计算模型1042利用FLUENT进行模拟计算1143模拟结果分析15431对攻角为36时的模拟结果分析15432对不同的攻角时的模拟结果分析17433对相同大小的正负攻角的模拟结果进行分析19结论23参考文献24致谢26热能与动力工程毕业设计11绪论11研究背景及意义风机是一种装有多个叶片的通过轴旋转推动气流的机械。叶片将施加于轴上旋转的机械能，转变为推动气体流动的压力，从而实现气体的流动。风机广泛应用于发电厂、锅炉和工业炉窑的通风和引风，矿井、隧道、冷却塔、车辆、船舶和建筑物的通风、排尘和冷却等1。尤其是在电站，随着机组向大容量、高转速、高效率、自动化方向的
8、发展，电站也对风机的安全可靠性提出了越来越高的要求，锅炉风机在运行中常发生烧坏电机、窜轴、叶轮飞车、轴承损坏等事故，严重危害设备、人身安全，也给电厂造成巨大的经济损失2。此外，风机一直是电站的耗电大户，电站配备的送风机、引风机和冷烟风机是锅炉的重要辅机，降低其耗电率是节能的一项重要措施。气体经过风机叶轮后能够获得相应的动能，但是，由于结构、工艺及流体黏性的影响，气体流经风机时不可避免地要产生各种能量损失，而使其实际可利用的能量降低。因此，尽可能地减少气体在风机内部的能量损失，对提高风机的效率，降低能耗，保证风机的经济性、安全性有着十分重要的意义。气体流经风机时的损失，按其能量损失的形式不同可分
9、为三种机械损失、容积损失和流动损失3。当风机工作时，气体流道的几何形状改变会使流体运动速度的大小和方向发生改变，从而产生流动分离。流动分离产生的冲击会造成流动损失。流体运动速度的大小和方向的改变，也会使得气体在进入叶片入口和从叶轮出来进入压出室时，流动角不等于叶片的安装角，从而产生冲击损失，影响风机的效率和性能。由于气体进入叶片入口时存在着冲击速度，使气体在风机叶片的吸力面上形成旋涡，造成边界层分离现象而会导致能量损失4。现在，全球学者都达成了优化叶片的设计是提高电厂风机效率，从而节省能源的一个有效途径这个共识5。风机的流动损失不仅仅影响到风机的效率，在流动损失过大时，它还会影响到风机的安全运
10、行，引发事故，造成更大的经济损失。2004年大唐唐山热电有限责任公司2300M机组锅炉，风机叶片背面流动恶化，层流边界受到破坏，在叶片背面尾端出现涡流区，此时，风机全压急剧降低，保护系统开关动作，风机停运，发生事故6。离心风机是通过降低二次流涡，涡舌和喷气攻角造成的能量损失来改善风机的气动性能的7，叶轮叶片的气动性能是决定风机性能优劣的主要因素，而叶轮叶片的剖面形状翼型又是决定风机性能的关键因素。从局部流动特性来看，机翼型叶片风机的气流匹配能力更强，气动损失更小，因此，其稳定工作范围也较宽，具有优良的气动和变工况性能，尤其是电站锅炉负荷受各方面的影响经常发生变化，与之匹配的风机风量也要随之改变
11、，为了适应电站锅炉阻力变化小，而风量变化要求较大的特点，在选用离心通风机时，一般首先选用机翼型叶轮。翼型的气动性能参数的确定是风机叶片设计的重要内容，通过实验来获取风机叶片设计所需翼型的所有性能参数将要花费太多的人力和时间，因此翼型数值模热能与动力工程毕业设计2拟准确性成为了风机叶片的设计的重要课题8。运用FLUENT数值计算软件，对翼型流动进行二维数值模拟，对不同冲角下的流动情况进行详细的研究，找出冲角与分离点位置的关系，对预测风机安全经济运行范围和风机的高效可靠运行具有重要的指导意义。12国内外研究现状和趋势我国风机拥有量约230万台以上，年耗电量约占全部发电量的10左右，因此风机的节电有
12、着十分重要的意义。2007年我国风机新增装机29617万KW，累计达到55617万KW，分别同比增长121、114，预计2008年市场需求还将远远超出预期，国产设备的新装机容量年增长速度为60709。但中国仍有多种低效旧风机需要更新换代，新推广的风机也有待于进一步完善。因此，中国应该不断提高风机产品质量、稳定市场需求，还要积极引进先进技术，提高技术开发能力。但是，据统计，风机的电能利用率超过50的仅占总数的542，而电能利用率超过60的只有36，如果将风机运行效率提高10，全国就可以节电150亿千瓦时了10。因此，如何能以科技为基础，发展、优化风机，从而提高其性能，降低经济损失，并将其转化为效
13、益成为一个十分重要的课题。2001年，山西原平化学工业集团有限责任公司的刘天灵，咸高创11就通过对风机的轴受力进行了数值计算分析，知道了该厂风机经常出现故障的原因。现代风机特点是转速高，压力大，叶轮流道窄，线速度高，叶轮所受传动扭的矩大，受力状态复杂且大，这要求叶轮制造有很高的精度12，因此对叶轮叶片的研究和设计是风力发电技术研究和开发的重要任务。传统风机的叶片多采用固定的翼型，但由于应用环境和应用目的不同，风机翼型的叶片并不能高效地进行能量转换。发达国家从20世纪80年代中期开始研究风机新翼型，并发展了各自的翼型系列。我国对风机翼型的研究主要在于测绘和仿制上，并且仅限于进行一些风机试验，由于
14、商业因素和技术保密等原因，我们不容易得到国外风机专用翼型相关的气动实验数据。一些设计和制机专利都是从国外引进的，严重制约了我国风机产业的发展。开发具有我国自主知识产权的风机翼型系列，研制我国新型高效的风机叶片，对促进我国风机事业的发展至关重要13。西华大学能源与环境学院的黄华，张礼达14基于翼型理论和线性动量理论对叶片翼型截面升力公式的计算，导出对非设计工况来流角计算的迭代式。应用牛顿拉普森迭代法对来流角进行计算，根据结果再计算叶片截面的升力、推力、切向力、功率等气动参数，提出一种风力机叶片翼型气动性能的计算和校核设计方法。传统风机设计是以实验为基础的设计，通过反复的设计计算和实验来确定最终设
15、计改进方案，设计周期长，费用也较高，对经验的依赖性较强，而USED技术已经改变了工程设计方法，它是一个用于分析流体现象和减少设计时间的有力工具15。2008年西华大学风电技术研究所的毛金铎，张礼达16应用USED流体力学软件对风力机叶片常用翼型THATHOW211进行数值分析，得出了其升力系数、阻力系数、升阻比以及翼型表面压力随来流攻角变化关系，并依据计算结果对FFAW3211翼型的气动性能进行分析。风机叶片翼型设计理论是决定风机功率特性和载荷特性的根本因素，一直是各国学者研究的热热能与动力工程毕业设计3点。现有翼型的表达都是通过离散的点来实现的，并不存在函数的具体表达形式。新翼型的设计也是基
16、于原有的翼型坐标，对其进行局部的调整，以获得性能更为优越的翼型。利用FLUENT有限元软件能很好地模拟离心风机流场，计算出风机的性能参数，可以节约成本，减短设计周期，并且能得到极具实际指导意义的结论。2008年辽宁工程技术大学机械工程学院李文华，范兴文17采用CFD商用软件FLUENT61对离心风机内部流场进行了三维数值模拟。计算中采用了标准K湍流模型与非结构化网格。通过模拟发现了蜗舌对叶轮中流动的影响和部分空气在叶轮中的螺旋状流动，捕捉到了离心通风机内部许多重要的流动现象，同时对计算结果进行了分析，对该类风机的性能改进提供了一定的依据。2009年重庆大学机械传动国家重点实验室的陈进，张晓，王
17、旭东18对某翼型扰流流动建立了二维可压缩湍流模型，利用商业软件FLUENT对翼型不同来流攻角下的气动特性进行了相应的数值模拟计算。湍流黏度采用基于RANS的SPALAALLMARAS湍流模型处理，得出了雷诺数在32106时，某翼型的升力系数、阻力系数和压力分布随来流攻角的变化关系。数字仿真能比真实实验提供更多结果，而且可以用于核实和完善实验结论19，可深化了解风机翼型的气动性能，对不同冲角下的流动情况进行详细的研究，找出冲角与分离点位置的关系，为风机叶片翼型选型和叶片翼型改型设计和研发工作提供技术参数和指导意见。13研究方法及主要内容由于叶轮机械内部流场非常复杂，并带有强烈的非定常特征，进行细
18、致的实验测量非常困难，目前尚没有完善的流体力学理论解释诸如流动分离、失速和喘振等流动现象，这就迫切需要可靠详细的流动实验和数值模拟工作来了解机械内部流动本质。本文将利用FLUENT软件对风机翼型叶片进行二维的数值模拟，研究空气以不同的方向流入翼型叶片入口所造成的流动分离。根据数值模拟的一般步骤利用GAMBIT创建二维模型，进行网格划分，设定边界条件和区域，输出网格，再利用FLUENT求解器求解，对不同空气来流攻角角下的流动进行二维数值模拟。在得到模拟结果后,对不同攻角下模拟所得到的速度矢量图进行比较分析，得出风机翼型边界层分离和攻角的关系。热能与动力工程毕业设计42翼型基本知识21几何参数翼型
19、的气动性能直接与翼型外形有关。通常，翼型外形由下列几何参数决定201翼弦CBO气动弦几何弦图21翼型的气动弦与几何弦翼型前缘点O与尾缘点B之间的连线称翼弦，翼弦OB的长度称作弦长，以C表示，它是翼型的基准长度，也称为几何弦。除几何弦外，翼型还有气动弦。当气流方向与气动弦一致时，作用在翼型上的升力为零，如图21所示。对称翼型的几何弦与气动弦重合，气动弦又称零升力线。2前缘半径和前缘角翼型前缘点的内切圆半径称为翼型前缘半径，亚音速翼型前缘是圆的，超音速翼型前缘是尖的。前缘点上下翼面切线的夹角就是前缘角。3厚度和厚度分布在计算翼型时通常采用如图22所示的直角坐标，X轴与翼弦重合，Y轴过前缘点。且垂直
20、向上。这样在X轴上方的弧线称为上翼面以XYU表示，下方的弧线称为下翼面以XYL表示。对应同一X坐标的上下翼面点距为翼型的厚度，以T表示，见图22。厚度随X的变化称厚度分布，以TX表示ULTXXXYY当XXC时，TTMAX称最大厚度。/CXTTC称为最大相对厚度，XC为最大厚度位置，其无因次量为CXXCC。通常，翼型的相对厚度即指最大相对厚度，以T表示。4中弧线气动弦几何弦BO热能与动力工程毕业设计5翼型内切圆圆心的连线叫做中弧线。只有对称翼型时中弧线与翼弦重合。CXCTMAXXTX上翼面下翼面XY图22翼型的厚度分布5弯度和弯度分布翼型中弧线和翼弦间的高度称为翼型的弯度，弧高沿翼弦的变化称为弯
21、度分布，以XYF表示XYXY/XYLUF21当XXF时，YXYFFFMAX称为最大弯度，以F表示。CFF称为最大相对弯度，XF为最大弯度位置，其无因次量为CXXFF。同样，通常翼型的相对弯度指最大相对弯度，用F表示。6尾缘半径和尾缘角翼型尾缘点B的内切圆半径称为翼型尾缘半径。若尾缘为尖的，则以尾缘点上下翼面的切线夹角表示，称为尾缘角。有的翼型尾缘是平的，则用尾缘厚度表示。22气动特性翼型所受的力是作用在上下表面的分布力之合力。表面力有两种，一种是法向力，即压力；一种是切向力，即摩擦力。这里定义和远前方来流相垂直的合力为升力，而与远方来流方向一致的合力为阻力。也像压强通常表示为无量纲的压强系数一
22、样，升力和阻力通常也表示为无量纲的升力系数CL和阻力系数CD，二者定义如下CBVLCL22121CBVDCD22122式中的L和D分别代表升力和阻力，单位为N；来流动压头为1/2V2，单位是N/M2，C是弦长，单位是M；B是垂直于纸面的尺寸，单位是M，为空气密度。热能与动力工程毕业设计6翼型是用来产生升力的，也就是说，要产生一个垂直于与翼型几何弦成小角度的入射流的力。翼型的几何形状和作用在翼型上的力如图23所示。相对速度V与翼型几何弦的夹角叫攻角。翼型的分布压力有个合力即升力，这个合力和翼弦的交点称为压力中心。压力中心的位置和翼面上的压力具体分布情况有关系。当攻角增大时未出现大分离以前，不仅上
23、翼面的吸力和下翼面的压力都增强了，而且吸力峰前移，结果压力中心前移。C旋转平面升力阻力合力V翼弦图23作用在翼型上的力从理论力学知道，一个平面力系是可以合成作用在某个指定点上的一个力和一个力矩。翼型上的分布压力也可以合成一个力升力和一个力矩，这个力矩名为俯仰力矩。这个指定点是一个特殊的点，称为气动中心，或者焦点。不论攻角多大，压力增大，压力中心前移，压力中心至气动中心的距离缩短，结果力乘以力臂的积，即俯仰力矩保持不变。翼型上确有这样的一个力矩值不变的点，是理论上证明了的。这一点的理论位置，薄翼型在距前缘1/4弦长处。实验测得的略有出入，大多数普通翼型的气动中心位于023024弦长处，而层流翼型
24、的则在026027弦长处。俯仰力矩系数记为CM，定义式如23所示。BCVMCM222123规定抬头力矩为正，低头力矩为负。俯仰力矩系数是翼型的重要气动参数之一，计算全机的平衡时必须用到它。23影响气动特性的主要因素1雷诺数V热能与动力工程毕业设计7影响低速翼型特性的最重要的流体因素是流体的粘性，它间接产生升力而直接产生阻力和造成流体分离。这种影响用翼型和流体组合的雷诺数来表示。现代风力机上翼型代表性的弦长典型地在3/4叶展处范围是从小型风力机的03米到兆瓦级风力机的2米。尖端速度通常是从45M/S到90M/S，因此水平轴风力机叶片3/4叶展处的切向速度的范围大概是从34M/S到68M/S。那么
25、，对于风力机翼型，雷诺数的范围是从10106一直到10706。这表明风力机翼型通常都不运行在敏感的低雷诺数范围一般低于10506，在这个敏感范围中，入射流湍流的变化、翼型自身的振动或翼型表面的粗糙度都会引起翼型性能的很大变化。雷诺数对翼型的升力特性和阻力特性有着重要的影响。随着雷诺数增加，升力曲线斜率增加，最大升力系数增加，失速攻角增加；随着雷诺数增加，最小阻力系数减小；同时雷诺数增加，翼型升阻比也增加。在低雷诺数RE1036情况下，翼型表面从层流边界层经过转捩发展为湍流边界层。不同的边界层发展情况对翼型空气动力特性，特别是阻力特性有较大的影响。层流翼型有较低的阻力系数和较高的升阻比。2粗糙度
26、翼型表面由于材料、加工能力以及环境的影响，使表面不可能绝对光滑，而总是凹凸不平。这些凹凸不平的波峰与波谷之间高度的平均值称为粗糙度。翼型表面的粗糙度对翼型气动特性有直接影响。通常粗糙的型面和光滑的型面相比，翼型的升力系数降低，阻力系数增加。当然其影响程度还和雷诺数、翼型形状等有关。通常翼型前缘向后到2030弦长处的上下表面对翼型气动特性影响尤为显著。3湍流度湍流度对翼型气动特性也密切相关。通常情况下，湍流度增加，翼型的阻力系数和最大升力系数增加，最大升阻比减小。4攻角上面叙述的气动特性大多是在中小攻角范围内的情况，在大攻角情况下其变化要复杂得多。风力机叶片的工况是很宽的，不仅涉及小攻角情况，而
27、且涉及到失速和大攻角范围的升力和阻力特性。由于大攻角范围的气动特性变化较复杂，纯理论计算很困难，因而大多依靠相应的实验和数值模拟得到较可靠的结果。热能与动力工程毕业设计83数值模拟理论31边界条件的确定FLUENT提供了十余种类型的进口、出口边界条件21，下面将本文涉及到的边界条件介绍如下1速度入口VELOCITYINLET给定入口边界上的速度及其他相关标量值。该边界条件适用于不可压缩流动问题，对可压缩问题不适合，否则该入口边界条件会使入口处的总温或总压有一定波动。输入量包括速度大小、方向或各速度分量、周向速度轴对称有旋流动、静温考虑能量等。2压力出口PRESSUREOUTLET对于有回流的出
28、口，压力出口比自由出流更容易收敛。给定出口边界上的静压强表压强。该边界条件只能用于模拟亚音速流动。如果当地速度已经超过音速，该压力在计算过程就不采用了。压力根据内部流动计算结果给定。其他量都根据内部流动外推出边界条件。该边界条件可以处理出口有回流的问题，合理的给定出口回流条件，有利于解决有回流出口问题的收敛困难问题。出口回流条件需要给定回流总温如果有能量方程、湍流参数湍流计算、回流组分质量分数有限速率模型模拟组分输运、混合物质量分数及其方差PDF计算燃烧。如果有回流出现，给定的表压将视为总压，所以不必给出回流压力。回流流动方向与出口边界垂直。3固壁边界WALL对于黏性流动问题，FLUENT默认
29、设置是壁面无滑移条件。对于壁面有平移运动或者旋转运动时，可以指定壁面切应力和与流体换热情况。壁面热边界条件包括固定热通量、固定温度、对流换热系数、外部辐射换热与对流换热等。32K模型K模型是两方程湍流模型中最具代表性的，同时也是工程中应用最为普遍的模式。湍流被称为经典力学的最后难题，原因在于湍流场通常是一个复杂的非定常、非线性动力学系统，流场中充满着各种大小不同的涡结构。整个湍流场的特征取决于这些涡结构的不断产生、发展和消亡，同时，这些涡结构之间又不断发生着复杂的相互作用，这就使得对湍流现象的理解、描述和控制变得十分困难。对于单相流动，科学界已经有较为成熟的湍流封闭模型。K模型包括标准的K模型
30、，RNGK模型和可实现的K模型，下面简单介绍一下221标准的K模型最简单的完整湍流模型是两个方程的模型，要解两个变量，速度和长度尺度。在FLUENT中，标准K模型自从被LAUNDER和SPALDING提出之后，就变成工程流场计算中热能与动力工程毕业设计9主要的工具了。适用范围广、经济、合理的精度。它是个半经验的公式，是从实验现象中总结出来的。湍动能输运方程是通过精确的方程推导得到，耗散率方程是通过物理推理，数学上模拟相似原型方程得到的。应用范围该模型假设流动为完全湍流，分子粘性的影响可以忽略，此标准K模型只适合完全湍流的流动过程模拟。2RNGK模型RNGK模型来源于严格的统计技术。它和标准K模
31、型很相似，但是有以下改进A、RNG模型在方程中加了一个条件，有效的改善了精度。B、考虑到了湍流旋涡，提高了在这方面的精度。C、RNG理论为湍流PRANDTL数提供了一个解析公式，然而标准K模型使用的是用户提供的常数。D、标准K模型是一种高雷诺数的模型，RNG理论提供了一个考虑低雷诺数流动粘性的解析公式。这些公式的作用取决于正确的对待近壁区域。这些特点使得RNGK模型比标准K模型在更广泛的流动中有更高的可信度和精度。3可实现的K模型可实现的K模型是近期才出现的，比起标准K模型来有两个主要的不同点可实现的K模型为湍流粘性增加了一个公式，为耗散率增加了新的传输方程，这个方程来源于一个为层流速度波动而
32、作的精确方程。术语“REALIZABLE”，意味着模型要确保在雷诺压力中要有数学约束，湍流的连续性。应用范围可实现的K模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测，而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。可实现的K模型和RNGK模型都显现出比标准K模型在强流线弯曲、旋涡和旋转有更好的表现。由于带旋流修正的K模型是新出现的模型，所以现在还没有确凿的证据表明它比RNGK模型有更好的表现。但是最初的研究表明可实现的K模型在所有K模型中流动分离和复杂二次流有很好的作用。该模型适合的流动类型比较广泛，包括有旋均匀剪切流，自由流射流和混合层，腔道流动和边界层
33、流动。对以上流动过程模拟结果都比标准K模型的结果好，特别是可再现K模型对圆口射流和平板射流模拟中，能给出较好的射流扩张。模型评价可实现的K模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然的湍流粘度，这是因为可实现的K模型在定义湍流粘度时考虑了平均旋度的影响。这种额外的旋转影响已经在单一旋转参考系中得到证实，而且表现要好于标准K模型。由于这些修改，把它应用于多重参考系统中需要注意。热能与动力工程毕业设计104数值模拟结果及分析整个研究过程是空气来流速度为50M/S，攻角在36到8之间，每两度进行一次模拟计算的过程。41利用GAMBIT建立计算模型下面选取0为例，讲述在整个研究过程中GA
34、MBIT的使用。1建立翼型轮廓和设定流动区域设定流动区域，如图41所示，其中，翼型弦长为224CM，流动区域左边为半径为672CM的半圆，右边为1124134CM2的矩形；翼型局部放大图如图42所示。图41流动区域图图42翼型局部放大图2划分网格和翼型边界层各条边上的节点数如表41所示，得到的整体网格效果图如图43所示表41各条边的上网格划分的节点数边ABBCCDDEEHHFAFFIIJJEAD节点数200545200200353520070570470ABCDFEGHFEIJHIJ热能与动力工程毕业设计11图43整体网格效果图3定义边界类型各边界的类型和所包含的边线如表42所示。表42各边界
35、的类型和所包含的边线组名包含的边线类型INLETADVELECITYINLETOUTLETBCPRESSUREOUTLETBODYSHANGHF，FIWALLBODYXIAHE，EJWALLBODYWEIIJWALLOPENSHANGABWALLOPENXIACDWALL4保存文件并输出网格42利用FLUENT进行模拟计算下面选取0为例，讲述在整个研究过程中FLUENT软件的使用。热能与动力工程毕业设计121读取翼型的MESH文件。2检查网格，最小面积不能出现负值。3确定长度单位为CM。4求解器参数设置设置参数如图44所示。5确定紊流模型选择KEPSILON2EQN模型，设置如图45所示。图4
36、4求解器设置对话框图45紊流模型选择对话框6不选用能量方程。7确认流体的物理属性确认选择流体为无黏空气，密度为1225KG/M3。8确定工作压强设置工作压强为默认的101325PA。热能与动力工程毕业设计139定义边界条件INLET和OUTLET打开边界条件设置对话框，OUTLET的设置如图46所示；INLET的设置如图47所示。其中，K，是由K计算器根据来流速度，特征长度和流体的动力粘度求出。图46压力出口设置对话框图47速度入口设置对话框10设置求解控制参数选择二阶求解。11求解初始化流场初始值设为入口流动的参考值。12设置求解过程残差监视器收敛准则设为0001。13求解过程升力监视器设置
37、打开FORCEMONITORS对话框，在COEFFICIENT项选择LIFT，设置如图48所示。其中X项是由SIN00，Y项由COS01确定。14求解过程阻力监视器设置在COEFFICIENT项选择DRAG，设置如图49所示。其中X项是由COS01，Y项由SIN00确定。热能与动力工程毕业设计14图48升力监测设置对话框图49阻力监测设置对话框15为迭代计算设置基本参考值在COMPUTEFROM项选择INLET，在REFERENCEZONE项选择FLUID。16保存CASE文件文件名为“FENGJIYIXING0”。17迭代求解计算先设置迭代1000次，计算开始。经过296次迭代后，计算收敛，
38、残差曲线、升力曲线，阻力曲线如下A残差监测曲线如图410所示；B升力监测曲线如图411所示；C阻力监测曲线如图412所示。图410迭代过程残差监测曲线热能与动力工程毕业设计15图411迭代过程升力监测曲线图412迭代过程阻力监测曲线18保存DATE文件文件名为“FENGJIYIXING0”。至此，来流攻角为0情况下数值模拟的过程结束了。可按照上述的步骤对其它攻角情况进行相同的模拟。43模拟结果分析431对攻角为36时的模拟结果分析攻角为36时模拟得到速度矢量图如图413所示，压强分布云图如图414所示；对速度矢量图进行放大，得到翼型附近的流场局部放大图如图415所示热能与动力工程毕业设计16图
39、413攻角为36时的速度矢量图图414攻角为36时的压力分布云图图415攻角为36时翼型附近的流场局部放大图热能与动力工程毕业设计17由图413，图414和图415可以看出，当风机翼型来流速度攻角较大的时候，会在翼型背面形成旋涡，破坏流场的线型，同时边界层分离，影响到风机翼型的气动性能。边界层分离出现的原因是，在分离点以后，壁面附近被黏性和逆压梯度滞止的流体质点逐渐增多，压强的进一步升高使被滞止的质点发生回流，而排挤上游来流边界层使其与壁面分离。流场中速度等于零的流体质点成为顺流和回流的分界面，该分界面极不稳定，稍经扰动便破裂形成旋涡被主流带走。这样，分离点后的旋涡不断地产生，又不断地被主流带
40、走，就在翼型的背面形成涡流区。432对不同的攻角时的模拟结果分析攻角为30图416，24图417，18图418，12图419，16图420，0图421时翼型附近的速度矢量图，如下图416攻角为30时翼型附近的速度矢量图图417攻角为24时翼型附近的速度矢量图热能与动力工程毕业设计18图418攻角为18时翼型附近的速度矢量图图419攻角为12时翼型附近的速度矢量图图420攻角为6时翼型附近的速度矢量图热能与动力工程毕业设计19图421攻角为0时翼型附近的速度矢量图比较图415到图421，可以看出，当来流攻角越大时，流场中形成的旋涡越大，风机翼型的边界层分离现象越明显，即是流场的线型受到的破坏越大
41、。而当来流攻角很小的时候，风机翼型就没有出现边界层分离现象，则风机翼型的气动性能较好。在很小角度攻角来流冲击的时候，翼型附近的流场都是顺压梯度DP/DX0区，壁面附近被黏性和逆压梯度滞止的流体质点逐渐增多，压强的进一步升高使被滞止的质点发生回流，在翼型的背面形成涡流区。另外，逆压梯度区足够长的话，逆压差和层外势流的减速使得边界层中流动减速，而近壁出流动的动能也愈来愈小，故在黏性和逆压梯度的双重作用下，流体质点会在壁面某处被停滞，因此，可以说攻角的角度越大，逆压梯度越大，边界层分离现象越容易出现，现象越明显。433对相同大小的正负攻角的模拟结果进行分析为了对大小相同，正负相反的攻角下模拟得到的速
42、度矢量图作一对一比较分析，继续做了攻角为10到36情况下的数值模拟正常的运行工况下是不会出现这么大的正角度攻角的，这里只是为了得到更为明显的模拟结果来进行比较而做的附加数值模拟攻角为36图422，30图423，24图424，18图425，12图426，6图427时翼型附近的速度矢量图，如下热能与动力工程毕业设计20图422攻角为36时翼型附近的速度矢量图图423攻角为30时翼型附近的速度矢量图图424攻角为24时翼型附近的速度矢量图热能与动力工程毕业设计21图425攻角为18时翼型附近的速度矢量图图426攻角为12时翼型附近的速度矢量图图427攻角为6时翼型附近的速度矢量图热能与动力工程毕业设
43、计22对图415和图422，图416和图423，图417和图424，图418和图425，图419和图426，图420和图427做分组比较分析，得出当来流攻角的大小一样时，风机翼型在来流攻角为正时出现的旋涡更大，即风机翼型在受到正攻角的来流冲击时，风机翼型非工作面更容易出现旋涡，更容易出现边界层分离现象。因为在G473型风机翼型的上翼面的压强梯度DP/DX0点后的逆压梯度区内的翼型壁面变化更快，边界层外势流通道更宽，则流速下降得更快，边界层外缘处的流速减少得更快，压力增加得更快，即逆压梯度更大，所以，攻角大小相同的话，正角度时风机翼型非工作面更容易出现旋涡，边界层分离现象更明显。热能与动力工程毕
44、业设计23结论本文用FLUENT软件对G473风机翼型受到速度大小为50M/S的空气来流在36到8攻角下为方便比较而增加了10到36攻角下的数值模拟冲击的情况进行了二维数值模拟，通过对模拟的结果进行分析，得出下面的结论1风机翼型在受到较大攻角的来流冲击时，就会使流场出现旋涡，边界层分离，而且攻角越大，现象越明显，即风机气动性能受到的影响越大；2当攻角大小一样时，攻角为正时风机翼型流场里出现的旋涡更大，边界层分离现象更明显。由于所做的只是二维的数值模拟，与三维的现实现象有出入，所以模拟的结果可能会有偏差，另外本人所学知识的有限，所做的分析还不够全面，不够详尽，存在很多不足，希望能在以后的学习工作
45、中能够对其进一步改善。热能与动力工程毕业设计24参考文献1安连锁泵与风机M北京中国电力出版社,20012袁春杭锅炉引风机事故的预防J中国锅炉压力容器安全,2005,14638393蔡兆林,吴克启,颖达离心风机损失的计算J工程热物理学报,1993,14153564王松岭流体力学M北京中国电力出版社,20045MSBHOWMICK,SCBERASTUDYTHEPERFORMANCESOFINDUCEDFANSANDDESIGNOFNEWINDUCEDFANFORTHEEFFICIENCYIMPROVEMENTOFATHERMALPOWERPLANTJIEEEXPLORE,2008,479166田小
46、冬,刘福君锅炉送风机失速故障一例J热力发电,2007,1079807SONGLINGWANG,LEIZHANG,ZHENGRENWU,HONGWEIQIANOPTIMIZATIONRESEARCHOFCENTRIFUGALFANWITHDIFFERENTBLADENUMBERANDOUTLETBLADEANGLEJIEEEXPLORE,2009,978148巫发明,杨从新,张玉良,赵伟国风力机翼型摩擦阻力数值计算中不同湍流模型的比较研究J流体机械,2008,123611149风机制造业可看好至2012年JGM通用机械,2008,55010任仁良风机效率测试仪的研制J仪器仪表学报,2001,32
47、230931111刘天灵,咸高创D3OO12煤气风机的故障分析及改进措施J,山西机械,2001,3484912蔡育新离心鼓风机叶轮制造工艺的改进J机电工程技术,2001,4586013张果宇,冯卫民,刘长陆,俞剑锋6种风力机叶片翼型的气动性能数值模拟研究J可再生能源,2009,227111514黄华,张礼达风力机叶片翼型气动性能设计计算方法的分析与研究J能源工程,2007,3454715ZHANGBO,HUANGSHUJUAN,SUNZHONGWEI,HUISHIENPERFORMANCEMPROVEMENTOFCENTRIFUGALFANBYMEANSOFNUMERICALSIMULATIO
48、NJCHINESEJOURNALOFMECHANICALENGINEERING,2006,191555816毛金铎,张礼达风轮叶片翼型气动特性分析J机械设计与制造,2008,11,394117陈进,张晓,王旭东通用风力机翼型气动特性数值模拟J重庆工学院学报,2009,723515418李文华,范兴文630型风机流场的三维模拟与分析J科技创新导报,2008,31161719GUANYUPU,ZHAOZHENHUA,CHENWEI,GAODEPINGFOREIGNOBJECTDAMAGETOFANROTORBLADESOFAEROENGINEPARTIINUMERICALSIMULATIONOFB
49、IRDIMPACTJCHINESEJOURNALOFAERONAUTICS,200821328334热能与动力工程毕业设计2520李杰风力机叶片翼型气动特性数值研究D北京华北电力大学,200921韩占忠FLUENT流体工程仿真计算实例与分析M北京北京理工大学出版社,200922韩占忠,王敬,兰小平FLUENT流体工程仿真计算实例与应用M北京北京理工大学出版社,2004热能与动力工程毕业设计26致谢经过三个多月的毕业设计，终于迎来了论文完稿的这一刻，在这个过程中，有过成功的喜悦，有过失败的沮丧，也有过彷徨和疑惑，这一切都让我感觉到求知过程的挑战和艰辛，同时也让我深刻明白学无止境，我将在以后的学习工作中继续努力，探索真知。本文是在导师王松岭老师的悉心指导和热情鼓励下顺利完成的。王老师渊博的知识、严谨的治学态度、诲人不倦的作风和兢兢业业的工作精神，让我受益匪浅。在做论文期间，王老师关心每一个环节，及时提出建议，悉心指导，直到我的论文顺利完成，而我在王老师身上学到的不只是学问，更是对科研孜孜不倦的追求精神。在此非常感谢王老师对我的谆谆教诲
