混合缩比连续翼型设计方法

对于飞机的结冰风洞试验,采用多控制点的遗传优化算法进行混合缩比翼型设计。以混合翼型与全尺寸机翼在前缘范围的压力系数相同为目标,通过数值计算分析翼型气动特性,并对比多种优化策略,改进优化算法,生成不同迎角下的混合翼型。以铝合金为材料考察了混合翼型蒙皮变形能力,计算得到蒙皮最大应力为500 MPa,超过其抗拉强度。最终采用分段蒙皮的方式实现混合翼型变形功能。     

适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷。方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混合翼型和原始翼型前缘表面的压力系数在不同迎角下能够保持一致,继而保证前缘结冰外形的一致性。设计的混合翼型较原始翼型弦长减小接近40%,在冰风洞中对混合翼型和原始翼型在目标结冰条件下进行试验,对比结果显示,在选取的多个目标迎角下,混合翼型和原始翼型二维截面的结冰外形基本一致,提出的混合翼型设计方法能够有效的模拟原始翼型在不同迎角下的结冰外形。     

结冰风洞过冷大水滴试验中混合翼设计

开展结冰风洞过冷大水滴结冰试验验证了现有混合翼设计准则。结冰风洞过冷大水滴结冰试验条件采用大、小水滴两种喷嘴组合喷雾实现了与冻毛毛雨水滴质量“双峰”分布曲线非常符合的过冷大水滴结冰环境模拟。以0.50 m弦长的NACA0012翼型为原始翼型，基于面向工程、面向适航的混合翼设计准则通过求解Navier-Stokes(NS)方程计算翼面压力分布，以混合翼与原始翼型驻点位置和前缘吸力峰值尽可能重合为目标，在保证前缘附近与原始翼型相同的基础上混合翼的设计弦长缩短50%。在结冰风洞中开展冻毛毛雨环境下的结冰试验验证混合翼设计效果，结果显示混合翼型与原始翼型模型前缘表面的冰形特征基本一致，局部冰形存在少量差异，混合翼型模型表面冰形冰高略偏大。所得结论表明现有混合翼设计准则在冻毛毛雨环境下依旧适用，可供后续开展相关民机型号的过冷大水滴结冰风洞试验参考。     

运输机翼型结冰的计算和实验

采用数值计算和结冰风洞实验相结合的手段,对某运输机机翼剖面缩比翼型的结冰特性进行了研究.结冰实验在气动中心0.3m×0.2m结冰风洞中进行,实验模型为弦长0.18m的层流翼型.翼型的绕流流场通过求解低速黏流的时均Navier-Stokes(N-S)方程得到,采用拉格朗日法计算水滴撞击特性,在此基础上,求解结冰热力学模型...     

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

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

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

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

数值模拟模型尺度与来流条件对实验数据的影响

采用商业计算软件CFD-Fastran,数值模拟了高超声速一体化飞行器在真实飞行以及风洞实验来流条件下的流场,通过对壁面压力与摩擦力进行积分,得到全机的阻力系数、升力系数以及俯仰力矩系数。对不同来流条件以及模型缩比尺度的气动数据相关性进行研究,其中风洞来流条件的马赫数与静压与真实飞行条件相同,但总焓以及模型尺度不同。数值模拟结果揭示了风洞来流条件以及模型缩比尺度对飞行器整体气动性能的影响,表明发动机内流道摩阻系数计算结果的不同是发动机通流状态下全尺寸模型与缩比模型的阻力系数不同的主要原因。该结果为风洞实验数据的修正提供了参考。     

无尾布局翼型的DISC设计研究

采用DISC方法耦合N-S方程求解程序研究了无尾布局翼型的设计。通过NLF1015,RAE2822翼型的反设计表明,该方法可用于亚、跨声速翼型设计,具有简单、计算量小的优点。完成了翼型改进和新翼型设计,对某BWB飞机内翼段翼型进行了改进设计,采用人工干预对现有翼型压力分布进行局部修正,达到快速削弱激波、改善翼型气动性能的目的;针对小展弦比翼身融合布局内翼段新翼型设计,采用基于多目标多约束压力分布优化的设计思想,得到满足气动和几何约束的目标压力分布。设计结果表明,数值优化目标压力分布方法较易满足总体与气动的设计目标和约束。     

基于粗糙度敏感性研究的风力机专用翼型设计

对几个具有代表性的风力机常用翼型的几何特性和气动特性进行分析研究,详细探讨了各种翼型的前缘粗糙度敏感特性情况。基于一种新的翼型几何形状表征形式,以设计攻角工况下光滑条件和粗糙条件下的升阻比加权值为目标,优化设计得到了一种相对厚度为24.7%的风力机专用新翼型。对新翼型在主要工作攻角范围的升力特性、阻力特性及升阻比特性进行了研究,并详细分析了新翼型在不同工况下的气动特性情况,研究结果表明:新翼型具有低前缘粗糙度敏感性和低湍流敏感性;在设计和非设计运行工况下都具有很好的气动性能;雷诺数的升高和三维旋转效应对风轮叶片新翼型的气动特性具有促进作用。     

结冰对翼型流场影响的研究

本文以NACA23012翼型为研究对象,采用CFD软件模拟了结冰对翼型流场的影响。结果表明：结冰对翼型的气动性能有着明显的影响,尤其是在大攻角时,升力系数显著下降,而阻力系数则明显增加;对于光滑翼型和结冰翼型．计算值与试验值在‘定攻角范围内吻合良好：马赫数对两种翼型的升力系数和阻力系数随攻角的变化影响均比较小,雷诺数也一样。     

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

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

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

高升力翼型的气动优化设计和实验研究

利用遗传算法进行了低雷诺数高升力翼型的气动优化设计,并利用风洞实验检查了设计的正确性。优化后的翼型其气动特性有明显改善,这说明了利用遗传算法进行低雷诺数翼型气动优化的可行性。风洞实验结果验证了优化设计方法的正确性,同时也表明新设计翼型有较高的升力系数和相对大的升阻比,其升阻比提高了10%。     

Robust design of NLF airfoils

 A robust optimization design approach of natural laminar airfoils is developed in this paper. First, the non-uniform rational B-splines (NURBS) free form deformation method based on NURBS basis function is introduced to the airfoil parameterization. Second, aerodynamic characteristics are evaluated by solving Navier-Stokes equations, and the γ-Re_θt transition model coupling with shear-stress transport (SST) turbulent model is introduced to simulate boundary layer transition. A numerical simulation of transition flow around NLF0416 airfoil is conducted to test the code. The comparison between numerical simulation results and wind tunnel test data approves the validity and applicability of the present transition model. Third, the optimization system is set up, which uses the separated particle swarm optimization (SPSO) as search algorithm and combines the Kriging models as surrogate model during optimization. The system is applied to carry out robust design about the uncertainty of lift coefficient and Mach number for NASA NLF-0115 airfoil. The data of optimized airfoil aerodynamic characteristics indicates that the optimized airfoil can maintain laminar flow stably in an uncertain range and has a wider range of low drag.     

Rotor airfoil aerodynamic design method and wind tunnel test verification

Well-designed airfoil is very important for high-performance rotor. This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis, and verified the performance of the optimized airfoil. Using CRA09-A as the baseline rotor airfoil, the CRA09-B optimized rotor airfoil was designed successfully. Combined with the foundation of high-precision rotor airfoil stationary test technology, the CRA09-B and CRA09-A rotor airfoils were tested in the S3MA high-speed wind tunnel of ONERA. In order to correct the aerodynamic data, a single parameter linear wall pressure method is used to consider the tunnel effects. The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable, and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.     

不同热边界对超音速机翼气动热数值分析的影响

飞行器在超音速飞行时受到的气动加热效应给结构强度及热防护设计带来极大影响，且真实状态下的气动热环境需要考虑外流场与结构的耦合及内壁面边界条件的影响。采用S-A湍流模型求解Navier-Stokes方程，通过流场与固体壁面交界处的信息传递，实现外流场与结构场的耦合数值分析。针对三种不同翼型的超音速绕流气动加热进行耦合数值研究，对比翼型内壁面在不同热边界条件下的气动热效应，结果表明：不同翼型具有与气动力相似的的气动热效应；内壁面考虑对流换热的边界条件最接近真实；考虑机翼燃油箱满油时，三维机翼前缘驻点处热流密度最高可达4200w/m2。     

翼型结冰过程数值模拟验算与分析

应用FENSAP-ICE结冰计算软件,对NACA0012翼型进行了流动特性、水滴撞击特性以及冰型生成过程的计算;同时,对结冰前后的翼型进行气动力特性计算对比分析,其中包括升力特性对比、阻力特性对比、流场细节分析以及压力系数分布对比。计算结果表明：翼型前缘结冰后,导致翼型前缘气流提前分离,最大升力系数、失速攻角大幅减小,...     

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

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

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

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