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

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

翼型后缘厚度对气动性能影响的数值分析

在研究DU系列翼型的空气动力特性时,发现DU97-W-300翼型的升阻力随时间变化。通过对比多种攻角和多种后缘厚度下的升阻力随时间变化曲线,发现该现象与边界层分离有关,而且后缘厚度对边界层分离有重要影响。最后得出结论:翼型后缘厚度存在一个临界值,当后缘厚度大于临界值时,受力不稳定,小于它时受力稳定。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

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

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

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

翼型襟翼受力的稳定性分析

采用数值求解N-S方程的方法,用DES湍流模型求解襟翼高度对风力机专用翼型DU95-W-180气动性能的影响。结果表明,添加一定高度的襟翼后,翼型的受力出现波动现象,而且波动幅度随着襟翼高度的增加而增大,进一步研究引起该现象的原因发现,添加襟翼的翼型其失速旋涡的脱落过程会引起翼型周围压力的周期性变化,进而导致翼型受力出现波动现象。     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

可连续光滑偏转后缘的变弯度翼型气动特性分析

以变弯度翼型为研究对象,计算了其六种外形的绕流流场,分析了不同的连续光滑变形翼型与传统偏转翼型的气动特性,讨论了变形参数对气动特性的影响,研究了气动特性的产生机理;与此同时,以形状记忆聚合物柔性蒙皮和机械结构实现了可连续光滑偏转后缘的变弯度翼型,并在风洞实验中测试了其气动特性。计算和实验结果表明:可连续光滑偏转后缘的变弯度翼型能改进传统主翼-简单襟翼翼型的气动特性和流场分离特性;可变形段范围、转轴位置、后缘偏转角度、后缘高度等变形参数对气动特性具有显著影响。     

多段翼型低速流场与气动特性数值计算研究

在结构化网格求解RANS方程的基础上,将嵌套网格技术应用在低速绕流气动特性计算中,并以NLR-7301两段翼型和某运输机5段翼型在起飞构型下对低速绕流流场进行数值模拟,将气动特性及流场特性的数值模拟数据与对接网格计算以及风洞实验数据进行了比较.结果表明,该气动特性计算方法有着较高的精度和可靠性.     

NPU-WA系列风力机翼型设计与风洞实验

针对兆瓦级大型风力机,研究发展了以具有更优良高雷诺数和高升力气动性能为特点的NPU-WA翼型族,风洞实验表明,该翼型族达到了在高雷诺数、高升力条件下实现高升阻比和外侧翼型对粗糙度不敏感的主要设计要求,为我国自主研发大型风力机提供了可以实际使用的翼型几何数据和雷诺数范围内1.0×106～5.0×106的风洞实验数据。     

扫描模式航管监视雷达风电场杂波检测与抑制

近年来,全球风力发电装机容量呈指数增长。研究表明,风轮机对其附近的航管(ATC)监视雷达会产生严重影响。风轮机杂波的有效检测及抑制,对于保证空中交通安全具有重要意义。首先提出了基于回波谱宽特征的航管监视雷达的风轮机杂波检测方法。针对扫描模式下航管监视雷达频谱分辨率较低的问题,将基于自回归(AR)模型的超分辨率方法和质量中心的概念应用于雷达回波的快速谱宽和谱中心估计算法中,提高谱宽估计的精度。其次针对扫描模式下的风轮机回波数据不是一个完整周期数据的问题,基于缺省数据幅度和相位估计(GAPES)算法实现了扫描模式下风轮机雷达回波缺省数据的估计,而后利用风轮机杂波的周期性抑制风轮机杂波。该算法实现了风轮机杂波的有效检测和抑制,并且其不受限于风轮机与飞机目标在同一个距离单元的情况。仿真结果验证了所提方法的有效性。     

密度比对涡轮叶片表面气膜冷却换热系数的影响

 采用放大的叶片模型,利用大尺寸低速线性叶栅风洞进行试验,测量了涡轮工作叶片表面不同位置处气膜孔下游叶片表面的换热系数,研究了不同吹风比、密度比和雷诺数的影响。风洞试验段由3个叶片组成,其中中间的叶片为试验叶片。试验叶片表面上开有6排气膜孔,其中吸力面1排,前缘区3排,压力面2排。试验结果表明：密度比对叶片表面气膜孔下游换热系数有影响,以往采用空气作为主流及二次流,在低温差下进行试验,所获得的叶片表面气膜孔下游的换热系数在用于涡轮叶片气膜冷却的实际设计时,必须进行修正。     

一种大型风力机叶片的气动优化设计方法

探索研究了一种大型风力机叶片优化设计的方法,采用有限叶片变环量气动计算模型,同时结合了PSO优化算法。以一3MW大型水平轴风力机叶片设计为例,并同传统的方法进行比较论证,显示了本文方法的可行性和优越性。     

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。     

低雷诺数涡轮流动损失控制技术

针对高空、低雷诺数条件下涡轮损失显著增大、效率显著降低的问题，应用理论分析与数值计算相结合的方法，从流动控制原理、设计方法、流场数值计算等多个方面研究了“微球窝”(Dimple)流动控制技术在航空燃气涡轮中应用的问题，分析了“微球窝”表面处理控制流动分离的基本原理，数值计算表明，围绕球窝周围的气流参数具有明显的变化，球窝对叶片表面流动过程的扰动，使得流动分离点向后推移，雷诺数在12208～15611范围的流动情况，分离点分别向后推迟3％～5％的叶片弦长。     

低雷诺数涡轮叶栅损失的实验与数值模拟

发动机涡轮部件在高空低速飞行条件下工作雷诺数降低,其损失显著增大、效率显著降低。应用实验分析与数值模拟相结合的方法,深入认识高空低雷诺数条件下涡轮流动损失的特征和规律,数值计算是基于Jame son中心差分和Runge Kutta时间推进的N S方程计算的有限体积方法。研究表明,随着雷诺数降低,涡轮叶栅流动损失增大,当雷诺数小于42000之后,涡轮叶栅流动损失呈明显增大的趋势。数值计算结果表明在低雷诺数条件下,涡轮叶栅吸力面后部流动产生了分离,这是流动损失增大的主要原因。数值预测的结果与实验测量结果的趋势吻合得相当好。     

涡轮叶片表面气膜冷却效率的实验研究

采用放大的叶片模型，利用大尺寸低速线性叶栅风洞进行实验，测量了涡轮导向叶片表面不同位置单排气膜孔的气膜冷却效率，研究了孔排位置、吹风比及来流雷诺数的影响。风洞实验段由3个叶片组成，其中中间的叶片为试验叶片，由优质木材制成。试验叶片表面上开有15排气膜孔，其中吸力面3排，前缘区6排，压力面6排。实验的参数变化范围是：基于叶片弦长的来流雷诺数250000-450000，吹风比0．5-2．5。结果表明，由于气膜孔排位置的不同，其下游冷却效率受来流雷诺数及吹风比影响的变化趋势也有所不同。     

Dynamics of Wind Generators on Electric Utility Networks

Dynamic interaction of wind-turbine-driven generators on electric utility networks was studied by computer simulation. Nonlinear representations of wind-turbine and various drive train elements and Park equation representations of synchronous and induction generators were implemented. An infinite capacity network was assumed. Time history responses for various system configurations were computed using as the input function severe wind gust data added to cyclical torque variations occurring at turbine blade frequency. Results indicated that severe transient mechanical and electrical stresses can be induced for certain system configurations. Best results were obtained by interposing rate or damped compliant couplings between the wind turbine and a synchronous generator. The induction generator did not appear to require such means. Blade pitch control (or equivalent) was required to limit output above rated wind velocities for wind turbines configured to produce maximum mum specific power. The blade pitch control loop must exhibit high performance to limit transient overshoots. An aerodynamically limited turbine driving an induction generator exhibited good response without the need for blade pitch control, but at the cost of increased turbine rotor diameter. Further work is indicated, taking into account wind-turbine aeroelastic effects, finite capacity networks, and other factors.