旋转冲压压缩转子试验系统变工况数值研究

采用FLUENT软件对旋转冲压压缩转子试验系统通流部分进行了不同出口反压和转速下的数值模拟.结果表明:旋转冲压压缩转子特性线较为陡峭,其稳定工作范围小,相对流量宽度由26.57%降到1.96%.反压大或转速低时,导叶受旋转冲压压缩转子低压区影响减弱,反压和转速对导叶段气流损失影响不大.旋转冲压压缩转子在反压增加时,其内部气流相继经历斜激波、λ激波/激波串、λ激波/正激波、λ激波4个压缩阶段.转速增加时,压缩面区域形成的激波在机匣壁面的入射点逐渐向下游发展,旋转冲压压缩转子压比增大的同时损失增加.出口支板段随反压增加存在最小损失工况,转速增加时,其总压恢复系数由0.944逐渐降低到0.875,气流损失增加.      

考虑激波串的超声速燃烧流场分析模型

发展了一种考虑激波串结构的超声速燃烧流场分析模型,以应用于吸气式高超声速飞行器的设计优化过程。模型通过求解耦合有限速率化学反应的刚性常微分控制方程组来描述燃烧室内点火、燃烧等气动热力现象,采用Billig激波串模型模拟燃烧高压前传过程。采用该模型对"等直段+扩张段"、"后向台阶"和"支板喷射"三种典型构型燃烧室流场进行了计算。结果表明:所建模型正确地反映了超声速燃烧室流动中的物理化学过程;与实验数据比较,模型计算出的壁面压力分布比没有考虑激波串结构的模型符合的更好。     

高背压进气道中内外流耦合作用的大涡模拟

采用大涡模拟研究了出口堵塞比为50.8%的轴对称进气道流动,重点考察了内外流耦合作用下流动的非定常特性。采用国家数值风洞（NNW）工程仿真软件进行数值模拟,得到的壁面平均压力、瞬时压力分布与试验数据符合良好。分析表明：为匹配出口背压,进气道在喉道区域形成激波串结构,使内流道流场分为上游超声速区、中部激波串区以及下游亚声速区;在激波串区,剧烈的逆压梯度产生了分离激波、激波串、分离区及分离剪切层等复杂结构;伴随着激波串运动和边界层大尺度分离,进气道壁面压力出现宽频脉动特征。脉动压力的时空分布表明：内流道脉动压力以扰动波的形式传播,为此建立的声反馈模型能较好地预测亚声速区的主导频率。相关性分析表明：激波串运动受上下游流动耦合作用,其中,频率为St=0.7的运动主要受上游流动影响,频率为St=0.9的运动主要受下游压力扰动波影响。     

不同反压下椭圆形隔离段流场特征与气动性能

在反压条件下,对椭圆形等直隔离段进行了全三维数值仿真研究.结果表明:在对称的几何条件和边界条件以及均匀来流条件下,当来流马赫数较小时,流场与激波串结构基本对称;当来流马赫数较大(大于2)时,流场与激波串结构出现明显不对称现象,并且随着来流马赫数增大,这种不对称现象加剧.椭圆形隔离段抗反压能力随着长短半轴比减小逐渐增强,出口马赫数和温升比基本不受长短半轴比影响.不同来流马赫数下,当隔离段承受最大反压时,出口马赫数变化不大.      

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

设计马赫数为5的一级二元混压式进气道典型状态下气动特性

对一种设计马赫数为5的一级二元混压式进气道再入大气层过程典型状态进行了仿真和风洞试验,得到了该进气道典型状态下的气动特性.结果表明:当来流马赫数高于设计马赫数（为5）时,进气道外压斜激波系提前汇合,与唇罩入射斜激波相互作用,产生了波-波干扰;尽管发生了流动分离,但当来流马赫数为7和6时进气道出口上游气流紊流度分别不超过3.337和3.256,且流道内动态压力信号的功率谱密度呈现白噪声特征,不会对发动机造成结构损伤.因此,对于宽来流马赫数工作范围的进气道来讲,为了提供足够的流量,可以适当降低进气道的设计马赫数.      

双模态燃烧室隔离段激波串位置的测量及控制

隔离段是双模态超燃冲压发动机隔离进气道和燃烧室相互干扰、实现亚燃-超燃双模态的重要部件.在发动机实际工作过程中,燃烧室反压引起的进气道不起动在飞行器加速爬升阶段是需要极力避免和预防的.针对双模态超燃冲压发动机整机模型和燃烧室模型进行了数值模拟研究,分析了激波串前沿位置与隔离段压力分布的关系,在此基础上介绍了三种通过隔离段壁面压力实时测量和监控隔离段激波串前沿位置的方法,并完成了验证实验.结果表明,所使用的计算方法有效可行;隔离段壁面压力分布能够很好地反应隔离段的激波串前沿位置,通过监控隔离段壁面压力分布,控制隔离段激波串前沿位置,能够有效避免和预防燃烧室反压过高引起的进气道不起动问题.     

管内激波串振荡和壁面脉动压力特性

分析了管内激波串振荡引起的壁面压力脉动特性。采用直联风洞实验方式结合动态压力测量技术获得了矩形管道内不同激波串位置下的壁面压力分布数据,对所得压力脉动幅值、功率谱密度曲线等结果进行分析,探讨激波串流动非定常特点与引起压力脉动的机制。     

中心线偏置隔离段内激波串振荡特性

采用非定常数值方法研究了一种中心线偏置隔离段内的激波串振荡特性,分析了隔离段出口反压及管道扩张比对激波串振荡特性的影响规律.结果表明:出口总压最大值对应激波串前缘最小位移,出口总压最小值对应激波串前缘最大位移.在一定大小反压条件下隔离段内出现激波串的自激振荡,在较高反压条件下激波串稳定在管道内部.在来流马赫数为2.0条件下,中心线偏置隔离段较等直隔离段更易出现激波串自激振荡.在出现激波串自激振荡时,在同一反压条件下,大扩张比（如37%）隔离段流场的振荡主频、幅度和出口平均总压均小于小扩张比隔离段.      

顶部间隙对超声速膨胀器流动特性的影响

隔板与机匣之间留有间隙,间隙的存在势必会对超声速膨胀器的内部流场和总体性能产生影响,为了获得超声速膨胀器内部间隙流动的流动细节,采用三维雷诺平均Navier-Stokes方程和标准k-ε湍流模型,就顶部间隙对超声速膨胀器流动特性的影响进行了数值研究。结果表明:膨胀流道出口斜激波导致吸力面压力高于压力面,隔板尾缘附近部分泄漏流体经间隙流回压力面侧;间隙的存在导致吸力面进口及中、后部近下端壁压力上升,而压力面前缘附近压力下降,对比同一隔板位置,间隙高度每增加1%喉部高度,超声速膨胀器隔板载荷系数最高下降2.6%;端壁损失和斜激波损失降低,但产生了泄漏损失,三维流道内总的流动损失增加,膨胀器效率降低,本文研究范围内效率最多下降8.8%;马蹄涡、泄漏涡及二者之间的相互作用是顶部区域的主要涡系结构;前缘附近气流经间隙流到吸力面侧和尾缘附近泄漏流体越过间隙重新流回压力面侧是间隙内气流的主要运动形式。     

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

凹腔前缘角对超声速燃烧室性能的影响

针对带有不同前缘角的凹腔内流动和燃烧过程,分别在冷态和燃烧条件下探讨了前缘角对凹腔内流动损失及阻力特性的影响.研究表明:在壁面垂直喷射的喷口上游和凹腔内部均会形成低速、高温回流区,有利于点火及火焰稳定,燃烧反压通过边界层的亚声速区域上传,形成激波/边界层干扰结构.减小前缘角,可使剪切层分离位置提前,更偏向凹腔内部,导致凹腔后壁面再附激波增强,进而增大了总压损失,降低了总压恢复系数;亦可导致凹腔前、后壁面压差阻力增大,阻力系数上升.进一步认识了凹腔内部流场及稳焰增混机理,进而为优化凹腔结构设计提供依据.      

MHD Flow Control of Oblique Shock Waves Around Ramps in Low-temperature Supersonic Flows

This article is devoted to experimental study on the control of the oblique shock wave around the ramp in a low-temperature supersonic flow by means of the magnetohydrodynamic(MHD) flow control technique. The purpose of the experiments is to take advantage of MHD interaction to weaken the oblique shock wave strength by changing the boundary flow characteristics around the ramp. Plasma columns are generated by pulsed direct current(DC) discharge, the magnetic fields are generated by Nd-Fe-B rare-earth permanent magnets and the oblique shock waves in supersonic flow are generated by the ramp. The Lorentz body force effect of MHD interaction on the plasma-induced airflow velocity is verified through particle image velocimetry(PIV) measurements. The experimental results from the supersonic wind tunnel indicate that the MHD flow control can drastically change the flow characteristics of the airflow around the ramp and decrease the ratio of the Pitot pressure after shock wave to that before it by up to 19. 66%, which leads to the decline in oblique shock wave strength. The oblique shock waves in front of the ramp move upstream by the action of the Lorentz body force. The discharge characteristics are analyzed and the MHD interaction time and consumed energy are determined with the help of the pulsed DC discharge images. The interaction parameter corresponding to the boundary layer velocity can reach 1. 3 from the momentum conservation equation. The velocity of the plasma column in the magnetic field is much faster than that in the absence of magnetic field force. The plasma can strike the neutral gas molecules to transfer momentum and accelerate the flow around the ramp.     

变截面隔离段流场数值分析

对一种由矩形转变为圆形的变截面隔离段流场进行了数值模拟,并将其与圆形和矩形等截面隔离段内流场进行了对比。在相同来流条件和不同进出口压比条件下,分析了隔离段内激波串区域流动特性、激波串起始位置、激波串长度、出口处流量平均参数以及隔离段出口流场流动均匀性等隔离段性能。结果表明:相对于矩形截面隔离段,变截面隔离段中的分离区域有变短、变宽的趋势。这不仅能够有效降低激波串长度和内阻,而且产生了主流和壁面之间的横向流动。最后,对变截面隔离段中激波串长度与压比之间的关系公式进行了修正。     

二维弯曲等截面管道中的激波串特性研究

利用Carroll的试验数据验证了数值方法的有效性之后,对二维等截面弯曲管道中激波串的特性进行了数值模拟试验,研究了管道弯曲对激波串的结构与特征长度、壁面沿程静压分布、出口截面马赫数与总压恢复、反压特性等的影响,研究中考虑了不同的进口马赫数和边界层厚度。结果表明,管道弯曲对流动的对称性有着明显影响,当马赫数较高时(如Ma₀=2.45)合适程度的管道弯曲有利改善直管道已有的流动不对称,使激波串长度缩短。管道弯曲能够有效抑制出口压力变化所导致的出口截面马赫数的大幅波动,考虑到低压比时(出口Ma_e>1)直、弯管道之间总压恢复系数存在明显差距,而当压比较高时(出口Ma_e<1)两者相当接近,因此亚燃发动机的超声速扩压器可适当使用大曲率以缩短管道的轴向长度。另外,鉴于弯曲管道与直管道内激波串长度之间的明显差异,已有的基于直管道的激波串长度经验公式不能很好地适用于弯曲管道。     

Stagnation temperature effect on the conical shock with application for air

The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high temperature, in the context of correcting the perfect gas model. In this case, the specific heat at constant pressure does not remain constant and varies with the increase of temperature. The stagnation temperature becomes an important parameter in the calculation. The mathematical model is presented by the numerical resolution of a system of first-order nonlinear differential equations with three coupled unknowns for initial conditions. The numerical resolution is made by adapting the higher order Runge Kutta method. The parameters through the conical shock can be determined by considering a new model of an oblique shock at high temperature. All isentropic parameters of after the shock flow depend on the deviation of the flow from the transverse direction. The comparison of the results is done with the perfect gas model for low stagnation temperatures, upstream Mach number and cone deviation angle. A calculation of the error is made between our high temperature model and the perfect gas model. The application is made for air.     

旋转爆震发动机喷管非定常流动特性的数值模拟研究

为分析旋转爆震喷管内非定常流动特性,本文基于特征线方法和旋转爆震燃烧室出口时均参数设计了一种塞式喷管模型,通过改变塞锥式喷管的进出口压比（压比分别为15、30和45）,分析了喷管内非定常流场结构与喷管工作性能之间的相互关系,并探究了喷管进出口压比变化对喷管工作性能的影响。研究结果表明：旋转爆震燃烧室下游斜激波进入喷管后,在喷管内形成一道绕喷管壁面螺旋分布的激波,激波结构主要由波前气流角和激波前气流速度决定;爆震喷管局部工作状态依据相对螺旋激波位置分别存在自由膨胀状态、欠膨胀状态、欠膨胀与过膨胀过渡状态以及过膨胀状态四种情况,喷管处于第二种工作状态（欠膨胀）时,喷管内存在激波,并随着激波的传播不断增强,在第四种工作状态达到最大（喷管处于过膨胀状态）,且喷管在第二到第四种工作状态内工作性能变差;在低压比（压比15）工作条件下喷管性能优于高压比（压比30）,其最佳推力和燃料比冲分别为992.05N和1769.52s,因而旋转爆震喷管设计应选择较高的进出口压比作为设计点,使旋转爆震喷管更多的工作在低压比工作状态,以提高旋转爆震喷管工作效率。     

超燃冲压发动机隔离段非对称来流下激波串受迫振荡流动研究

用非对称马赫1.98的隔离段直连实验研究了不同波形不同频率的脉动反压对隔离段内流动的影响。实验结果表明:隔离段能有效地隔离不同频率不同波形周期性反压脉动对上游进气道的影响;反压脉动以第二特征波速向上游传播,但在激波振荡区域,压力脉动主要受激波振荡的影响;厚附面层一侧的下壁面激波振荡区域内压力脉动前传时以指数规律衰减,衰减的程度随着频率和波形的改变而改变。而薄附面层的上壁面激波振荡区域内压力脉动前传时不是单调下降,而是波动的。     

波转子非定常泄漏流动机理

针对制约波转子性能的泄漏问题,通过提取波转子中与泄漏相关的主要流动现象并建立波转子单通道泄漏模型,对非定常泄漏流动机制进行了详细的数值研究,对本文给出的非定常泄漏损失预测模型进行了数值验证。结果表明:转子通道中存在不同程度的压力波动,波动幅值与间隙宽度有关;连续反射膨胀波、周期性出现的弓形激波及其反射激波是压力波动的根本原因;间隙内部泄漏过程存在3个主要的流动阶段;泄漏过程中通道激波传播速度不变、波后时均压力不变;在一定间隙宽度范围内,激波马赫数、激波静增压比与无量纲间隙宽度均呈线性关系,当间隙宽度从0增大到0.08时,激波马赫数衰减7.3%,激波静增压比衰减10.1%;泄漏流动通过泄漏产生的主膨胀波对激波传播过程施加影响,通道激波衰减本质上是理想激波与主膨胀波叠加效应的结果;泄漏损失预测模型与数值结果吻合良好。