HyTRV标模下表面边界层稳定性特征

刘姝怡, 王傲, 段茂昌, 徐洋, 黄刚雷, 陈坚强

刘姝怡, 王傲, 段茂昌, 等. HyTRV标模下表面边界层稳定性特征[J]. 实验流体力学, doi: 10.11729/syltlx20240014.
引用本文: 刘姝怡, 王傲, 段茂昌, 等. HyTRV标模下表面边界层稳定性特征[J]. 实验流体力学, doi: 10.11729/syltlx20240014.
LIU S Y, WANG A, DUAN M C, et al. Characteristics of boundary layer stability of HyTRV model bottom[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20240014.
Citation: LIU S Y, WANG A, DUAN M C, et al. Characteristics of boundary layer stability of HyTRV model bottom[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20240014.

HyTRV标模下表面边界层稳定性特征

详细信息
    作者简介:

    刘姝怡: (1994—),女,云南宣威人,博士,助理研究员。研究方向:边界层转捩。E-mail:LSYNWPU@126.com

    通讯作者:

    黄刚雷: E-mail:huangganglei2021@163.com

  • 中图分类号: V211.7

Characteristics of boundary layer stability of HyTRV model bottom

  • 摘要: 高超声速转捩研究飞行器(Hypersonic Transition Research Vehicle, HyTRV)是为高超声速复杂三维边界层转捩研究设计的升力体标模,更接近真实飞行器外形。为深入认识HyTRV边界层稳定性特征及转捩机制,验证数值模拟和理论分析结果,本文在马赫6常规风洞中,采用红外热成像技术和高频压力传感器获得了0°和2°迎角时HyTRV标模下表面(迎风面)边界层转捩阵面及不稳定波演化特征,并与一维稳定性分析结果进行了对比。研究结果表明:迎角为2°时,标模下表面横流失稳区较0°迎角时缩小且后移,标模中心线处转捩先于两侧横流区。另外,横流区存在宽频的高频信号,信号幅值在转捩完成前达到饱和。该信号包含多频段扰动,多种扰动之间相互作用,使得频谱变宽,边界层最终转捩变为湍流。
    Abstract: The HyTRV (Hypersonic Transition Research Vehicle) is designed for the study of hypersonic complex three-dimensional boundary layer transition, the shape of which is close to that of the real aircraft. The principal objective of this research is to gain deeper understanding of the stability characteristics and transition mechanism of the boundary layer over the HyTRV model and verify the results of numerical simulation and theoretical analysis. In the Mach 6 wind tunnel, the transition front and unstable wave evolution characteristics at the bottom of the HyTRV model, obtained by the infrared thermal imaging technology and high-frequency pressure sensors respectively, were compared with the results of one-dimensional stability analysis. At the angle of attack of 2 degrees, the instability crossflow region on the windward side, which is located after the front of transition at the centerline, shrinks and moves backwards compared to the results with zero angle of attack. Broadband high-frequency signals, containing multiple disturbances in different frequency bands, were found in the crossflow region. The interaction between various disturbances widens the frequency spectrum. Then the boundary layer appears to breakdown. The signal amplitude reaches saturation before the transition is complete.
  • 作为汽车的主要动力源之一, 汽油发动机具有动力响应性好、升功率高、成本低等优点, 但是其燃油经济性受到汽油特性的限制, 相比柴油机要低20%~ 30%.针对该问题, 学界提出火花点火-可控自燃复合燃烧(Spark ignition-Controlled auto-ignition Hybrid Combution, SCHC)这一概念, 将汽油和柴油发动机可控自燃(Controlled Auto-Ignition, CAI)的燃烧方式加以组合, 不仅可以实现汽油机均质压燃, 还可以通过火花点火的火焰传播对CAI燃烧进行有效调控, 减轻中高负荷下CAI燃烧的爆震问题, 拓宽CAI燃烧的运行极限.但是, 随着研究的深入, SCHC燃烧暴露出高废气率条件下的平均指示有效压力(Indicated Mean Effective Pressure, IMEP)循环变动以及可控性差等问题[1-2].研究表明, 缸内湍流对初期的火核形成及火焰传播起到重要作用, 被认为是引起循环变动的重要原因之一[3-4]。为解决该问题, 本课题组提出二甲醚(Dimethyl ether, DME)微火源引燃复合燃烧概念[5], 利用高活性燃料二甲醚自燃形成高活性燃料分层, 形成微火源并引燃油气混合气进而引发后期自燃。微火源可以增强高废气稀释条件下火核形成及前期火焰传播过程, 提高燃烧稳定性, 从而实现发动机热效率的提高。

    韩国汉阳大学Junepyo Cha等[6]对比了DME-汽油双燃料与DME单燃料下可控自燃(CAI)燃烧过程。结果表明, DME喷射时刻可以直接控制着火时刻, 当着火时刻在上止点之前时, 采用DME-汽油双燃料具有更好的经济性。Brunel大学的Zhang等[7]在一台单缸汽油机上同样利用DME-汽油双燃料喷射技术实现了CAI燃烧。研究结果表明, 不同的DME喷射策略下的燃烧特征呈现较大差异。当DME单次早喷时, DME在气缸内与汽油混合气混合较为均匀, 因此整缸混合气活性分布较为一致, 燃烧过程呈现CAI燃烧特征; 而当DME单次晚喷时, 燃烧过程呈现火焰传播的特征。徐康[8]在单缸光学发动机上开展了DME微火源引燃复合燃烧的可视化实验。结果表明, 微火源的分布受DME分布影响。微火源引燃复合燃烧的初始反应区位置(微火源)呈现集中分布的特征。但是, 已有的实验研究对高活性燃料DME在发动机缸内这一受限空间受流动影响的扩散和分布机制尚不明确, 仅通过实验研究无法对缸内微火源产生及引燃汽油的作用机制进行解释。因此, 需要通过实验结合仿真手段对实验现象中的燃烧过程进行分析。

    本文基于光学发动机平台, 通过粒子图像测速(Particle Image Velocimetry, PIV)手段探究微火源引燃特征最为明显的单次晚喷工况下缸内流场分布, 并对其成因进行分析。通过Rayleigh散射测量缸内DME分布, 研究缸内复杂的流动条件对DME活性分布的影响。同时, 通过高速摄影结合三维计算流体力学(Computational Fluid Dynamics, CFD)仿真手段对DME分布及流动特征对燃烧放热过程影响的机理进行解释。

    实验均在单缸可视化光学发动机实验平台上完成。光学发动机实验平台系统如图 1所示, 其中可视化窗口的直径为71mm。在光学测量系统方面, 配备4倍频高能钇铝石榴石Nd:YAG (Neodymium-Yttrium-Aluminium Garnet)激光器、双脉冲532nm波长的Nd:YAG激光器、增强型光电耦合探测器(Intensified Charge-Coupled Device, ICCD)和高速摄影机。光学发动机主要参数如表 1所示。发动机采用的进气道及燃烧室设计可以使进气在缸内形成涡核, 从而增强缸内涡流运动。发动机配备了不同升程和持续期的凸轮轴系, 可以根据实验需求灵活调整进排气升程, 并通过机械错齿改变气门相位。由于光学发动机较低的压缩比, 采用辛烷值较高的汽油燃料难以在光学发动机上实现混合气自燃, 因此, 实验中采用的气道喷射燃料为辛烷值较低的基础燃料(Primary Reference Fuel, PRF)PRF40(PRF40表示在基础燃料中, 异辛烷的占比为40%, 余数为正庚烷的比重)。

    图  1  光学发动机台架系统示意图
    Fig.  1  Scheme of experimental setup for the optical engine
    表  1  光学发动机主要参数
    Table  1  Specifications of the optical engine
    发动机形式 四冲程单缸汽油机
    缸径 95mm
    冲程 95mm
    气门升程及相位 可变
    排量 0.67L
    压缩比 9.24
    燃烧室结构 棚顶室
    燃油喷射方式 气道喷射+缸内直喷
    气道喷射燃料 PRF40
    气道喷油压力 300kPa
    直喷燃料 DME
    直喷喷油压力 4MPa
    进气方式 自然吸气
    发动机转速 600r/min
    节气门 有节气门
    下载: 导出CSV 
    | 显示表格

    由于DME微火源引燃要求较低的DME循环喷射量, 因此DME油轨内压力不宜过高。综合考虑涉及工况点的DME喷射背压, 在实验中, DME共轨压力保持为4MPa。DME直喷喷油器采用的是6孔直喷汽油机喷油器, 单个喷孔直径为0.2mm, 该喷油器采用侧置式布置方案安装于进气道的正下方, 喷油器轴线与气缸横截面呈15°夹角。

    由于光学发动机较短的运行时间, 常规的外部废气再循环(Exhaust Gas Recirculation, EGR)系统难以从排气中将废气引入进气系统中。因此, 实验中使用氮气(N2)和二氧化碳(CO2)按照一定比例混合来模拟外部EGR。其中, N2和CO2的比例根据实验中所采用的燃料化学式, 在当量燃空比和完全燃烧的假设下计算得出。虽然在模拟的外部EGR中不包含燃烧产物H2O, 且完全燃烧的假设也忽略了外部EGR中微量活性小分子物质的存在(这会对结果产生一定的影响), 但是, 本课题组之前的燃烧采样分析结果显示[9], 外部EGR中小分子组分的变化对燃烧相位并无显著影响。因此, 实验中对模拟外部EGR组分的简化假设并不会对结果产生明显影响。

    使用由美国Convergent Science公司开发的计算流体力学软件Converge构建了发动机三维模型, 如图 2所示(图中A, B为截面位置)。其中, 发动机内流场包含燃烧室和进、排气道等区域。计算起始时间为进气门开启前10°曲轴转角(Crank Angle, 后文均以℃A表示), 计算结束时间为放热终止时刻(压缩上止点后40℃A)。在整个计算过程中, 最大网格尺寸为基础网格尺寸4mm, 最小网格尺寸为自适应网格加密极值0.125mm, 网格数量为120~440万。

    图  2  发动机几何模型
    Fig.  2  Engine geometry model

    在模拟缸内直喷DME过程时, 为了更好地模拟液滴撞壁现象, 采用Wall Film[10]模型, 该模型可模拟燃油撞壁后的黏附、反弹、摊布和飞溅等现象; 液滴之间的相互作用采用NTC Collision模型[11]进行模拟; 使用KH-RT Breakup Length模型[12]可以较好地模拟液滴的雾化破碎过程, 其优势在于可以较准确地预测油束的贯穿度; 而Frossling Drop Evaporation模型耦合Boiling模型[13]的使用可以准确模拟DME的蒸发过程。为验证所选用模型的准确性, 基于定容弹实验平台对喷雾模型进行了标定。实验与仿真贯穿距测量验证结果如图 3所示。结果表明:采用的DME喷雾和蒸发模型喷雾贯穿距、喷雾形态及锥角能够较好地模拟DME的发展情况, 满足本文研究的需求。

    图  3  喷雾模型贯穿距验证
    Fig.  3  Validation of spray penetration length between simulation and experiments

    在DME-PRF复合燃烧的化学反应动力学机理方面, 基于Cai和Pitsch等[14]构建的PRF机理和Burke等[15]构建的DME化学反应动力学机理, 构建了包含348种组分和1125步反应的DME-PRF详细化学反应动力学机理, 并基于雅可比辅助的直接关系图(Direct Relation Graph with Error Propagation, DRGEP)方法[16-17], 结合目标检索算法(Targeted Search Algorithm, TSA)、准稳态假设(Quasi-Steady State Assumption, QSSA)和氢离子敏感性分析对反应机理进行了简化, 构建包括了143种组分、729步反应的DME-PRF骨架反应机理。详细机理的点火延迟和火焰传播速度如图 4所示(相关信息请参考文献[21-26]), 详细机理与骨架机理的对比如图 5所示。结果表明, 构建的DME-PRF机理与实验验证结果吻合较好, 可以用于本文所进行的CFD仿真研究。

    图  4  详细化学反应动力学机理着火延迟与层流火焰传播速度验证
    Fig.  4  Validation of auto-ignition delay time and laminar flame speeds between the 348-species detailed mechanism and experiments
    图  5  骨架化学反应动力学机理与详细化学反应动力学机理着火延迟与火焰传播速度对比
    Fig.  5  Comparison of laminar flame speeds between 348-species detailed and 143-species skeletal mechanisms

    针对发动机燃烧室内的复杂流动条件, 采用LaVison开发的FlowMaster PIV系统测量缸内流动特征, 如图 6所示。该系统以2台SOLO120 Nd:YAG激光器作为光源, 以Laskin喷雾器将癸二酸二异辛酯溶液以300~500kPa的喷射压力喷入气流中。采用CCD照相机作为传感器进行图像采集, 并使用DaVis 7.2软件进行图像处理。根据之前的研究结果[18], 选取上止点前20°曲轴转角(Crank Angle Before Top Dead Center, ℃A BTDC)作为测量时刻, 此时缸内流场分布会对DME喷雾扩散及后续的燃烧过程产生显著影响, 且在该时刻活塞上行对测量没有影响。

    图  6  PIV测量系统示意图
    Fig.  6  Scheme of the PIV system

    除DME喷射外, 本研究所采用的发动机运行工况与PIV实验工况条件均一致, 如表 2所示。众所周知, 火花点火发动机的缸内流动主要受进气和压缩过程中缸内滚流的影响。为了更好地形成缸内涡流, 对采用的光学发动机进气道及燃烧室进行了增强缸内涡流的设计。随着活塞上行, 滚流在燃烧室及进气的影响下逐渐转化为涡流。实验测量了100个循环, 图 7展示了其中的8个随机循环的缸内流场及其涡心分布, 而图 8显示了该工况下PIV测量获得的100个循环的平均流场。缸内流场从进气侧开始, 向光学视窗的右上流动, 并沿缸壁从窗口左侧流出, 以窗口中心偏左侧为中心形成了涡流结构。

    表  2  发动机运行工况
    Table  2  Engine operation setup
    发动机转速/(r·min-1) 600
    排气门开启时刻[CA]/(°) 180
    排气门关闭时刻[CA]/(°) 328
    排气门升程/mm 2
    进气门开启时刻[CA]/(°) 396
    进气门开启时刻[CA]/(°) 508
    进气门升程/mm 1
    循环气道喷油量/(mg·cycle-1) 14.8
    气道喷油时刻[CA]/(°) 330
    循环DME直喷油量/(mg·cycle-1) 1.8
    DME直喷时刻(Start of Injection, SOI)[CA BTDC]/(°) 25
    过量空气系数(lambda) 1
    外部EGR率 27%
    进气温度/℃ 40±1
    冷却水温/℃ 85±2
    下载: 导出CSV 
    | 显示表格
    图  7  8个随机循环的缸内速度场结构和涡心分布特征
    Fig.  7  Random selected 8 cycles of in-cylinder velocity field structure and vortex center distribution
    图  8  上止点前20℃A缸内集总平均速度场(100个循环)
    Fig.  8  Ensemble averaged in-cylinder velocity field at 20℃A BTDC in the studied DME MFI case(100 cycles)

    三维CFD仿真所得到的缸内流场形成过程如图 9所示。进气结束时刻(200℃A BTDC)缸内主要流动形式为逆时针滚流, 其滚流比为1.1。同时, 在缸内进排气中轴线一侧出现约4mm的涡流核心, 该核心位置由进气道及燃烧室形状共同决定。随着活塞运行, 在燃烧室形状的影响下, 缸内滚流逐渐减弱, 滚流涡核逐渐消失, 滚流在壁面和初始涡核引导下逐渐转化为涡流, 缸内涡流比逐渐提升至0.8, 形成以涡流为主的流动特征。

    图  9  缸内流动形成过程仿真结果
    Fig.  9  The simulation of flow field development

    针对所测得的PIV流动图像及100个循环的平均流场, 对缸内流动的循环变动进行了分析。其中, 通过计算每循环缸内湍动能与平均流场的标准差来计算缸内流场的循环变动(COV)。针对缸内湍动能, 首先通过二维快速傅里叶变换(Fast Fourier Transform, FTT)将二维瞬时速度变换为空间频率域, 然后再将高于给定空间截止频率的所有傅里叶系数设置为0后执行逆FTT。截止频率通过每个测量周期的速度波动分量的功率谱密度(Power Spectral Density, PSD)确定[19-20]。缸内湍动能及其循环变动如表 3所示。结果表明, 低频湍动能约占总湍动能的90%, 此时缸内涡流以大涡为主, 同时伴随着部分小尺度涡流在废气中耗散。结合PIV分析结果, 高频湍动能的循环变动表明缸内湍流循环变动主要是由小尺度涡核位置的循环变动产生的。

    表  3  缸内湍动能及其循环变动
    Table  3  Turbulent kinetic energy and COV in cylinder
    湍动能/(m2·s-2) 循环变动(COV)/%
    总湍动能 2.12275 23.20492
    高频湍动能 0.21347 26.52597
    低频湍动能 1.86785 25.39271
    下载: 导出CSV 
    | 显示表格

    缸内DME分布及扩散过程受到流场的调控作用。由于实验条件限制, 无法在DME喷射时利用PIV手段测量缸内流场。故通过Mie散射及Rayleigh散射对DME喷雾进行测量, 并使用仿真手段解释流动对DME分布的影响。其中, Mie散射是入射光波长远小于测量介质粒子直径时出现的一种散射现象, 主要用于测量液相喷雾形态; 而Rayleigh散射在前期液相Mie散射较强时被Mie散射所覆盖, 因此主要用于测量小液滴及气相信号。图 10展示了喷雾发展历程中50个循环的平均值。在喷射初期DME未完全雾化时, 测量结果中会同时存在Rayleigh散射和Mie散射, 但由于Mie散射信号较强, 在直喷开始时刻后(After Start of Injection, ASOI)1.44℃A和2.16℃A的拍摄角度下, 可能存在Rayleigh散射信号被Mie散射所掩盖的情况。从图中只能观测到液相喷雾的形态, 而在2.88℃A ASOI时, 由于DME迅速雾化, 此时Mie散射信号大幅削弱, 已经能从图中观测到部分散射光信号由强至弱的过渡现象。在4.32℃A ASOI后(对应图 7中的结果), 散射光信号图像中已经观测不到明显的局部高亮度区, 而呈现连续的均匀过渡, 表明此时DME已经基本雾化完毕, Mie散射信号基本消失, 此时测量平面内采集到的是入射光在气相的DME中产生的Rayleigh散射信号。

    图  10  DME SOI 25℃A BTDC工况下喷雾发展历程(50循环平均值)
    Fig.  10  DME spray development with DME SOI 25℃A BTDC (averaged over 50 cycles)

    三维CFD仿真得到的5℃A ASOI缸内DME分布及流场分布如图 11所示。其中缸内整体流场形态仍为逆时针涡流, 同时, 在DME喷雾的影响下, 大涡核心逐渐向气缸的左侧移动, 使缸内左侧部分的流动增强。此时缸内DME喷雾受到缸内逆时针涡流影响而形成逆时针扩散, 而缸内涡核位置分布在DME喷雾外围, 影响DME边缘的扩散及其浓度梯度, 但不会对DME喷雾的整体形态产生影响。二者共同作用使DME移动到进排气中轴线位置形成集聚型分布, 此时集聚区域DME浓度与其他区域的比值为5:1。

    图  11  DME SOI 25℃A BTDC工况下喷雾发展历程及流场仿真结果
    Fig.  11  Simulation of DME spray development and flow field with DME SOI 25℃A BTDC

    在所选工况下, DME微火源引燃复合燃烧呈现多阶段放热的特征。图 12中展示了光学发动机实验的100个循环平均缸压、放热率以及三维仿真获得的缸压、放热率及缸内温度历程数据。主燃烧放热过程可分为3个主要阶段。第一阶段为剧烈放热阶段, 燃烧放热率在不到5℃A的时间内超过30J/℃A; 接下来的第二阶段放热率小幅下降并出现平台期, 持续时间略短于第一阶段的快速放热; 而在主燃烧放热的第三阶段, 放热率曲线斜率猛增, 缸内出现迅猛的燃烧放热反应, 此时放热达到峰值。通过多阶段放热, 可以对缸内燃烧的放热率及最大压升率进行调控, 从而抑制燃烧过程中的爆震现象, 提高缸内燃烧的稳定性。

    图  12  DME SOI 25℃A BTDC工况缸压、放热率及缸内温度历程
    Fig.  12  Pressure, heat release rate and temperature in-cylinder with DME SOI 25℃A BTDC

    为进一步了解DME、MFI复合燃烧的燃烧放热过程, 使用高速摄影采集的数据对放热过程进行分析, 如图 13所示。由于其循环变动特征十分微弱, 因此仅选择其中一个循环的放热过程进行分析。缸内首次观测到着火点的时刻为1.4℃A BTDC, 初始自燃点首先出现在可视区域的中部偏左区域, 但是该阶段的自燃反应强度较弱, 并伴随着较弱的自燃化学发光。随后, 自燃着火面积迅速扩大, 逐渐向视窗的右下方发展。从0.7℃A ATDC之后, 自燃着火的放热逐渐减缓, 此时从高速摄影的结果可以看出火焰锋面正在逐渐成型, 在2.2℃A ATDC时, 已经可以看到一个带着明亮反应锋面的火焰结构形成, 并在接下来的一段时间内, 火焰传播覆盖区域缓慢扩大, 火焰锋面的亮度也有一定程度的增加, 但是由于较慢的火焰传播速度和活塞的下行, 导致了此阶段燃烧放热率的小幅度下降, 放热率出现平台区域。在5.8℃A ATDC附近, 未燃区的PRF40燃料达到自燃条件, 开始自燃。此时, 图中火焰锋面外侧(视窗右上角)出现了十分明亮的着火点, 其亮度远高于火焰传播覆盖的区域。随后外围的混合气迅速自燃, 并覆盖几乎整个火焰锋面以外的区域。通过上述分析可知, DME微火源引燃复合燃烧的主放热阶段呈现自燃-火焰传播-自燃的3阶段顺序放热特征。

    图  13  DME SOI 25℃A BTDC工况下同步采集的燃烧放热高速摄影结果
    Fig.  13  Synchronized heat release and high-speed imaging results with DME SOI 25℃A BTDC

    三维CFD仿真获得的缸内DME分布及火焰锋面的结果如图 14所示。其中以1000K等温面来体现火焰锋面发展过程。结果表明, 在2℃A BTDC时, 初始火核在缸内高活性燃料DME集聚分布区域的中心产生, 随后向未燃区进行火焰传播, 其初始自燃区域与DME集聚区域吻合。可以看到, 在2℃A ATDC时刻放热量为总放热量的37.5%, 原高活性区域内的DME几乎已经被完全消耗, 高活性燃料DME自燃起到高能点火源的作用。随后火焰向燃烧室外围传播, 但传播速度减缓, 同时, 火焰传播使缸内的温度从850K提高到1250K。5℃A ATDC时刻, 较高的缸内温度诱发了已燃反应区外围的自燃, 形成如图 11所示的多阶段放热特征。

    图  14  DME SOI 25℃A BTDC工况下缸内DME分布与火焰锋面分布
    Fig.  14  Simulation of the distribution of DME and flame frontal results with DME SOI 25℃A BTDC

    根据前文分析可知, 缸内涡核位置会对DME集聚区域边缘的浓度梯度产生影响, 而随着涡核位置接近DME集聚区域位置, DME分布区域边缘流动增强, 扩散更为充分, DME分布区域边缘浓度梯度降低, 使DME自燃-火焰传播这一转换过程更加平顺, 从而稳定复合燃烧过程。

    基于光学可视化和计算流体力学仿真平台, 研究了缸内湍流运动对DME微火源引燃高活性燃料分层及燃烧过程的影响, 主要结论如下:

    (1) 进气结束时刻缸内初始流场为滚流, 滚流比为1.1, 同时燃烧室和进气道结构的改变在缸内产生了4mm的涡流涡核。随着活塞移动, 滚流在燃烧室形状及涡流涡核的引导下逐渐转化为涡流, 涡流比为0.8, 此时缸内主要循环变动形式为涡核位置的循环变动。

    (2) 在缸内涡流作用下, DME喷雾随流动向逆时针方向扩散, 在进排气轴线靠近进气侧分布, 该区域DME浓度与其他区域比约为5:1。此时缸内涡核位置对DME分布区域边缘的浓度梯度产生影响。

    (3) 随着活塞上行, DME活性浓区由于自燃需求较低而开始自燃, 其放热率在不到5℃A的时间内超过30J/℃A, 该自燃过程放热占总热值的37.5%, DME几乎被完全消耗。活性燃料DME的自燃及火焰传播过程使缸内温度提高了400K, 引发了缸内未燃混合气自燃。

  • 图  3   Re = 2.0 × 107 m−1,基本流

    Fig.  3   Re = 2.0 × 107 m−1,baseflow

    图  1   HyTRV标模

    Fig.  1   HyTRV model

    图  2   传感器位置示意图

    Fig.  2   The position of sensor

    图  4   不同Re下的温升分布

    Fig.  4   Temperature rise distribution under different Re

    图  5   传感器安装位置温升曲线

    Fig.  5   Temperature rise curve along the position of sensor

    图  6   Re = 1.5 × 107 m−1时,Kulite传感器测量信号功率谱密度分布

    Fig.  6   Re = 1.5 × 107 m−1, PSD of Kulite

    图  7   Re = 1.5 × 107 m−1,PCB传感器测量信号功率谱密度分布

    Fig.  7   Re = 1.5 × 107 m−1, PSD of PCB

    图  8   小波系数分布

    Fig.  8   Wavelet coefficient

    图  9   双谱分析结果

    Fig.  9   Bispectral analysis result

    图  10   不同Re下的高频信号幅值演化曲线

    Fig.  10   The amplitude evolution of high frequency signal at different Re

    图  11   α = 2°,不同Re下的标模表面温升分布

    Fig.  11   α = 2°, Temperature rise distribution under different Re

    图  12   α = 2°,Re = 2.0 × 107 m−1,基本流

    Fig.  12   α = 2°, Re = 2.0 × 107 m−1, baseflow

    图  13   Re = 2.0 × 107 m−1时,Kulite功率谱密度分布

    Fig.  13   Re = 2.0 × 107 m−1, PSD of Kulite

    图  14   Re = 2.0 × 107 m−1时,PCB传感器信号功率谱密度分布

    Fig.  14   Re = 2.0 × 107 m−1, PSD of PCB

    图  15   小波系数分布

    Fig.  15   Wavelet coefficient

    图  16   双谱分析结果(C4)

    Fig.  16   Bispectral analysis result (C4)

    图  17   N值分布和最不稳定频率分布

    Fig.  17   Distribution of N-factor and the most unstable frequency

    图  18   横流模态不稳定频域范围

    Fig.  18   The frequency range of unstable crossflow mode

    图  19   第二模态不稳定频率范围

    Fig.  19   The frequency range of unstable second mode

    表  1   实验工况

    Table  1   Experimental conditions

    马赫数Ma 迎角α/(°) Re/(m−1)
    6 0 1.0 × 107、1.5 × 107、2.0 × 107
    6 2 1.5 × 107、2.0 × 107
    下载: 导出CSV
  • [1] 余平, 段毅, 尘军. 高超声速飞行的若干气动问题[J]. 航空学报, 2015, 36(1): 7–23. DOI: 10.7527/S1000-6893.2014.0224

    YU P, DUAN Y, CHEN J. Some aerodynamic issues in hypersonic flight[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 7–23. doi: 10.7527/S1000-6893.2014.0224

    [2] 陈坚强, 涂国华, 张毅锋, 等. 高超声速边界层转捩研究现状与发展趋势[J]. 空气动力学学报, 2017, 35(3): 311–337. DOI: 10.7638/kqdlxxb-2017.0030

    CHEN J Q, TU G H, ZHANG Y F, et al. Hypersonic boundary layer transition: what we know, where shall we go[J]. Acta Aerodynamica Sinica, 2017, 35(3): 311–337. doi: 10.7638/kqdlxxb-2017.0030

    [3] 罗纪生. 高超声速边界层的转捩及预测[J]. 航空学报, 2015, 36(1): 357–372. DOI: 10.7527/S1000-6893.2014.0244

    LUO J S. Transition and prediction for hypersonic boundary layers[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 357–372. doi: 10.7527/S1000-6893.2014.0244

    [4] 解少飞, 杨武兵, 沈清. 高超声速边界层转捩机理及应用的若干进展回顾[J]. 航空学报, 2015, 36(3): 714–723. DOI: 10.7527/S1000-6893.2014.0245

    XIE S F, YANG W B, SHEN Q. Review of progresses in hypersonic boundary layer transition mechanism and its applications[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(3): 714–723. doi: 10.7527/S1000-6893.2014.0245

    [5]

    MORKOVIN M V. Transition in open flow systems-a reassessment[J]. Bulletin of the American Physical Society, 1994, 39: 1882.

    [6] 杨武兵, 沈清, 朱德华, 等. 高超声速边界层转捩研究现状与趋势[J]. 空气动力学学报, 2018, 36(2): 183–195. DOI: 10.7638/kqdlxxb-2018.0011

    YANG W B, SHEN Q, ZHU D H, et al. Tendency and current status of hypersonic boundary layer transition[J]. Acta Aerodynamica Sinica, 2018, 36(2): 183–195. doi: 10.7638/kqdlxxb-2018.0011

    [7]

    LI F, CHOUDHARI M, CHANG C L, et al. Transition analysis for the ascent phase of HIFiRE-1 flight experiment[J]. Journal of Spacecraft and Rockets, 2015, 52(5): 1283–1293. doi: 10.2514/1.a33258

    [8]

    MALIK M R. 40 years of transition research at NASA: a personal perspective[C]//SHERWIN S, SCHMID P, WU X. IUTAM Laminar-Turbulent Transition. Cham, Switzerland: Springer, 2022: 33-53. doi: 10.1007/978-3-030-67902-6_3.

    [9]

    CRAIG S A, SARIC W S. Crossflow instability in a hypersonic boundary layer[J]. Journal of Fluid Mechanics, 2016, 808: 224–244. doi: 10.1017/jfm.2016.643

    [10]

    LIU S S, YUAN X X, LIU Z Y, et al. Design and transition characteristics of a standard model for hypersonic boundary layer transition research[J]. Acta Mechanica Sinica, 2021, 37(11): 1637–1647. doi: 10.1007/s10409-021-01136-5

    [11] 陈坚强, 涂国华, 万兵兵, 等. HyTRV流场特征与边界层稳定性特征分析[J]. 航空学报, 2021, 42(4): 124317. DOI: 10.7527/S1000-6893.2020.24317

    CHEN J Q, TU G H, WAN B B, et al. Characteristics of flow field and boundary-layer stability of HyTRV[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(4): 124317. doi: 10.7527/S1000-6893.2020.24317

    [12] 郑文鹏. HyTRV外形的三维边界层转捩实验研究[D]. 长沙: 国防科技大学, 2021.

    ZHEN W P. Experimental study on three-dimensional boundary layer transition of HyTRV[D]. Changsha: National University of Defense Technology, 2021.

    [13] 陈久芬, 徐洋, 蒋万秋, 等. 升力体外形高超声速边界层转捩红外测量实验[J/OL]. 实验流体力学. doi: 10.11729/syltlx20220030.

    CHEN J F, XU Y, JIANG W Q, et al. Infrared thermogram measurement experiment of hypersonic boundary-layer transition of a lifting body[J/OL]. Journal of Experiments in Fluid Mechanics, 2022. doi: 10.11729/syltlx20220030.

    [14]

    MEN H Y, LI X L, LIU H W. Direct numerical simulations of hypersonic boundary layer transition over a hypersonic transition research vehicle model lifting body at different angles of attack[J]. Physics of Fluids, 2023, 35(4): 044111. doi: 10.1063/5.0146651

    [15]

    QI H, LI X L, YU C P, et al. Direct numerical simulation of hypersonic boundary layer transition over a lifting-body model HyTRV[J]. Advances in Aerodynamics, 2021, 3(1): 31. doi: 10.1186/s42774-021-00082-x

    [16]

    CHEN X, DONG S W, TU G H, et al. Boundary layer transition and linear modal instabilities of hypersonic flow over a lifting body[J]. Journal of Fluid Mechanics, 2022, 938: A8. doi: 10.1017/jfm.2021.1125

    [17]

    KIMMEL R L, ADAMCZAK D, HARTLEY D, et al. HIFiRE–5b flight overview[C]//Proc of the 47th AIAA Fluid Dynamics Conference. 2017. doi: 10.2514/6.2017-3131.

    [18]

    TUFTS M W, GOSSE R, KIMMEL R L. PSE analysis of crossflow instability on HIFiRE 5b flight test[C]//Proc of the 47th AIAA Fluid Dynamics Conference. 2017. doi: 10.2514/6.2017-3136.

    [19]

    TUFTS M W, GOSSE R C, KIMMEL R L. Parabolized stability equation analysis of crossflow instability on HIFiRE–5b flight test[J]. Journal of Spacecraft and Rockets, 2018, 55(6): 1369–1378. doi: 10.2514/1.A34149

    [20]

    BORG M P, KIMMEL R L. Ground test measurements of boundary-layer instabilities and transition for HIFiRE–5 at flight-relevant attitudes[C]//Proc of the 47th AIAA Fluid Dynamics Conference. 2017. doi: 10.2514/6.2017-3135.

    [21]

    BORG M P, KIMMEL R L. Ground test of transition for HIFiRE–5b at flight-relevant attitudes[J]. Journal of Spacecraft and Rockets, 2018, 55(6): 1329–1340. doi: 10.2514/1.A34163

    [22]

    JULIANO T J, PAQUIN L A, BORG M P. HIFiRE–5 boundary-layer transition measured in a Mach-6 quiet tunnel with infrared thermography[J]. AIAA Journal, 2019, 57(5): 2001–2010. doi: 10.2514/1.j056750

    [23]

    BORG M P, KIMMEL R L. Measurements of crossflow instability modes for HIFiRE–5 at angle of attack[C]//Proc of the 55th AIAA Aerospace Sciences Meeting. 2017. doi: 10.2514/6.2017-1681.

    [24]

    SHI M T, ZHU W K, LEE C B. Crossflow instability of hypersonic flow over an elliptic cone[C]//Proc of the AIAA Aviation 2020 Forum, Virtual Event. 2020. doi: 10.2514/6.2020-3073.

    [25]

    PAREDES P, GOSSE R, THEOFILIS V, et al. Linear modal instabilities of hypersonic flow over an elliptic cone[J]. Journal of Fluid Mechanics, 2016, 804: 442–466. doi: 10.1017/jfm.2016.536

    [26]

    RIHA A K, GROOT K J, MOYES A, et al. Secondary-instability-mode identification in hypersonic crossflow-dominated boundary layers[C]//Proc of the AIAA Scitech 2020 Forum. 2020. doi: 10.2514/6.2020-2242.

    [27]

    FISHER D, HORSTMANN K H, RIEDEL H. Flight test measurement techniques for laminar flow. Volume 23(Les techniques de mesure en vol des ecoulements laminaires)[R]. AGARDograph 300. Flight Test Techniques Series, Volume 23. Washington: RTO/NATO, 2003.

    [28]

    BORG M P, KIMMEL R L, STANFIELD S. Traveling crossflow instability for HIFiRE–5 in a quiet hypersonic wind tunnel[C]//Proc of the 43rd Fluid Dynamics Conference. 2013. doi: 10.2514/6.2013-2737.

图(19)  /  表(1)
计量
  • 文章访问数:  100
  • HTML全文浏览量:  39
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-02-01
  • 修回日期:  2024-02-20
  • 录用日期:  2024-02-26
  • 网络出版日期:  2024-05-06

目录

/

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