基于自回归模型的RBCC隔离段激波串位置识别与压力值预估

为了清楚客观地判断火箭基组合循环发动机(Rocket-based combined-cycle,RBCC)隔离段激波串位置,将Ma=6,4,3.5工况下直连试验中实测得到的RBCC隔离段测压点压力数据按照时间的先后顺序排列形成一时间序列,建立自回归(Auto-Regressive,AR)模型并计算赤池信息准则(Akaike information criterion,AIC)值,完成了不同工况下激波串前缘位置的识别。研究表明:当隔离段测压点没有受到激波串影响时,实时压力值仅存在微弱波动,模型AIC值变化较为平稳;当激波串运动至测压点处时,该点压力升高,振荡幅度明显增加,AIC值随之瞬时增大。取同一时间段内发动机沿程测压点中首个AIC值增加500以上,并在不改变工况的情况下始终保持较大值的测点位置为激波串前缘位置。与压比法相比,时间序列分析法能敏感监测到实时压力值的升高和振荡,激波串前缘位置识别更为准确。通过建立自回归模型还可以实现激波串内部压力值预估,记录连续160 ms内Ma=6,4,3.5工况下测压点压力数据,采样频率1 kHz,使用前80 ms数据建立自回归模型,完成后80 ms压力值预估及准确性检验,得到三个工况下预估平均误差分别为3.21%,7.68%,6.49%。

内凹通道中激波串流场特性的实验研究

三维内转式进气道以其高总压恢复系数和高压缩效率等优势,已成为高超声速进气道的发展趋势。然而,内转式进气道流场中更为复杂的激波串结构以及激波串/边界层相互作用成为其性能发挥的关键制约因素。针对内转式进气道中复杂的流动特性,文章开展了高超声速地面风洞实验研究,并对流场进行了可视化测量。首先,搭建了超声速直连式风洞实验台;然后,设计了方转圆段将内转式进气道中独特的内凹通道和矩形喷管平滑连接。结果表明:低来流马赫数下,内凹通道中的激波串前缘激波呈“λ”型,高来流马赫数下则呈“X”型,并且平直壁侧发生流动大分离,流场具有显著非对称性。此外,随着背压的增加,激波串继续向前移动,前缘激波扫掠时的测量点压力发生显著的振荡跳跃。

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

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

来流Mach数连续变化下的激波串运动特性

为了研究入射激波变化的隔离段内激波串的运动特性,设计并搭建了直连式变Mach数实验系统,捕捉了相同来流Mach数变化速率、不同背压变化速率下激波串的运动行为,揭示了入射激波与背压同时变化对激波串运动的影响机理。入射激波与背压同时变化时,共有3个方面的因素会影响激波串整体的上下游运动趋势,其一,Mach数变化,Mach数增大导致激波串向下游运动,该影响随Mach数增大逐渐减弱;其二,背景波系移动,背景激波反射点靠近激波串前缘时,可能引起激波串的突跳;其三,背压压比变化,背压压比增大时激波串向上游运动,该影响随背压压比增大逐渐增强。三方面因素共同作用下激波串表现出复杂的运动。

激波串的三维结构分析

针对二维数值模拟是否能够完全反映激波串真实结构的问题,研究了Mach数1.5和3.0来流条件下激波串的三维模拟结构与二维模拟结构。研究采用Reynolds平均算法的SST湍流模型对喷管及等直段分别进行三维和二维的模拟,网格量分别为1.8×107和2×105。首先,当喷管满流状态完全收敛时,通过在出口施加线性增长反压的方式诱导出激波串结构;然后,对三维模拟结构在各个视角的不同切面与二维模拟结构进行对比分析。研究发现:来流Mach数1.5下的激波串在xy与zx两个平面视角的结构都为对称结构,与二维模拟结果均高度相似;来流Mach数3.0下的激波串只有一个视角为对称结构,另一个视角为非对称结构,且非对称视角的结构与二维模拟结果高度相似。因此,对激波串的二维模拟结果仅可反映激波串三维结构中的一个“准二维区域”:Mach数1.5的激波串在两个视角的中心切面及两侧对称的一定范围内均可近似为二维结构(准二维区域占比40%左右);而对于Mach数3.0的激波串,只有对称视角中心切面及两侧一定范围内的区域可近似为二维结构(准二维区域占比60%左右)。

隔离段反压对激波串起始位置的影响

使用商业软件FLUENT对三维矩形截面隔离段流场进行了数值模拟。在来流马赫数2.3条件下,通过改变不同的反压来研究反压对不同长高比等截面隔离段和带扩张角隔离段激波串起始位置的影响。结果表明:隔离段长高比对反压与激波串起始位置之间的关系有非常大的影响,当长高比不是非常大时,激波串起始位置与反压成线性关系,当长高比超过此范围后,线性关系不再成立。

超声速扩压器中激波串结构的数值模拟

通过求解由BL湍流模型封闭的二维、轴对称及三维雷诺平均N S方程 ,数值模拟了等截面超声速扩压器中由激波 /附面层干扰诱导的复杂流场 ,比较了二维直管、圆截面直管及三维矩形截面直管中的流场特性、激波串长度及压强恢复程度。在来流马赫数为 3 0的二维直管计算中 ,采用四步Runge Kutta显式方法数值仿真了激波串自激振荡过程 。

Ma5斜激波串动态特性实验研究

为加深对激波串动态特性的认识,进行了简单管道模型在Ma风洞中的实验研究,探讨了两种不同背压变化条件下斜激波串的流动结构和动态特性。背压的生成和控制通过模型尾部两块斜板的闭合运动实现。为了更好观察实验现象,实验过程中动态压力传感器和高速纹影实行同步采集。研究发现:随着背压升高,激波串逐渐前移,由对称形态发展为非对称形态,并持续整个移动过程。背压增加的速度对激波串前缘的前移速度没有影响。激波串在管道内的前移过程并不是一个恒定运动过程,而是存在有稳定前移和急剧前移两种状态,其特性与管道内流场结构、壁面压强分布、背压大小以及距离管道出口的距离等因素相关。在不同的背压条件下,管道内的斜激波串表现出不同的频谱特征,但在同一工况中,激波串区域内不同位置处壁面压强的频谱特性相似。在Case 1状态中,压强振荡主频为f512Hz在Case 2状态中,振荡主频为f578Hz次级频率约为f260Hz。两次实验的主频均大于Pipon.niau模型理论计算值,而Case 2中频率f和声学振荡频率相近。

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

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

隔离段内正激波串受迫振荡特性研究

为了解隔离段在真实燃烧室反压作用下的工作特性,通过数值方法模拟燃烧室振荡燃烧引起的脉动反压,在来流Mai=2条件下,探讨了脉动反压引起的激波串受迫振荡特性。结果表明,在反压脉动条件下,激波串前缘位置周期振荡,其振荡频率与反压脉动频率一致,振幅与反压脉动频率负相关,反压脉动频率600Hz时,振幅仅为管道长2.97%;就抗反压能力而言,在反压脉动条件下,隔离段能承受比定常状态更大的反压峰值,且脉动频率越大,可承受瞬态反压峰值越大。激波串的受迫振荡包含两种形式:激波串整体位置的前后运动和激波节之间的相对运动。在一定程度内,反压脉动频率较低时(50Hz,100Hz),激波串整体位置的前后运动是主要的受迫振荡形式,激波节之间的相对运动很弱;随反压脉动频率升高(300Hz,600Hz),激波节之间的相对运动加剧,逐渐成为激波串受迫振荡的主要形式。另外,在激波串振荡过程中,激波串往上下游运动经过同一位置时,激波串形态会出现迟滞,这是导致总压恢复系数迟滞的原因。

等截面隔离段中激波串结构的数值模拟

采用二阶TVD格式的有限体积法偶合Baldwin Lomax代数湍流模型求解雷诺平均Navier stokes方程 ,数值模拟了二维矩形等截面隔离段中的流动现象 ,很好地模拟出由激波 /附面层干扰所形成的复杂的激波串流场结构。计算结果与国外有关文献的结果进行了比较。

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟。数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象。本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响。

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

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

中心线偏置隔离段内激波串迟滞特性研究

采用定常与非定常数值计算相结合的方法研究了一种中心线偏置的隔离段流场,分析了不同反压作用下的激波串特征。通过模拟隔离段反压升高和降低的过程,研究了隔离段内激波串的迟滞特性,并比较了扩张比为0,10%,37%隔离段内迟滞特性的差异。结果表明,在来流马赫数2.0条件下,所研究隔离段内存在两种类型的迟滞现象;在同一反压条件下,降压路径对应的激波串更靠近管道入口。当反压接近隔离段所能承受的最大反压时,流场迟滞现象消失。隔离段扩张比越大(如37%),激波串位置出现迟滞的反压范围越宽,迟滞量越小。最后利用流量匹配的观点从无粘角度解释了有粘流道内的激波串迟滞现象。

烧蚀后C-SiC复合材料隔离段内激波串长度预测

为了研究畸变气流影响下烧蚀后的C-SiC复合材料隔离段内激波串长度预测问题,设计了能够模拟进气道喉道流场畸变的隔离段直连台,并开展了马赫数2.5和3.5来流的两种不同烧蚀程度的C-SiC复合材料隔离段直连实验,获得了有效的纹影和压力数据,结合实验和理论分析,采用无量纲不平度hs/Dhy来表征烧蚀后C-SiC隔离段内壁面严重不规则粗糙特征,尝试在Waltrup公式中增加包含无量纲不平度的因子(1+4hs/Dhy),以此形成新的公式来预测烧蚀后C-SiC隔离段内激波串长度,预测精度较原来公式大幅提高,按照10%的误差带分析,修正后的Waltrup公式预测精度提高了33.3%。

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

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

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

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

超燃冲压发动机隔离段中激波串演化发展过程研究

隔离段内流场波系结构对超燃冲压发动机进气道的启动起着重要的作用,它是衔接进气道和燃烧室的桥梁。采用二维耦合隐式NS方程和RNGk-ε湍流模型对二维超燃冲压发动机隔离段内流场进行了数值模拟,研究了其内部激波串的演化发展过程。结果表明:出口反压的升高可使激波串沿隔离段内流道往上游移动;伴随着激波串的前移和马赫盘的出现,隔离段壁面附近涡亦成对出现;随着隔离段出口反压的逐步升高,激波串强度逐渐增强,由斜激波逐步转变为正激波串、λ正激波,直至将激波串推出隔离段,造成进气道的不启动;等直隔离段相比单壁扩张隔离段和双壁扩张隔离段,其抗反压能力更强,更能有效隔离燃烧室压力对进气道工作状态的影响。

超燃冲压发动机隔离段内附面层/激波串相互干扰研究

隔离段内激波串的产生和发展以及激波/附面层相互干扰现象是极为复杂的,有效地进行激波串的组织是研究隔离段的关键所在,而其性能的好坏直接影响着超燃冲压发动机的性能。采用数值模拟的方法对不同来流附面层厚度工况的二维轴对称隔离段内流场流动特性进行了数值计算,分析了附面层/激波相互作用机理和附面层对隔离段激波串及隔离段性能的影响。结构表明:压缩-膨胀-再压缩-再膨胀……的气流流动挤压过程导致激波串的形成,逆压梯度的存在构成了附面层分离;附面层厚度的增加影响着激波串起始位置和结构;随着附面层厚度的增加,出口总压恢复系数和质量平均马赫数降低,且随着反压增大,变化趋势趋于明显。

动态节流下激波串运动特性的模拟和分析

为研究动态反压下的激波串特性,针对一种带凹腔的二元进气道/隔离段构型,在马赫数为6的来流下模拟了堵塞比从0.20增长到0.32再保持不变的动态节流流动,分析了堵塞比增长时间(1~10 ms)对激波串运动的影响。结果表明:激波串在节流变化初期向下游运动,随后向上游运动并最终稳定在某一位置。当堵塞比增长时间在5 ms以内时,激波串向下游和向上游运动的幅度分别为3 mm以内和约18 mm,且激波串运动滞后于节流变化,滞后时间随着增长时间的延长而缩短。当增长时间大于等于6 ms时,激波串可向下游运动到凹腔中部,幅度可达31 mm,并伴随着流动振荡;向上游运动幅度仍约为18 mm,激波串运动与节流变化近似同步。分析表明:较短增长时间工况下,激波串运动滞后主要是因为节流引起反压升高、传播时间大于堵塞比增长时间;较长增长时间工况下,凹腔内流动振荡主要是堵塞比增长初期凹腔亚声速区排出流量增加,回流区横向尺度减小,导致凹腔超声速区膨胀并出现“壅塞”,产生分离激波与回流区相互作用、发生振荡。工程设计时应考虑激波串运动的滞后及其对流动性能的影响。