FAE装置参数对燃料抛撒与爆炸威力影响的实验研究

 采用高速运动分析系统对FAE实验装置爆炸抛撒过程进行观测,描述了燃料抛撒过程的不同阶段,实验研究了比药量、长径比、壳体材质等装置参数对燃料抛撒与爆炸威力的影响。结果表明,燃料抛撒过程可分为射流形成与扩散运动阶段、燃料两向膨胀运动阶段和气液融合运动阶段。不同阶段对应不同的流体动力学特征,对云团形成的贡献不同。在实验装置总体优化条件下,适当增大比药量可提高云团覆盖面积与体积;在实验范围内长径比不是影响云团状态的显著因素,但长径比较大时可使燃料抛撒均匀性更好;采用钢质壳体时云雾抛撒状态明显优于铝质壳体。实验证明,采用碳钢壳体、比药量3%左右、长径比为3~5且装置参数良好匹配时,可获得理想的云团状态和高威力爆炸波毁伤效应。     

燃料抛撒成雾及其燃烧爆炸的光滑离散颗粒流体动力学方法数值模拟研究

燃料在炸药爆炸驱动下形成燃料空气爆炸云团, 进而引燃爆炸, 对目标造成毁伤. 本文在前期提出的光滑离散颗粒流体动力学方法(SDPH)的基础上, 引入描述炸药由爆轰到膨胀整个过程的Jones-Wilkins-Lee状态方程及描述气体快速燃烧过程的EBU-Arrhenius燃烧模型, 建立了求解战斗部起爆、燃料抛撒和燃料二次引燃爆炸问题的新型SDPH方法. 设计了圆环形燃料颗粒在炸药爆炸驱动下运动抛撒的算例进行数值验证, 结果与理论相符; 对燃料空气炸药(FAE)云雾的形成和发展过程进行了数值模拟, 分析了云雾的形态, 并与实验结果进行对比, 符合较好, 同时分析了不同起爆方式对云雾团成型的影响; 最后, 在云雾团成型的基础上, 引入蒸发燃烧模型对FAE的燃烧爆炸过程进行了模拟研究. 结果表明, 本文建立的数学模型和计算方法可以较好的模拟燃料空气炸药抛撒成雾及云雾燃烧爆炸过程, 为该类武器装备的设计研究提供了较好的数值方法.     

炸药爆炸作用下液体破碎后颗粒尺寸分布的研究

 对液体抛撒的液滴尺寸进行研究在军事和民用上是很重要的,国内刚开始使用激光散射仪开展此项研究工作。利用R. A. Dobbins等人的液体颗粒测量技术,研制了一套既简单又实用的测量液体抛撒过程中液滴尺寸的实验装置——激光散射仪。对于激光与液体微粒的相互作用,当微粒的反射与折射和吸收效应可被忽略时,可导出液体微粒对激光散射的光强公式。只要测量激光被微粒散射的光强,就可推算出微粒的Sauter平均直径。在使用激光散射仪测量液体抛撒液滴尺寸的实验中,用水代替爆炸抛撒液体,测量结果表明：液体抛撒二次破碎中,在固定位置测量到的云雾区液滴Sauter平均直径随测量时间的增加呈现出减小的趋势;而云雾区的宽度则随着与抛撒中心距离的增大而呈现出增加的趋势;云雾区前沿的液滴Sauter平均直径随着与抛撒中心距离的增加而呈现出先逐渐增大然后迅速减小的趋势。为便于比较,对燃料抛撒二次破碎进行了回收法测量和数值模拟计算,其测量与计算结果与用激光散射仪测量的结果有较好的一致性。     

FAE燃料配方和工艺研究

设计了应用于二次点火体系的含铝燃料配方和工艺，通过5L燃料抛撒装置的静态爆轰试验，研究了铝粉含量和形状对爆轰威力的影响。试验结果表明：(1)相对于采用球状铝粉的燃料配方，采用片状铝粉的燃料配方具有引爆延时范围宽、起爆可靠性好、爆轰威力高等优点。根据爆轰过程的录像，我们对试验结果进行了分析。球状铝粉在抛撒分散过程中具有较高的初始速度，减速形成紊流云团的时间较短、直径较大，因而难于找到最佳爆轰浓度形成的时机；片状铝粉在抛撒分散过程中的初始速度相对较低，减速形成紊流云团的时间较长、直径较小，因而易于找到最佳爆轰浓度形成的时机。不过，采用球状铝粉的燃料配方，其抛撒分散形成云团的形状较规则、直径较大，估计在大装药量的装置中有一定的应用前景。(2)相对于采用片状铝粉的燃料配方，采用球状铝粉的燃料配方具有很好的工艺性能。其原因在于，球状铝粉的形状规则、对称，颗粒相互之间的接触面积较小，     

铜丝水中电爆炸能量沉积特性

对s脉冲电压作用下铜丝水中电爆炸的能量沉积过程进行了实验研究,利用自积分Rogowski线圈和电阻分压器分别测量铜丝电爆炸时的电流和电压。利用测量电压波形确定了熔融起始、熔融结束、汽化起始和击穿时刻点,将铜丝电爆炸划分成熔融、液态和汽化3个阶段。通过数学方法计算了3个阶段和击穿前的沉积总能量。通过实验和计算,分析了电路参数,包括放电电压和回路电感,以及铜丝特性,包括铜丝长度和直径,对铜丝电爆炸过程中3个阶段和击穿前沉积总能量的影响。结果表明:在s脉冲电压作用下,放电电压、回路电感、铜丝长度和直径对熔融阶段能量沉积影响较小,但对液态和汽化阶段能量沉积影响较大,通过调节电路参数提高电流上升速率,可以显著提高汽化和击穿前的沉积能量。     

发射过程中混合燃料应力状态的数值分析

 炸药装药在发射条件下的安全性问题已成为高能量炸药应用的一个技术障碍，引起了各国的广泛重视。燃料空气炸药（FAE）战斗部是近年来倍受国内外关注的新概念武器。目前关于FAE装药的发射安全性特别是高过载环境下的发射安全性研究还不系统、不够深入，也没有给予足够的重视。运用爆炸力学程序Object MMIC程序分析了FAE战斗部发射过程中的应力和应力率变化规律。分析结果表明，FAE战斗部发射过程中受到的应力和应力率不同于传统的发射过程，因此，有必要建立适用于FAE武器的全新的发射安全性理论。这些研究对提高FAE战斗部发射安全性提供了依据。     

柴油的抛撒成雾及燃爆特性研究

为了研究柴油、阻燃柴油、阻燃抑爆柴油的抛撒成雾及燃爆特性,使用高速摄像系统记录云雾形成过程,采用红外热成像系统记录火球演变情况。研究表明:在中心药柱爆轰驱动下,油料液滴与空气快速混合形成云雾,其尺寸随油料运动黏度的增大而减小;3种柴油形成的云雾均在二次引爆后产生火球,阻燃柴油的火球最高温度较柴油低35%,阻燃抑爆柴油较阻燃柴油低23%;阻燃柴油的火球最大直径和高度较柴油小54%、42%,阻燃抑爆柴油较阻燃柴油小46%、55%;火球持续时间,阻燃柴油较柴油短38%,阻燃抑爆柴油较阻燃柴油短47%。说明在抑制火球温度上升和尺寸增长的效果上,阻燃抑爆柴油优于阻燃柴油,阻燃柴油优于柴油。     

爆炸驱动动能杆定向抛撒散布试验

动能反导技术作为新型的反导技术，原理是在拦截弹和目标的相对速度方向上，抛射出大量动能杆条毁伤目标。其核心技术之一就是动能元素的抛撒，特别是定向抛撒。国外研究一般集中针对既定抛撒装置结构的杆条抛撒结果验证性报道，缺乏对杆条抛撒影响因素的详细研究。本文提出了采用周向弧形状药选择区域爆炸驱动中心杆条束定向侧向飞散的动能扦定向抛撒装置方案，设计不同结构参数的杆条定向抛撒装置的试验，研究抛撒结构参数对动能杆的驱动特性影响规律，为新概念反导战斗部设计提供技术帮助。     

真空中金属丝电爆炸数值模拟研究

利用金属丝电爆炸物理数学模型对电爆炸物理过程开展了数值模拟，研究了不同直径铝丝电爆炸特性，进一步分析了金属丝内沉积能量、电压击穿时刻、电压峰值随金属丝直径的变化规律，并与相关实验数据作了对比。     

云爆弹外壳体厚度对威力的影响规律研究

通过对不同厚度壳体的速度和动能进行理论分析,给出了壳体对燃料抛撒的影响规律;同时开展了不同厚度壳体的云爆装置试验,得到了装置威力与壳体厚度的关系.研究结果表明:采用薄壳体,并加强壳体与端盖的连接强度,可以提高中心抛撒药能量的利用率,进而提高一次起爆型燃料空气弹(SEFAE)的威力,试验结果与理论分析结果一致.     

超高速发射实验模型的数值计算

 采用三阶精度PPM（Parabolic Piecewise Method）方法和VOF（Volume of Fluid）相结合,运用Lagrange-Remapping算法,编制了多介质流体高精度欧拉计算程序MFPPM,可以对多介质复杂流场进行数值计算。对程序计算精度进行了测试,并应用于超高速发射实验模型数值计算,对Sandia实验室不带会聚作用的实验模型进行了一维计算以及带会聚作用的实验模型进行二维计算并与其实验结果和CTH计算结果进行对比,其中一维计算最大相对误差1%,二维计算相对误差4.7%,在此基础上对超高速发射设计模型进行了初步计算。     

Thermal explosion in a droplet-gas cloud

The effect of the presence of a spray of liquid fuel on thermal explosion in a combustible droplet-gas cloud is investigated. By ‘thermal explosion’ we refer exclusively to the initial stages of the behaviour of the combustible medium as its temperature begins to rise and various competing physical and chemical processes are called into play. A qualitative analysis of the system of governing equations is carried out using an advanced geometrical asymptotic technique (the integral manifold method). Possible types of dynamical behaviour of the system are classified and parametric regions of their existence are determined analytically. It is demonstrated that the original problem can be decomposed into two subproblems, due to the underlying hierarchical time scale structure. The first subproblem relates to the droplet heat up period, for which a relatively rapid time scale is applicable. The second subproblem begins at the saturation point. For the latter, more significant second stage, it is found that there are five main dynamical regimes: slow regimes, conventional fast explosive regimes, thermal explosion with freeze delay and two different types of thermal explosion with delay (the concentration of the combustible gas decreases or increases). Upper and lower bounds for the delay time are derived analytically and compared with results of numerical simulations, with rather satisfactory agreement.     

Experimental investigation and numerical validation of explosion suppression by inert particles in large-scale duct

A large-scale duct with an explosion suppressor was designed to investigate experimentally the explosion suppression by inert particles for a CH₄/O₂/N₂ mixture. The duct is 25 m long and has an internal diameter of 700 mm. Pressure and flame signals were recorded some distance away from ignitor in the duct. Pressure tracking lines of the shock front for the different inert particle cloud densities and the inert particle diameters were made. The measured results indicate that the shock front is decoupled from the flame front in the inert particle cloud, which leads to a suppression of explosion. Also, the experiments suggest that increasing the inert particle cloud density or decreasing the inert particle diameter can enhance the ability to suppress explosion. For the purpose of validation, a two-dimensional numerical model coupled with the element chemical reaction mechanism for the simulations of the CH₄/O₂/N₂ mixture explosion suppression by the inert particles has been developed. This model makes use of the second-order TVD scheme and the MacCormack scheme to calculate gas-phase and particle-phase equations, respectively. The Strang splitting technique is used to treat the stiffness due to the coupling of the governing equations, while the implicit Gear algorithm is used to treat the stiffness due to the chemical reactions. The effect of inert particle cloud density on explosion suppression was investigated using the model. The calculated results indicate that the accumulation of inert particles slows the propagation of the gas-phase shock front and results in explosion suppression. With increased inert particle cloud density, the explosion suppression is more prominent. The calculated results show a qualitative agreement with the measured results in the large-scale duct experiment.     

Investigation of contact line dynamics under a vapor bubble at boiling on the transparent heater

The paper presents the results of an experimental study of dynamics of vapor bubble growth and departure at pool boiling, obtained with the use of high-speed video recording and IR thermography. The study was carried out at saturated water boiling under the atmospheric pressure in the range of heat fluxes of 30?150 kW/m². To visualize the process and determine the growth rates of the outer bubble diameter, microlayer region and dry spot area, transpa-rent thin film heater with the thickness of 1 μm deposited on sapphire substrate was used in the experiments, and video recording was performed from the bottom side of the heating surface. To study integral heat transfer as well as local non-stationary thermal characteristics, high-speed infrared thermography with a frequency of up to 1000 FPS was used. High-speed video recording showed that after formation of vapor bubble and microlayer region, dry spot appears in a short time (up to 1 ms) under the vapor bubble. Various stages of contact line boundary propagation were ob-served. It was shown that at the initial stage before the development of small-scale perturbations, the dry spot propaga-tion rate is constant. It was also showed that the bubble departure stage begins after complete evaporation of liquid in the microlayer region.     

Quantitative, dynamic fuel distribution measurements in combustion-related devices using laser-induced fluorescence imaging of biacetyl in iso-octane

Knowledge of in-situ fuel distributions in practical combustion devices, such as internal combustion engines, is crucial for research and devlopment purposes. Numerous imaging techniques, mostly based on laser-induced fluorescence (LIF), have been developed and yield high levels of 2-D spatial information, but generally lack the temporal resolution (frame rates) necessary to resolve important timescales at sub-millisecond levels for sustained times. A planar LIF technique for quantitatively visualizing fuel distribution is presented which gives not only high spatial resolution, but also high temporal resolution. Using a high-speed CMOS camera, a lens-coupled image intensifier, and frequency-tripled diode-pumped Nd:YAG laser allows for capturing LIF images of biacetyl that is used as a fluorescence tracer at 12 kHz (one crank-angle resolution at 2000 RPM) for hundreds of consecutive engine cycles. The LIF signal strength of biacetyl doped in iso-octane is shown to vary substantially over a wide range of temperatures and pressures. The low absorption coefficient at 355 nm and a longpass filter in the detection path exclude bias errors due to laser beam attenuation and fluorescence trapping. An intensifier gate time of 350 ns is shown to suppress the detection of phosphorescence signals under practical conditions. An example for a quantitative high-speed measurement of fuel concentration at varying pressure and temperature conditions is presented. Quantitative equivalence ratio maps are shown for the fuel injection event within a single cycle in a spark-ignition direct-injected engine, showing the ability of the technique to not only reveal static fuel concentration maps, but also the motion of the fuel cloud along with very steep gradients. Spray velocities determined from the moving fuel cloud are in agreement with previous particle image velocimetry measurements.     

Evolution of temperature of a droplet of liquid composite fuel interacting with heated airflow

The macroscopic patterns of a temperature change at the center of a droplet of three-component (coal, water, petroleum) composite liquid fuel (CLF) were studied using a low-inertia thermoelectric converter and system of high-speed (up to 10⁵ frames per second) video recording during the induction period at different heating intensity by the air flow with variable parameters: temperature of 670?870 K and motion velocity of 1?4 m/s. The studies were carried out for two groups of CLF compositions: fuel based on brown coal and coal cleaning rejects (filter cake). To assess the effect of liquid combustible component of CLF on characteristics of the ignition process, the corresponding composition of two-component coal-water fuel (CWF) was studied. The stages of inert heating of CLF and CWF droplets with characteristic size corresponding to radius of 0.75?1.5 mm, evaporation of moisture and liquid oil (for CLF), thermal decomposition of the organic part of coal, gas mixture ignition, and carbon burnout were identified. Regularities of changes in the temperature of CLF and CWF droplets at each of identified stages were identified for the cooccurrence of phase transitions and chemical reactions. Comparative analysis of the times of ignition delay and complete combustion of the droplets of examined fuel compositions was performed with varying droplet dimensions, temperatures, and oxidant flow velocity.     

Combustion effects in confined explosions

Results of shock-dispersed-fuel (SDF) explosion experiments are presented. The SDF charge consisted of a spherical 0.5-g PETN booster surrounded by 1 g of fuel, either flake aluminum (Al) powder or TNT. The charge was placed at the center of a sealed chamber. Three cylindrical chambers (volumes of 6.6, 20, and 40 l with L/D = 1) and three tunnels (L/D = 3.8, 4.65, and 12.5) were used to explore the influence of chamber volume and geometry on completeness of combustion. Detonation of the SDF charge created an expanding cloud of explosion product gases and hot fuel (Al or TNT). When this fuel mixed with air, it formed a turbulent combustion cloud that consumed the fuel and liberated additional energy (31 kJ/g for Al or 15 kJ/g for TNT) over and above detonation of the booster (6 kJ/g) that created the explosion. Static pressure gauges were the main diagnostic. Pressure and impulse histories for explosions in air were much greater than those recorded for explosions in nitrogen—thereby demonstrating that combustion has a dramatic effect on the chamber pressure. This effect increases as the confinement volume decreases and the excess air ratio approaches values between 2 and 3.5.