气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明:气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。     

爆轰波在扩张管内临界传播特性的研究

由于稀疏波的作用，爆轰波在变截面扩张管内传播时会发生衰减。临界条件下，爆轰波恰能在扩张管内维持自持传播的特性，称为爆轰波的临界传播特性。针对Ｃ＿２Ｈ＿２－Ｏ＿２－Ａｒ不同的爆轰气体成份，在不同扩张半角的方形扩张管内，作者研究了临界初压ｐ＿（ｃｒ）、胞格尺寸Ｓ与扩张半角的关系，用Ｗｈｉｔｈａｍ理论与重点火准则分析了爆轰波在扩张管的传播，得到胞格尺寸的计算值与实验值基本相符。     

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。     

爆轰加载下金属锡层裂破碎数值模拟

对爆轰加载下低熔点金属锡的层裂破碎问题开展了数值模拟。在利用实验数据对所采用数值方法和材料模型开展对比验证的基础上，通过对样品内部物理量时间及空间分布演化对比分析，剖析了冲击加-卸载中样品内部应力波与材料相互作用过程。此外，通过对比分析不同厚度锡样品在爆轰加载下的动态行为特征，进一步认识了自由面反射稀疏波、边侧稀疏波和入射稀疏波共同作用下层裂破碎演化机制。结果表明，当样品较薄时，层裂破碎行为由反射稀疏波主导；随着样品厚度的增大，反射稀疏波主导区缩小，入射稀疏波和边侧稀疏波主导区逐渐增大。     

爆轰波拐角传播三维数值模拟

为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。     

气体非稳定爆轰波在具有声学吸收壁的直角三通弯管中传播特性研究

本文实验研究了丙烷、氧气和空气的预混气体非稳定爆轰波在具有声学吸收壁的直角三通弯管中传播特性，分析了声学吸收材料厚度对非稳定爆轰波传播特性的影响．实验结果表明，在直角三通弯管中安装声学吸收材料对气体非稳定爆轰波有明显的衰减作用，在弯管前后非稳定爆轰波的传播速度和压力都有很大程度的减小．而且随着声学吸收材料厚度的增加，气体非稳定爆轰波强度衰减幅度增大．这一实验结论对工业安全具有重要的参考价值．     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体，在初始压强20kPa的条件下， 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动 力学过程.     

气相爆轰波传播过程中的自点火效应

基于基元反应模型和单步反应模型,对直管道中H₂-air混合气体中爆轰波的传播过程进行了数值模拟,揭示了气相爆轰波传播过程中的自点火效应。利用数值模拟方法计算了不同爆轰模型的点火延迟时间,并得到了爆轰波三波点的传播过程以及所形成胞格结构的尺寸。结果表明,胞格宽度与点火延迟时间成正比;爆轰波诱导区内气体的点火延迟时间与三波点的运动周期基本一致。进一步对结果分析可知,爆轰波的自维持传播取决于点火延迟时间(表征化学反应的特征时间)和三波点的运动周期(表征流动的特征时间)的匹配;当二者相匹配时,经过前导激波压缩后形成的高温高压爆轰气体,在短时间内实现了自点火,同时释放出大量的能量推动了爆轰波的前进,即爆轰波的稳定自维持传播依靠其自点火机制。     

不同角度分叉管道内氢气-空气爆轰传播特性

在3种角度分叉管道内开展化学计量比氢气-空气爆轰实验,采用自制的火焰传感器和烟迹法分别获得了爆轰波传播速度和胞格结构,探究了不同角度管道分叉对爆轰传播的影响。结果表明:氢气-空气爆轰在经过分叉三通时受分叉口稀疏波影响导致爆轰波衰减解耦,但随着入射激波与下游管道壁面碰撞,逐渐由规则反射向马赫反射转变,最终完成重起爆过程。其中,直通支管内爆轰衰减主要受支管入口面积的影响,随着分叉角度增大,入口面积减小,爆轰衰减程度和重起爆距离也随之减小;而分叉支管内,爆轰衰减受支管入口面积与入口渐扩程度共同影响,但随着分叉角度的增大,入口面积变为主要影响因素。不同角度分叉管内的实验结果均表明,初始压力升高能显著提高爆轰稳定性,从而削弱分叉几何结构的影响。     

激波与爆轰波对撞的数值模拟研究

用二阶精度NND差分格式和改进的二阶段化学反应模型模拟了爆轰波与激波的对撞过程,研究了不同强度入射激波对爆轰过渡区域的影响. 当对撞激波较弱时,透射爆轰波演变主要受流动膨胀作用的影响,可划分为对撞影响区、爆轰恢复区和稳定发展区3个阶段. 在爆轰恢复区和稳定发展区,前导激波压力经历一个过冲、然后向稳定爆轰过渡的过程,表现了爆轰波熄爆和再起爆的物理特征. 当对撞激波较强时,可燃混合气体的高热力学参数导致了更高的化学反应活化程度,形成了弱爆轰向稳定爆轰的直接转变.      

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项，用基元反应来描述爆轰化学反应过程，获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响，爆轰波阵面发生畸变. 但由于弯管曲率半径较大，未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧，且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少，爆轰波衰减. 在出口直段，受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

一维与二维爆轰传播的时空关联特性数值研究

一维爆轰传播的理论完备、计算准确, 二维斜爆轰传播由于壁面与黏性效应, 大尺度、高精度预测还有一定难度. 利用Euler方程和H₂-Air基元反应模型, 对二维有限长楔面诱导的斜爆轰和活塞驱动一维非定常正爆轰进行计算比较研究, 从时空两个维度方面, 分析了两者在起爆过程、稀疏波传播、爆轰波面演化中的关联特性. 研究结果表明: 在过驱动度相同的条件下, 经过时空变换的活塞驱动一维爆轰传播与二维驻定斜爆轰在起爆区波系结构、波面演化特征和主要参数分布规律方面无论定性或者定量对比均符合较好, 所以, 一维非定常爆轰和二维驻定斜爆轰具有时空相关性. 两者的差异主要体现在过驱动斜爆轰受稀疏波影响过渡到近Chapman-Jouguet (C-J)爆轰状态所需的弛豫时间不同, 原因可能是起源于活塞和壁面稀疏波强度的差异. 本文提出的一维与二维爆轰传播的时空关联方法不仅有助于认知斜爆轰起爆、过驱爆轰产生、胞格爆轰演化的三阶段规律, 还可以对比揭示壁面、边界层和黏性效应的影响, 应用在斜爆轰发动机燃烧室设计中能够有效节约计算时间和成本, 并降低复杂度.      

气相爆轰波绕射流场显示研究

采用基于红宝石激光器（波长694.3 nm）的纹影系统,对气相爆轰波绕射进行了初步的流场显示研究。采用单色激光和合适半带宽（15 nm）的滤光片,有效地消除了爆轰波自发光对流场显示的影响。合理设置激光器同步控制系统的触发延时,得到了序列的爆轰波阵面纹影照片。结果表明:图像清晰地显示了爆轰波阵面的诱导激波、横波及化学反应区。当爆轰波在左尖点处绕射,受稀疏波作用,诱导激波与化学反应区明显分离,导致爆轰波衰减为爆燃。分离的诱导激波和折皱的化学反应区在纹影图上清晰可见。诱导激波在垂直支管右壁面反射,诱导二次起爆。畸变爆轰波在水平和垂直支管中均发生马赫反射。提高初压,爆轰波受分叉口几何属性的影响减小,畸变爆轰波在水平和垂直支管下游较易恢复为自持爆轰波。     

可爆性气体爆炸极限和爆燃转变成爆轰的研究

我们成功地研制了长为6.6m、内径为100mm的柱形爆炸激波管。利用此激波管我们研究了氢气和氧气混合物的爆炸特性。研究表明:氢气和氧气混合物的可爆性极限为25％H2至84％H2（体积比）;混合物起爆的临界初始压力P0c近似与混合物的浓度无关。我们还用石英传感器技术测量了混合物从爆燃转变为爆轰的过程（DDT）,确定了爆炸速度、爆炸压力与混合物初始压力P0的关系。爆燃波速度D和压力P快速增加,并且DG0（声速）时,爆燃波能转变成爆轰。爆燃波速度D和压力P衰减时,爆燃波将熄灭。氢氧混合物浓度接近上限时,爆燃向爆轰转变时的爆轰呈过激励状态,然后逐渐趋于正常爆轰。C-J理论可近似预估气体爆轰参数。     

前向爆轰驱动变截面激波管特性的数值模拟

讨论的激波管由氢氧前向爆轰产生高压驱动气体、并配有收缩截面段,它能产生高焓超高速气流.对此进行了模拟.既对变截面角度的作用,两端截面积比的作用又对比正向反向爆轰有何不同效果做了探讨.正向爆轰时主膜处的收缩段产生的汇聚作用既加强主激波又产生反向激波,缓解了爆轰波阵面后紧跟着的稀疏波导致主激波衰减偏快的不利影响.当收缩角度在30°和45°之间时,主激波的强度较高,衰减较小;当β接近90°时在主激波传播过一小段距离之后,主激波后高温高速气流较均匀,且主激波的衰减最小,具有实用价值.驱动段与被驱动段的截面积比越大,主激波的强度越高,但是最初阶段衰减也越快.反向爆轰时产生的主激波衰减最缓,但是同样的主激波强度需要的驱动段与被驱动段初始压力比前向爆轰高1个量级.     

氢氧爆轰波在变截面扩张管道中传播的数值模拟

采用二阶精度频散控制耗散格式(DCD)和8组分20个方程的基元反应模型,对轴对称变截面管道中氢氧爆轰波传播进行数值模拟。结果表明,爆轰波传播至突变截面扩张管道时,由于稀疏波的作用可能会使爆轰波局部熄爆甚至完全熄爆,对于某些敏感度高的反应气体爆轰波可以二次起爆。而在渐变截面扩张管道爆轰波相对不易熄爆。     

惰性颗粒抑爆过程的数值模拟

对高温火球诱导的爆炸及惰性颗粒的抑爆过程进行了数值研究。由于扰动特征时间 ,两相弛豫特征时间和化学平衡特征时间的差异 ,这种带激波的两相化学反应流守恒方程具有刚性。利用处理包括组元守恒方程的全耦合TVD格式和可处理方程源项的Lax Wendroff Rubin格式分别求解气相和颗粒相 ,并根据各特征时间的比值用分步法处理方程的刚性。计算结果反映了均温燃烧火球在可燃环境中形成加速火焰和激波 ,进而成长为爆轰波的过程 ,以及惰性颗粒对爆轰波的抑制过程 ,反映了激波、化学反应和惰性颗粒之间的相互作用和耦合。计算结果表明 ,仅当颗粒浓度大于某值时才能有效抑制爆轰 ,否则无论多大的颗粒散布区域 ,皆不能使爆轰波完全抑制 ,爆轰波或是在其中发展为两相爆轰波 ,或是穿越该区域后重新成长为气相爆轰。     

乳化炸药铁板凹槽弯曲装药的实验研究

用铁板凹槽装填乳化炸药,通过离子探针法测得不同装药截面积和不同装药曲率条件下乳化炸药 爆炸传播速度,分析弯曲装药条件下装药直径和装药曲率对爆轰波在乳化炸药中传播速度的影响。结果表 明,爆轰波传播速度随着装药直径的二次方和曲率半径的减小而线性衰减,当装药直径或曲率半径小到某一 临界值时其爆轰反应中止。通过大量实验数据计算得到乳化炸药特定曲率半径和装药直径条件下的爆轰波 传播速度计算的经验公式。     

气相爆轰波在半圆形弯管中传播现象的实验研究

对气相(2H2/O2/Ar系统)爆轰波在半圆形弯管中的传播现象进行实验研究。用烟迹膜记录了弯管中爆轰波的胞格结构,采用压电传感器测量了沿弯管内外母线指定点的压力时间曲线,得到了爆轰波沿弯管内、外母线的平均速度和胞格尺寸的变化。结果表明:当平面爆轰波进入弯管后,受壁面的几何形状作用,诱导激波阵面发生弯曲。沿诱导激波阵面,自内母线到外母线方向,激波强度逐渐增大。同时,爆轰波后的化学反应区也受到影响,胞格尺寸发生较明显的变化。在本文条件下,当初压p08.00kPa,受扰动的爆轰波在弯管出口下游仍恢复为强度不变的稳定爆轰。胞格记录的三波点迹线表明:受扰动的爆轰波在出口段发生了马赫反射。实验结果还表明:当p0降至5.33kPa,平面稳定爆轰波经过半圆形弯管后,其强度发生衰减并直至出现熄灭。     

锚固洞室中不同方向爆炸应力波传播规律及裂纹形成机理研究

基于抗爆模型试验结果,利用数值分析软件研究锚固洞室中不同方向爆炸应力波传播规律及裂纹形成机理。通过分析锚固洞室爆炸压应力时程曲线规律,发现压应力时程曲线特征符合应力波的一般传播作用规律,这说明数值分析的结果比较合理。集中装药爆炸后,应力波以球面向四周扩散传播,随着时间推移,快速衰减。当应力波向上传播到地表时,会在地表附近发生多次“层裂”形成大面积裂纹,爆源从拱顶至侧墙,地表面附近受拉破坏越来越轻,其中直墙侧爆不发生受拉破坏;当应力波向下传播至地下洞室时,由于反射拉伸,会在锚固区及其末端发生“层裂”。随着应力波继续传播,经由上地表反射的拉伸波会与由洞室表面反射的拉伸波发生相遇,形成加载波,一旦加载波的强度大于围岩的动态抗拉强度,会在相应的位置形成裂纹。