背景波系下的隔离段激波串运动特性及其流动机理研究进展

徐珂靖, 常军涛, 李楠, 鲍文, 于达仁

徐珂靖, 常军涛, 李楠, 鲍文, 于达仁. 背景波系下的隔离段激波串运动特性及其流动机理研究进展[J]. 实验流体力学, 2019, 33(3): 31-42. DOI: 10.11729/syltlx20180196
引用本文: 徐珂靖, 常军涛, 李楠, 鲍文, 于达仁. 背景波系下的隔离段激波串运动特性及其流动机理研究进展[J]. 实验流体力学, 2019, 33(3): 31-42. DOI: 10.11729/syltlx20180196
Xu Kejing, Chang Juntao, Li Nan, Bao Wen, Yu Daren. Recent research progress on motion characteristics and flow mechanism of shock train in an isolator with background waves[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(3): 31-42. DOI: 10.11729/syltlx20180196
Citation: Xu Kejing, Chang Juntao, Li Nan, Bao Wen, Yu Daren. Recent research progress on motion characteristics and flow mechanism of shock train in an isolator with background waves[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(3): 31-42. DOI: 10.11729/syltlx20180196

背景波系下的隔离段激波串运动特性及其流动机理研究进展

基金项目: 

国家自然科学基金项目 51722601

详细信息
    作者简介:

    徐珂靖(1987-), 男, 浙江宁波人, 博士研究生.研究方向:激波串运动特性.通信地址:黑龙江省哈尔滨市南岗区一匡街2号哈工大科学园知源楼(150001).E-mail:648460636@qq.com

    通讯作者:

    常军涛, E-mail:changjuntao@hit.edu.cn

  • 中图分类号: V211.48

Recent research progress on motion characteristics and flow mechanism of shock train in an isolator with background waves

  • 摘要: 对高超声速进气道-隔离段激波串在复杂背景波系下的突跳运动特性及其流动机理的最新研究进展进行了综述,涵盖了背景波系作用下的激波串运动特性、突跳机理和突跳运动特性的数学描述方法,以期对高超声速进气道相关研究工作提供一定的参考。首先,对固定背景波系和变化背景波系下的激波串运动特性和突跳机制进行阐述,指出隔离段壁面压力顺压力梯度和逆压力梯度的交替变化是激波串突跳特性产生的内在物理机制。其次,对背景流场下隔离段激波串突跳运动的触发机理和触发条件进行了讨论。最后,基于对运动特性和突跳机制的认识,尝试给出了背景波系作用下的隔离段激波串运动特性的数学模型,为激波串前缘位置控制提供参考。
    Abstract: The present paper aims to provide a summary report on recent research progress about motion characteristics and flow mechanism of a shock train in a hypersonic inlet-isolator with complex background waves to help the researchers working on hypersonic inlet-isolator easily with their further work. It covers shock train motion characteristics, mechanism of a shock train jumps, and the method for the model of a shock train jumps with complex background waves. At first, the investigations for the motion characteristics of shock train with fixed or variable background waves are described, which point out that the boundary layer separation caused by the alternating favorable or adverse pressure gradient in the isolator is the physical mechanism of the shock train motion characteristics. Followed, the triggering mechanism and condition of shock train jumps in an isolator with background waves are discussed. At last, based on the understanding of the motion characteristics and jump mechanism, the method for the mathematical model of the shock train motion in an isolator with background waves is given to provide a reference for the control of the shock train leading edge.
  • 超燃冲压发动机利用进气道-隔离段的激波增压原理对来流进行压缩。隔离段是位于进气道和燃烧室之间的重要部件,它在宽工况下通过激波串流场结构来平衡隔离段进口与燃烧室入口的巨大压比,隔离了燃烧过程和压缩过程之间的相互干扰,使进气道捕获流量在宽燃烧工况下不受燃烧反压的影响,保证超燃冲压发动机能在宽工况下稳定可靠地工作[1]

    隔离段激波串的长度受多方面因素的共同影响,其中,燃烧室燃烧释热引起的高反压对其形成有重要作用。当隔离段入口为均匀来流,随着当量比的提高,激波串前移,其长度与隔离段出、入口压比呈二次增长关系[2-3]。但对于二元压缩超燃冲压发动机进气道构型,隔离段入口大多存在唇罩入射激波,导致入口流场非均匀,并导致隔离段存在复杂的背景波系。背景波系在隔离段内的反射形成了交替变化的顺、逆压力梯度[4],在流向和横向形成多次参数间断,改变了近壁面低能流的流动特性,影响了激波串的运动状态[5]。最近的研究发现,激波串在背景波系下的前移过程中存在突跳运动特性[6]。激波串突跳运动会影响发动机可靠工作、稳定控制,因此,针对背景波系下的隔离段激波串运动特性、流动机理和数学描述方法等开展研究,总结了前人研究成果,并对该主题的研究进展进行了阐述,最后对该主题未来的研究工作做出展望。

    对于均匀入口条件下的等直管道,Waltrup和Billig[2]通过试验研究发现,激波串长度受多方面因素的影响,包括来流马赫数、上游边界层动量厚度、管道高度和管道出入口压比。对于给定的来流工况和隔离段构型,激波串长度与上下游压比呈单调的二次增长关系[3]。Waltrup和Billig的激波长度经验公式被广泛地应用于隔离段设计等工程应用,但由于该公式是在均匀来流等直管道试验基础上得到,对具有唇罩激波和喉部分离激波入射的进气道-隔离段来说,该公式并不完全适用于描述激波长度与上下游压比的关系。

    对于具有激波入射背景流场的隔离段构型,激波串的前移过程就复杂多了。Wagner等[6]在试验中初步研究了1个马赫数为5的进气道-隔离段构型,发现激波串在该构型中不同位置时的前移速度有显著的差别,如图 1所示。首先,激波串以一个相对较高的速度(35m/s)在隔离段中段区间内(6.0<x/h < 10.0)前移,越过一个背景波系在上壁面的反射点(x/h=6.8,见图 2(a));然后在隔离段前段区间(4.0<x/h < 6.0)的下壁面局部顺压力梯度区内减速至19m/s;最后,当激波串被推出隔离段时,其在唇罩激波入射位置形成的逆压力梯度区内被加速至74m/s。图 1(a)中上半部分给出进气道-隔离段的构型和压力测点位置;下半部给出激波串运动到隔离段各位置时的速度,速度通过相邻压力测点检测到激波引起的压力上升信号的时间间隔得到;图 1(b)给出下壁面在背景波系作用下的压力分布。Wagner[7]随后又继续通过PIV方法对该进气道-隔离段构型的速度分布进行测量,还发现激波反射引起的流动分离对激波串前移过程起到关键作用。图 2(a)展示了隔离段纹影图,其中箭头A标记了激波串前缘,箭头B标记了背景波系在上壁面的1个反射点(x/h=6.8)。图 2(b)展示了流场速度分布图, 在t=2.0ms时刻,下壁面x/h=10.0处激波串前缘的下缘存在明显的流动分离现象,并形成局部的低速回流。在t=4.0ms,激波串前缘的下缘移动至x/h=7.0附近,并在下壁面形成流动分离区域。t=4.0~5.0ms时,在背景波系诱导下,背景波系在上壁面反射点处的流动分离急剧加强(箭头B处,图 2(b)),并形成了图 2(a)中的激波串前缘分离激波。这种激波诱导形成的激波串前缘流动分离加剧的现象也在张航等[8]的CFD仿真计算中被发现。田旭昂等[9]Ma=5.0侧压式进气道试验中也发现激波串的突跳现象与背景波系形成的压强分布有关。

    图  1  激波串在隔离段不同位置前移速度变化[6]
    Fig.  1  Velocities of a shock train moved forward in an isolator[6]
    图  2  隔离段内的速度分布[7]
    Fig.  2  Velocity contour of the flowfield in an isolator[7]

    隔离段入口的非均匀流场和激波入射主要由进气道唇罩激波和喉部分离激波引起。这些波系入射到隔离段内部,在隔离段壁面上发生多次反射,并形成复杂的波系结构。这种波系结构及来流条件和进气道构型有关,作为一种背景流场结构而存在。谭慧俊等[10]在2012年首次提出背景波系的概念,并通过试验对这种复杂背景波系作用下的激波串长度变化过程进行了研究,发现激波串前缘呈现走走停停的非线性特性。黄河峡等[11]发现超声速气流在背景波系的激波和膨胀波作用下经历反复的压缩膨胀过程,在隔离内形成了上升下降的压力梯度。

    背景波系强度,及其形成的局部流动分离强度对激波串过程中的突跳运动起决定作用。Koo和Raman[12]采用与Wagner试验相同的进气道-隔离段构型,使用大涡模拟方法进行计算和对照,复现了激波串前缘在靠近隔离段入口时发生的前移加速现象。Wagner使用的试验构型背景波系较强,其大涡模拟准确预测了激波串突跳过程中的加速运动,但该方法却高估了流动分离的强度,导致激波串前移速度的计算结果比试验结果更高。这是因为大涡模拟的滤波模型在近壁区不是很准,很可能使壁面分离预测过高,导致回流流速比试验高,最后造成激波串前移速度超过试验所得结果。但以上计算结果也间接表明, 激波串急剧前移过程由强背景波系引起的流动分离主导。

    徐珂靖等[13]使用数值模拟的方法研究了隔离段激波串突跳特性,通过分析发现背景波系引起的壁面压力梯度变化对激波串突跳特性起着重要作用。图 3(a)展示了1次激波串突跳过程的数值纹影,图 3(b)以相同的沿隔离段的x坐标给出了突跳前后的壁面压力分布。当激波串前缘从下游逐渐接近隔离段下壁面的某一反射点时,激波串前缘下缘进入一个局部逆压力梯度区域。随后,激波串前缘的流动分离在背景波系反射作用下急剧加强并向上游快速扩展。李楠等[14]比较了不同背景波系下的激波串前移路径,并对背景波系下的激波串前移路径做了数学建模。研究发现,当背景波系较强时,激波串前移过程中会发生突跳;但是,对于背景波系较弱的进气道构型,激波串突跳现象则不太明显。

    图  3  激波串过程突跳纹影序列和壁面压力分布[13]
    Fig.  3  Schlierens and pressure distributions when a rapid movement of shock train occured[13]

    上述研究均采用固定的来流条件,通过提高下游反压的方式使得激波串向前运动,但是对变化背景波系下的激波串运动特性研究比较少。徐珂靖等[15]通过数值模拟对反压固定条件下由背景波系变化驱动激波串前移的过程进行了试验,其中背景波系的变化通过改变迎角实现。随着迎角降低,背景波系在隔离段内发生流向方向的伸展变化。如图 4所示,激波串前缘的下缘向上游运动并靠近背景波系在下壁面的1个反射点(t=40ms);随后,激波串发生突跳并快速越过该反射点(t=42ms)。

    图  4  背景波系伸展过程中激波串突跳现象的数值纹影序列[16]
    Fig.  4  Sequence of schlierens when a rapid movement of shock train occurred[16]

    激波串突跳运动由背景波系流场下的局部逆压力梯度区的移动引起。随着背景波系的结构发生伸展,上游压力降低,激波串缓慢地向上游前移。背景波系形成的局部逆压力梯度区域也在此伸展过程中向下游运动。这就使激波串前缘和局部逆压力梯度区的下游边界有了接触机会。它们一旦接触,由于流动分离会快速在逆压力梯度区里向上游发展,便会形成激波串突跳现象。图 5展示了图 4激波串突跳过程中的壁面压力细节,其中绿色箭头区域标注了局部逆压力梯度区。当t=40ms,激波串前缘移动到x=0.342m处,并接触逆压力梯度区的下游侧的边缘,激波串突跳运动发生。

    图  5  激波串下缘突跳前后的壁面压力分布[15]
    Fig.  5  Pressure distributions when a rapid movement of shock train occured[15]

    前两小节中的研究都是在冷流条件下用各种方法模拟下游反压。在真实的发动机工作条件下,反压是由超燃燃烧室的燃烧反应区释热引起的,并对隔离段激波串运动特性产生重要影响。Dessorness等[16]通过地面试验和数值模拟相结合的方法研究了燃料当量比对燃烧模态的影响,发现随着燃料当量比的增加,燃烧室由超燃模态进入亚燃模态。Denis等[17]通过支板喷射燃料的双模态超燃冲压发动机燃烧室试验,发现模态转换取决于氢燃料的喷注压力。因此,如果不讨论燃烧释热引起的激波串前移过程,则不能全面认识激波串前移运动的流动机理。

    国内外很多研究发现,随着当量比提高到一定阈值,燃烧模态会由超燃模态突变为亚燃模态。Chun等[18]采用试验方法对燃烧模态转换边界做了研究,发现随着来流总压升高或燃烧室扩张角增大,从超燃模态转为亚燃模态所需要的当量比也越高。Le等[19]通过调节燃料当量比来实现模态转换,结果发现随着燃料当量比的增加(从0.08到0.34),燃烧室峰值压力与进口静压之比也逐渐增加。当压力比值超过一定阈值(2.8)时,发动机从超燃模态进入亚燃模态。模态转换时,燃烧室压力分布出现突变,在超燃模态下,峰值压力在物理斜坡之后获得;而在亚燃模态下,峰值压力在物理斜坡之前获得。于达仁等[20]通过机理模型分析了燃烧模态转换过程的突变与滞后等非线性现象,发现燃烧模态转换边界在扰动变量空间分布存在拓扑不变性。

    以上试验均未关注激波串前移过程的运动细节。使用直连试验台,隔离段入口的气流品质较好,因此在隔离段和燃烧室内并没有比较明显的背景波系射入。激波串的突跳运动特性因为试验条件过于理想而被忽略。比如,Kobayashi[21]的直连发动机试验中,发现发动机推力增量与有效当量比呈线性关系。但是,还是有少数燃烧试验,发现了激波串突跳现象。秦斌等[22]使用带有燃烧室的隔离段进行试验,试验隔离段入口存在激波,隔离段内存在背景流场,随着当量比线性上升,激波串前缘(Shock Train Leading Edge, STLE)位置呈现走走停停的非线性特性。在某些区域,激波串发生向前突跳的动作,如图 6所示。

    图  6  激波串前缘位置随当量比变化过程[22]
    Fig.  6  Climbing path of a shock train[22]

    Gnani[23]总结了隔离段激波串前缘的激波附面层交互结构模型。在激波串前缘初始的增压过程中,核心流面积收缩,这种向流场中部传播的压缩扰动形成了激波串前缘的第1道斜激波。这道斜激波在隔离段壁面多次反射,以压缩膨胀交替的形式向下游传播,最后形成激波串结构,并匹配了上下游压力。Tu等[24]给出了高来流马赫数条件下的X型激波串结构(如图 7所示)。区域1的超声速核心流因为分离流(区域3)造成的面积收缩而压力上升,通过一系列激波和膨胀波,减速到亚声速状态。

    图  7  激波串流场结构示意图[24]
    Fig.  7  Schematic of the flowfield with a shock train[24]

    徐珂靖等[25]对激波串前缘压缩结构进行了数值研究,并分析了激波串前缘的波系结构。如图 8所示, 隔离段内因流动分离而增厚的边界层推挤并压缩主流形成激波串前缘最初的分离激波。图 8(a)流线图显示,核心流的气流角度在分离附面层的推挤作用下向内偏折,核心流的流通面积降低,在激波串前缘局部流场中存在1个由分离附面层组成的类似喉道的结构。与进气道外收缩段和内收缩段的固壁压缩过程不同,激波串前缘的收缩结构由分离附面层形成。不过,它也同样形成了等效的主流偏折和压力上升效果,所以,为了更准确地描述这种由附面层推挤而形成的主流收缩结构,可以将其称为“等效喉道”[25]

    图  8  等效喉道压缩结构示意图[25]
    Fig.  8  Sketch of the throat-like structure at shock train head[25]

    Kantrowitz收缩极限[26]被广泛应用于进气道流场状态的分析,用于判断进气道外收缩段和内收缩段的斜激波是否能稳定挂载在压缩斜面前缘。徐珂靖等[25]参考流道面积收缩极限思想,猜想等效喉道可能存在一种最小收缩面积限制,可以作为激波串突跳运动的触发条件。

    激波串前缘的初始激波角、气流转折角和2道初始激波相交后形成的第2道激波的激波角、第2道激波形成的转折角以及第2道激波后方的等效喉道面积θ图 9所示。假设在局部喉道和隔离段初始流通面积之比达到最小收缩比时,分离边界层形成的第1道分离激波是一个斜激波,然后用Kantrowitz收缩极限来分析第2道激波的压缩条件。第1道激波的激波角是个变量,它的值必须大于马赫角;但需要保证波后依然是超声速,这样第2道激波才可以正常产生,然后应用Kantrowitz极限对下游流场进行分析。局部喉道和隔离段初始流通面积之比可以表示为如图 10所示的三维曲面。这个曲面在任意隔离段马赫数下都有一个最小值。于是,任意隔离段马赫数下的触发突跳的最小收缩比如图 11中的实线所示,同时,图中虚线表示等熵压缩至Ma=1的最小收缩比,虚线以下的区间是不可能达到的。通过理论推导,可以猜想,当收缩比从实线以上到越过实线进入实线和虚线之间的区间时,激波串突跳会被触发。

    图  9  激波串前缘初始激波和附面层流场结构示意图[25]
    Fig.  9  Sketch of flowfield at shock train head[25]
    图  10  激波串前缘初始激波角与收缩极限的关系[25]
    Fig.  10  Relationship between angle of initial shock to contraction ratio[25]
    图  11  隔离段马赫数与收缩极限的关系[25]
    Fig.  11  Relationship between Mach number to the minimum value of contraction ratio[25]

    可以观察到,激波串前缘等效喉道的收缩比会在1次突跳前后发生明显的变化,如图 12所示。图 13的黑色实线给出试验结果中等效喉道收缩比的测量值,当该值降到最低点时,激波串突跳发生。随后,通过系数修正,使最小收缩比值(红色虚线)与突跳运动触发时的实际收缩比值尽量靠近,初步改善了该分析方法的鲁棒性,其中修正误差在10%左右。

    图  12  激波串突跳前后纹影图[25]
    Fig.  12  Schlierens of the flowfield when a rapid movement of shock train occurred[25]
    图  13  突跳前后等效喉道收缩比变化[25]
    Fig.  13  Contraction ratio of the throat-like structure when a rapid movement of shock train occurred[25]

    上述分析方法中10%的误差主要来自分析模型的局限性,没有考虑壁面摩擦下初始激波下方附面层的分离条件。于是,Xu等[27]通过超声速自由相互作用理论对等效喉道的最小收缩比分析方法做了改进。通过超声速自由相互作用理论得到第1道激波能形成的最大压比,然后通过激波极线图分析确定激波角和气流折转角。第2道激波的参数也可以通过激波极线图分析得到解析解,然后这个等效喉道的收缩比便可以通过几何关系求解,如图 14(a)所示。另一方面,当第1道激波的参数确定时,参考Kantrowitz极限思想,假设第2道激波形成了正激波,然后等熵收缩加速到局部喉道恢复Ma=1;这个收缩比值可以用于分析激波串是否能从突跳状态恢复至缓慢前移,如图 14(b)所示。可以推理得到,当收缩比小于0.43时,激波串突跳被触发;当收缩比大于0.48时,激波串从突跳状态恢复至缓慢前移。与试验数据对照发现,该推理结果不需要系数修正就能与试验结果相符合,如表 1所示。

    表  1  激波串运动过程的激波串前缘收缩比测量值[27]
    Table  1  Measured values of contraction ratio of the throat-like structure at shock train head[27]
    Data No. Motion state Actual CR Explain Correct
    1 Rapid motion 0.41 < 0.43 Yes
    2 Back to slow 0.50 >0.48 Yes
    3 Rapid motion 0.41 < 0.43 Yes
    4 Back to slow 0.56 >0.48 Yes
    5 Rapid motion 0.41 < 0.43 Yes
    6 Back to slow 0.50 >0.48 Yes
    7 Rapid motion 0.41 < 0.43 Yes
    8 Back to slow 0.51 >0.48 Yes
    9 Rapid motion 0.41 < 0.43 Yes
    下载: 导出CSV 
    | 显示表格
    图  14  激波串前缘最小收缩极限分析示意图[27]
    Fig.  14  Sketch of the analysis method for the minimum value of the contraction ratio[27]

    文献[27]对等效喉道最小收缩比的分析方法是在模型简化为二维压缩的前提下得到的,只适用于具有较强背景波系且带有等直隔离段的二元压缩进气道。这种进气道是一种高超声速推进系统中常见的进气道。研究中发现的各种现象和由此推出的结论揭示的是二维等直管道内激波-附面层交互引起的激波串流动分离机理。其误差主要来自激波-附面层交互的三维效应。其中,激波串向上游急剧前移的触发条件由激波串前缘等效收缩比分析得出;其中,最小收缩比又是由附面层可承受最大压比进行推导,与当地马赫数和壁面摩擦系数有关[27]

    通过第2.2节超声速自由相互作用理论分析,激波串前缘上游的壁面摩擦系数对激波串前缘位置可承受的最大反压起决定作用,可以此计算等效喉道的最小收缩比[27]。然而壁面摩擦系数又与当地近壁面的附面层流动状态相关,并受入射激波的影响。

    近年来,从三维角度研究隔离段内背景激波入射诱导的流动分离现象,发现了更多激波附面层相互作用的细节。对于矩形等直隔离段,由于激波附面层干涉,壁面对称面附近形成中部分离涡,而在侧壁与底面的交汇处形成角涡。黄河峡等[28]从三维角度研究了矩形隔离段内的角涡流动。唇罩激波与管道侧壁面的冲击与交互作用形成了扫掠激波诱导涡(Glancing Shock Induced Vortex, GSIV),侧壁的旋涡沿隔离段横向下扫,下扫气流与唇罩激波形成的分离包相互作用,使侧壁低能气流窜入角部分离区域,形成角涡流动,如图 15所示。Funderburk等[29]研究了背景波系反射处的非定常流动分离,发现激波附面层干涉形成的角涡与中部分离涡流动现象,如图 16所示。

    图  15  角涡流线图[28]
    Fig.  15  Sketch of streamlines of corner vortex[28]
    图  16  角涡和中部分离涡的流场可视化[29]
    Fig.  16  Visualization of corner vortex and central vortex[29]

    Pirozzoli等[30]在直接数值模拟中研究了湍流结构在由激波诱导的分离流中的发展过程。从分离点开始发展,湍流结构最初完全与壁面分离,然后在下游与壁面再附,其流场结构如图 17(a)所示,壁面摩擦系数如图 17(b)所示,灰色区域代表壁面摩擦系数为负值。结果表明,流动分离区域具有高度的三维特征,并在激波附面层相互作用区域存在回流的散点。

    图  17  激波入射引起的激波附面层交互与壁面摩擦系数分布[30]
    Fig.  17  Shock-boundary layer interactions and friction coefficient distribution induced by incident shock[30]

    对于实际的进气道,因为侧壁的存在,流场呈方形截面。Bermejo-Moreno等[31]通过大涡模拟方法研究了激波入射诱导的激波附面层相互作用现象,在激波附面层交互区域,流场结构存在逆流,并且分离包主要集中在壁面的中部,如图 18所示。由于壁面摩擦系数受近壁面流速影响,通过超声速自由相互作用理论可知,这些区域无法承受激波串前缘的压力梯度,因此,激波串将会在这些区域发生突跳运动。

    图  18  三维效应下的流场速度分布[31]
    Fig.  18  Velocity contours of flowfield with sidewall effect[31]

    激波串对扰动的响应过程非常复杂,并且其非线性特征具有累加性及局部性[32]。Bruce和Babinsky[33]将激波动态解释为:激波通过运动对背压变化做出响应,所以其相对运动马赫数可以与波后压升相匹配。在崔涛等[34]的研究中,激波的动态特性也同样和其相对运动马赫数建立了关系。除此之外,激波的动态、分离激波下游亚声速区域体积和背压之间相互耦合,具有很强的相互依赖性。

    苏伟仪等[35]通过数值仿真也发现激波串在前移过程中会呈现振荡特性。谭慧俊等[10]通过试验发现,当激波串非常靠近背景激波的反射点附近时,上游湍流边界层内或者激波串下游的一个小扰动就可能诱导激波串前移,使得激波串内的分离区和背景激波诱导形成的分离区合并,导致分离包显著增加,激波串快速前移。激波串前移过程中分离区尺度增加,激波串内的逆压梯度不断降低,使得分离区开始逐渐趋于再附,激波串开始后移。激波串在后移过程中,内部的流体不断积蓄,逆压梯度逐渐增加,当激波串退回到其初始位置时,边界层再次分离,迫使激波串前移,至此完成一个完整的振荡周期。可见,在激波串的周期性振荡过程中,背景激波入射点附近的边界层起到了扰动放大器的作用。

    为进一步证实激波串的振荡机理,需要进行额外的试验。如上所述,这种剧烈的振荡是由背景波系入射点附近边界层分离引起的扰动、扰动向下游的传播过程以及反压之间的强烈相互作用引起的。减弱这部分相互作用,使压力恢复到稳定位置,可以减小激波串及其内部压力的剧烈波动。考虑到在有入射激波的情况时,激波附面层的作用不可避免,因此可以通过调节背压来减弱背压与上游扰动之间的相互作用。当上游出现扰动时反压上升足够快,激波串区域内的压力分布可以立即被推至稳定状态。同时,激波串与激波附面层作用区域之间的相互作用时间也会缩短。

    李楠等[36]对快速变化背压下的激波串运动进行了研究,试验段构型和压力测点布置如图 19所示,其中TC1~10为上壁面压力传感器,TB1~10为下壁面压力传感器。在文献[31]研究中, 工况D、E和F中下游堵块的旋转时间分别为0.1338、0.0808和0.0675s。不同条件下背压对输入信号的响应如图 20所示,然而由于执行机构特性及节流效应的限制,背压上升的时间稍长。从背压变化时序中可以辨别出,工况C、D、E和F中的背压实际上升时间分别约为5.30、0.20、0.12和0.10s。与工况C类似,激波串首先被推入试验段,然后执行试验时序。为了得到相同的背压变化范围,挡板的初始位置和最终位置均与工况C相同。

    图  19  试验段构型和压力测点布置[36]
    Fig.  19  Configuration of test sections and the pressure measuring points[36]
    图  20  背压和挡板角度变化时序:(a)工况D,(b)工况E,(c)工况F。时间根据试验时长进行了无量化处理[36]
    Fig.  20  Time histories of the backpressure ratio and the flap angle: (a) case D, (b) case E, (c) case F. The time is non-dimensionalized by the total duration[36]

    对于工况D,典型的壁面压力时序如图 21(a)所示。从图中可以发现,TC2 (TB2)和TC6 (TB6)压力时序图中出现突然的上升和下降,表明这种剧烈振荡仍然存在。在工况E中,背压上升的速度变快,TC2(TB2)和TC6(TB6)压力的波动变弱。在工况F中,背压上升速率进一步提高,压力波动变得不再明显,特别是在TC2(TB2)处,工况C和F中激波串前缘轨迹如图 22所示。由于来流条件相同,两种情况下激波串前缘的初始和终止位置基本一致。在工况C中,当激波串前缘通过时,激波串振荡是显著的。随后,当反压继续上升时,激波串前缘稳定在激波附面层作用区域上游。然而对于工况F,由于背压变化速率增加,扰动作用的时间变短,并且激波串内部的压力被快速推至稳定位置,使得振荡明显变小。虽然TC10的压力信号在结束时出现了大幅上升,但激波串运动过程中沿程的压力变化要平缓的多,尤其是在TB2位置, 甚至没有出现压力的突变。

    图  21  试验过程中典型压力时序[36]
    Fig.  21  Comparisons of typical pressure histories during the test conditions[36]
    图  22  工况C和F中激波串前缘轨迹对比:(a)上壁面,(b)下壁面。时间根据试验时长进行了无量化处理[36]
    Fig.  22  Comparisons of the trajectories of the shock train leading edge with case C and case F conditions at the top wall (a) and bottom wall (b). The time is non-dimensionalized by the total duration[36]

    激波串内部的压力分布很难通过解析表达式获得。一方面,关系式要能很好地描述来流静压通过激波串之后的压升;另一方面,关系式要能描述升高到此压力所需要的长度。因此,大部分是通过经验公式、半经验公式(如Billig关系式及其修正公式)近似描述激波串内部的压升。另外还需要寻找合适的表达式来描述第1道分离激波及整个激波串对上下游扰动的响应,由此可见,激波串的动力学模型既要能描述其时间尺度的特性又要保证其空间尺度上的分布参数特性。但是,目前激波串的数学描述一直停留在以数据建模的方式,这种形式的模型不能很好地反应出激波串运动与各个因素之间的联系。

    李楠等[37]对激波串的一般动力学建模开展了研究。和亚燃冲压发动机内的正激波相似,激波串会对上游、下游的扰动做出响应。图 23所示的物理过程可以很好地描述各个因素对激波串运动的影响。当下游背压pb发生变化时(扰动的产生),扰动波会在激波串内的亚声速区域向上游传播,传播时间为τ;当扰动传播至激波串的第1道分离激波时,其波后压力会发生改变。和正激波类似,当波后压力降低时,激波会向下游移动;当波后压力增加时,激波会向上游移动。当第1道激波发生位移时,其后部容积发生相应变化,其内部压力分布也发生变化,进而导致第1道分离激波和背压之间存在着相互耦合作用。

    图  23  激波串运动物理过程简图[37]
    Fig.  23  Physical process of the shock train movement[37]

    激波串的形成起始于边界层内大尺度分离区的形成。边界层内充斥着低能流体,如果其不能抑制较大的逆压力梯度,其中的流动将会分离,随着分离流向上游移动,激波串也向前移动。因此,激波串的动力学特征与近壁面分离流的流动特性密切相关。激波串被定义为承受一定压力梯度的分离流以及核心流的流体结构。这里激波串的位置被定义为其第1道分离激波的前缘,因此激波串的整体运动特性被视为第1道分离激波的运动特性,并且分离激波角度的非定常特性被忽略。考虑到第1道分离激波对附面层流动严重的依赖,平衡状态下特定位置处分离激波的波前、波后压力可以通过自由干涉理论建立起联系,如下式所示[37]

    (1)

    根据上述的关系式可知,附面层内流动所能承受的最大压力p2和来流动压q1,壁面静压p1,来流马赫数Ma1和壁面摩擦系数有关,其中下标“1”和“2”分别表示分离激波上游和下游的条件。为了简化计算,式(1)中(Ma12-1)0.5Ma1q1 =0.5γp1Ma12F函数和雷诺数有关,化简得到[37]

    (2)

    结合平衡流型、扰动波传播及容积效应,一般性激波串动力学模型被构建如下[37]

    (3)

    其中,

    (4)

    通过地面试验在Ma=1.85及2.72条件下对模型进行定性验证,其结果如图 24所示。考虑建模过程中的假设,尽管存在偏差(可被接受),此低阶模型可以定性地对激波串稳态特性进行描述。为进一步验证激波串的动态响应,展开了额外的非定常数值模拟研究,模拟结果和模型预测结果对比如图 25所示。从图中可以明显发现,激波串对台阶信号的响应接近数值模拟的结果。考虑到扰动波的传播速度是通过壁面滞止温度计算而得,模型预测的响应滞后要小于数值模拟结果。并且由于数值模拟中考虑了附面层内的粘性影响,阻尼弱化了分离激波的振荡特性。通过数值模拟和模型预测的对比可以发现,低阶模型可以定性地对激波串的动态特性进行描述,并且具有的偏差可被接受,可适用于二元压缩的高超进气道的等直隔离段。这些偏差来源于三方面:首先,分离区内的稳态压力分布关系在激波串充分发展的情况下具有良好的适用性;其次,激波串动力学模型严格依赖于当地流场参数,并且文中的参数是通过数值计算而得,试验结果与数值计算的偏差也构成了最终结果的差异;最后,通过一维分析方法去近似三维复杂流动,并且简化模型所提出的假设条件也造成结果必然的偏差。

    图  24  无激波入射下低阶模型激波串前缘预估与试验结果比较[37]
    Fig.  24  Trajectories of the shock train leading edge estimated by the low-order model compared with the experimental values[37]
    图  25  模型预测激波串动态特性与数值模拟结果对比,背压压比变化为:(a) 2.83~3.05, (b) 3.05~3.39, (c) 3.39~3.62[37]
    Fig.  25  Dynamic features of the shock train compared with the CFD results, where backpressure ratio is from (a) 2.83 to 3.05, (b) 3.05 to 3.39, (c) 3.39 to 3.62[37]

    由分析可知,激波串前缘的运动受到隔离段沿程流场参数的影响。在激波附面层作用区域,附面层饱受一定的逆压力梯度进而导致其对背压引起的逆压力梯度抵抗能力下降。这种影响可以体现在隔离段沿程壁面摩擦系数Cf上,因此当以壁面静压及摩擦系数为模型输入时,其预测结果如图 26所示。

    图  26  有激波入射下低阶模型激波串前缘预估与试验结果比较[37]
    Fig.  26  Trajectories of the shock train leading edge estimated by the low-order model compared with the experimental values under incident shocks[37]

    背景波系下的激波串运动特性对超燃冲压发动机进气系统保护控制非常重要,是一个必须面对的重要问题。对背景波系下的激波串突跳运动特性和流动机理的最新研究进展进行了阐述,给出了国内外对背景波系作用下的隔离段激波串运动特性、突跳机理和突跳运动特性的数学描述方法等方面的研究进展。

    当前的研究现状对背景波系下的激波串突跳运动、触发机制和突跳运动数学模型已有较深入的认识,但是对激波串突跳运动抑制的流动控制手段还很少,未来的工作包括:

    (1) 宽工况下进气道-隔离段激波串突跳位置预测方法。

    (2) 宽工况下抑制激波串突跳运动特性的分布式流动控制方法研究。

    (3) 考虑背景波系的隔离段-燃烧室耦合作用下的激波串突跳特性、机理和抑制方法。

  • 图  1   激波串在隔离段不同位置前移速度变化[6]

    Fig.  1   Velocities of a shock train moved forward in an isolator[6]

    图  2   隔离段内的速度分布[7]

    Fig.  2   Velocity contour of the flowfield in an isolator[7]

    图  3   激波串过程突跳纹影序列和壁面压力分布[13]

    Fig.  3   Schlierens and pressure distributions when a rapid movement of shock train occured[13]

    图  4   背景波系伸展过程中激波串突跳现象的数值纹影序列[16]

    Fig.  4   Sequence of schlierens when a rapid movement of shock train occurred[16]

    图  5   激波串下缘突跳前后的壁面压力分布[15]

    Fig.  5   Pressure distributions when a rapid movement of shock train occured[15]

    图  6   激波串前缘位置随当量比变化过程[22]

    Fig.  6   Climbing path of a shock train[22]

    图  7   激波串流场结构示意图[24]

    Fig.  7   Schematic of the flowfield with a shock train[24]

    图  8   等效喉道压缩结构示意图[25]

    Fig.  8   Sketch of the throat-like structure at shock train head[25]

    图  9   激波串前缘初始激波和附面层流场结构示意图[25]

    Fig.  9   Sketch of flowfield at shock train head[25]

    图  10   激波串前缘初始激波角与收缩极限的关系[25]

    Fig.  10   Relationship between angle of initial shock to contraction ratio[25]

    图  11   隔离段马赫数与收缩极限的关系[25]

    Fig.  11   Relationship between Mach number to the minimum value of contraction ratio[25]

    图  12   激波串突跳前后纹影图[25]

    Fig.  12   Schlierens of the flowfield when a rapid movement of shock train occurred[25]

    图  13   突跳前后等效喉道收缩比变化[25]

    Fig.  13   Contraction ratio of the throat-like structure when a rapid movement of shock train occurred[25]

    图  14   激波串前缘最小收缩极限分析示意图[27]

    Fig.  14   Sketch of the analysis method for the minimum value of the contraction ratio[27]

    图  15   角涡流线图[28]

    Fig.  15   Sketch of streamlines of corner vortex[28]

    图  16   角涡和中部分离涡的流场可视化[29]

    Fig.  16   Visualization of corner vortex and central vortex[29]

    图  17   激波入射引起的激波附面层交互与壁面摩擦系数分布[30]

    Fig.  17   Shock-boundary layer interactions and friction coefficient distribution induced by incident shock[30]

    图  18   三维效应下的流场速度分布[31]

    Fig.  18   Velocity contours of flowfield with sidewall effect[31]

    图  19   试验段构型和压力测点布置[36]

    Fig.  19   Configuration of test sections and the pressure measuring points[36]

    图  20   背压和挡板角度变化时序:(a)工况D,(b)工况E,(c)工况F。时间根据试验时长进行了无量化处理[36]

    Fig.  20   Time histories of the backpressure ratio and the flap angle: (a) case D, (b) case E, (c) case F. The time is non-dimensionalized by the total duration[36]

    图  21   试验过程中典型压力时序[36]

    Fig.  21   Comparisons of typical pressure histories during the test conditions[36]

    图  22   工况C和F中激波串前缘轨迹对比:(a)上壁面,(b)下壁面。时间根据试验时长进行了无量化处理[36]

    Fig.  22   Comparisons of the trajectories of the shock train leading edge with case C and case F conditions at the top wall (a) and bottom wall (b). The time is non-dimensionalized by the total duration[36]

    图  23   激波串运动物理过程简图[37]

    Fig.  23   Physical process of the shock train movement[37]

    图  24   无激波入射下低阶模型激波串前缘预估与试验结果比较[37]

    Fig.  24   Trajectories of the shock train leading edge estimated by the low-order model compared with the experimental values[37]

    图  25   模型预测激波串动态特性与数值模拟结果对比,背压压比变化为:(a) 2.83~3.05, (b) 3.05~3.39, (c) 3.39~3.62[37]

    Fig.  25   Dynamic features of the shock train compared with the CFD results, where backpressure ratio is from (a) 2.83 to 3.05, (b) 3.05 to 3.39, (c) 3.39 to 3.62[37]

    图  26   有激波入射下低阶模型激波串前缘预估与试验结果比较[37]

    Fig.  26   Trajectories of the shock train leading edge estimated by the low-order model compared with the experimental values under incident shocks[37]

    表  1   激波串运动过程的激波串前缘收缩比测量值[27]

    Table  1   Measured values of contraction ratio of the throat-like structure at shock train head[27]

    Data No. Motion state Actual CR Explain Correct
    1 Rapid motion 0.41 < 0.43 Yes
    2 Back to slow 0.50 >0.48 Yes
    3 Rapid motion 0.41 < 0.43 Yes
    4 Back to slow 0.56 >0.48 Yes
    5 Rapid motion 0.41 < 0.43 Yes
    6 Back to slow 0.50 >0.48 Yes
    7 Rapid motion 0.41 < 0.43 Yes
    8 Back to slow 0.51 >0.48 Yes
    9 Rapid motion 0.41 < 0.43 Yes
    下载: 导出CSV
  • [1]

    Curran E T, Heiser W H, Pratt D T. Fluid phenomena in scramjet combustion systems[J]. Annual review of Fluid Mechanics, 1996, 28(1):323-360 DOI: 10.1146/annurev.fl.28.010196.001543

    [2]

    Waltrup P J, Billig F S. Structure of shock waves in cylindrical ducts[J]. AIAA Journal, 1973, 11(10):1404-1408. DOI: 10.2514/3.50600

    [3]

    Billig F S. Research on supersonic combustion[J]. Journal of Propulsion and Power, 1993, 9(4):499-514. DOI: 10.2514/3.23652

    [4] 张堃元, 王成鹏, 杨建军, 等.带高超进气道的隔离段流动特性[J].推进技术, 2002, 23(4):311-314. DOI: 10.3321/j.issn:1001-4055.2002.04.012

    Zhang K Y, Wang C P, Yang J J, et al. Investigation of flow in isolator of hypersonic inlet[J]. Journal of Propulsion Technology, 2002, 23(4):311-314. DOI: 10.3321/j.issn:1001-4055.2002.04.012

    [5]

    Wang C P, Zhang K Y, Yang J J. Analysis of flows in scramjet isolator combined with hypersonic inlet[R]. AIAA-2005-24, 2005.

    [6]

    Wagner J L, Yuceil K B, Valdivia A, et al. Experimental investigation of unstart in an inlet/isolator or model in Mach 5 flow[J]. AIAA Journal, 2009, 47(6):1528-1542. DOI: 10.2514/1.40966

    [7]

    Wagner J L, Yuceil K B, Clemens N T. Velocimetry measurements of unstart of an inlet-isolator model in Mach 5 flow[J]. AIAA Journal, 2010, 48(9):1875-1888. DOI: 10.2514/1.J050037

    [8] 张航, 谭慧俊, 孙姝.进口斜激波、膨胀波干扰下等直隔离段内的激波串特性[J].航空学报, 2010, 31(9):1733-1739. http://d.old.wanfangdata.com.cn/Periodical/hkxb201009005

    Zhang H, Tan H J, Sun S. Characteristics of shock train in a straight isolator with interference of incident shock waves and corner expansion waves[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(9):1733-1739. http://d.old.wanfangdata.com.cn/Periodical/hkxb201009005

    [9] 田旭昂, 王成鹏, 程克明. Ma5斜激波串动态特性试验研究[J].推进技术, 2014, 35(8):1030-1039.

    Tian X A, Wang C P, Cheng K M. Experimental investigation of dynamic characteristics of oblique shock train in mach 5 flow[J]. Journal of Propulsion Technology, 2014, 35(8):1030-1039.

    [10]

    Tan H J, Sun S, Huang H X. Behavior of shock trains in a hypersonic inlet/isolator model with complex background waves[J]. Experiments in Fluids, 2012, 53(6):1647-1661. DOI: 10.1007/s00348-012-1386-1

    [11]

    Huang H X, Sun S, Tan H J, et al. Characterization of two typical unthrottled flows in hypersonic inlet/isolator models[J]. Journal of Aircraft, 2015, 52(5):1715-1721. DOI: 10.2514/1.C033190

    [12]

    Koo H, Raman V. Large-eddy simulation of a supersonic inlet-isolator[J]. AIAA Journal, 2012, 50(7):1596-1613. DOI: 10.2514/1.J051568

    [13]

    Xu K J, Chang J T, Zhou W X, et al. Mechanism and prediction for occurrence of shock-train sharp forward movement[J]. AIAA Journal, 2016, 54(1):1403-1412. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=260709d7c2fe50b58d7bdb3ff4625e54

    [14]

    Li N, Chang J T, Yu D R, et al. Mathematical model of shock-train path with complex background waves[J]. Journal of Propulsion and Power, 2017, 33(2):468-478. DOI: 10.2514/1.B36234

    [15]

    Xu K J, Chang J T, Zhou W X, et al. Mechanism of shock train rapid motion induced by variation of attack angle[J]. Acta Astronautica, 2017, 140:18-26. DOI: 10.1016/j.actaastro.2017.08.009

    [16]

    Dessornes O, Scherrer D, Novelli P. Tests of japhar dual mode ramjet engine[R]. AIAA-2001-1886, 2001.

    [17]

    Denis S R, Kau H P, Brandstetter A. Experimental study on transition between ramjet and scramjet modes in a dual-mode combustor[R]. AIAA 2003-7048, 2003.

    [18]

    Chun J, Scheurermann T, von Wolfersdorf J, et al. Experimental study on combustion mode transition in a scramjet with parallel injection[R]. AIAA-2006-8063, 2006.

    [19]

    Le D B, Goyne C P, Krauss R H, et al. Experimental study of a dual-mode scramjet isolator[J]. Journal of Propulsion and Power, 2008, 24(5):1050-1057. DOI: 10.2514/1.32591

    [20]

    Yu D R, Cui T, Bao W. Catastrophe, hysteresis and bifurcation of mode transition in scramjet engines and its model[J]. Science in China Series E:Technological Sciences, 2009, 52(6):1543-1550. DOI: 10.1007/s11431-009-0181-6

    [21]

    Kobayashi K, Tomioka S, Kato K, Performance of a dual-mode combustor with multistaged fuel injection[J]. Journal of Propulsion and Power, 2006, 22(3):518-526. DOI: 10.2514/1.19294

    [22]

    Qin B, Chang J T, Jiao X L, et al. Unstart margin characterization method of scramjet considering isolator-combustor interactions[J]. AIAA Journal, 2015, 53(2):493-500. DOI: 10.2514/1.J053547

    [23]

    Gnani F, Zare-Behtash H, Kontis K. Pseudo-shock waves and their interactions in high-speed intakes[J]. Progress in Aerospace Science, 2016, 82:36-56. DOI: 10.1016/j.paerosci.2016.02.001

    [24]

    Tu Q Y, Segal C. Isolator/combustion chamber interactions during supersonic combustion[J]. Journal of Propulsion and Power, 2010, 26(1):182-186. DOI: 10.2514/1.46156

    [25]

    Xu K J, Chang J T, Li N, et al. Preliminary investigation of limits of shock train jumps in a hypersonic inlet-isolator[J]. European Journal of Mechanics-B/Fluids, 2018, 72:664-675. DOI: 10.1016/j.euromechflu.2018.07.015

    [26]

    VeillardX, Tahir R, Timofeev E, et al. Limiting contractions for starting simple ramp-type scramjet intakes with overboard spillage[J]. Journal of Propulsion and Power, 2008, 24(5):1042-1049. DOI: 10.2514/1.34547

    [27]

    Xu K J, Chang J T, Li N, et al. Experimental investigation of mechanism and limits for shock train rapid forward movement[J]. Experimental Thermal and Fluid Science, 2019, 98:336-345. https://www.sciencedirect.com/science/article/pii/S0894177718303959

    [28]

    Huang H X, Tan H J, Sun S, et al. Evolution of supersonic corner vortex in a hypersonic inlet/isolator model[J]. Physics of Fluids, 2016, 28(12):126101. DOI: 10.1063/1.4971448

    [29]

    Funderburk M, Narayanaswamy V. Experimental investigation of primary and corner shock boundary layer interactions at mild back pressure ratios[J]. Physics of Fluids, 2016, 28(8):086102. DOI: 10.1063/1.4960963

    [30]

    Pirozzoli S, Grasso F. Direct numerical simulation of impinging shock wave/turbulent boundary layer interaction at Ma=2.25[J]. Physics of Fluids, 2006, 18(6):065113. DOI: 10.1063/1.2216989

    [31]

    Bermejo-Moreno I, Campo L, Larsson J, et al. Confinement effects in shock wave/turbulent boundary layer interactions through wall-modelled large-eddy simulations[J]. Journal of Fluid Mechanics, 2014, 758:5-62. DOI: 10.1017/jfm.2014.505

    [32]

    Chakravarthy R V K, Nair V, Muruganandam T M, et al. Analytical and numerical study of normal shock response in a uniform duct[J]. Physics of Fluids, 2018, 30(8):086101. DOI: 10.1063/1.5027903

    [33]

    Bruce P J K, Babinsky H. Unsteady shock wave dynamics[J]. Journal of Fluid Mechcanics, 2008, 603:463-473. DOI: 10.1017/S0022112008001195

    [34]

    Cui T, Wang Y, Yu D R. Bistability and hysteresis in a nonlinear dynamic model of shock motion[J]. Journal of Aircraft, 2014, 51(5):1373-1379. DOI: 10.2514/1.C032175

    [35]

    Su W Y, Zhang K Y. Back-pressure effects on the hypersonic inlet-isolator pseudoshock motion[J]. Journal of Propulsion and Power, 2013, 29(6):1391-1399. DOI: 10.2514/1.B34803

    [36]

    Li N, Chang J T, Xu K J, Yu D R, et al. Oscillation of the shock train in an isolator with incident shocks[J]. Physics of Fluids, 2018, 30(11):116102. DOI: 10.1063/1.5053451

    [37]

    Li N, Chang J T, Xu K J, et al. Prediction dynamic model of shock train with complex background waves[J]. Physics of Fluids, 2017, 29(11):116103. DOI: 10.1063/1.5000876

  • 期刊类型引用(4)

    1. 赵照,熊建军,冉林,易贤. 大型结冰风洞热气供气防除冰试验技术. 航空动力学报. 2024(05): 16-22 . 百度学术
    2. 陈旷,章荣平,晋荣超. 动力短舱校准装置的供气流量控制策略. 航空动力学报. 2023(01): 240-249 . 百度学术
    3. 王万波,姜裕标,黄勇,张鑫,魏然. 大型飞机襟翼吹气增升风洞试验. 航空学报. 2023(13): 73-83 . 百度学术
    4. 赵照,易贤,冉林,熊建军. 多路热气流量控制及其在防除冰试验中的应用. 测控技术. 2023(10): 18-23+37 . 百度学术

    其他类型引用(1)

图(26)  /  表(1)
计量
  • 文章访问数:  419
  • HTML全文浏览量:  469
  • PDF下载量:  14
  • 被引次数: 5
出版历程
  • 收稿日期:  2018-12-12
  • 修回日期:  2019-04-28
  • 刊出日期:  2019-06-24

目录

/

返回文章
返回
x 关闭 永久关闭