Integral inequalities in the dynamics of a gravitating gas

In the framework of Newtonian mechanics, a study is made of the spherically symmetric problem of the adiabatic motion of a gravitating perfect gas in the presence of a shock wave produced by inhomogeneous gravitational collapse or a point explosion. The method of Golubyatnikov [1, 2] is used to construct a system of integro—differential inequalities that determine, in particular, the law of motion of the shock wave if the initial state of the gas is known. The investigated examples include some self-similar and nearly self-similar solutions to the problem of the gravitational contraction of dust with the formation of a strong shock wave, possibly with the release of energy; the self-similar problem of a point explosion in a gas at rest; and also the problem of the equilibrium of a gas sphere for =4/3 and arbitrary distribution of the entropy. In these cases, the inequalities reduce to algebraic relations and can be solved numerically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 169–173, September–October, 1980.     

The detonation burning of a fuel mixture resulting from a double explosion

A study is made of the propagation of a multifront detonation burning in a fuel mixture consisting of a gaseous fuel and an oxidant with additions of combustible solid or liquid particles arising as a result of a double point explosion. In such combustible media it is possible for there to be propagation of several detonation or burning fronts following one after the other. The easily igniting gaseous fuel burns in the first detonation wave, which propagates in the gaseous mixture with particles which are heated by the products of the explosion, ignite and burn in the second detonation wave or in the flame front. Self-similar regimes of propagation of such waves in an idealized formulation were studied in [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–131, March–April, 1985.     

Strong explosion in a combustible gas mixture

Assume that a planar, cylindrical, or spherical point explosion takes place in a combustible mixture of gases. As a result of the explosion a strong shock wave develops and triggers chemical reactions with the release of heat. The solution of the problem for the case in which the thickness of the heat release zone is neglected (the infinitely thin detonation wave model) was obtained in [1–3].It was emphasized in [4] that these solutions can be considered only as asymptotic solutions for time and distance scales which are large in comparison with the scales which are characteristic for the chemical reactions, and under the assumption that as the overdriven detonation wave which is formed in the explosion is weakened by the rarefaction waves it does not degenerate into an ordinary compression shock. Here the question remains open of the possibility of obtaining such asymptotic solutions with account for finite chemical-reaction rates.In conclusion the authors wish to thank E. Bishimov for carrying out most of the computations for this study.     

On regular reflection of detonation waves

The boundaries of regular reflection of detonation waves by a rigid wall are calculated. It is assumed that detonation is initiated at the point of reflection when a shock wave is incident on the wall at a finite angle in a gas fuel mixture, the detonation propagating instantaneously along the reflected front.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 178–180, March–April, 1983.     

Breakdown of a discontinuity,taking account of the real equation of state of the detonation products from a condensed explosive

The problem is considered concerning the breakdown of an arbitrary discontinuity, due to the emergence of a detonation wave at the boundary of a condensed explosive charge. The real equation of state of the detonation products of Hexogen was used in the numerical calculations. A u vs p diagram is constructed, which allows graphical calculations to be carried out of the discontinuity breakdown for different media. A comparison is carried out of the calculated values of the initial shock-wave velocities with the experimental data obtained at a certain distance from the explosive charge. It is shown that an increase of the pressure of the gas in which dispersion of the explosion products occurs leads to a reduction of the initial shock wave velocity and to an enhancement of its attenuation during its further motion in the shock tube.     

切缝药包爆炸作用机理数值模拟

为研究切缝药包定向控制爆炸机理,基于双缝/耦合切缝药包爆轰行为高速纹影实验,建立"炸药-切缝管-空气"模型,采用数值模拟方法研究切缝药包爆炸过程中冲击波相互作用、爆炸流场压力时空分布和切缝管形态变化。研究结果表明:切缝药包爆炸过程中切缝管能够有效控制爆炸能量释放和爆生气体动力学行为;切缝方向压力超前且高于非切缝方向;在爆炸冲击波及爆生气体共同作用下,切缝管曲率不断变小且从起爆点处以相同的变形特征沿切缝管轴向发展;数值模拟结果与实验结果吻合。     

Numerical analysis of oscillations of a gaseous cavity in a compressible liquid

The radial oscillations of a gas-filled spherical cavity in an infinite compressible ideal liquid are investigated. The numerical solution of the equation of motion of a gaseous sphere is used for the comparison of acoustic and energy characteristics of underwater trotyl explosion and small-amplitude explosion sound sources (detonation of explosive gaseous mixture, expansion of a compressed inert gas).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 176–182, March–April, 1971.     

一种水下爆炸冲击波压力调控方法

起爆位置和装药形状对水下爆炸冲击波压力有较为显著的影响,这使得利用小当量装药在局部方向形成与大当量装药一定程度等效的冲击波成为可能。为了能够在小当量装药条件下开展舰船结构及设备抗水下爆炸冲击实验,基于细长装药结构和参数优化设计,设计了一种冲击波压力幅值和持续时间可调的装药方法。首先,基于简单波理论给出了水下爆炸冲击波压力调控的原理,以及装药参数优化设计的目标函数和约束条件;然后,采用自主数值模拟软件研究了细长装药的水下爆炸能量输出规律,通过实验验证了数值模拟的置信度,研究发现起爆位置和装药形状对水下爆炸冲击波压力峰值和持续时间的影响是显著的,在炸药爆速一定的情况下,长药柱水下爆炸冲击波压力的持续时间可通过几何近似确定;最后,为了进一步考察该方法的有效性,以1 000 kg TNT和100 m爆距的水下爆炸冲击波压力-时间曲线作为原型,设计了2种与该原型冲击波压力等效的装药方案,并通过数值模拟进行了验证。研究结果表明:设计的装药能够在预定的持续时间内,在装药起爆端一侧形成与原型等效的冲击波压力-时间曲线。由于没有考虑对气泡载荷的等效,因此该调控方法仅适用于中远场爆炸冲击问题。     

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟。研究结果表明，标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波，爆轰波与激波和接触间断相互作用产生了复杂的波系结构；爆轰波爆点位置在对称轴上并不是固定的点，而是随着初始激波马赫数的变化而发生移动；可燃气体初始温度和压力对起爆临界马赫数都有影响，但是初始温度的影响大得多。     

锥形长药柱水下爆炸冲击波参数计算方法

为了计算锥形长药柱水下爆炸冲击波压力,以及研究长脉宽冲击波的传输特性,基于叠加原理建立了冲击波压力-时间曲线的计算方法,通过实验验证了该方法的有效性,在此基础上分析了锥形长药柱不同方位冲击波压力的分布规律。研究结果表明:锥形长药柱产生的冲击波压力具有各向异性,在起爆端一侧形成的是具有厚波头特征的低幅值长脉宽冲击波,在装药径向形成的是接近指数衰减的高幅值冲击波,而在远离起爆端的冲击波压力幅值和脉宽则介于前两者之间。锥形长药柱与球形装药冲击波分布的差异是由于装药形状和起爆方式的改变所导致的,由于装药不同部位起爆的时间差,导致水下爆炸冲击波在不同位置的叠加效果存在明显差异,药柱周围流场中形成的冲击波压力具有方向性。利用提出的计算方法得到的计算结果与实验结果和数值模拟结果吻合较好,研究结果可为锥形长药柱水下爆炸冲击波威力场和毁伤评估提供参考和依据。     

Detonation Wave Propagation in Rotational Gas Flows

This paper studies the propagation of detonation and shock waves in vortex gas flows, in which the initial pressure, density, and velocity are generally functions of the coordinate — the distance from the symmetry axis. Rotational axisymmetric flow having a transverse velocity component in addition to a nonuniform longitudinal velocity is considered. The possibility of propagation of Chapman–Jouguet detonation waves in rotating flows is analyzed. A necessary conditions for the existence of a Chapman–Jouguet wave is obtained.     

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

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

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

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

铝粉含量和粒度对CL-20含铝炸药水中爆炸反应特性的影响

为了研究CL-20基含铝炸药的爆炸反应机理,利用水中爆炸实验,测量了不同铝粉含量和粒度的CL-20炸药水中爆炸的冲击波参数、二次压力波参数,计算了冲击波能和气泡能。结果表明,水中爆炸的冲击波能和气泡能表征了爆轰和二次反应两个阶段的炸药爆炸能量分配,CL-20炸药中的铝粉主要在二次反应阶段发生反应,只有少部分的铝粉参与了早期的爆轰反应。气泡脉动形成的二次压力波能描述铝粉含量和粒度对二次反应过程的影响,铝粉含量对炸药的二次反应有显著的影响;铝粉粒度对炸药的水下爆炸的初始冲击波参数、冲击波能和气泡能的影响很小,对铝粉与爆轰产物的二次反应速率影响较大。     

小不耦合系数装药爆破孔壁压力峰值计算方法

基于轮廓爆破孔壁压力峰值计算方法的相关研究,充分考虑空气冲击波的传播与爆轰产物膨胀的过程,理论分析了小不耦合系数装药爆破过程中空气冲击波与炮孔壁的相互作用,建立了三维空气介质径向不耦合装药单孔爆破有限元模型,研究了工程爆破中常用的多种小不耦合系数装药组合工况下,炸药单点起爆后的炮孔壁压力峰值,并获得了相应工况下的孔壁压力峰值较爆生气体准静态等熵膨胀压力的压力增大倍数。结果表明:小不耦合系数装药爆破过程中,爆轰产物参数会对空气冲击波波后物质参数产生显著影响,揭示了小不耦合系数装药爆破与轮廓爆破在孔壁压力峰值计算方法上的本质差异;柱状装药结构爆轰波沿轴向传播使得空气冲击波撞击炮孔壁时存在叠加效应,孔壁压力峰值也相应增大,通过统计分析不同炸药类型、不同岩石类型工况下压力增大倍数与不耦合系数的关系,发现压力增大倍数随不耦合系数的增大近似呈线性增长;基于理论推导结果及常用爆破孔壁压力峰值计算形式,综合考虑炸药性能、孔壁岩石介质条件、不耦合装药系数对空气冲击波撞击炮孔壁后压力增大倍数的影响,提出了不耦合系数较小时爆破孔壁压力峰值计算方法。     

Population inversion behind a detonation wave propagating in a variable-density medium

The creation of an active medium by means of detonation has been investigated on a number of occasions. It was suggested that one could use the expansion of the detonation products of an acetylene-air mixture in vacuum [1] or the cooling of the detonation products of a mixture of hydrocarbons and air through a nozzle [2, 3]. In [4], the detonation of a solid high explosive was used to produce population inversion in the gas mixture CO₂-N₂-He(H₂O). Stimulated emission from HF molecules was observed in [5] behind the front of an overdriven detonation wave propagating in an F₂-H₂-Ar mixture in a shock tube. Population inversion behind a detonation wave was studied in H₂-F₂-He mixtures in [6–8] and in H₂-Cl₂-He mixtures in [9] with energy release on a plane and on a straight line in a medium with constant density. Similar problems were solved for shock waves propagating in both a homogeneous gaseous medium [7, 10] and in the supersonic part of an expanding nozzle. In the present paper, we study theoretically population inversion behind an overdriven detonation wave propagating in a mixture (fine carbon particles + acetylene + air) which flows through a hypersonic nozzle. The propagation of detonation in media with variable density and initial velocity was considered, for example, in [11, 12]. Analysis of the gas parameters behind a detonation wave propagating in a medium with constant density (for a given fuel) showed that the temperature difference across the detonation front is insufficient to produce population inversion of the vibrational levels of the CO₂ molecule.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 65–71, January–February, 1980.I am grateful to V. P. Korobeinikov for a helpful discussion of the results.     

Waves from an underground explosion

The problem of the propagation of a spherical detonation wave in water-saturated soil was solved in [1, 2] by using a model of a liquid porous multicomponent medium with bulk viscosity. Experiments show that soils which are not water saturated are solid porous multicomponent media having a viscosity, nonlinear bulk compression limit diagrams, and irreversible deformations. Taking account of these properties, and using the model in [2], we have solved the problem of the propagation of a spherical detonation wave from an underground explosion. The solution was obtained by computer, using the finite difference method [3]. The basic wave parameters were determined at various distances from the site of the explosion. The values obtained are in good agreement with experiment. Models of soils as viscous media which take account of the dependence of deformations on the rate of loading were proposed in [4–7] also. In [8] a model was proposed corresponding to a liquid multicomponent medium with a variable viscosity.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 34–41, May–June, 1984.     

点火能对丙烷-空气预混气体爆炸过程及管壁动态响应的影响

在长12 m的无缝不锈钢直管中,通过改变初始点火能量,探究了点火能对封闭管道内丙烷-空气混合气体爆炸传播特性和激波对管壁动态加载的影响。结果表明,初始点火能对预混气体爆炸火焰传播规律以及管壁的动态响应有显著影响:点火能越大,爆炸越剧烈,爆炸压力峰值压力和管壁最大应变就越大,且压力波和管壁应变的发展一致。火焰在传播过程中受到管道末端反射波的作用会发生短暂熄灭和复燃;管壁承受冲击波加载,应变信号主要分布在0~781.25 Hz,管壁最大应变率大于10-3 s-1,实验工况下管壁应变属动态响应。     

Coupling Relation Between Air Shockwave and High-Temperature Flow from Explosion of Methane in Air

After a methane-in-air explosion in a coal mine tunnel, a secondary explosion of coal dust is prone to happen. The shockwave in the gas explosion produces a coal dust suspension, and the peak temperature band may detonate that suspension. This secondary detonation depends on the space-time relation between the shockwave and the peak temperature band. This paper presents a methodology to estimate the coupling relation between the air shockwave and high-temperature flow from the explosion of methane in air. The commercial software package AutoReaGas was used to carry out the numerical simulation for the explosion processes of methane in air in the tunnel. Based on the numerical simulation and its analysis, the coupling relation between the leading shock wave and high-temperature flow was demonstrated for a methane-in-air explosion in a tunnel. In the near field of the ignition point, the deflagration wave transmits energy by heat, and the temperature load is in the front of the pressure wave. With development of deflagration and deflagration-to-detonation transition, the corresponding mechanism of energy transmission is changed from heat conduction to shock compression, and a precursor pressure wave is formed gradually. The time interval between the precursor pressure wave and high-temperature flow behind the wave increases with distance. Attenuation of the precursor shock wave and high-temperature flow depends on the length of the methane-in-air space in a tunnel. Beyond the methane-in-air space, the quantitative relation of the time interval between the precursor shock wave and high-temperature flow with axial distance from ignition and the length of methane-in-air space was proposed.     

On the initial stage of a point explosion in a radiating gas

A study is made of the initial stage of a point explosion in a radiating gray gas whose absorption coefficient is approximated by the dependenceK=x()e ^–n ,where is the density and e is the internal energy of the gas. It is shown that for n > —1/3 the initial stage of the process differs significantly from the solution of the problem in not only the classical adiabatic case [1, 2] but also in the case of a medium with nonlinear thermal conductivity [2–4]. The supply of energy to the medium at a point leads to instantaneous heating of the complete medium. The form of this heating is found analytically. The method of matched asymptotic expansions is used to investigate the behavior of the solution in the neighborhood of the center. It is found that for definite conditions at the center of the perturbed region there are formed a shock wave and a region of reverse flow of the gas.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 75–82, May–June, 1980.I should like to thank V. P. Korobeinikov for interest in the work and a helpful discussion of it.