前体涡发生器对轴对称高超声速进气道激波振荡流动的影响实验1)

引言

进气道作为吸气式高超声速飞行器助燃气流的捕获与压缩部件[1]，必须处于起动状态以保证飞行器的高效可靠运行.相比起动状态，高超声速进气道不起动时流场存在大范围的低速区，气流捕获量减少、压缩效率急剧下降[1-2]，甚至可能引起飞行器失控与飞行试验失败[3-4].高超声速进气道不起动问题的影响因素众多.当飞行器处于加速/巡航阶段时，燃烧室内部的高温高压气流会向上游传播并可能导致进气道不起动.这种由燃烧室反压引起的不起动流动通常伴随着强非定常性，进气道流场中可能出现周期性的分离区扩大/缩小和分离激波沿流向的前后移动[5-6]，常被称为激波振荡或喘振.由于激波振荡伴随大范围流动分离以及剧烈的非定常气动力载荷[5-7]，会严重危害飞行器控制与飞行安全，因此必须给予足够重视.

那一天，我们正在上写字课。教室里静悄悄的，大家都在认真练字。忽然，教室角落里传来一阵窃窃私语声，声音虽然细如蚊蝇，但在安静的教室里还是显得格外刺耳。我看向坐在讲台上的班长，只见她脸色大变，怒气腾腾地快步向发出声音的地方走过去。她阴着脸，手里紧握着黑板擦。她对着刚才窃窃私语的两个人打量了一番，然后把黑板擦重重地拍在桌子上，大声冲着那二人吼道：“你们两个说什么话！没看到别人正在写字吗？”这声音如同天上的惊雷，教室里立刻炸开了锅，同学们纷纷往这边看。那二人立刻涨红了脸，低下头去，不敢再言语。班长见自己的“狮吼功”起了作用，便微微地笑了起来，像一个得胜的将军那样走开了。

目前针对高超声速进气道激波振荡问题的研究多采用地面风洞实验[5-8]与数值模拟[9-10].研究时燃烧室反压多采用出口节流[5-7]、冷态喷流[11]、热态燃烧[12]等方法产生，进而形成不起动与激波振荡流动.Tan等[5,8,13-14]对于多种进气道构型的不起动与激波振荡问题进行了细致深入研究，观测到“小喘”与“大喘”两种模式[5]；并在激波振荡流动机理、频率预测[5,14]以及不起动预防[8]等方面进行细致的分析.Chang等[15]实验观测到“小喘”与“大喘”的混合模式.Li等[6]通过激波风洞实验研究发现激波振荡的产生机理与不起动时进气道入口附近的分离激波以及再起动时隔离段内部运动激波密切相关.Wagner等[7,16]实验观测了出口节流产生的进气道非定常不起动流动，发现流场的非定常振荡与节流区域的三维分离区相关.

根据现有研究，高超声速进气道激波振荡与入口/节流区附近分离区以及进气道内部的激波串密切相关[5-6,16].在此基础上，已有针对激波振荡的预防或控制研究.Valdavia等[17]使用安装在隔离段内的涡发生器与射流涡发生器对不起动/流场振荡现象进行控制.结果表明两者结合使用可降低激波串与壁面压强的振荡幅度，并明显提高进气道不起动时的反压阈值.Im等[18]试图通过等离子激励以减缓不起动问题，激励后边界层抵抗分离能力增强，不起动得到延迟.

根据前期研究，涡发生器或射流涡发生器产生的流向涡会影响下游分离区流动[17,19-20]，进而影响激波振荡流动.等离子体激励则通过改变边界层参数进而提升其对反压的抵抗力延缓进气道不起动.但在实际应用中，射流涡发生器、等离子体激励均需要附加设备，不利于飞行器的实际操控.在主动控制方法外，高超声速飞行器通常在前体的离散粗糙带实现边界层的强制转捩[21-23].其作用机理与涡发生器的作用机制相似，是通过尾流中的流向涡来增强边界层扰动而促进转捩[24-25].据此，采用安装在前体转捩带位置的涡发生器来控制激波振荡流动可能是一种更加有利的方案.一方面前体位置的涡发生器属于被动控制，不会增大飞行器结构与控制的复杂程度.另一方面，该方案可能实现“强制转捩”与“激波振荡控制”的双重作用，利于飞行器的设计优化.但迄今为止，虽然已有针对前体位置的粗糙带或涡发生器对进气道分离的控制研究[26-27]，但鲜见有关前体涡发生器对激波振荡流动的影响研究.

因此，本文开展前体转捩带位置的涡发生器对高超声速进气道激波振荡流动的影响研究，以期为实际工程应用提供指导.为避免二元/侧压缩进气道前体流场中角涡的干扰[28-29]，选择一高超声速轴对称进气道构型为研究对象.

1 实验模型与方法

1.1 实验模型

根据前期开发的设计方法[30-31]，设计出本文的实验模型如图1(a)所示.进气道设计来流马赫数Ma=6.5，捕获半径64mm.外压缩面总偏转角为19.7°，由10°锥面与偏转角为9.7°的弯曲压缩面组成.进气道鼻部采用可拆卸设计，并进行R=0.8mm的圆弧钝化.进气道唇口气流偏转角9.1°，内收缩比1.58，总收缩比为6.41.进气道内压缩段采用圆弧曲线过渡连接高4.72mm、长50mm的水平隔离段.隔离段下游长为60mm、单侧扩张角1.9°的简化燃烧室.燃烧室的下游至模型出口为水平等截面流道，在此安装环形堵块以产生燃烧室反压，激发激波振荡流动.

采用节流的面积比TR衡量堵塞程度，对应计算公式为

  

图1 实验模型Fig.1 Experimental model

 

式中，At,b对应堵塞后出口面积，Ae为无堵塞时出口面积.

实验时采用的涡发生器如图1(b)所示，采用聚四氟乙烯加工成边长3mm，内角60°的菱形，厚度分0.5mm，1.0mm两种，并粘附在鼻锥构件末端.同层流数值模拟结果相比，0.5mm和1mm涡发生器高度与当地位移边界层厚度比分别为2.1与4.2.

1.2 实验方法

实验研究在中国科学技术大学KDJB300激波风洞中进行[32].风洞来流马赫数5.9，采用平衡接触面运行，典型实验来流总温880K，总压1.3MPa，对应的单位雷诺数约4.67×106m−1.图2为典型实验工况驻室压强曲线，从图中可判断总压平衡后时间在20ms以上.若考虑总压平衡前的时间，可用的实验气流约35ms.

  

图2 典型工况驻室压强曲线Fig.2 Chamber pressure history of a typical run

为考察前体涡流发生器对进气道激波振荡流动的影响，选择无涡流发生器、0.5mm和1mm涡流发生器3种工况进行实验对比研究.为保证激波风洞实验结果的可靠性，各实验工况均进行了两次重复性实验.

实验观测采用高速纹影拍摄结合瞬态压强测量.其中纹影拍摄范围为图1(a)中直径200mm的圆形区域，拍摄帧率为15000帧/秒.压强测量点位于中心体表面，如图1(a)中CH3-CH15共13个测点.由于实验均采用无攻角来流，模型测压仅分布在中心体一侧.压强测量采用上海天沐NS-2型传感器，其最大响应频率约20kHz.实验时压强信号的采样率为1MHz.纹影拍摄与压强测量采用DG645信号延时器同步触发.

2 实验结果

2.1 TR=0%，起动流场

对于涡发生器对激波振荡流动的影响机理，这里进行定性的分析.涡发生器会在尾流中产生流向涡，促进主流区与近壁面低能流区域的动量交换.这将提高壁面附近边界层对进气道内部反压的抵抗能力，使得分离激波强度更高，在图6(b)和图7(b)纹影照片中体现更加“陡峭”.由于流向涡的三维流动特性，在同分离区相互作用会产生较明显的三维流动现象，从而使得分离激波波面、分离区内部结构比较“紊乱”.尤其在分离区增大过程中，流向涡与分离区运动方向相对，两者相互作用较剧烈；而在分离区减小过程中，两者运动方向相同，相互作用较弱.这导致涡发生器对再起动过程中(图6(d)和图7(d))流场结构的影响不如外部分离区向上游运动(图6(b)和图7(b))过程影响明显.

分析图3纹影照片中的流场结构，涡发生器作用下的流场激波波系(图3(b)和图3(c))同自然转捩(图3(a))时的激波波系大体相同.同图3(a)相比，图3(b)和图3(c)中唇缘处前缘激波略向外偏移，并且随着涡发生器高度增大偏移程度增加.图3(b)和图3(c)中存在涡发生器诱导的弱激波(图3(b)和图3(c)中红色线段标识的弱激波).对于下游弯曲压缩面波系，涡发生器的影响不显著.

 

  

图3 TR=0%不同工况下实验纹影照片(竖直刀口)Fig.3 Schlieren images as TR=0%(vertical knife)

在流场结构外，将3种实验工况下的中心体压强分布进行对比分析.如图4所示，3种工况下中心体压强变化趋势、数值基本一致，其中1mm涡发生器工况因传感器损坏缺少1个测点信号.同自然转捩结果相比，1mm涡发生器会使得外压缩面压强略有增大，在进气道内部压强存在一定差异.0.5mm涡发生器对中心体压强的影响可忽略.

在流场结构的定性分析外，为考察涡发生器对表面压强的影响，选择前体CH8，喉道CH12和隔离段CH14三个测点压强数据进行分析.实验过程中三个测点的动态压强数据如图8所示.

在激波波系外，图3(a)中壁面附近存在较明显灰度较亮的区域，但在图3(b)与图3(c)中该区域并不明显.图3(a)中近壁面灰度较亮区域对应壁面附近的边界层.而在图3(b)与图3(c)中是由于涡发生器产生的尾流涡会增强下游边界层与主流区的动量混合，使得边界层与主流区的分界变得模糊[21].

  

图4 TR=0%时中心体沿程压强分布Fig.4 Centerbody pressure distribution as TR=0%

总体而言，在进气道起动状态下，前体涡发生器会诱导弱激波，使得前缘激波位置向外偏移.涡发生器会模糊壁面附近边界层同主流区的分界.根据纹影照片与中心体压强分布，1mm高度内的涡发生器对进气道流场结构、壁面压强的影响较小.

2.2 TR=66%，激波振荡流动

2.2.1 纹影照片结果

当出口堵塞比(TR)增大到66%时，进气道出现明显的激波振荡流动.图5～图7分别列出无涡发生器、0.5mm涡发生器以及1.0mm涡发生器3种工况下某振荡周期内典型流场结构的纹影照片.

如图5(a)所示，以进气道不起动瞬间为振荡周期的开始，记t=0ms.图5(a)中进气道外压缩面存在明显的分离激波，唇缘位置存在溢流.此时进气道内部存在较高的反压，驱动分离区向上游扩大，分离激波向上游移动.当t=0.27ms，图5(b)中分离区已经达到纹影观测的边界，分离区及分离激波仍向上游移动.当t=1.07ms时，图5(c)分离区已经超出纹影图片的观测范围.根据纹影照片中气流运动方向，推测此时分离区范围达到最大.在分离激波向上游运动过程中，进气道溢流量逐步增大，进气道内部高压气体从唇口与出口逐渐排出，内部反压逐步降低.由于内部高压气体的排出，进气道流场会再次起动.如图5(d)与图5(e)所示，再起动过程中外压缩面的分离激波向下游移动，并且在t=2.0ms进气道再次起动.由于再起动过程进气道的气流捕获能力逐步恢复，气流捕获量超出节流口的气流排泄能力，进气道腔体内部会再次聚集高压气体.当t=4.0ms时，如图5(f)所示，进气道再次不起动，新的振荡周期开始.

 

  

图5 无涡发生器时激波振荡周期内典型时刻纹影照片(竖直刀口)(1∶前缘激波,2∶分离激波)Fig.5 Typical schlieren images in an oscillatory cycle of no VG(vertical knife)(1∶leading edge shock waves,2∶separated shock waves)

将图6、图7同图5中纹影照片对比，有涡发生器时激波振荡流动同样经历“进气道不起动→分离区增大，分离激波向上游移动→分离区缩小，分离激波向下游移动→进气道再起动→进气道不起动”这一过程.但添加涡发生器后，激波振荡周期明显缩短.相比自然转捩时4ms的振荡周期，0.5mm和1mm涡发生器对应的振荡周期分别为3.33ms与3.13ms.

  

图6 0.5mm涡发生器时激波振荡周期内典型时刻纹影照片(竖直刀口)(1∶前缘激波,2∶VG诱导激波,3∶分离激波)Fig.6 Typical schlieren images in an oscillatory cycle of 0.5mm VG(vertical knife)(1∶leading edge shock waves,2∶VG induced shock waves,3∶separated shock waves)

  

图6 0.5mm涡发生器时激波振荡周期内典型时刻纹影照片(竖直刀口)(1∶前缘激波,2∶VG诱导激波,3∶分离激波)(续)Fig.6 Typical schlieren images in an oscillatory cycle of 0.5mm VG(vertical knife)(1∶leading edge shock waves,2∶VG induced shock waves,3∶separated shock waves)(continued)

在振荡周期外，为考察振荡周期内进气道沿程压强分布规律，选择图9中各阶段起始、结束时刻±0.02ms内时间平均的压强分布进行分析.图10展现了Phase 1～Phase 5各阶段内沿程压强的变化趋势.

（1）冠状动脉造影技术。临床冠脉检查以通过冠状动脉造影（DSA）作为“金标准”，可以有效明确冠状动脉解剖畸形以及阻塞性病变的具体程度、具体位置以及范围等。临床运用冠状动脉造影技术可以有效诊断不典型心绞痛症状（如中老年人群心脏扩大等）、心电图异常、严重心律失常、心力衰竭、疑似冠状动脉病变、疑似冠状动脉畸形等，从而给予临床诊断及治疗提供有效参考依据。但由于冠状动脉造影技术属于有创性检查，且该检查方法的费用高，检查技术水平、检查设备水平、检查操作者水平要求均较高等，因此受到较大应用限制。

  

图7 1mm涡发生器时激波振荡周期内典型时刻纹影照片(竖直刀口)(1∶前缘激波,2∶VG诱导激波,3∶分离激波)Fig.7 Typical schlieren images in an oscillatory cycle of 1mm VG(vertical knife)(1∶leading edge shock waves,2∶VG induced shock waves,3∶separated shock waves)

在激波振荡研究前，首先考察出口无堵塞、进气道起动状态下前体涡发生器对进气道流场的影响.图3列出3种工况的外部流场纹影照片.如图3，进气道外部流场无明显分离，均为起动状态.

2.2.2 表面压强数据

通过图5可以发现，椒江流域土壤侵蚀量与流域内年均降雨能够保持一致的趋势变化，说明降雨是椒江流域土壤侵蚀的主导因素之一.

  

图8 TR=66%，实验时间内典型测点压强--时间变化曲线Fig.8 Pressure-time data of typical survey points as TR equal to 66%

图8中实验时间内各测点压强信号均表现出典型的周期振荡现象.在图8三个测点压强信号中，进气道内部的CH12,CH14压强信号的周期性较明显，前体表面的 CH8压强信号周期性较差.对于无、0.5mm、1mm涡发生器工况下压强信号的振荡周期，选择周期性明显的CH12与CH14数据进行快速傅里叶分析(FFT)，并选择第一主频及其对应的振荡幅值进行对比分析.如表1，无、0.5mm、1mm涡发生器对应的压强振荡第一主频率分别为260Hz,288Hz和304Hz.图5～图7中对应纹影照片得到振荡频率分别为250Hz,300Hz,319Hz.纹影照片与压强信号的振荡频率基本一致.此外，以来流静压(P0)无量纲化，上述3种工况CH12测点第一主频对应的振荡幅值为12.1P0,9.7P0,12.5P0,CH14测点第一主频对应的振荡幅值为40.1P0,30.8P0,37.6P0.这说明对于喉道附近的CH12测点，相比无涡发生器工况，0.5mm涡发生器会降低压强信号的振荡幅度，1mm涡发生器使得压强振荡幅度略微升高.对于隔离段中的CH14测点，添加涡发生器后压强振荡幅度均有降低，其中0.5mm涡发生器降低幅度约23%.但如图8所示，有涡发生器时外压缩面CH8测点的压强振荡幅值、压强均值均有较明显增大.对于涡发生器在振荡频率以及压强幅度变化的影响机理，需要结合纹影与压强信号进行分析.

发改委11月1日公布2018年前三季度天然气运行情况。前三季度，天然气产量1156亿立方米，同比增长6.3%。

 

表1 CH12，CH14测点FFT分析结果Table 1 FFT results of CH12 and CH14 survey points

  

Cases Primary frequency/Hz Amplitude P0 CH12 CH14 CH12 CH14 no VG 260 260 12.1 40.1 0.5mm VG 288 288 9.7 30.8 1.0mm VG 304 304 12.5 37.6

2.3 涡发生器对激波振荡流动的影响分析

2.3.1 振荡周期

根据上述分析，涡发生器会缩短激波振荡周期，提高振荡频率.为分析涡发生器的影响机理，根据纹影照片与各测点的压强信号，将典型的激波振荡周期分成5个阶段(phase).以无涡流发生器工况为例，各阶段的划分如图9所示.

由于我国现有的森林资源有限，能源面临着短缺的局面。因此，在进行森林资源的合理开发利用过程中，在保证森林资源经济效益的基础上，要注重森林资源开发利用的生态效益，只有做到生态效益与经济效益的协调发展，才能够实现森林资源的合理利用。此外，在森林资源的开发利用过程，要基于长远发展的角度考虑开发问题，制定科学合理的开发方案，体现出森林资源开发的合理性。基于我国可持续发展战略目标下，必须结合目前森林资源的实际现状，适当的调整森林资源开发利用战略目标，不但要保证其更好的适应我国经济发展的需要，而且还要充分发挥其生态功能，调节环境气候，为人类的生存发展提供重要的环境基础，使其为我国社会的发展发挥其作用。

结合图5中纹影照片,与图9中t=0ms时刻对应图5(a)中纹影照片，即进气道不起动瞬间.从该时刻开始，在图5(c)中分离区达到最大时结束对应Phase 1阶段.Phase 2对应分离区减小、分离激波向下游移动，进气道再起动过程.当t=2.0ms时，图5(d)中进气道再起动，Phase 2结束，Phase 3开始.

Phase 3对应从进气道再起动到节流反压传递到隔离段的历程.考虑CH15测点可能受到模型装配垫块附近流动干扰，以CH14测点压强上升的时刻作为Phase 3与Phase 4的界限.在Phase 4阶段，下游反压沿着隔离段向上游传播，当图9(a)中CH9测点压强上升时，表明下游反压已通过边界层传递到进气道入口，此时Phase 4结束，Phase 5开始.当进气道再次不起动时，Phase 5结束，新的振荡周期开始.

  

图9 典型振荡周期内测点的动态压强信号Fig.9 Dynamic pressure signals of transducers in a typical oscillatory cycle

同样分析0.5mm、1mm涡发生器工况的激波振荡周期，得到Phase 1到Phase 5的时间间隔如表2所示.相对无涡发生器工况，0.5mm和1.0mm涡发生器工况的振荡周期在Phase 1,Phase 3和Phase 5出现较明显的缩短.这表明激波振荡周期的缩短主要是外压缩面分离激波向上游运动过程、燃烧室反压形成并前传到隔离段过程以及反压在边界层与主流区前传的时间差的缩短引起的.此外，0.5mm涡发生器工况Phase 4略有增加，即反压从隔离段末端传递到进气道入口要耗费更多的时间.其他阶段的时间间隔变化不明显.

 

表2 激波振荡周期内各阶段时间间隔分析Table 2 Time intervals of various phases in a typical oscillatory period

  

Cases Time intervals/ms Phase 1 Phase 2 Phase 3 Phase 4 Phase 5 No VG 1.07 0.93 0.63 0.87 0.50 0.5mm VG 0.73 0.94 0.39 1.04 0.23 1.0mm VG 0.73 0.94 0.44 0.88 0.14

2.3.2 沿程压强分布

在喘振周期外，添加涡发生器后外部分离区的流动结构出现较显著的变化.一方面，添加涡流发生器后，分离区运动范围明显缩小.对比图5(c)、图6(c)和图7(c)，可以发现无涡发生器时分离区运动范围远超出纹影系统的观测范围；添加0.5mm厚涡发生器后，图6(c)中的最大分离区范围略超过纹影观测范围；图7(c)中1.0mm厚涡发生器工况分离区范围在纹影的观测范围内.另一方面，如图6(b)、图7(b)所示，分离区增大过程中，有涡发生器工况分离激波更加“陡峭”、激波波面较紊乱，并且分离激波的波后流场更加紊乱.而在图6(d)、图7(d)对应的分离区减小过程中，涡发生器也存在类似效应，但影响程度较弱.

所谓信息化支撑：借助离退休人员信息管理系统，建立分门别类、种类齐全的基础信息数据库，实现动态化管理；构建了一个天上有“云”（云计算中心）、地上有格（党员管理网格）、中间有网的新型服务管理信息化支撑体系，规范服务管理工作流程。

如上文分析，图10(a)中Phase 1对应自进气道不起动到外部分离区达到最大的过程.如图5～图7中纹影照片，Phase 1过程中分离区范围增大，分离激波与腔体内的高压气流沿前体向上游运动，进气道溢流量逐渐增大.该时间历程中，图10(a)中进气道外压缩面(X＜0.26m)压强升高、内部压强降低.其中，外压缩面压强的升高是由于内部高压气流向上游移动引起的.与此同时，在内部高压气流向上游扩张、唇罩处溢流以及出口泄流的综合作用下，进气道内部压强降低.对比图10(a)中的3种工况，无涡发生器时进气道内部压强下降幅度较大，有涡发生器后进气道内部压强下降幅度明显减小，1mm涡发生器工况内部压强的下降幅度略小于0.5mm涡发生器工况.而对于外压缩面，1mm涡发生器时外压缩面测点压强上升幅度最高、无涡发生器时压强上升幅度最小.结合图5(b)、图6(b)与图7(b)，有涡流发生器时分离激波较强、分离区范围较小.较强的分离激波使得外压缩面测点上升幅度较高.而较小的分离区范围导致进气道溢流量相对无涡流发生器时减少，内部高压气流的泄流量减少，因此内部测点压强下降幅度较小.较小的分离区范围和压强下降幅度也使得分离区向上游运动时间缩短.

在Phase 2进气道再起动过程中，图10(b)沿程测点压强整体呈下降趋势.这是由于再起动过程中进气道的气流捕获能力尚未恢复，进气道腔体中高压气流的排泄量仍高于气流的捕获流量，整体压强分布保持下降.在进气道内部，三种工况的压强变化幅度相近.由于Phase 2初始时刻有涡发生器工况压强较高，因此Phase 2结束时有涡发生器工况压强分布仍明显高于无涡发生器工况.但在外压缩面，无涡发生器时进气道再起动后压强恢复到起动数值，但有涡发生器工况前体压强仍偏高.这与图4中起动状态下涡发生器的影响规律相矛盾.究其原因，涡发生器产生的尾流涡增强了边界层与主流区的能量交换，从而“阻碍”传感器腔体内高压气体的排泄，使得外压缩面测点响应变低、压强偏高.

从高校办学职能出发，高校投入成本及功能产出的形式是多样的，尤其在建立指标体系时，各种功能、成本指标量化的数据具有不同的量纲，因此，对功能指标、成本指标无量纲处理是价值管理理论与方法应用于高校绩效考核中必须解决的重要环节。在实际应用中，可引入功能绩点与成本级点的概念，功能绩点反映的是功能产出的水平，成本级点反映的是资源投入的情况，制定统一的定量标准，并对照标准进行定量计算，实现功能与成本数据的无量纲化。

 

  

图10 振荡周期不同阶段始末时刻沿程压强分布Fig.10 Pressure distributions at the start and the end moment of various phases in an oscillatory period

进气道再起动后在风洞高速来流作用下进气道内部可能形成短暂的超声速流动，这将使得沿程测点压强降低.但再起动后进气道气流捕获量恢复最大，在出口节流效应下进气道内部会再次聚集高压气体，并且自节流位置向上游传播.由于图10(c)中Phase 3的结束时刻以反压传递到 CH14测点为依据，因而其下游的CH15测点已经受到反压影响、数值升高，CH14及其上游测点在起动的高速气流作用下压强降低.在phase 3阶段，有涡发生器时沿程测点压强仍高于无涡发生器结果.

图10(d)中Phase 4进气道捕获气流在腔体内持续聚集，使得出口反压持续增高，影响范围从CH14测点传递到CH9测点.在此过程中，受出口反压影响，进气道内部测点压强显著提高.对于Phase 4结束时刻进气道内部测点压强分布，CH15测点3个工况压强数值相近，其他的内部测点无涡发生器工况明显偏低，有涡发生器时测点压强偏高，1mm高度的涡发生器工况最高.由于进气道仍处于起动状态，外压缩面流态不变.除因涡发生器引起的传感器响应滞后，外压缩面测点压强不变.

构音障碍又称运动性构音障碍,有研究表明[2]构音障碍主要的发病机制是由于运动障碍造成的。西医在临床治疗方面多以构音障碍训练为主。

如上文分析，图10(e)中Phase 5的起始时刻对应反压通过边界层传递到CH9测点，此时进气道仍处于起动状态.当反压对主流区影响传递到进气道入口及其上游时，进气道再次不起动，Phase 5结束.Phase 5过程内部测点压强逐渐增高，外部测点压强基本不变，沿程压强恢复到Phase 1初始时刻数值.涡发生器主要影响该阶段的时间间隔，从无涡发生器时的0.5ms缩短到1mm涡发生器工况时的0.14ms.Phase 5时间间隔的缩短主要是由于有涡发生器时尾流涡促进了边界层与主流区的动量交换，这使得反压在主流区与边界层传递的时间滞后缩短.

兰州新区的西部药谷规划，已经明确提出要构建生产、研发、物流、营销、检验、新药申报、医疗健康于一体的产业园，这将是兰州生物医药产业基地科研、服务体系新的核心载体，为兰州生物医药产业在更大空间、更高层次上的发展创新提供强有力的支持。

2.3.3 激波振荡流动的控制效果分析

通过上述分析，前体涡发生器的影响主要通过其尾流中流向涡的作用施加.流向涡促进了边界层同主流区的动量交换，使得起动时外压缩面边界层同主流区的分界变模糊，并提升了边界层抵抗反压的能力.因此在激波振荡流动时，有涡流发生器工况对应的分离激波较强，激波角度较大.并且在分离区向上游扩张时，流向涡同相向运动的分离区相互干扰，进一步限制了外压缩面分离区运动范围.这一方面缩短了分离区向上游运动时间，另一方面也使得激波振荡过程中沿程压强下降幅度减小.从而进一步影响后续振荡过程中下游反压形成以及前传的时间，并综合导致振荡周期的缩短、振荡周期内进气道内部压强时均值上升、振荡幅值降低.对于激波振荡过程中外压缩压强的上升，一方面由于分离激波强度增大，另一方面涡发生产生的尾流涡会“阻碍”传感器腔体内高压气流的排泄，影响传感器的响应.

上述前体涡发生器对流动参数的改变对实际激波振荡的控制与不起动预防的影响是复杂的.激波振荡周期的缩短意味着更高频率的非定常气动荷载，这对飞行器控制与材料强度提出更高要求.对于振荡周期内各阶段的时间变化，有涡流发生器时反压传递到进气道入口至进气道不起动(Phase 5)时间间隔的缩短使得以进气道入口处压强测点作为不起动监控的方案[8,33]的响应时间缩短.0.5mm涡发生器延长了反压从隔离段向上游传播(Phase 4)的时间，这有利于实际流动控制与不起动预防.从振荡周期的变化，本文中前体涡发生器对激波振荡的控制或对这种较“剧烈”的不起动流动的预防是不利的.

霍尼韦尔UOP(www.uop.com)是国际领先的供应商和授权商，专门面向炼油、石化和天然气加工行业提供加工技术、催化剂、吸附剂、加工设备和咨询服务。霍尼韦尔UOP隶属于霍尼韦尔特性材料和技术集团。集团旗下还包括霍尼韦尔过程控制部(www.honeywellprocess.com)，是为石油天然气、炼油、石化、化工等各领域提供自动化控制、仪器设备和服务的先驱者。

前体涡发生器显著缩小了激波振荡过程中外压缩面的分离区范围.从流动控制的角度，虽然前体涡发生器未能完全抑制本文中的激波振荡流动，但前体涡发生器能够对激波振荡流动的范围产生有效的控制.此外，本文中0.5mm涡发生器使得进气道内部测点压强的振荡幅值有较明显的下降，尤其是隔离段内部的CH14测点，相比无涡发生器时下降幅度可达23%，这有效降低了非定常气动力载荷，对于飞行控制与结构安全是有利的.1mm涡发生器也使得部分测点的振荡幅值降低，但下降幅度较小.

对比 0.5mm与 1mm涡发生器的作用效果，本文研究工况中，同当地位移边界层厚度比为 2.1的0.5mm涡发生器的综合效果较好.0.5mm涡发生器对激波振荡过程中外部分离区的限制效果略小于1mm涡发生器，控制效果偏弱.但0.5mm涡发生器的振荡周期以及不起动过程中时间间隔要大于1mm涡发生器结果，更利于不起动控制.此外，0.5mm涡发生器沿程压强时均值较低，并且显著较低了进气道内部测点压强的振荡幅值，这更有利于飞行控制与安全.由于本文中激波振荡流动为剧烈的非定常不起动流动状态，今后还需针对起动/不起动临界状态下涡发生器的影响以及不起动控制开展进一步研究.

3 结论

综上所述，本文结论如下.

(1)前体0.5mm与1mm涡发生器对出口无堵塞，进气道起动状态的主流流场与壁面压强影响较小.纹影照片与壁面测压数据表明，涡发生器会产生弱激波，造成进气道唇缘处激波略向外偏移，并使得外压缩面压强略有提高.此外，无涡发生器时外压缩面边界层与主流区存在较清晰的分界，添加涡发生器后两者分界变模糊.推测是由于涡发生器产生的流向涡促进了边界层与主流区的动量交换，弱化了边界层同主流区的参数差异.

(2)当出口堵塞比TR=66%时，本文设计的轴对称进气道构型出现激波振荡现象.无涡发生器时振荡周期约4ms，外压缩面分离区范围远超出纹影系统的观测范围.0.5mm和1mm高度涡发生器分别将激波振荡周期缩短为3.33ms和3.13ms，并显著限制分离区范围，使得1mm工况最大分离区在纹影观测范围内.相比无涡发生器工况，涡发生器明显提升了激波振荡周期内的壁面压强时均数值，降低了进气道内部测点压强振荡幅度.其中0.5mm工况隔离段CH14测点压强振荡幅度下降可达23%.

(3)分析前体涡发生器的影响机理，涡发生器产生的流向涡促进了边界层同主流区的动量交换，提升了边界层抵抗反压的能力.一方面使得激波振荡过程中外压缩面分离激波更强，激波角度增大.另一方面减弱了反压沿主流区和边界层传播的时间迟滞.当分离区向上游扩张时，流向涡同相向运动的分离区的相互干扰进一步限制了外压缩面分离区运动范围，缩短了分离区向上游运动时间，降低了沿程壁面压强的下降幅度.这些综合影响振荡过程中下游反压形成以及前传的时间，使得振荡周期缩短、振荡周期内进气道内部压强时均值上升、振荡幅值降低.

20例OSAHS患者的手术前后的相关数据，采用用Windows SPSS 13.0进行统计学分析，本临床实验中涉及的数据以均数±标准差形式表示，采用t检验分析手术前后相关参数变化是否具有统计学意义，本实验中t检验结果以P＜0.05认为差异具有统计学意义。

(4)实际应用中应充分重视前体位置的类涡发生器(如强制转捩带)对高超声速进气道不起动流动的影响.本文中激波振荡流动为剧烈的非定常不起动流动状态，今后还需开展起动/不起动临界状态下涡发生器的影响以及不起动控制研究.此外，对于进气道构型的影响，尾流涡同分离区的作用过程以及相应的流动机理仍需要更加深入的观测与分析.

参考文献

1 Van Wie DM.Scramjet inlets.Scramjet Propulsion,2000,189∶447-511

2 Van Wie DM,Kwok FT,Walsh RF.Starting characteristics of supersonic inlets.AIAA Paper,1996-2914,1996

3 Voland RT,Auslender AH,Smart MK.CIAM/NASA Mach 6.5 scramjet fligh and ground Test.AIAA Paper,1999-4848,1999

4 Rozenberg Z.Second X-51 hypersonicfligh endsprematurely. Flight International. Flightglobal. 15 June 2011,http∶//www.flightglobal.co

5 Tan HJ,Sun S,Yin ZL.Oscillatory flws of rectangular hypersonic inlet unstart caused by downstream mass-flw choking.Journal of Propulsion and Power,2009,25(1)∶138-147

6 Li ZF,Gao WZ,Jiang HL,Yang JM.Unsteady behaviors of a hypersonic inlet caused by throttling in shock tunnel.AIAA Journal,2013,51(10)∶2485-2492

7 Wagner JL,Valdivia A,Clemens NT,et al.Experimental investigation of unstart in an inlet/isolator model in Mach 5 flw.AIAA Journal,2009,47(6)∶1528-1542

8 Tan HJ,Li LG,Wen YF,et al.Experimental investigation of the unstart process of a generic hypersonic inlet.AIAA Journal,2011,49(2)∶279-288

9 Trapier S,Deck S,Duveau P.Delayed detached-eddy simulation and analysis of supersonic inlet buzz.AIAA Journal,2008,46(1)∶118-131

10 Hong W,Kim C.Computational study on hysteretic inlet buzz characteristics under varying mass flw conditions.AIAA Journal,2014,52(7)∶1351-1373

11 Do H,Im S,Mungal MG,et al.The influenc of boundary layers on supersonic inlet flw unstart induced by mass injection.Exp Fluids,2011,51(3)∶679-691

12 Im S,Baccarella D,McGann B,et al.Unstart phenomena induced by mass addition and heat release in a model scramjet.Journal of Fluid Mechanics,2016,797∶604-629

13 Huang HX,Tan HJ,Sun S,et al.Behavior of shock train in curved isolators with complex background waves.AIAA Journal,Published online.DOI∶10.2514/1.J056166(check)

14 Zhang QF,Tan HJ,Sun S,et al.Unstart of a hypersonic inlet with side compression caused by downstream choking.AIAA Journal,2016,54(1)∶28-38

15 Chang JL,Wang L,Bao W,et al.Novel oscillatory patterns of hypersonic inlet buzz.Journal of Propulsion and Power,2012,28(6)∶1214-1221

16 Wagner JL,Yuceil KB,Clemens NT.Velocimetry measurements of unstart in an inlet-isolator model in Mach 5 flw.AIAA Journal,2010,48(9)∶1875-1888

17 Valdivia A,Yuceil KB,Wagner JL,et al.Control of supersonic inletisolator unstart using active and passive vortex generators.AIAA Journal,2014,52(6)∶1207-1218

18 Im S,Do H,Cappelli MA.The manipulation of an unstarting supersonic flw by plasma actuator.J Phys D:Appl Phys,2012,45(48)∶485-502

19 Lu FK,Li Q,Liu CQ.Microvortex generators in high-speed flw.Progress in Aerospace Sciences,2012,53∶30-45

20 童福林，李新亮，唐志共.激波与转捩边界层干扰非定常特性数值分析.力学学报,2017,49(1)∶93-104(Tong Fulin,Li Xinliang,Tang Zhigong.Numerical analysis of unsteady motion in shock wave/transitional boundary layer interaction.Chinese Journal of Theoretical and Applied Mechanics,2017,49(1)∶93-104(in Chinese))

21 Berry S,Daryabeigi K,Wurster K.Boundary-layer transition on X-43A.Journal of Spacecraft and Rockets,2010,47(6)∶922-934

22 Schneider S.Effects of roughness on hypersonic boundary-layer transition.Journal of Spacecraft and Rockets,2008,45(2)∶193-209

23 宋威,蒋增辉,贾区耀.细长锥边界层绊线转捩风洞自由飞试验.力学学报,2016,48(6)∶1301-1307(Song Wei,Jiang Zenghui,Jia Quyao.Wind-tunnel free-fligh test of boundary layer transition induced by tripped thread for slender cone.Chinese Journal of Theoretical and Applied Mechanics,2016,48(6)∶1301-1307(in Chinese))

24 Williams OJH,Smits AJ.Effect of tripping on hypersonic turbulent boundary-layer statistics.AIAA Journal,2017,55(9)∶3051-3058

25 朱德华，袁湘江，沈清等.高超声速粗糙元诱导转捩的数值模拟及机理分析.力学学报,2015,47(3)∶381-388(Zhu Dehua,Yuan Xiangjiang,Shen Qing,et al.Numerical simulation and mechanism analysis of hypersonic roughness induced transition.Chinese Journal of Theoretical and Applied Mechanics,2015,47(3)∶381-388(in Chinese))

26 高文智，李祝飞，杨基明.一种鼻锥钝化高超声速轴对称进气道流动特性实验.航空学报，2015,36(1)∶302-310(Gao Wenzhi,Li Zhufei,Yang Jiming.Flow characteristics experiments of a hypersonic axisymmetric inlet with nose bluntness.Acta Aeronautica et Astronautic Sinica,2015,36(1)∶302-310(in Chinese))

27 Fan XQ,Li H,Jia D,et al.Experimental investigation on forced boundary-layer transition of axisymmetric inlet and numerical verification Journal of Aircraft,2006,43(5)∶1544-1549

28 Eagle WE,Driscoll JF.Shock wave–boundary layer interactions in rectangular inlets∶three-dimensional separation topology and critical points.Journal of Fluid Mechanics,2014,756∶328-353

29 Panaras AG.Review of the physics of swept-shock/boundary layer interactions.Progress in Aerospace Sciences,1996,32(2)∶173-244

30 高文智,李祝飞,杨基明等.一种用于前缘钝化高超声速进气道设计评估的 CFD/特征线组合方法及其应用.推进技术,2013,34(12)∶1585-1592(Gao Wenzhi,Li Zhufei,Yang Jiming,et al.A combined CFD/characteristic method for design and performance prediciton of hypersonic inlet with leading edge bluntness.Journal of Propulsion Technology,2013,34(12)∶1585-1592(in Chinese))

31 Gao WZ,Li ZF,Yang JM.A hybrid CFD/characteristics method for fast characterization of hypersonic blunt forebody/inlet flw.Science China G:Physics,Astronomy&Mechanics,2015,58(10)∶38-45

32 李祝飞,高文智,李鹏等.二元高超声速进气道激波振荡特性实验.推进技术,2012,33(5)∶676-682(Li Zhufei,Gao Wenzhi,Li Peng,et al.Experimental investigation on the shock wave oscillation behaviors in a two-dimensional hypersonic inlet flw.Journal of Propulsion Technology,2012,33(5)∶676-682(in Chinese))

33 Srikant S,Wagner JL,Valdivia A,et al. Unstart detection in a simplified-geometrhypersonic inlet-isolator flw.Journal of Propulsion and Power,2010,26(5)∶1059-1071