NPU翼型的气动力分析和改进设计

 在飞行器设计中用计算方法设计超临面翼型已完全取代了选用现成翼型的设计方法。为考察已设计出的NPU翼型是否满足飞行器设计要求我们对其进行了全面气动分析,发现这些翼型尚有不足之处,有必要进行改进设计。     

跨声速叶栅抽吸流、激波以及分离流相干效应

以某高负荷、跨声速压气机叶栅为研究对象,应用数值模拟手段探讨通过抽吸控制激波从而控制附面层发展的可行方法。研究结果表明：随着抽吸量的增加吸力面马赫数峰值提高,激波损失增加,同时使得吸力面马赫数峰值点位置后移,附面层分离减弱,分离的减弱所导致的总压恢复系数增加量要远大于激波强度增加所导致的总压恢复系数减小量;抽吸对叶栅性能改善存在一个最佳抽吸量1.2%;在保证叶栅静压压升不变的前提下,相对于未抽吸条件1.2%抽吸使得叶栅总压恢复系数提高10%,扩散因子降低18%,落后角减小5°;通道激波后实施附面层小流量抽吸不能有效改善附面层内部流动参数,当实现前缘入射斜激波投射点位于通道激波上游时,叶表附面层流动得到较大改善。     

电弧放电等离子体激励控制超声速压气机叶栅激波/边界层干扰仿真研究

为研究电弧放电等离子体激励对超声速压气机叶栅激波/边界层干扰的控制作用,建立了模拟等离子体激励作用效果的唯象学模型,进一步以ARL-SL19超声速叶栅为对象,通过数值仿真研究了电弧放电等离子体与叶栅通道内部流动的相互作用及其对叶栅流动损失的影响。结果表明:等离子体唯象学模型能够较好模拟电弧放电等离子体诱导产生冲击波的气动特性。电弧放电等离子体激励对叶栅通道内部流动主要具有三种作用效果:在放电区,注入的热量会产生阻塞效应,增加近壁面气流的流动损失;在激波/边界层相互作用区,能够改变激波系结构,减小激波损失;在尾迹区,冲击波会诱导产生脱落涡。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

基于线性稳定性理论的转捩控制试验

 使用线化稳定性理论研究了某后掠机翼试验模型在风洞试验条件下可能的试验结果。在当前试验条件下,三维边界层中横流驻波的不同波长扰动发展趋势不同,根据理论分析和数值计算给出了最大放大率波长为4.0 mm。通过升华法试验验证了理论分析与数值计算结果。在试验中使用丝网印刷技术在后掠机翼前缘添加粗糙带,成功地引入了理论计算的最大放大率波长。根据预期影响了三维边界层的转捩。     

高超声速边界层流动稳定性实验研究

高超声速条件下边界层转捩相关研究是近年来空气动力学领域的研究热点。通过采用基于纳米粒子示踪平面激光散射(Nano-tracer-based planar laser scattering,NPLS)技术以及温敏漆等测试技术,对高超声速边界层的流动稳定性开展了实验研究。通过NPLS技术,对圆锥边界层中第二模态波精细结构进行了测量,并基于时间相关的测量结果,对第二模态波的波长和频率进行了分析。针对裙锥边界层NPLS结果,计算了特定位置上功率谱空间分布结果,测量得到了高次谐波成分。通过温敏漆和NPLS结果,发现主导三角翼前缘边界层转捩的模态为行进横流模态,分析了该模态的特性,并且与Kulite传感器测量得到的频率进行比较。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

端壁吸/吹气对超声速压气机叶栅流场影响机理的对比研究

附面层吸/吹气是抑制流动分离、提高压气机叶片负荷的有效技术途径。针对超声速压气机叶栅内激波诱导的角区分离,分别采用多种不同的端壁吸/吹气方案对其进行流动控制,旨在探索端壁吸/吹气对激波干涉下角区分离的控制机理,并对比分析端壁吸/吹气对超声速压气机叶栅角区分离的控制效果。结果表明:在激波/端壁附面层干涉下,该超声速压气机叶栅内存在大范围的激波诱导角区分离,角区分离使得该超声速叶栅存在强三维效应,二维叶栅中的单正激波变为"斜激波+正激波"结构,叶中吸力面尾缘开式分离变为闭式分离;端壁吸气可有效抑制该超声速叶栅的角区分离,吸气后近端壁区损失系数大幅降低,最优端壁吸气缝方案的起始点与亚声速压气机叶栅相同,但端壁吸气后叶中的双激波结构变为单正激波结构,叶中流动分离增大;端壁吹气也可有效抑制角区分离,其控制效果略优于端壁吸气,其原因是吹气缝处的静压高于吸气缝,对激波的增强作用弱于端壁吸气;与端壁吸气方案不同的是,最优端壁吹气缝方案的起始点位于叶片前缘。     

高超声速双模块内转式进气道的流动特性研究-Part I： 设计状态

本文通过特征线法设计了基于多道激波+等熵压缩波的基准流场,在此基础上通过流线追踪法设计了双模块内转式进气道。通过数值仿真和风洞试验相结合的方法,获得了内转式进气道的内外流特性。研究结果表明：在内转进气道最大半径对应的角区位置存在大量的边界层堆积,受第二道激波/边界层干扰,在激波根部卷起锥形旋涡;在内转式进气道内部,唇罩激波和管道边界层干扰显著,管道内存在自唇罩指向压缩面的强周向压力梯度,从而诱导管道内边界层均往一处汇聚,卷起大尺度流向涡。仿真和试验结果表明在来流马赫数5.74,攻角0度状态下,进气道气动性能优良,出口总压恢复系数系数达到0.58,最大抗反压为112倍。     

可压紊流边界层掺混方法力学模型的改进

 本文给出了二维可压紊流边界层的一种改进的掺混算法,方法简单,计算机时短,可用于计算激波后的紊流边界层的发展。六个算例的计算结果表明,本方法具有良好的数值精度。同P.D.史密斯的方法相比,结果有明显的改进。     

Steady and unsteady numerical investigation of transitional shock-boundary-layer-interactions on a fan blade

The present numerical study was performed to investigate the impact of both the Reynolds number variation and the used turbulence model to capture the boundary layer development on the characteristic of a BR710 fan blade. A one-equation model of Spalart Allmaras with an optional semi-empirical transition model of Abu-Ghanam Shaw has been applied. The transition model allows the boundary layer development from a laminar to a turbulent behaviour to be taken into account. This is of particular importance at low Reynolds numbers and thus high operating altitudes of airplanes when the transition location moves further downstream to the trailing edge. Therefore the interaction between the shock waves of the transonic fan blade row and the boundary layer leads to a significant change of the characteristic. The steady numerical studies of the 3-D blade passage have been carried out with the commercial 3-D Navier–Stokes solver NUMECA. Additionally an unsteady calculation has been applied to explain the penalty in efficiency on high operating altitudes.     

局部扰动对平板边界层流动稳定性影响的研究

在壁面上构建局部扰动,数值研究不同特征的局部扰动对平板边界层流动稳定性影响的作用机理,进一步探讨壁面扰动的分布、类型、强度大小对三维扰动波在平板边界层流动中的非线性演化规律以及失稳机制的问题.确定什么样的局部扰动能使三维扰动波快速增长或者是什么样的局部扰动能对流动起抑制的作用,同时寻求在三维扰动波向下游发展过程中的传播方向会发生怎样的变化等.对该问题的深入研究,将使人们能深入理解认识平板边界层流动过程中的转捩发生、流动失稳及湍流形成的理论机制.     

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

局部抽吸对边界层流动稳定性的影响和控制研究

在边界层壁面上,设计局部抽吸结构,采用直接数值模拟的方法,获得稳定的三维基本流.在此基础上,研究稳定及最不稳定的二维扰动T-S波的时、空演化机制;进一步探讨了局部抽吸的形式、强度大小及分布结构对二维T-S波的非线性演化影响及其对增长率的贡献大小.结果表明,局部抽吸结构诱导产生的三维基本流是扰动波得以快速增长的一个关键性因素,这是由于平均流剖面的改变及展向速度的出现,增强了流体运动中的不稳定性、扩大了中性曲线的不稳定区域范围.在最不稳定的二维扰动T-S波的非线性演化过程中,由于非线性作用的不断增强,逐渐激发产生出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率相同;同时展向速度的大小对二维扰动波的增长、流动的失稳、流向涡的形成等方面都起着激励的作用.随着时、空的不断发展和非线性作用的迅速加强,正、负相间的流向涡逐渐形成,强度逐渐增大,流向涡的影响区域也在不断扩大,涡的形状逐渐拉伸变长,并出现强的剪切层,流动开始失稳等其它机制;这些结论与文献[4、5]的结果相吻合.     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

跨音弯掠风扇转子叶片的气动弹性剪裁

跨音风扇进口级增压比的进一步提高,主要受两方面的制约,一是效率,二是气动弹性失稳。这两个因素与转子流场中的激波和激波诱导的大尺度分离紧密相关。激波结构又和转子叶片前缘空间曲线形状直接联系。因此,如何精心设计转子叶片前缘空间曲线形状来控制激波结构,就成了当代风扇气动力学的前沿,并导致风扇弯掠空气动力学概念的出现。在具体设计中,此问题可称之为气动与气动弹性综合剪裁。本文简述综合剪裁中的一个组成部分,即气动弹性剪裁问题,并对一个单级风扇转子叶片给出实例。     

数值模拟中的边界处理方法研究

为适应于模拟超声速流场中的一些特殊物理现象,研究了不同物理现象对边界处理的要求。针对激波、旋涡和声波等流场中典型的流动现象,采用常用的外推、特征和缓冲三类边界处理方法分别进行数值模拟。计算结果表明超声速流场边界采用外推技术,只适于计算带有激波的边界。特征方法还能适于旋涡运动的计算。声波的计算只能采用缓冲层技术。     

Design method with controllable velocity direction at throat for inward-turning inlets

In the design of a hypersonic inward-turning inlet by applying the traditional basic flowfield, a reflected shock-wave is formed in the isolator due to the continuous reflection of the cowlreflected shock wave in the basic flow-field, which interacts with the boundary layer to produce a considerable influence on the performance of the inlet. Here, a basic flow-field design method that can control the velocity direction at the throat section is developed, and numerical simulations are conducted to demonstrate the effectiveness of this method. The method presented in this paper can achieve the absorption of the reflected waves at the shoulder of the basic flow-field by adjusting the variation law of the center radius in the basic flow-field, and a smooth transition between the compression surface and the isolator can also be produced. The Mach number and total pressure recovery coefficient of the inlet designed according to this method are 3.00 and 0.657, respectively, at design point(the incoming flow Mach number Ma1= 6.0). The results show that with this method, the inlet can efficiently weaken both the reflection of the shock wave and the interaction between the boundary layer and the reflected shock waves, which improves the aerodynamic performance of the inlet.     

不同转速下跨声速轴流压气机内部流动失稳的机理

以跨声速轴流压气机转子NASA Rotor 67为研究对象,采用数值模拟方法,开展100%、80%及60%转速下跨声速轴流压气机内部流动失稳触发机制的机理研究。数值结果与实验数据的对比分析表明:在3个转速下,数值总性能曲线的变化趋势与实验数据符合一致。通过压气机内部流场的详细分析,得出其基本流动机理。在3个转速下,随着压气机节流,叶顶泄漏涡(TLV)的起始位置逐渐向叶片前缘移动,叶顶泄漏涡也逐渐向相邻叶片压力面偏转,相比近峰值效率点,近失速点时在100%、80%以及60%转速下叶顶泄漏涡的偏转角度分别为3°、6°和9°。在100%和80%转速下,叶顶泄漏涡与激波相互作用所导致的堵塞是触发压气机内部流动失稳的机制,并且在80%转速下,叶顶泄漏涡发生破碎;而在60%转速下,泄漏涡在相邻叶片出现的叶顶前缘溢流(LESF)是触发压气机内部流动失稳的主要机制,叶片吸力面尾缘出现的小尺度附面层气流分离(BLFS)不是主要机制。