Spherically imploding detonation waves initiated by two-step divergent detonation

The detonation chamber developed by K. Terao and H. G. Wagner in Göttingen has been improved, in order to concentrate the combustion energy more effectively to a focus, so that imploding detonation waves are initiated by two-step divergent detonation waves in a hemispherical space having an effective diameter of 500 mm. The imploding detonation waves are investigated by measuring their propagation velocity using ion probes and pressure variations at different positions in the space by a piezoelectric pressure transducer, while the temperature in the implosion center is measured by a laser light scattering method. The results suggest that the peak pressure at the detonation front increases with the propagation to the center more rapidly than that in the Göttingen apparatus, while the propagation velocity is almost the same. A temperature from 10⁷ K to 10⁸ K is also observed in the focus of the imploding detonation waves.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

球面聚心气相爆轰波传播过程的数值研究

应用频散可控耗散格式对球面聚心气相爆轰波的传播过程进行了数值模拟 研究. 通过跟踪波阵面上压力和温度变化，分析这些参数在爆轰波传播过程中的演变规律， 及其与几何尺度和初始条件之间的依赖关系. 研究结果表明，仅在远场波面压力的变化近似 只依赖于r/R，对称中心附近则需要考虑初始半径$R$的影响；波面压力与初始压力的变化呈 线性关系；汇聚过程中温度升高比压力慢得多.     

Numerical studies of shock focusing induced by reflection of detonation waves within a hemispherical implosion chamber

The initiation and the propagation of detonation waves in a hemispherical chamber and the imploding shock waves that are the reflected detonation waves at the chamber wall are numerically investigated. The effects of the boundary layer and the non-uniformity of the flow field induced by the detonation wave on the imploding shock stability are examined. It is found that the effect of the boundary layer separation on the chamber wall has the strongest effect on the implosion focus.     

多障碍物通道中激波诱导气相爆轰的数值研究

应用数值模拟方法，研究了直通道中激波经过多块矩形障碍物时诱导 H2/O2混合气体起爆的物理机制. 研究表明：在前导激波强度不足以诱导波后气 体直接起爆的条件下，经过激波压缩的可燃气体也可能在远离激波的障碍物之间的凹槽部位 起爆；障碍物表面产生的压缩波、膨胀波和气流滑移面对可燃气体的起爆、爆轰波的形成和 传播过程有重要的影响；添加不同稀释比的氮气可以影响爆轰波后流场的温度分布；增加障 碍物的间距可以改变可燃气起爆位置.     

以DSD理论和LS方法为基础的程序燃烧法

给出了二维弯曲爆轰波后产物流场计算方法。爆轰波阵面传播规律满足Detonation Shock Dynamics (DSD)理论并用level set (LS)方法计算,波阵面传播规律与波后流场的耦合通过程序燃烧法实现,反应进程变量可作为LS函数的函数给出。爆轰波从刚性细管道向粗管道传播产生绕射的二维计算结果表明,化学反应速率不影响波后流场分布,只影响反应区结构。此方法可用于钝感炸药的驱动计算问题。     

Behavior of detonation waves at low pressures

With respect to stability of gaseous detonations, unsteady behavior of galloping detonations and re-initiation process of hydrogen-oxygen mixtures are studied using a detonation tube of 14 m in length and 45 mm i.d. The arrival of the shock wave and the reaction front is detected individually by a double probe combining of a pressure and an ion probe. The experimental results show that there are two different types of the re-initiation mechanism. One is essentially the same as that of deflagration to detonation transition in the sense that a shock wave generated by flame acceleration causes a local explosion. From calculated values of ignition delay behind the shock wave decoupled from the reaction front, the other is found to be closely related with spontaneous ignition. In this case, the fundamental propagation mode shows a spinning detonation. Received 10 March 1997 / Accepted 8 June 1997     

Three-dimensional numerical simulation of detonations in coaxial tubes

Three-dimensional numerical simulation of detonations in both a circular tube and a coaxial tube are simulated to reveal characteristics of single spinning and two-headed detonations. The numerical results show a feature of a single spinning detonation which was discovered in 1926. Transverse detonations are observed in both tubes, however, the single spinning mode maintains the complex Mach reflection whereas the two-headed mode develops periodically from the single Mach reflection to the complex one. The calculated cell aspect ratio for the two-headed mode changes from 1.09 to 1.34 as the radius of axial insert increases from r ₁/R = 0.1 to 0.9. The calculated cell aspect ratio for r ₁/R = 0.1 is close to the experimental results without an axial insert. The formation of an unreacted gas pocket behind the detonation front was not observed in the single spinning mode; however, the two-headed mode has unreacted gas pocket behind the front near the axial insert.      

氢氧燃烧及爆轰驱动激波管

分析并观察了沿驱动段轴向分布多火塞燃烧驱动段的性能．提出主膜处同一管截面均匀分布三火花塞引燃的点火方法．用这种点火方法驱动产生的入射激波强度重复性较高,激波后气流速度、温度和压力的定常性亦大大改善,可满足气动试验实际要求．提出在驱动段尾端串接卸爆段来消除爆轰波反射高压,从而可使反向爆轰驱动段用来产生高焓高密度试验气流．这种反向爆轰驱动产生的入射激波重复性高,激波衰减弱．在主膜处的收缩段产生的反射波可缓解爆轰波后跟随的稀疏波的不利影响,从而使前向爆轰驱动具有实用性．在产生的入射激波强度相同条件下,前向爆轰驱动所需的爆轰驱动段可爆混合气初始压力可较反向爆轰低近一个量级．     

Numerical simulation of spinning detonation in square tube

A single spinning detonation wave propagating in a square tube is simulated three-dimensionally with the detailed chemical reaction mechanism for hydrogen/air mixture proposed by Petersen and Hanson. The spinning detonation is composed of a transverse detonation rotating around the wall normal to the tube axis, triple lines propagating partially out of phase, and a short pressure trail. The formation of an unburned gas pocket behind the detonation front was not observed in the present simulations because the rotating transverse detonation completely consumed the unburned gas. The calculated profiles of instantaneous OH mass fraction have a keystone shape behind the detonation front. The numerical results for the pitch and track angle on the tube wall agree well with the experimental results. This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems at Montreal, Canada, from July 31 to August 5, 2005.     

低温下小尺度钝感炸药非理想爆轰实验研究

钝感炸药爆轰波法向速度Dn与当地曲率密切相关,为研究这种非理想行为,测量了－30 ℃下直径为10、12.5、15、30 mm的钝感炸药的定态爆轰速度和爆轰波形,并与24 ℃下的实验结果进行了比较。研究表明,低温下爆速随药柱直径的增加而增加;直径相同时,温度越高,爆轰波形越平坦;温度相同时,直径越大,爆轰波形越平坦;直径较小时,低温下的爆速小于常温下的爆速,相应的法向爆速与当地曲率的关系曲线Dn()位于常温下曲线的下方;直径较大时,低温下的爆速大于常温下的爆速,相应的Dn()关系曲线位于常温下曲线的上方。     

氢气/丙烷/空气预混气体爆轰性能的实验研究

通过采用压力传感器和烟灰板两种测试设备,开展了常温常压下氢气/丙烷和空气混合气体爆轰性能的实验研究。实验过程中观察到自持爆轰波,爆轰速度比值在0.99~1之间,爆轰压力比值在0.8~1.2之间。爆轰胞格尺寸在10~50 mm范围内,建立了爆轰胞格尺寸和化学诱导长度的关系式。随着丙烷不断添加,爆轰速度减小,而爆轰压力和胞格尺寸增加。这种变化趋势起初较快,而后变缓。因为起初氢气摩尔分数较大,混合气体趋向于氢气/空气的爆轰性能;而后因丙烷摩尔质量较大,丙烷逐渐起主要作用,混合气体表现出丙烷/空气的爆轰性能。     

激波管双爆轰驱动段性能的数值模拟研究

采用频散可控的耗散格式(DCD)，求解Euler方程和一种改进的二阶段化学反应模型， 对氢氧反向-正向双爆轰驱动段激波管进行了数值模拟. 计算结果表明：当辅驱动段与主驱动 段初始压力比小于临界值时，Taylor波仍会出现，但波扇夹角较单一前向爆轰驱动段小，入 射激波马赫数衰减率变小；当初始压力比等于临界值时，主驱动段中的Taylor波完全被消除， 入射激波马赫数不再衰减. 当初始压力比大于临界值时，在主驱动段中能产生过驱动爆轰波， 不仅Taylor波被完全消除，而且驱动能力较单一前向爆轰驱动段强.     

气相爆轰波正反射激波加速研究

实验采用压力传感器测量了指定点压力时间曲线。数值模拟基于二维反应欧拉方程和基元反应模型,采用二阶附加半隐的龙格-库塔法和5阶WENO格式分别离散时间和空间导数项,获得了指定点数值压力时间曲线。理论分析基于爆轰理论和激波动力学,分析了气相爆轰波反射过程所涉及的复杂波系演变并获得了反射激波速度。结果表明:本文数值模拟和理论计算定性上重复并解释了实验现象。气相爆轰波在右壁面反射后,右行稀疏波加速反射激波。其加速原因是:尽管激波波前声速减小,但激波马赫数增大,波前气流速度减小。在低初压下,可能还由于爆轰波后未反应或部分反应气体的作用,导致反射激波加速幅度比高初压下大。     

T型管内流动气体中爆轰绕射过程的数值模拟

基于带化学反应的二维Euler方程,对H₂、O₂、Ar体积比为2:1:1的混合气体系统在T型管内的爆轰绕射进行了数值模拟。用二阶附加半隐的Runge-Kutta法和五阶WENO格式分别离散欧拉方程的时间和空间导数项,采用9组分48步基元反应简化模型描述爆轰波在静止系统和流动系统中的传播过程,得到了温度、压力、典型组元H质量分数的分布及数值胞格结构。结果表明:在流动系统中,迎风面上波阵面为斜爆轰结构,静止系统两侧和顺风面上的波阵面为完全解耦的前导激波;在水平管中,波阵面与上下壁面经历一系列马赫碰撞后,最终形成正爆轰;在流动系统中,胞格结构明显向下游偏移;横向爆轰波的产生对爆轰波的再生起到了关键作用。     

Oxyhydrogen combustion and detonation driven shock tube

The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incidence shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation. The project supported by State Science and Technology Committee, National Natural Foundation of Science of China (19082012), Chinese Academy of Sciences and Project of National High Technology of China. In memory of academician Kuo Yonghuai's 90th anniversary.     

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

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

Interaction of detonation with inert gas zone

The results of experimental study on detonation interaction with the regions of low reactivity, generated by the injection of an inert gas into reactive mixture, are reported. A square cross-section 60×60 mm, 3.6 m long detonation channel was used. The experiments were done for stoichiometric H₂−O₂ mixture at 0.3 bar and 0.5 bar initial pressure and room temperature. The inert gas (Ar, He, N₂ or CO₂) was injected from 0.523 dm³ container into the main channel 1 s before ignition. The size of the inert zone was controlled by inert initial pressure. The behavior of detonation was studied using direct streak photography and pressure transducers. The study has shown that at low pressure of Ar, N₂ and CO₂ injection only a slight decrease of detonation velocity occurs. At higher injection pressures complete damping of detonation and flame extinguishment occur, followed by flame reigniton and DDT outside the inert zone. For low He injection pressures an increase in detonation velocity was recorded. For higher injection pressures, detonation damping occurred, followed by DDT process. The results have shown that CO₂ has the strongest effect on damping 2H₂+O₂ detonation, with N₂ and Ar in the next places, and He very far behind. The effectiveness of inert gas in detonation damping was found proportional to its molecular weight and reciprocal to its specific heat ratio. The numerical simulations of detonation propagation through inert gas zone were also performed using the one- dimensional Detonation Lagrangean code with simple energy release model. The results of simulations are in good qualitative agreement with experimental tendencies. In particular, the model has shown that the re-initiation of detonation is enhanced by smooth concentration gradients at inert/reactive interface. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

Wave dynamic processes in cellular detonation reflection from wedges

When the cell width of the incident detonation wave (IDW) is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H ₂/O ₂ diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes. The English text was polished by Yumming Chen.     

高-低爆速圆板炸药串联爆轰引起平面爆轰波的变凸现象

为分析高-低爆速圆板炸药于空中串联爆轰时平面爆轰波在传播过程中演变成凸面波的现象,分别对平面爆轰波阵面后爆轰产物状态参数的分布、有效药量与柱面装药几何尺寸的关系、爆轰产物的状态方程及强爆轰关系式进行了讨论。并以100 mm50 mm的TNT与RC炸药串联爆轰为例,描述了平面爆轰波演变为凸面爆轰波的过程,预估了爆轰波的平面范围和波形差。预估结果与实验结果基本符合。