Gasdynamic characteristics of toroidal shock and detonation wave converging

The modified CCW relation is applied to analyzing the shock, detonation wave converging and the role of chemical reactions in the process. Results indicate that the shock wave is strengthened faster than the detonation wave in the converging at the same initial Mach number. Euler equations implemented with a detailed chemical reaction model are solved to simulate toroidal shock and detonation wave converging. Gasdynamic characteristics of the converging are investigated, including wave interaction patterns, observable discrepancies and physical phenomena behind them. By comparing wave diffractions, converging processes and pressure evolutions in the focusing area, the different effects of chemical reactions on diffracting and converging processes are discussed and the analytic conclusion is demonstrated through the observation of numerical simulations.     

Double detonation drivers for a shock tube/tunnel

~~Double detonation drivers for a shock tube/tunnel~~     

爆轰波在炸药-金属界面上的折射分析

针对爆轰波在炸药-金属界面上折射时由实验获得的金属折射冲击波压力与经典爆轰波极曲线理论预测的压力存在显著差异这一问题, 本文展开了进一步的理论和数值模拟分析研究. 首先通过分析指出经典爆轰波极曲线理论的缺陷, 并对爆轰波极曲线理论进行了改进, 改进爆轰波极曲线理论给出了炸药爆轰波折射类型以及折射冲击作用点处的压力值. 然后发展了一个基于次特征理论来数值求解爆轰反应流动控制方程的二阶中心型Lagrange方法, 并数值模拟了一个典型的炸药爆轰波折射实验. 改进爆轰波极曲线理论和数值模拟分析结果表明, 爆轰波折射类型有三种:反射冲击波的正规折射、带Mach反射的非正规折射、无反射波的正规折射, 并且金属折射冲击波压力值随入射角增大而单调减小.     

Double detonation drivers for a shock tube/ltunnel

Recent progress on detonation drivers is reviewed. Performances of the forward detonation driver and backward detonation driver have been obsenred. To eliminate occurrenoe of a Taylor wave following the detonation wave in the primary driver and to improve the performance of the detonation driver, an additional backward detonation driver was proposed to attach to the end of the forward detonation driver. When the ratio of the initial pressures between the additional and the primary drivers becomes larger than or equal to a critical value, the Taylor wave will disappear, and thus a homogeneous driving gas with high pressure and high temperature can be generated. Furthermore, an over-driving detonation wave will be also obtained, which can increase the driving capability.     

Experimental study of the influence of annular nozzle on acoustic characteristics of detonation sound wave generated by pulse detonation engine

Acoustic characteristics of the detonation sound wave generated by a pulse detonation engine with an annular nozzle, including peak sound pressure, directivity, and A duration, are experimentally investigated while utilizing gasoline as fuel and oxygen-enriched air as oxidizer. Three annular nozzle geometries are evaluated by varying the ratio of inner cone diameter to detonation tube exit diameter from 0.36 to 0.68. The experimental results show that the annular nozzles have a significant effect on the acoustic characteristics of the detonation sound wave. The annular nozzles can amplify the peak sound pressure of the detonation sound wave at 90° while reducing it at 0° and 30°. The directivity angle of the detonation sound wave is changed by annular nozzles from 30° to 90°. The A duration of the detonation sound wave at 90° is also increased by the annular nozzles. These changes indicate that the annular nozzles have an important influence on the acoustic energy distribution of the detonation sound wave, which amplify the acoustic energy in a direction perpendicular to the tube axis and weaken it along the direction of the tube axis.     

Modelling detonation of heterogeneous explosives with embedded inert particles using detonation shock dynamics: Normal and divergent propagation in regular and simplified microstructure

This paper discusses the mathematical formulation of Detonation Shock Dynamics (DSD) regarding a detonation shock wave passing over a series of inert spherical particles embedded in a high-explosive material. DSD provides an efficient method for studying detonation front propagation in such materials without the necessity of simulating the combustion equations for the entire system. We derive a series of partial differential equations in a cylindrical coordinate system and a moving shock-attached coordinate system which describes the propagation of detonation about a single particle, where the detonation obeys a linear shock normal velocity-curvature (D_n–κ) DSD relation. We solve these equations numerically and observe the short-term and long-term behaviour of the detonation shock wave as it passes over the particles. We discuss the shape of the perturbed shock wave and demonstrate the periodic and convergent behaviour obtained when detonation passes over a regular, periodic array of inert spherical particles.     

Cell-like structure of unstable oblique detonation wave from high-resolution numerical simulation

A comprehensive numerical study was carried out to investigate the unsteady cell-like structures of oblique detonation waves (ODWs) for a fixed Mach 7 inlet flow over a wedge of 30° turning angle. The effects of grid resolution and activation energy were examined systematically at a dimensionless heat addition of 10. The ODW front remains stable for a low activation energy regardless of grid resolution, but becomes unstable for a high activation energy featuring a cell-like wave front structure. Similar to the situation with an ordinary normal detonation wave (NDW), a continuous increase in the activation energy eventually causes the wave-front oscillation to transit from a regular to an irregular pattern. The wave structure of an unstable ODW, however, differs considerably from that of a NDW. Under the present flow condition, triple points and transverse waves propagate downstream, and the numerical smoke-foil record exhibits traces of triple points that rarely intersect with each other. Several instability-driving mechanisms were conjectured from the highly refined results. Since the reaction front behind a shock wave can be easily destabilized by disturbance inherent in the flowfield, the ODW front becomes unstable and displays cell-like structures due to the local pressure oscillations and/or the reflected shock waves originating from the triple points. The combined effects of various instability sources give rise to a highly unstable and complex flow structure behind an unstable ODW front.     

三硝基甲苯起爆机理实验研究

 用单脉冲化学激波管实验方法,使用自由基清扫剂和对比速率,测定了邻位和对位硝基甲苯在高温下,裂解瞬间（500 μs） 的化学反应机理,并测定了化学反应速率常数。作为邻位硝基甲苯的同分异构体,对位硝基甲苯的主要裂解通道与其不同。通过实验发现了邻位硝基甲苯的裂解重要通道,测到它的产1氧-2氮-3,4-环丁二稀基异嚙唑（Anthranil） 在瞬间随温度变化生成和很快消失的过程。由此,测得这一化学性质极不稳定的产物的消失速率常数为：k(Anthranil)＝3.7×10¹⁵·exp(－25 800/T) s^-1。分析这一过程的机理,认为第一步是硝基甲苯的裂解,第二步是Anthranil的生成,第三步是Anthranil中的N－O键的断裂。     

Effect of thermal pressure in converging detonation waves

Propagation of converging detonation waves in various explosives is studied using the equation of state, which considers both the thermal and elastic pressures. It is seen that the rate of increase of thermal pressure is higher than that of the elastic pressure during convergence. The present equation of state is better since it also gives the variation of temperature, whereas the polytropic form of the equation of state is independent of temperature. It is seen that the total detonation pressure is slightly greater than the elastic pressure. Results are compared with those reported elsewhere.     

激光引致激波演化数值研究

在高重复频率激光推进的研究中,激波的合并是发生在激波演化后期的,同时由于脉冲间隔短,脉冲宽度对流场演化的影响也需要详细研究。考虑了激光辐照过程对流场演化的影响,通过数值计算对激波演化特性进行了研究。结果表明,初期波阵面的椭球形状逐渐转化为一个球形,球心与击穿点的距离随时间逐渐减小并最终趋于稳定。基于激波合并的应用,当激波马赫数在1~2之间时,给出了激波波阵面半径随时间的变化规律,以及激波高压区长度和波峰压强随激波波阵面半径变化规律的经验公式。     

会聚球面波圆孔衍射的一般积分式

分别从菲涅尔-基尔霍夫和瑞利-索莫非衍射积分公式出发,导出了考虑倾斜因子后会聚球面波经圆孔衍射的一般衍射积分表达式。所得公式不仅可计算衍射场的焦移,而且对弱聚焦和强聚焦情况均适用。并与Li和Wolf公式的结果进行了比较,结果证实所得公式适用广泛。     

基于压电探测判定激光支持燃烧波和爆轰波的点燃阈值

为了获得激光支持燃烧波和爆轰波的点燃阈值,采用压电探测器检测波长为1 064 nm的Nd:YAG激光作用在铝靶表面所产生的应变和冲压。从实验结果观察到压电信号的变化分为3个阶段,分别为光热弹性应变阶段、等离子体增强耦合阶段和激光支持爆轰波对靶表面的压力阶段,并从理论上研究了这3个阶段的激光与靶材料相互作用的机理,从而可以从压电信号是否发生跃变判断出激光支持燃烧波和激光支持爆轰波的点燃阈值,与其它方法所得到结果基本吻合。     

凝聚炸药爆轰波在高声速材料界面上的折射现象分析

凝聚炸药爆轰在边界高声速材料约束下传播时,爆轰波会在约束材料界面上产生复杂的折射现象.本文针对凝聚炸药爆轰波在高声速材料界面上的折射现象展开理论和数值模拟分析.首先通过建立在爆轰ZND模型上的改进爆轰波极曲线理论给出爆轰波折射类型,然后发展一种求解爆轰反应流动方程的基于特征理论的二阶单元中心型Lagrange计算方法来数值模拟典型的爆轰波折射过程.从改进爆轰波极曲线理论和二阶Lagrange方法数值模拟给出的结果看出,凝聚炸药爆轰波在高声速材料界面上的折射类型有四种:反射冲击波的正规折射、带束缚前驱波的非正规折射、带双Mach反射的非正规折射、带λ波结构的非正规折射.     

Unstable combustion induced by oblique shock waves at the non-attaching condition of the oblique detonation wave

The instability of oblique shock wave (OSW) induced combustion is examined for a wedge with a flow turning angle greater than the maximum attach angle of the oblique detonation wave (ODW), where archival results rarely exist for this case in previous literatures. Numerical simulations were carried out for wedges of different length scales to account for the ratio of the chemical and fluid dynamic time scales. The results reveal three different regimes of combustion. (1) No ignition or decoupled combustion was observed if a fluid dynamic time is shorter than a chemical time behind an OSW. (2) Oscillatory combustion was observed behind an OSW if a fluid dynamic time is longer than a chemical time behind an OSW and the fluid dynamic time is shorter than the chemical time behind a normal shock wave (NSW) at the same Mach number. (3) Detached bow shock-induced combustion (or detached overdriven detonation wave) was observed if a fluid dynamic time is longer than a chemical time behind a NSW. Since no ignition or decoupled combustion occurs as a very slow reaction and the detached wave occurs as an infinitely fast reaction, the finite rate chemistry is considered to be the key for the oscillating combustion induced by an OSW over a wedge of a finite length with a flow turning angle greater than the maximum attach angle for an ODW. Since this case has not been previously reported, grid independency was tested intensively to account for the interaction between the shock and reaction waves and to determine the critical time scale where the oscillating combustion can be observed.     

来流边界层效应下斜坡诱导的斜爆轰波

以超声速预混气中的斜爆轰波为研究对象,对其在来流边界层效应下的特性进行了实验研究.在马赫数为3的超声速预混风洞中,通过斜坡诱导产生了斜爆轰波.当来流的当量比较低时,预混气中产生的是化学反应锋面与激波面非耦合的激波诱导燃烧现象.此时边界层分离区中的化学反应放热将使分离区尺度显著增大,流场非定常性显著增强,激波位置剧烈振荡.当来流的当量比较高时,预混气将产生斜爆轰波.此时边界层分离区会影响到斜爆轰波起爆时的形态.在小尺度分离区下,斜爆轰波起爆时呈突跃结构(有横波);在中等尺度分离区下,流场固有的非定常性使斜爆轰波呈间歇突跃结构;在大尺度分离区下,斜爆轰波起爆则呈完全的平滑结构(无横波).     

基于压电探测判定激光支持燃烧波和爆轰波的点燃阈值

为了获得激光支持燃烧波和爆轰波的点燃阈值,采用压电探测器检测波长为1 064nm的Nd:YAG激光作用在铝靶表面所产生的应变和冲压。从实验结果观察到压电信号的变化分为3个阶段,分别为光热弹性应变阶段、等离子体增强耦合阶段和激光支持爆轰波对靶表面的压力阶段,并从理论上研究了这3个阶段的激光与靶材料相互作用的机理,从而可以从压电信号是否发生跃变判断出激光支持燃烧波和激光支持爆轰波的点燃阈值,与其它方法所得到结果基本吻合。     

Combustion wave propagation and detonation initiation in the vicinity of closed-tube end walls

There are not many studies on DDT with no obstacles and the initiation of DDT near the end of a closed tube. Therefore in the present study we experimentally investigate the mechanism of the combustion wave transition to a detonation wave when there are no obstacles. In particular, we show that a local explosion near the tube wall is necessary for the initiation of a detonation. Parameters that we varied are the wall configuration, distance between the ignition point and the wall, and initial filling pressure. The combustion waves and the compression waves are visualized using the Schlieren optical system. From the results, we found it is necessary for the combustion wave to reach four walls so that the detonation could be initiated by the local explosion. In the conditions of the present experiment, we exhibited that the local explosion did not occur in the vicinity of a single wall and four orthogonal walls; instead, the local explosion occurred in a situation with five orthogonal walls. The time of the local explosion and the detonation initiation is 2.6 ± 1.1 and 2.0 ± 0.1 times the characteristic time for the combustion wave to propagate hemispherically from an ignitor and reach the four walls.     

冲击波物理与爆轰物理研究进展

在中国工程物理研究院(以下简称中物院)建院50周年之际,文章作者对近年来在中物院开展的冲击波物理与爆轰物理研究工作进展情况作了介绍.首先,对该项研究工作的意义作了说明.其次,简单回顾了过去完成的几个代表性研究工作.然后,对近期主要的研究方向,其中包括冲击波温度测量和高压熔化规律研究、高密度气体物态方程研究、固体物质的宽区物态方程实验和理论研究、高压本构关系研究、材料的损伤演化及动态破坏研究、爆轰波传播规律与驱动飞层研究,以及高能炸药起爆性能与热点形成机理研究在实验技术、规律性探索到理论建模方面取得的近期主要进展作了重点介绍.文章最后对未来研究作了展望.     

液相放电等离子体特性与激波特性的数值分析

以能量平衡方程为基础,采用不同的电导率唯象模型描述了液相放电等离子体圆柱形通道特性,得到了通道内半径、温度、电阻、电流和耗散能量随时间的变化关系,还给出了距离放电间隙中心一定距离处的冲击波压力变化,并与前人利用等离子体通道球状模型计算得到的结果进行了比较。结果表明：把等离子体通道看成球状和看成圆柱状在描述通道压力和通道半径时差异显著,而在描述其他物理特性时差别不大;三种电导率模型在描述等离子体通道物理特性时,变化趋势大体相同,而在描述激波特性时,电导率模型σ₂更符合实际;通过对比电学参数与压力参数的变化,就可以在实验中根据实验数据以及具体的研究问题进行模型的适用性选择。