激波串与进气道肩部分离泡相互作用

针对高超声速进气道激波串与肩部分离泡相互作用时的流动振荡问题，在来流马赫数为6的激波风洞中，采用高速纹影拍摄结合壁面动态压力测量，研究了有/无抽吸情况下激波串与肩部分离泡的相互作用过程。结果表明：当激波串前移至肩部附近时，有抽吸进气道也会产生大尺度的分离泡，进而有/无抽吸进气道内的激波串均会与肩部分离泡形成耦合振荡，并造成严重的脉动压力。在激波串的推动下，分离泡能够自由地越过肩部凸拐角，使得其自身的低频振荡特性能够显现。激波串内的压力波动会显著改变分离泡的形态，而分离泡形态的变化又会影响激波串内的压力，两者相互耦合从而维持这种低频振荡。无抽吸进气道具有相同的低频耦合振荡特性；而抽吸缝阻碍了上下游的信息传递，使得有抽吸进气道的分离泡低频振荡显著，而激波串振荡具有一定的宽频特性。经分离激波振荡范围和进气道入口速度无量纲后，有/无抽吸进气道低频耦合振荡的St均处于0.011～0.021，与经典分离泡的低频振荡特性相当。     

二元高超声速进气道激波振荡特性实验

通过在进气道/隔离段模型出口附近设置固定的堵塞楔块提高反压,并采用高速纹影拍摄同步壁面压强测量的手段,在马赫数5.9的激波风洞中研究了二元高超声速进气道在不同堵塞度下的流动特征。研究结果表明,在较低的堵塞度下,进气道仍然能够保持起动状态,而在较高的堵塞度下,进气道出现激波振荡。上游产生的压缩波/激波在节流段的反射是出现激波振荡的重要原因之一。随着堵塞度的增加,激波振荡的频率有所升高。     

周期节流扰动下激波串振荡流动的数值模拟

基于动网格方法数值模拟并分析来流马赫数为6，二元进气道/隔离段构型在频率为50～500 Hz周期节流下的激波串振荡流动。结果表明：当节流比在0.2～0.32范围内周期变化时，隔离段出现与节流频率相同的激波串振荡现象。节流频率会影响激波串振荡幅度和壁面压强波动特性。50 Hz与100 Hz工况的激波串流向振幅相近，100～500 Hz范围内随频率增加，流向振幅从15.5 mm减小至10.8 mm。壁面压强随频率的变化规律更加复杂，以凹腔中部为界，其上游壁面压强时均值、均方差峰值整体随频率增加而降低，其中50 Hz工况唇口侧壁面压强均方差峰值可达21倍来流静压，但其下游壁面压强无明显规律。分析表明节流频率对激波串振荡的影响与节流扰动的传播时间相关，工程设计中需综合考虑构型与反压参数对激波串振荡的影响。     

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

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

隔离段抽吸引起的激波迟滞现象研究

为了研究壁面抽吸条件下隔离段流动迟滞现象,采用数值模拟和理论分析相结合的方法,模拟反压升高再降低过程,对隔离段的激波形态进行了研究。基于Zhukoski提出的中等雷诺数下(3×10~4Re_δ1.2×10~6)分离区压力与马赫数的量化关系,发展了无控制措施条件下激波串首道激波的理论模型,发现来流马赫数大于2.0时激波串首道激波反射类型为规则反射,且不会出现激波反射迟滞现象。而壁面抽吸使首道激波固定在抽吸缝位置,导致激波串首道激波强度随反压升高不断增强,边界层分离角和激波角不断增大,从而进入Von Neumann准则的双解区甚至马赫反射区,在升高及降低反压的过程中隔离段出现流动迟滞现象。研究结果进一步揭示了壁面抽吸引起的流动迟滞现象不仅包含常规RR?MR激波反射迟滞,而且包含了一种新的迟滞现象——边界层分离迟滞。     

弯曲扩张段内激波串自激振荡特性

为了研究一种典型超声速进气道弯曲扩张段内流场结构及激波串自激振荡特性,采用数值模拟方法分析了扩张段出口反压对激波串振荡的影响规律及气流分离和激波串振荡的关系。结果表明:随着反压增大,激波串的大幅振荡和小幅振荡交错出现。当激波串头激波与背景激波相交时,激波串出现高频小幅振荡,振动频率大于900Hz;当激波串前缘点的位置处在扩张段上壁面背景激波反射点时,激波串出现低频大幅振荡,振动频率为200~500Hz;而处在扩张段下壁面背景激波反射点时,不会出现低频大幅振荡。激波串振荡频率与扩张段内部分离包变化频率相近,壁面形状和背景激波的逆压力梯度是激波串大幅自激振荡的影响因素。     

壁温对隔离段激波串振荡的影响

为了探索存在换热条件下超燃冲压发动机隔离段内部的复杂流动机理，采用2阶精度的非定常数值模拟程序研究了壁温对隔离段内部激波串振荡特性的影响。结果表明:压力振荡曲线与激波串前缘位置之间存在紧密关系。此外，随着壁温升高，隔离段壁面压力振荡的频率减小、振幅增加。研究还发现，隔离段压力振荡标准差最大值位于激波串前缘激波及其与上下壁面边界层相互作的位置，而湍动能最大值则分布在分离区后侧剪切层区域内。      

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

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

一种高超声速二元混压式进气道的研究

 针对飞行马赫数为6.00的二元进气道模型开展了高焓脉冲风洞试验研究,分析了进气道在不设置反压和设置反压两种情况下的激波结构、内通道皮托压分布及隔离段出口的性能,并结合数值仿真分析了通道内的流场特性。研究结果表明:在无反压情况下,进气道内通道激波反射明显,靠近下壁面的皮托压值均低于其他测点,在隔离段出口截面,靠近侧壁皮托压有所降低;在一定反压条件下,结尾激波系上传至隔离段内,结尾激波位置不对称;堵塞度为62%的反压条件下,结尾激波系位于喉道位置,隔离段出口截面下半部分已经是亚声速流动;在来流马赫数Ma=6.07,迎角α=4.5°无反压情况下,隔离段出口总压恢复系数为0.477,平均马赫数为2.72,增压比为44,流量系数为0.81,表明进气道性能良好。     

边界层抽吸和脉动反压作用下隔离段内流动特性研究

为了解上游边界层抽吸控制和下游周期脉动反压作用下隔离段内流动特性,采用非定常数值模拟和理论分析相结合的方法,对来流Ma=2情况下的隔离段内激波串动态演化特性、激波串形态结构变化以及激波串演化迟滞现象进行了研究。结果表明,在脉动反压和边界层抽吸作用下,激波串在上游抽吸狭缝与下游隔离段出口之间周期振荡,振荡频率与脉动反压一致。在振荡过程中,首道激波串形态在规则反射与马赫反射以及马赫反射与弧形激波(包含正激波)之间相互转换。边界层抽吸将激波串固定在抽吸狭缝位置,有效提高了隔离段抗反压性能,脉动频率越大,可承受的瞬态反压峰值越大。在一个振荡周期内,激波串向上移动速度较向下移动更快,且在上下移动过程中形态变化存在迟滞现象。     

反压诱导方式对超燃冲压发动机隔离段流动特性影响

为了探究传统机械限流方式同燃烧释热对隔离段流动特性影响的异同,采用实验结合数值计算的方法研究了限流方式与实际燃烧诱导压升的差异.使用氢燃料进行不同当量比的燃烧实验,并使用燃烧室出口安装楔块的方式进行限流实验,对比了来流条件相同,且隔离段出口压比相同时两种实验下的压力分布.使用经过验证的数值方法模拟了不同隔离段出口压力下燃烧状态及对应的限流状态,对比了两种状态下隔离段流场细节.实验结果表明:激波链即将进入隔离段时,两者的压力分布大致相同;激波链进入隔离段后,隔离段出口压比2.3,两者的壁面压力分布有明显差别.此状态下的计算结果表明:燃烧状态下隔离段内分离区首先出现于下壁面,激波链向上偏折;而限流状态分离区出现于上壁面,两者的流场会有一定差异.反压继续增大的计算结果表明:隔离段出口压比达3.0时,两者的隔离段内流场差别逐渐减小并最终趋于一致.燃烧反压场与限流实验的模拟反压场一致时,才能直接采用限流实验的结果评估反压对进气道的扰动风险.      

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

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

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

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

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

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

一种CARET进气道喉部附面层抽吸数值研究

以一种CARET(后掠双斜切双压缩面)进气道为研究对象,设计喉道附面层抽吸槽以控制流动分离。采用CFD数值计算软件对进气道在设计点工况下（马赫数为2.0)下内、外流场进行计算,以总压恢复系数和进气道出口总压畸变为评价指标分析不同抽吸方案的效果。结果表明:喉道附面层抽吸能够稳定结尾正激波,削弱激波/附面层干扰,抑制流动分离,显著改善流场,提高总压恢复系数,减小出口畸变;喉道段抽吸槽位置靠前能够明显降低出口畸变;随着抽吸量的增大,附面层抽吸对进气道内特性性能提升的贡献越来越小。      

来流总温对双模态燃烧室模态转换边界的影响

针对煤油燃料双模态超声速燃烧室，开展了来流总温对燃烧室模态转换边界影响的试验研究。试验采用甲烷燃烧加热直连式试验系统，隔离段进口马赫数保持2.0不变，总压为1.05 MPa，来流总温分别为885、1 085、1 285 K。试验中采集了燃烧室沿程壁面压力，并采用一维分析方法得到了燃烧室的工作模态。试验结果表明:来流总温不同时，燃烧室壁面峰值压力位置相同，同时压力峰值与隔离段壁面压力分布和激波串起始位置存在一一对应关系；来流总温上升导致燃烧室超燃-亚燃模态转换时的当量油气比上升；在燃烧室当量油气比不变的条件下，来流总温上升能够导致燃烧室壁面压力下降，隔离段内激波串长度缩短。      

Unsteady supercritical/critical dual flowpath inlet flow and its control methods

The characteristics of unsteady flow in a dual-flowpath inlet, which was designed for a Turbine Based Combined Cycle(TBCC) propulsion system, and the control methods of unsteady flow were investigated experimentally and numerically. It was characterized by large-amplitude pressure oscillations and traveling shock waves. As the inlet operated in supercritical condition,namely the terminal shock located in the throat, the shock oscillated, and the period of oscillation was about 50 ms, while the amplitude was 6 mm. The shock oscillation was caused by separation in the diffuser. This shock oscillation can be controlled by extending the length of diffuser which reduces pressure gradient along the flowpath. As the inlet operated in critical condition, namely the terminal shock located at the shoulder of the third compression ramp, the shock oscillated,and the period of oscillation was about 7.5 ms, while the amplitude was 12 mm. At this condition,the shock oscillation was caused by an incompatible backpressure in the bleed region. It can be controlled by increasing the backpressure of the bleed region.