高超声速马赫反射

通过实验和数值计算模拟了激波（马赫数７．６）绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项、粘性项和化学源项分别采用了迎风ＴＶＤ格式、中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

用热流探针测量激波速度

激波速度测量是激波管和激波风洞运行状态的一个重要参数,压电传感器或光学方法测速系统成本高,而传统电离探针在激波马赫数较低、波后温度达不到空气电离程度的情况下无法满足实验要求。提出了一种使用同轴热电偶作为测速探针来测激波速度的方法,弥补了电离探针在激波马赫数较低时的不足。介绍了同轴热电偶探针测速原理,并设计了测量激波速度的系统电路。通过信号放大电路锁定激波冲激信号,触发脉冲信号发生电路,实现了一种单通道、多测点的激波风洞测速系统。分别开展以温度与热流为触发信号的风洞实验,结果表明只有使用热流信号才能满足激波测速的时间要求。     

扰动激波冲击界面不稳定性:反射激波效应

通过平面激波绕刚体圆柱的方法形成扰动激波,采用无膜技术形成N₂/SF₆均匀界面,在竖式激波管中开展了扰动激波冲击界面Richtmyer-Meshkov (RM)不稳定性实验研究。针对3种不同的无量纲距离η(圆柱到界面距离与圆柱直径之比)情形,利用高速纹影技术及平面Mie散射技术,获得了反射激波二次冲击作用下的界面演化图像。前期工作(邹立勇等,2017)显示,入射激波冲击后,界面发展为包括中心气腔和两侧台阶的Λ形结构。研究结果表明:反射激波二次冲击后,Λ形界面首先经历相位反转,然后扰动逐渐发展增强。在η=2.0情形,界面演化为气泡,而当η=3.3和4.0时,在整体的气泡结构之外,界面中心发展为尖钉结构。获得了反射激波作用后的混合区宽度,并与理论模型结果进行了比较。在界面演化线性阶段,Meyer-Blewett(MB)线性模型结果和实验结果吻合较好。在界面演化非线性阶段,Dimonte-Ramaprabhu (DR)模型结果和实验结果吻合较好。特别地,当η=4.0时,理论与实验结果差别最小。     

激波在通道口部绕射反射研究

在地下工程研究和设计中,地面运动的高速空气冲击波进入通道后的波动问题是非常重要和急待解决的问题。为了弄清激波进入通道后的传播规律,我们利用新近研制的光学激波管和其他激波管设备,对侧向激波进入通道后的各种模型进行了系统试验研究,观察了激波进入通道后的波动图像,获得了一些很有价值的成果。本文介绍了试验方法和设备,描述了地面激波进入通道后在口部的绕射、反射波动图像和波后多旋涡运动规律,定量给出了通道口部绕射区压力变化曲线。     

近轴对称内收缩流场中的激波干扰

针对高超声速三维内转式进气道中的复杂曲面激波干扰问题,提炼出不同来流迎角下的内收缩直锥流场,采用激波风洞实验观测结合数值模拟的方法,高效揭示小幅偏离轴对称状态时,内收缩直锥流场中的激波汇聚效应及干扰机理。结果表明:在轴对称情况下,即使内收缩直锥的前缘压缩角很小,由于汇聚效应,激波逐渐增强,在轴线上必然发生马赫反射并形成马赫盘,继而终止了激波的进一步汇聚增强;而在来流有迎角情况下,流场小幅度偏离轴对称后,表现出复杂的三维特征;沿流向迎风侧激波的汇聚增强比背风侧更快,在对称面上迎风侧和背风侧激波干扰的位置偏离轴线,并且可以发生规则反射;在对称面发生规则反射的情况下,激波反射后局部能够达到的压力比发生马赫反射的情况下马赫盘后的压力更高;随着内收缩直锥前缘压缩角的增大,对称面上发生规则反射时的临界迎角也呈现增大趋势。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

二维非定常化学非平衡流动高精度数值模拟及其与实验的比较

本文成功地将二阶迎风 TVD 数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中。气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向,文中详细给出了气体力学方程和组元质量守恒方程组全耦合求解过程,最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论,本文所得结果与文献〔4〕提供的实验结果和文献〔1〕提供的数值结果相符合。     

激波管中生物体的荷载机制和过程及效应研究

本文介绍了国内外利用激波管对生物冲击伤的实验研究进展。通过空中爆炸产生的激波对生物的荷载机制分析,选用62只成年雄性杂种狗,采用 BST-I 型生物激波管产兆的激波作为政伤源,实验研究了对暴露于规则反射区和马赫反射区中的动物造成冲击伤的特征和分布规律及荷载过程。结果表明:1.在相同入射超压值148.3±18.5kPa条件下,两个区域中造成动物的伤情相差3～4级;2.两个区域中造成动物各脏器损伤的分布规律显著地不同;3.两个区域中机体的荷裁机制和过程不同;4.在规则反射区中,狗体表不同部位的荷载有显著的差异。文中的一些结果能提供冲击伤、工程效应、激波管气动力和结构设计作为进一步研究的参考。     

激波风洞高焓流动及其驻点对流和辐射热流测量

在激波风洞中用氢氧燃烧驱动方法获得了总压１４ＭＰａ,总温高达７２００Ｋ的高超声速高焓平衡流,可以模拟再入飞行速度４至５ｋｍ／ｓ的真实气体效应。本文还介绍了高温气流中驻点对流和辐射传热测量技术及其测量结果。     

激波绕楔形物体的流场分析

使用数值计算方法和数值干涉技术并结合实验干涉图片,分析了激波过拐角绕楔形物体运动的流场.并采用数值彩色干涉技术,将数值计算结果转换为彩色条纹编码的数值干涉图,并与实验干涉图比较,可以辨识过激波发生跳跃后的干涉条纹的位移数,从而确定实验干涉照片的密度场.数值计算采用显式二阶迎风TVD格式.采用数值干涉技术将计算密度场转化成数值干涉图与实验干涉图进行了比较,检验了数值结果的可靠性.数值计算结果显示了包括激波、接触间断和旋涡在内的详细的流场结构,与实验相比,更为细致地刻画了流场演变过程的细节.通过数值计算结果观察到激波传播到旋涡中心附近并没有发生明显耗散,旋涡上半部的激波得到减速,下半部分的激波得到加速.     

激波穿过孔板系统的传播及其诱导的流场

基于Ｅｕｌｅｒ方程，采用高精度、高分辨率的ＧＲＰ有限差分法及算子分裂技术，对激波通过孔板系统的传播、绕射和反射，以及诱导的流场进行了数值研究。计算结果清楚地给出了激波传播、绕射、反射、及相互干扰的复杂流场结构。与实验结果比较，早期两者吻合较好，后期由于湍流发展的影响，必需考虑粘性修正。     

含灰气体平板气动加热实验研究

利用灰尘激波管入射激波后平衡区两相流开展平板传热实验。在低灰尘负荷率、来流分别为亚声速和超声速流条件下测量平板在零攻角及小攻角（α＝３°）情况沿轴向的热流分布。     

用数值模型研究树脂传递成型工艺中的树脂流动和温度变化

用控制体／有限元方法建立了树脂传递成型工艺过程中树脂流动和温度变化的数学模型。考虑到热传导效应和树脂放热反应对温度变化的影响，对树脂传递成型工艺过程中树脂流动、温度变化的数值求解方法及相应的精度和收敛性进行了分析，最后对两种不同树脂传递成型工艺过程中的温度变化进行了计算，并与相应实验结果作了对比。结果表明，本建立的数学模型能有效地确定树脂传递成型工艺过程中树脂流动和温度变化，且与试验结果吻合的较好。     

高超声速压缩拐角再附峰值热流率

给出三维和二维压缩拐角诱导激波与高超声速层流和湍流边界层相互作用的流场结构和前向压缩表面的热流率分布。实验气流M数为7．8，单位长度Re数为3．5×107／m。结果表明：激波与边界层相互作用能促进边界层转捩。只有基于再附峰值热流率位置长度的参考温度条件下的雷诺数足够低，相互作用过程才能是纯层流的。在全再附高超声速分离流中，再附峰值热流率可用Simeonides通用相关式进行估算。