悬浮RDX炸药粉尘爆轰的数值模拟

用两相流模型对悬浮RDX炸药粉尘爆轰波进行了数值模拟。RDX炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,颗粒表面发生熔化。参考液滴在高速气流作用下剥离的效应,假设炸药熔化部分在高速气流的作用下发生剥离,破碎成极小的颗粒,瞬时发生分解反应,释放出能量支持爆轰波传播。数值模拟了在不同粒径和浓度的悬浮RDX炸药粉尘中爆轰波的发展与传播过程,得到了爆轰波流场中气-固两相的物理量分布,并确定了爆轰波参数。在较低的RDX粉尘浓度条件下,爆轰波阵面压力的峰值曲线出现振荡。当RDX粉尘浓度在80~150 g/m3时,数值模拟得到的爆轰波阵面压力峰值曲线的振荡是规则的;当RDX粉尘浓度为70 g/m3时,爆轰波阵面压力峰值曲线出现不规则振荡。     

Numerical Simulation of Explosive Dust Detonation with CE/SE Method

采用CE/SE方法数值模拟悬浮在空气中的RDX炸药粉尘的两相爆轰过程.炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,释放能量支持爆轰波传播.数值模拟爆轰波管中的粉尘爆轰,得到爆轰波流场中的物理量分布,确定爆轰参数,数值结果与文献符合较好.数值模拟复杂通道中的炸药粉尘爆轰,预测了爆轰波的发展和传播过程以及爆轰波后的流场演化.数值结果表明CE/SE方法能成功模拟气体-固体两相爆轰,为粉尘爆轰的研究提供了新的数值预测手段.     

CE/SE方法数值模拟炸药粉尘爆轰

采用CE/SE方法数值模拟悬浮在空气中的RDX炸药粉尘的两相爆轰过程.炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,释放能量支持爆轰波传播.数值模拟爆轰波管中的粉尘爆轰,得到爆轰波流场中的物理量分布,确定爆轰参数,数值结果与文献符合较好.数值模拟复杂通道中的炸药粉尘爆轰,预测了爆轰波的发展和传播过程以及爆轰波后的流场演化.数值结果表明CE/SE方法能成功模拟气体-固体两相爆轰,为粉尘爆轰的研究提供了新的数值预测手段.     

化学反应和球面聚心爆轰波的相互作用

考察球面爆轰波聚心传播过程中,波阵面附近压力和温度不断升高引起化学反应进程的改变;对比氢氧可燃气体与氮气的数值模拟结果,分析化学反应对波面温度和压力的影响,从而考察Zeldovich理论预测聚心爆轰波后参数的精确性.数值结果表明,爆轰波聚心传播初期,放热的燃烧反应对波后热力学参数起主导作用;传播后期,波阵面趋近于对称中心时,吸热的气体解离反应变得非常活跃,解离反应对后期的汇聚压力影响不大,但会在很大程度上限制汇聚温度的升高.     

来流温度影响驻定爆轰波结构和性能的数值研究

 采用数值模拟的方法研究了超音速来流撞击圆锥体诱导的驻定爆轰波的结构和性能,讨论了不同初始来流温度下,爆轰波结构和波后压力的变化。结果表明,随着来流初始温度的降低,爆轰波的增压比逐渐增加,有利于提高驻定爆轰发动机的推进性能。当初始温度较高时,爆轰波阵面呈现光滑平直的形态,而初始温度较低时,爆轰波呈现出明显的三波点结构,且温度越低,三波点数量越多排列越紧密,这些具有较高压力的三波点对爆轰波波后压力的变化起到了重要作用。     

激光支持等离子体爆轰波温度的实验测量

 高能量激光聚焦空气产生等离子体,等离子体进一步吸收激光能量会形成激光支持等离子体爆轰波。等离子体爆轰波温度是表征爆轰波的一个重要参数,研究等离子体爆轰波温度对于深入了解激光支持等离子体爆轰波形成机理有重要意义。分析了激光聚焦空气形成等离子体爆轰波过程和影响等离子体爆轰波温度的主要因素。采用多通道瞬态光学高温计,测量了不同激光发射能量下空气中形成的激光支持等离子体爆轰波的辐射强度,获得了一系列等离子体爆轰波温度动态变化曲线。测量结果表明：等离子爆轰波温度在随时间演化过程中出现3个峰,最高温度在7 000～10 000 K范围内;激光能量与等离子体爆轰波温度没有明显的相关性。     

空气中激光支持爆轰波实验及理论分析

为了研究激光击穿空气产生的等离子体爆轰波形成机制和传播规律,利用高能量CO2激光器产生强激光,进行了空气中产生激光支持等离子体爆轰波实验。实验中:设置了诱导靶板,用于诱发和定位空气中的激光支持爆轰波;以激光器升压过程球隙放电产生的光信号作为触发源,触发高时间分辨率(纳秒级)的高速相机,记录了激光支持爆轰波的成长和传播全过程。分析了激光支持爆轰波的形成机理和传播规律。采用C-J爆轰理论,计算了激光支持爆轰波的压力和温度。研究结果表明:激光支持等离子体爆轰波形成初期,等离子体爆轰波发光体为球形;随着时间增加,等离子体爆轰波发光体的形状类似流星,且头部为等离子体前沿吸收层,亮度较高,而尾部等离子体温度较低,亮度较弱。等离子体爆轰波高速向激光源的方向移动,爆轰波速度高达18 km/s,温度约为107K。随着激光强度的减弱,爆轰波速度迅速按指数规律衰减,当爆轰波吸收的激光能量不能有效支持爆轰波传播时,爆轰波转变为冲击波。     

起爆药爆轰场和燃烧场的干涉法显示及其图像处理

本文介绍了利用脉冲YAG激光引爆起爆药产生爆轰流场或燃烧温度场,用时间序列干涉法首次获得了爆轰流场和起始燃烧过程的时间序列干涉图的实验结果.在此基础上,采用条纹曲线检测法和双三次多结点插值样条函数从干涉图上提取爆轰场冲击波的波阵面,并通过序列波阵面配准算法,定量地计算了冲击波的传播速度.同时,还定量计算了燃烧场温度随时间变化的分布.     

猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射

 采用高速转镜分幅相机和电探针技术研究了猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射图像和不同位置上的爆轰波传播时间。从研究得出，两种炸药都在拐角顶点附近绕射，爆轰波传播时间增长，爆速变小。但是两种炸药绕射爆轰波的状态不一样，钝感炸药IHE-2中爆轰波绕过直角时，在拐角顶点附近约10 mm范围内炸药未完全反应，猛炸药RHT-901中爆轰波绕过直角时未出现类似现象。两者相比，钝感炸药中绕射爆轰波速度变化大，波阵面曲率半径小，而猛炸药的绕射爆轰波速度变化小，波阵面曲率半径大。这说明炸药的爆轰波绕射与炸药的冲击感度、反应区宽度有关。     

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

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

激光驱动冲击波在铝-金阻抗匹配靶中的传播稳定性

 针对神光Ⅱ第九路激光条件，利用1维JB程序和多台阶靶技术对冲击波在铝 金阻抗匹配靶中的传播稳定性进行了理论和实验研究，实验结果与理论结果基本吻合。结果表明，激光驱动冲击波在铝 金阻抗匹配靶中传播时，高阻抗待测材料金中的冲击波稳定传播最大距离急剧减小。因此，在进行铝-金阻抗匹配靶物理参数设计时，应保证高阻抗材料金台阶厚度满足冲击波传播稳定性，然后再按照阻抗匹配实验中两种材料的冲击波速度比来确定低阻抗标准材料铝的台阶厚度。根据神光Ⅱ第九路激光条件，铝-金阻抗匹配靶中铝基底厚度选取为30 μm左右较好，金台阶和铝台阶厚度应分别小于10 μm和17 μm。     

Deflagration-to-detonation transition in mixtures of finely divided ammonium perchlorate with submicron aluminum powder

Deflagration-to-detonation transition in binary mixtures of fine ammonium perchlorate (20-μm grains) with submicron ALEX-L aluminum powder (0.2-μm particles) is studied using high-speed photography and pressure recording with quartz crystal sensors. The test mixtures were loaded in thin-walled quartz tubes of inner diameter 10 mm. The charges had a porosity of ~50%. It has been shown that, even under very mild conditions (low-strength shell and a weak source of initiation), the deflagration mode of mixture combustion easily transforms into the detonation mode. The shortest length of the region of transition from deflagration to normal detonation (not more than 30 mm) was observed for a lean mixture, with an aluminum content of ~5%. The mechanism of transition to detonation involves the stage of convective combustion, resulting in the formation of a brightly luminescent crescent-shaped area behind the primary flame front, which, in turn, generates a forward (in the direction of propagation) and a backward wave. The forward wave gives rise to low-speed detonation, which later transforms into normal detonation. The pressure profile within the region of low speed detonation is measured. A comparison with similar experiments in which ALEX-L alu- minum powder was replaced by ASD-4 aluminum (4 μm particles) shows that ALEX-L sensitizes the mixture, resulting in a dramatic reduction of the length of the transition region, making it possible to produce normal detonation in low-strength shells.     

Nonequilibrium Zeldovich-von Neumann-Doring theory and reactive flow modeling of detonation

This paper discusses the Nonequilibrium Zeldovich-von Neumann-Doring (NEZND) theory of self-sustaining detonation waves and the Ignition and Growth reactive flow model of shock initiation and detonation wave propagation in solid explosives. The NEZND theory identified the nonequilibrium excitation processes that precede and follow the exothermic decomposition of a large high explosive molecule into several small reaction product molecules. The thermal energy deposited by the leading shock wave must be distributed to the vibrational modes of the explosive molecule before chemical reactions can occur. The induction time for the onset of the initial endothermic reactions can be calculated using high pressure-high temperature transition state theory. Since the chemical energy is released well behind the leading shock front of a detonation wave, a physical mechanism is required for this chemical energy to reinforce the leading shock front and maintain its overall constant velocity. This mechanism is the amplification of pressure wavelets in the reaction zone by the process of de-excitation of the initially highly vibrationally excited reaction product molecules. This process leads to the development of the three-dimensional structure of detonation waves observed for all explosives. For practical predictions of shock initiation and detonation in hydrodynamic codes, phenomenological reactive flow models have been developed. The Ignition and Growth reactive flow model of shock initiation and detonation in solid explosives has been very successful in describing the overall flow measured by embedded gauges and laser interferometry. This reactive flow model uses pressure and compression dependent reaction rates, because time-resolved experimental temperature data is not yet available. Since all chemical reaction rates are ultimately controlled by temperature, the next generation of reactive flow models will use temperature dependent reaction rates. Progress on a statistical hot spot ignition and growth reactive flow model with multistep Arrhenius chemical reaction pathways is discussed. The text was submitted by the authors in English.     

表面放电等离子体产生冲击波特性

采用压力传感器,分析XeF(C→A)激光器高电压表面放电光泵浦方式不同位置、不同放电电压条件下冲击波特性。研究结果表明,表面放电等离子体膨胀产生的冲击波持续时间为ms量级,在冲击波过后伴随有相对较弱的压力波。随着充电电压增加,放电过程释放能量增多,冲击波峰值压强也变大,最大压强可达1 MPa,同时冲击波峰值延后,冲击波持续时间变长。     

Time resolved laser-induced plasma dynamics

The structure and evolution of the laser-induced vapor plume and shockwave were measured from femtosecond time resolved shadowgraph images. By changing the wavelength of the probe beam (400 and 800 nm), differences in the opacity of the vapor plume were measured as a function of delay time from the ablation laser pulse. The evolution of the temperature and electron number density during and after the ablation laser pulse were determined and compared for ablation in argon and helium background gases. A laser supported detonation wave (LSD) observed for ablation in argon, blocks the incoming laser energy and generates a high-pressure region above the vapor plume.     

Shadowgraph investigation of plasma shock wave evolution from Al target under 355-nm laser ablation

The propagation of a plasma shock wave generated from an Al target surface ablated by a nanosecond Nd:YAG laser operating at 355 nm in air is investigated at the different focusing positions of the laser beam by using a time-resolved shadowgraph imaging technique. The results show that in the case of a target surface set at the off-focus position, the condition of the focal point behind or in front of the target surface greatly influences the evolution of an Al plasma shock wave, and an ionization channel forms in the case of the focal point set in front of the target surface. Moreover, it is found that the shadowgraph with the evolution time around 100 ns shows that a protrusion appears at the front tip of the shock wave if the focal point is at the target surface. In addition, the calculated results of the expanding velocity of the shock wave front, the mass density, and pressure just behind the shock wave front are presented based on the shadowgraphs.     

Numerical studies on autoignition and detonation development from a hot spot in hydrogen/air mixtures

Detonation development inside spark ignition engines can result in the so called super-knock with extremely high pressure oscillation above 200?atm. In this study, numerical simulations of autoignitive reaction front propagation in hydrogen/air mixtures are conducted and the detonation development regime is investigated. A hot spot with linear temperature distribution is used to induce autoignitive reaction front propagation. With the change of temperature gradient or hot spot size, three typical autoignition reaction front modes are identified: supersonic reaction front; detonation development and subsonic reaction front. The effects of initial pressure, initial temperature, fuel type and equivalence ratio on detonation development regime are examined. It is found that the detonation development regime strongly depends on mixture composition (fuel and equivalence ratio) and thermal conditions (initial pressure and temperature). Therefore, to achieve the quantitative prediction of super-knock in engines, we need use the detonation development regime for specific fuel at specific initial temperature, initial pressure, and equivalence ratio.     

Blast waves in a cylindrical channel generated by the nonideal detonation of aluminum-Teflon-RDX high-density formulations

The pressure at the front and the pressure impulse of blast waves generated in a cylindrical tube by the expanding products of the nonideal detonation of low-porosity charges prepared by pressing of fine-grained powders of aluminum, Teflon, and RDX were measured. The measured parameters are compared to the same parameters of blast waves produced by the detonation of TNT charges of identical mass. The relative quantities were used to evaluate the effectiveness of blast waves with respect to those generated by TNT. Mixed compositions differing in the shape (brand) of the aluminum powder particles and the ratio between the components at 30% RDX are studied. It is shown that, for the investigated compositions, the pressure at the leading front of the wave exceeds the pressure achieved during TNT explosion on average by 10–30%, almost independently of the distance traveled along the tube in the range from 0.8 to 3.8 m. The dependence of the wave amplitude on the particle shape and aluminum content was weak. In the same range of distances, the relative pulse pressure increases strongly, from 0.5 to 2.1 and higher, mainly due to an increase in the width of the wave. This result is of interest from the point of view of achieving a high pressure impulse of the blast wave in an area remote from the charge. The obtained data suggest that RDX mainly reacts in the detonation wave, with the chemical transformation of Teflon and aluminum in the detonation wave and near-to-charge zone occurring, if at all, to a small extent. On the contrary, as the blast wave front moves through the channel, the burning of aluminum in the fluoride formed during the decomposition of Teflon provides an appreciable support to the blast wave, causing a significant increase in the pressure impulse.