风力机叶片翼型的结冰数值模拟研究

研究了风力机叶片翼型的结冰数值模拟方法,及翼型结冰后对其气动性能的影响。求解雷诺平均N-S方程,引入k-ωSST湍流模型封闭方程,获得风力机叶片流场;采用拉格朗日法计算风力机叶片翼型周围的水滴运动轨迹,同时考虑了多尺寸水滴分布的影响,获得翼型表面的局部水收集系数分布;根据质量守恒和能量守恒原理,计算翼型表面的各项热流,获得翼型表面的结冰速率和结冰冰形;考虑到翼型结冰是时间的动态函数,采用多时间步长法完成结冰数值模拟。计算了风力机叶片翼型在不同环境条件和气象条件下的结冰冰形,同时模拟了风力机翼型结冰后周围流场的变化,并与干净翼型的气动特性进行了对比。结果表明,环境温度较高时形成的明冰对翼型气动性能的影响较大,结冰导致翼型升力下降,升阻比减小,最大减小幅度达到61%,同时结冰后的翼型会提前进入失速区,导致桨叶气动性能恶化。     

基于N-S方程和自由转捩预测耦合求解的钝后缘翼型气动性能计算

大厚度钝后缘翼型由于结构和气动性能上的优点,被用于作为大型风力机叶片设计的内侧翼型。而由于其大厚度的特点和风洞实验阻塞度的限制,大雷诺数的风洞实验数据很少,给此类翼型的设计和使用带来困难。本文基于雷诺平均Navier-Stokes(RANS)方程和自由转捩预测耦合求解方法,进行了大厚度钝后缘翼型的气动性能计算研究;针对由于钝后缘后的涡脱落而造成的翼型表面压力分布的脉动,发展了有时均效应的转捩模型以考虑这种周期脉动的时均效应,并对传统的耦合求解方法进行了修改。通过对风力机翼型DU97-W-300的钝后缘改形DU97-Flat翼型的气动性能计算,分析了网格数对计算结果和计算效率的影响。并将计算的DU97-W-300翼型的气动性能和实验结果进行了比较,表明通过本文发展的耦合求解方法可以在更少的计算网格时得到比参考文献更吻合风洞实验结果的气动性能,为此类翼型的设计和使用提供数值计算基础。     

后缘锯齿襟翼对翼型气动特性的影响研究

针对在翼型后缘增加锯齿襟翼与格尼襟翼的差异性,以及二者对翼型的空气动力特性有何改善等问题,在某民用运输机机翼翼型下表面后缘处增加锯齿襟翼,对其进行了空气动力特性试验和计算。结果表明,增加一定高度的后缘锯齿襟翼,能够减弱和改变翼型下表面后缘处气流的流动,从而减少和改变翼型下表面尾涡的涡量和涡量变化率,达到增加翼型升力、减小翼型阻力的目的。     

转捩对风力机翼型和叶片失速特性影响的数值模拟

采用基于k-ωSST湍流模型的Gamma-Theta转捩模型对S809翼型和NREL Phase VI叶片进行了考虑转捩的气动力数值模拟,研究了转捩对其气动特性特别是失速特性的影响。首先对S809翼型在迎角0°～30°范围内开展了数值模拟,比较了转捩模拟和全湍流模拟获得的翼型升力系数和流场特征,发现前缘层流分离泡的存在显著影响了翼型的失速特性;然后对NREL Phase VI叶片开展了类似的数值模拟,结果显示转捩对叶片失速特性和翼型失速特性产生影响的作用方式是相似的。经过对以上数值模拟结果加以分析认为,转捩对翼型和叶片失速特性的影响主要通过前缘层流分离泡的作用体现出来,前缘层流分离泡的存在使翼型和叶片更早地进入深失速。     

小型垂直轴风力机叶片结冰风洞试验与数值计算

利用风洞试验和数值模拟相结合的手段,研究了小型垂直轴风力机叶片在旋转状态下的结冰特性以及结冰后翼型与风力机的气动特性变化,以期为建立较复杂的大型水平轴风力机叶片旋转试验系统和研究其结冰机理、防除冰技术提供参考。试验在东北农业大学自行设计的利用自然低温的结冰风洞中进行,获得了采用NACA0018翼型的小型2叶片垂直轴风力机风轮在5种尖速比下的结冰分布：风力机叶片结冰遍布叶片整个表面,随着尖速比的增大,结冰形状出现不对称性。同时,数值模拟结果表明：叶片结冰后,随着尖速比的增加和结冰量的增多,升力系数降低阻力系数增大的趋势明显,风力机的功率系数也随之下降。分析发现,叶片结冰导致不同旋转角下叶片翼型周围的压力场和速度场发生了不同程度的变化,从而气动特性发生变化,影响了风力机性能。     

后缘隔板对充气式机翼翼型气动力影响分析

以充气式机翼为背景,研究了安装后缘隔板充气式翼型的气动特性。以翼型E387为基础生成的钝后缘翼型为研究对象,采用大涡模拟方法对安装中部隔板的充气式翼型进行了流场数值计算,并与未安装隔板的情况进行了对比。总结了后缘中部隔板对充气式翼型气动特性影响的变化规律,为充气式机翼设计提供参考。     

前缘带光滑霜冰的NACA0012翼型表面声学特性计算

结冰将改变飞机空气动力表面形状,不仅使飞机空气动力性能下降,还会导致气动噪声的变化。为研究结冰对翼型气动噪声的影响,采用计算流体力学方法对前缘带光滑霜冰的NACA0012翼型表面声学特性进行了数值计算。采用C型网格拓扑结构对结冰翼型的计算区域进行了划分,采用不可压缩雷诺平均N-S方程对结冰翼型周围黏性流场进行了数值计算,采用基于Proudman理论的宽频噪声模型和Curle的表面积分方法预测了结冰翼型的表面声学参数,获得了沿结冰翼型弦向分布的表面声功率和表面声功率级。研究表明,0°或小攻角时,靠近前缘霜冰区域的流动转捩或流动分离使结冰翼型的表面声功率更高;较大攻角时,靠近后缘的区域发生流动分离,使后缘的表面声功率增加,进一步增加了结冰翼型的表面声功率。前缘霜冰产生的流动转捩和流动分离是结冰翼型气动噪声增加的主要原因。     

零质量射流对翼型增升影响的数值模拟研究

针对安装了激振器的NACA0012二维翼型,采用数值模拟的方式对零质量射流技术进行研究.流场解算器采用了中国航空工业空气动力研究院计算所自主研发的ENSMB V1.2计算平台.计算结果给出了湍流模型对零质量射流数值模拟技术的影响,并且分析了激振器参数对翼型气动力特性的影响规律.     

基于致动线方法的风力机气动数值模拟

致动线方法是通过引入体积力代替叶片的致动线技术与三维Navier-Stokes方程相结合来获得风力机周围流场信息的一种方法。该方法避免了花费大量网格与计算资源去求解风力机叶片的附面层,从而可以把更多的网格与计算资源用于风力机尾流流场的模拟,非常适合用于风力机尾流流场的研究。以NH1500叶片为计算模型,从叶片载荷分布和功率系数两个方面,将引入高斯分布的致动线方法的数值模拟结果与BEM理论、CFD方法以及风洞实验进行了系统的比较,验证致动线方法用于风力机气动数值模拟的可行性,并且对致动线方法计算出的尾流流场进行简要的分析。     

大厚度钝后缘翼型气动性能计算研究

大厚度钝后缘翼型由于其气动和结构的优点,近来被作为大型风力机叶片设计时的内侧翼型成为风力机翼型设计的热点之一EULL由于此类翼型厚度和风洞阻塞度的限制,大雷诺数的风洞实验数据很少。耦合RANS方程和基于线性稳定性分析EULL捩预测方法,进行了大厚度钝后缘翼型气动性能计算研究。使用Gridgen程序生成计算网格,湍流模型为S-A模型。对常用的钝后缘修形方式如直接截断、对称增加厚度等并且有风洞实验结果的几种翼型进行了CFD计算并和实验结果及文献计算结果进行了比较。结果表明升力系数在线性段,计算结果和实验结果吻合很好,但是对失速迎角的捕捉能力差;阻力系数计算结果和实验值吻合较好。不同修形方式的计算结果为钝后缘修形设计提供了参考。     

桨叶结冰对尾桨气动性能的影响

基于某直升机尾桨桨叶地面结冰试验数据,建立了结冰桨叶的气动计算模型。对NACA0012翼型的气动特性计算表明,翼型升力系数和阻力系数的计算结果与试验数据吻合良好。然后根据所建模型,利用Fluent软件分别计算了某尾桨桨叶翼型结冰前后剖面的气动特性,发现结冰使桨叶翼型升力系数降低,阻力系数增大。最后采用动量一叶素理论结合...     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

低雷诺数下螺旋桨翼型非定常气动性能的比较

为得到适用于平流层螺旋桨桨叶叶素的翼型,采用数值方法,通过改变气流速度和攻角,比较研究了低雷诺数下8种典型翼型的非定常气动性能.结果显示,大攻角工况下翼型容易发生流动分离现象,翼型吸力面的分离越严重,气动力的振荡幅度越大.通过比较不同攻角和不同速度下翼型的升阻比大小和气动稳定性得出,SD 8000 - PT翼型具有较高的升阻比和较稳定的气动力,适合作为螺旋桨的翼型.     

GAW-1翼型前后缘变弯度气动性能研究

传统增升装置主要用于提高飞机起降气动性能。利用计算流体力学(CFD)的方法,引入了通用飞机翼型的前后缘变弯装置的概念,数值模拟了GAW-1翼型在爬升状态时,前缘变弯装置、后缘襟翼/副翼偏转以及前后缘装置综合偏转对翼型气动特性的影响。研究表明,前缘变弯装置可以有效地改善翼型的失速特性,失速迎角提高了3°左右,最大升力系数提高了4.56%;同时提高升阻比50%~120%;但在设计升力系数下,升力系数和阻力系数都略微减小。另一方面,后缘变弯装置可以改变最大升阻比所对应的迎角,以及在小迎角时,提高升力系数6%左右。翼型综合偏转可以在小迎角时增加升力系数,在大迎角时增加升阻比。     

基于Isight的倾转旋翼飞行器前飞状态翼型优化

倾转旋翼飞行器被认为是下一代旋翼类飞行器的主要发展方向,研究其机翼的气动优化设计,对于提高该类飞行器的飞行性能具有重要意义。以NACA2412为原始翼型,首先,采用Hicks-Henne方法进行翼型参数化,并确定设计变量;其次,采用Isight集成翼型生成、网格划分、流场求解等软件,建立翼型自动优化平台;然后,采用基于最优拉丁超立方设计(Opt LHD)和径向基函数(RBF)的代理模型,并用多岛遗传算法(MIGA)进行机翼优化;最后,将优化后的翼型生成三维机翼,进行气动特性计算。优化过程中,对比两种边界条件的优化结果,以证明所用优化方法的有效性;为了减少计算量,使用动量源方法用作用盘代替旋翼。结果表明:根据机翼展向来流速度分布进行翼型优化,在前飞状态下优化后的机翼的升阻比提高了66.03%。     

翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响

为了研究翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响,对带气室的、展弦比为3的不同特征几何参数翼伞模型的流场进行了三维、定常数值模拟。运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes (RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω两方程湍流模型对湍流进行模拟。数值模拟得出的原始翼伞的气动性能参数与试验数据在总趋势上符合很好,不同几何参数翼伞模型计算结果表明:翼伞弧面下反角越大,升力及诱导阻力越小,升阻比变化不大;前缘半径、厚度小的翼伞翼型,阻力更小,升阻比大;前缘切口对翼伞影响区域限于前缘附近,压力分布同干净翼类似,降低了其失速迎角,对升力影响不大,但明显增大阻力。该数值方法可为进一步研究更多不同几何参数的翼伞模型提供参考。     

自适应翼型的气动外形优化设计

 二维翼型自适应的研究是设计自适应机翼的基础。提出了在不同 Ma数、迎角下,用 Powell法优化二维翼型前、后缘襟翼的偏转角,以获得比常规翼型在亚声速时升阻比大而在超声速时阻力系数小的自适应翼型的研究方案。并与原始翼型以及气动双目标 (亚声速时,大升阻比;超声速时,小阻力系数 )的优化翼型进行了比较,证明了自适应翼型比气动双目标的优化翼型气动效率更高。初步探讨了二维翼型前、后缘襟翼的偏转位置对气动效率的影响。     

后缘连续变弯度对跨声速翼型气动特性的影响

针对后缘连续变弯度对跨声速翼型气动特性的影响进行了研究。首先不考虑翼型后缘连续变弯度,基于搭建的优化设计系统对跨声速翼型进行气动减阻优化设计,通过添加不同的约束优化得到两种跨声速翼型:无激波翼型和超临界翼型。然后在这两种翼型的基础上,以后缘偏转角度为设计变量、以阻力系数最小为目标,针对不同的升力系数分别进行优化设计,并根据优化结果深入分析后缘连续变弯度对这两种翼型极曲线特性的影响机理。优化结果表明:无激波翼型与超临界翼型相比,其设计点处的气动特性较好,但鲁棒性较差;升力系数小于设计升力系数时,应用后缘连续变弯度后,无激波翼型的极曲线特性明显提高,减阻最高达到3.9%,而超临界翼型的极曲线特性提高不明显;升力系数大于设计升力系数时,应用后缘连续变弯度后,无激波翼型和超临界翼型的极曲线特性都明显提高,减阻分别达到2.4%~18.1%和1.7%~13.2%。     

大型飞机铰链下垂前缘多段翼型气动优化设计与研究

随着增升装置的发展,新型商用飞机机翼内侧采用了结构简单的铰链下垂前缘,与使用最多的前缘缝翼相比,它具有减阻降噪,提高升阻比等优点。针对大型飞机铰链下垂前缘翼型进行了气动优化设计及数值研究。首先设计出铰链下垂前缘的二维翼型,然后通过iSIGHT平台对翼型进行气动优化,分析铰链下垂前缘翼型的低速气动特性,并且得出最佳翼型,最后与前缘缝翼翼型的气动性能进行对比,说明使用铰链下垂前缘的可行性。     

等离子体热效应对NACA0012翼型增升减阻的研究

基于等离子体热效应机理,在来流速度为34m/s和攻角0～12°内,对NACA0012翼型在等离子体激励下的流场特性进行数值模拟。通过研究等离子体激励的位置和数量对翼型的升阻力特性的影响,得出翼型增升减阻的最佳位置和数量。为保证计算模型的准确性,将未激励的翼型流场参数与NASA实验数据进行对比验证。结果表明:未激励翼型的流场计算参数与实验结果吻合度较高;在等离子体单激励下,最佳减阻位置位于翼型下表面的前缘,最佳增升位置位于翼型下表面的后缘,且二者受攻角的影响较大;在翼型下表面的前缘和后缘同时施加激励时,翼型的减阻比约为20%,最大增升比为52%。