水下爆炸冲击波数值仿真研究

由于在水下爆炸冲击波的数值仿真研究中,水的状态方程、人工黏性系数和网格尺寸对数值计算结果影响很大,采用常规TNT炸药的水下爆炸为例,以冲击波的峰值压力和比冲量为衡量指标,研究了这3个主要影响因素对数值仿真结果的影响。首先,通过采用常用的5种水的状态方程进行系列仿真,给出了各种状态方程的适用范围;其次,讨论了人工黏性系数对计算结果的影响,并给出了一次与二次人工黏性系数的建议取值范围;最后,通过对不同炸药当量及不同网格尺寸开展系列运算,从而得到不同炸药当量在满足工程计算精度要求下所对应的建议网格尺寸,并得到了不同炸药当量所对应的建议网格尺寸的表达式。     

炸药混凝土中爆炸能量释放规律的研究

 通过对炸药在混凝土中爆炸的并行计算,分析其数值计算结果,并对爆炸力学基础理论进行研究后认为：炸药在土岩介质中爆炸,炸药能量将通过爆炸空腔中爆炸产物内冲击波的多次往复反射的方式进行释放,同时破坏并推动土岩介质移动,土岩介质与爆炸气体的界面在爆腔内稀疏波的作用下进行持续减速扩张运动。     

Three-dimensional detonation structure and its response to confinement

Three-dimensional (3D) detonation simulations solving the compressible Navier-Stokes equations with detailed chemistry are performed in both square channel and round tube geometries. The simulations are compared with each other and with two-dimensional (2D) channel simulations and round tube experiments of identical mixture and conditions (stoichiometric hydrogen-oxygen with 3000 PPMv ozone at 300 K and 15 kPa) with the goal of understanding the effect of confinement and boundaries on detonation structure. Results show that 3D detonations propagate with highly inhomogeneous blast dynamics, where blasts emerge not only from intersections of two transverse waves (similar to 2D propagation) but also from intersections of many transverse waves (unique to 3D detonations in the confinements tested). Intersections of many transverse waves lead to extreme thermodynamic states and highly overdriven wave velocities, well in excess of those seen in the ZND model and in 2D simulations. 3D simulations in the square tube show highly regular blast latticing, smaller detonation cells, and highly oscillatory velocities when compared to the round tube simulations. Round tube simulations show more spatially non-uniform blast dynamics. The conclusions reached in the current work are found irrespective of numerical grid resolution.     

含铝炸药水下爆炸性能的实验研究

对RS211、HL-1、HL-2这3种含铝炸药和TNT炸药进行了水下爆炸实验,测量了4种炸药水下爆炸冲击波压力剖面和气泡脉动周期,获得了4种炸药水下爆炸冲击波的峰值超压、冲量、能量和气泡能。实验结果表明,在所研究的范围内,与TNT炸药相比,含铝炸药的冲击波能量有明显增加,约为TNT炸药的1.20～1.35倍,气泡能有显著增加,约为TNT炸药的1.50～2.30倍,表明在炸药中加入铝粉对于提高炸药水中爆炸威力是有益的。     

High strain rate response of an elastomer

Abstract

Plate impact experiments have been performed to examine the high strain rate response of an elastomer (explosive binder) shocked to 8 kbar stress. Particle velocity data have been obtained using in-material electromagnetic velocity gauges. Normal impact experiments have provided compression and release wave velocities and stress-strain results; ramp wave experiments indicate that the elastomer response is independent of loading rate; inclined plate impact measurements show that the material does not support any shear even at these strain rates. Additionally, numerical simulations have been performed to show that shock propagation in viscous fluids canbe simulated without artificial viscosity but by using the material viscosity in the constitutive relations. Such simulations will be useful in obtaining material viscosities from the measured wave profiles.     

氢化镁储氢型乳化炸药的爆炸特性研究

 通过理论计算和水下爆炸实验,初步研究了MgH₂敏化储氢型乳化炸药的爆炸特性和爆轰反应机理。结果表明：与玻璃微球敏化的乳化炸药相比,MgH₂敏化的乳化炸药水下爆炸的冲击波超压、比冲量、比冲击波能、比气泡能及水下爆炸比总能量显著增加,其中冲击波超压和水下爆炸总能量分别增加了20.5%和31.0%。MgH₂储氢型乳化炸药的爆轰机理与玻璃微球敏化乳化炸药不同,MgH₂在乳化炸药中起到了敏化剂和含能材料的双重作用,即MgH₂在乳化基质中水解产生均匀分布的氢气泡,起到了敏化作用,同时氢气参与爆炸反应,提高了炸药的爆炸能量和做功能力。     

Non-invasive determination of shock wave pressure generated by optical breakdown

Shock waves generated by a laser-induced plasma were investigated using a pump-and-probe technique. Both 7-ns and 40-ps laser pulses at 1.06 m were employed to initiate breakdown in water. Two He-Ne laser beams were used as a velocity probe, allowing the accurate measurement of the shock velocity around the plasma. The maximum shock pressure was determined from the measured shock velocities, the jump condition and the equation of state for water. The conservation of the total momentum of the shock front was used to derive expressions for the shock velocity, particle velocity and shock pressure vs. the distance (r) from the center of the plasma. For a shock wave of spherical symmetry, the shock pressure is proportional to 1/r ². Our work shows that the expanding plasma initially induces a shock wave; the shock wave dissipates rapidly becoming an acoustic wave within 300–500 m.     

Blast waves generated in an unconfined space by a nonideal detonation of high-density aluminum-enriched formulations

Experiments on the detonation of high-density (1.8 g/cm³) aluminum-ammonium perchlorate-paraffin-RDX formulations in an unconfined space demonstrated their high efficiency at pressure amplitudes within 0.3–7.0 atm. The relative pressure amplitude and impulse of the blast waves with respect to the analogous characteristics of TNT charges of the same mass were found to be 1.3–2.4. The TNT equivalents in pressure and impulse vary with the distance nonmonotonically, ranging within 1.4–2.8. The blast wave produced by an infield explosion of a 1.42-kg composite charge demonstrated its high performance characteristics. Measurements at blast wave amplitudes of 1 to 20 atm gave a TNT equivalent in pressure of up to 3 and a TNT equivalent in impulse of 1.3 to 1.8. The high parameters of blast waves in an unconfined space originate from both the high-energy characteristics of the systems themselves and the afterburning of excess metal fuel in air. To estimate the extent of participation of the reaction of excess metal fuel with air in supporting the blast wave, numerical simulations of the generation of blast waves for various rates of mixing of detonation products with air at the contact surface were conducted. The main elements of the mechanisms of the processes that determine the efficiency of explosive systems with a heat release spread in space and time were considered. It was concluded that an optimal regime of blast wave generation, capable of ensuring a prolonged attenuation of the wave with the distance, could be realized for low-velocity detonation.     

爆炸信号中气泡脉动去除方法及其应用

水下爆炸声源激发的声信号包括冲击波和气泡脉动,气泡脉动严重干扰冲击波的传播特性。当气泡脉动的幅度小于冲击波的幅度时,通过对此类混合爆炸信号进行复倒谱分析,提出了基于卷积型的爆炸信号模型的同态解卷积气泡脉动去除方法,其中,针对同态滤波系统中通常采用的“梳状”滤波器会在倒谱的尖峰处产生不连续点的问题,采用了对尖峰附近的采样点进行多项式函数拟合,并根据拟合函数对尖峰采样点插值的改进型滤波方法。数值仿真和实验信号的处理及分析表明,应用该方法可以明显消除爆炸信号自相关曲线中由气泡脉动造成的对称次尖峰,并在信号的时频分布图中再现冲击波的简正波特性。这为进一步利用爆炸声源研究海洋信道的声传递函数和反演海洋环境参数提供了技术途径。      

炸药水下爆炸能量输出特性试验研究

为了解炸药水中爆炸能量的输出特性,特别是冲击波能随距离的衰减规律,对有效比冲能进行了研究。经过理论推导发现,传统的有效比冲能计算公式仍然满足相似律。TNT炸药水下爆炸试验结果表明,炸药的有效比冲能、比气泡能和总能量的测量结果与经验公式计算结果吻合较好。通过对试验结果进行分析,拟合得到了有效比冲能随距离变化的趋势线。对比经验公式结果,发现在小药量试验中冲击波能随距离的衰减更明显。     

非均匀炸药冲击起爆和起爆后的行为

 介绍并分析了Campbell 等人及其他作者研究非均匀炸药冲击起爆和起爆后行为所获得的实验结果，但不涉及其冲击起爆条件。足够强的冲击波进入非均匀炸药后，爆轰将瞬时（指不经过感应时间）且直接（指不经过其他过程，如爆燃）被引发；非均匀炸药起爆后，其中传播的自始至终是一个不断增长的爆轰波，直至发展为正常爆轰，整个过程都是爆轰的增长（新定义）过程。不存在由反应冲击波不断增长并转变为爆轰波的所谓向爆轰的增长。所谓向爆轰的增长，实际上是爆轰的增长（按新定义）的初期；Craig原定义的爆轰的增长，实际上是爆轰的增长（按新定义）的后期；而所谓反应冲击波，实际上是增长中的初期爆轰波。爆轰的增长（按新定义）是所有猛炸药的特性，炸药反应不充分并逐渐趋于充分是爆轰的增长的化学机制。     

二氧化硅气凝胶的防爆震性能及机理研究

利用聚偏二氟乙烯压电传感器研究了爆炸加载下冲击波在二氧化硅气凝胶中的传播特性,并对冲击波在二氧化硅气凝胶和泡沫铝中的传播特性进行了比较.结果表明:在二氧化硅气凝胶中冲击波的强度随传播距离的增加呈现指数衰减的趋势.冲击波在二氧化硅气凝胶中衰减比在泡沫铝中衰减明显.由于二氧化硅气凝胶内部特殊的纳米多孔网状结构,导致冲击波在二氧化硅气凝胶中的衰减效果较好.冲击波在二氧化硅气凝胶中的传播速度极低,因此冲击波在二氧化硅气凝胶中传播时卸载波的追赶卸载效应非常明显,这又进一步促进了冲击波的衰减.     

A Composite Scheme for Gas Dynamics in Lagrangian Coordinates

One cycle of a composite finite difference scheme is defined as several time steps of an oscillatory scheme such as Lax–Wendroff followed by one step of a diffusive scheme such as Lax–Friedrichs. We apply this idea to gas dynamics in Lagrangian coordinates. We show numerical results in two dimensions for Noh's infinite strength shock problem and the Sedov blast wave problem, and for several one-dimensional problems including a Riemann problem with a contact discontinuity. For Noh's problem the composite scheme produces a better result than that obtained with a more conventional Lagrangian code.     

Numerical Study of Air Flow Induced by Shock Impact on an Array of Perforated Plates

This study is focused on the propagation behavior and attenuation characteristics of a planar incident shock wave when propagating through an array of perforated plates. Based on a density-based coupled explicit algorithm, combined with a third-order MUSCL scheme and the Roe averaged flux difference splitting method, the Navier–Stokes equations and the realizable k-ε turbulence model equations describing the air flow are numerically solved. The evolution of the dynamic wave and ring vortex systems is effectively captured and analyzed. The influence of incident shock Mach number, perforated-plate porosity, and plate number on the propagation and attenuation of the shock wave was studied by using pressure- and entropy-based attenuation rates. The results indicate that the reflection, diffraction, transmission, and interference behaviors of the leading shock wave and the superimposed effects due to the trailing secondary shock wave are the main reasons that cause the intensity of the leading shock wave to experience a complex process consisting of attenuation, local enhancement, attenuation, enhancement, and attenuation. The reflected shock interactions with transmitted shock induced ring vortices and jets lead to the deformation and local intensification of the shock wave. The formation of nearly steady jets following the array of perforated plates is attributed to the generation of an oscillation chamber for the inside dynamic wave system between two perforated plates. The vorticity diffusion, merging and splitting of vortex cores dissipate the wave energy. Furthermore, the leading transmitted shock wave attenuates more significantly whereas the reflected shock wave from the first plate of the array attenuates less significantly as the shock Mach number increases. The increase in the porosity weakens the suppression effects on the leading shock wave while increases the attenuation rate of the reflected shock wave. The first perforated plate in the array plays a major role in the attenuation of the shock wave.     

An adaptive moving mesh method for two-dimensional ideal magnetohydrodynamics

This paper presents an adaptive moving mesh algorithm for two-dimensional (2D) ideal magnetohydrodynamics (MHD) that utilizes a staggered constrained transport technique to keep the magnetic field divergence-free. The algorithm consists of two independent parts: MHD evolution and mesh-redistribution. The first part is a high-resolution, divergence-free, shock-capturing scheme on a fixed quadrangular mesh, while the second part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative-interpolation formula is used to calculate the remapped cell-averages of the mass, momentum, and total energy on the resulting new mesh; the magnetic potential is remapped to the new mesh in a non-conservative way and is reconstructed to give a divergence-free magnetic field on the new mesh. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy, track and resolve strong shock waves in ideal MHD problems, and preserve divergence-free property of the magnetic field. Numerical examples include the smooth Alfvén wave problem, 2D and 2.5D shock tube problems, two rotor problems, the stringent blast problem, and the cloud–shock interaction problem.     

A study for sound wave scattering by corrugated ground with complex trench structures

Several trench structures in corrugated ground are investigated for the possibility of mitigating gun blast noise by numerical simulations. The blast noise usually includes large explosive energy with nonlinearity in the near field and exhibits a very low-frequency spectrum. In this study, the linearity approximation for the noise is taken because the nonlinearity of the wave reaching the scatterer is not serious for many proved guns and the low-frequency characteristic is concentrated. The structures are designed based on the surface impedance design approach proposed in our previous work and arbitrary three-dimensional (3D) geometries within a truncated ground are now assumed. The acoustic characteristic of the structures is evaluated by using a fast numerical solver. The solver employs the multilevel fast multipole algorithm (MLFMA) as an accelerator and can solve very large acoustic wave scattering problems with millions of unknowns on workstations within several days. This tool allows us to truncate the ground as large as needed for accurate modeling. Four structures are mainly considered in the design, namely, concentric trenches, sectorial trapezoidal trenches, interlaced arc trenches and parabolic reflectors. Some of them may have a sloped inner wall or tilted surface as a means of adjustment. Numerical simulations show that the concentric trench design has a very good mitigation behavior for linear and continuous noise sources and the structure is further studied for mitigating real-world gun blast noise.     

炸药爆轰驱动超高速击波管的数值摸拟

 高能炸药（HE）柱壳沿轴向传播的稳定态爆轰，在置于炸药柱壳内部的泡沫塑料（工作流体）中产生很强的轴向击波，达到极高的压力和能量密度。适当设计，可使工作流体的流动接近一稳定态，与炸药爆轰以相同的相速度前进。击波阵面后的工作流体可用Nozzle流动方程很好地描述^[1]，Nozzle收缩段中的流动不受散开段流动的影响。利用这一性质，可以设计出达到极高压力和能量密度的超高速击波管，并且不受管壁物质强度的限制。数值模拟计算给出了上述稳定态物理图象，并显示出开始阶段轴向击波的形成过程及其后对稳定态的逼近，计算结果及物理图象与理论分析符合得很好。     

An experimental investigation of the propagation mechanism of critical deflagration waves that lead to the onset of detonation

The reflection of a CJ detonation from a perforated plate is used to generate high speed deflagrations downstream in order to investigate the critical conditions that lead to the onset of detonation. Different perforated plates were used to control the turbulence in the downstream deflagration waves. Streak Schlieren photography, ionization probes and pressure transducers are used to monitor the flow field and the transition to detonation. Stoichiometric mixtures of acetylene–oxygen and propane–oxygen were tested at low initial pressures. In some cases, acetylene–oxygen was diluted with 80% argon in order to render the mixture more “stable” (i.e., more regular detonation cell structure). The results show that prior to successful detonation initiation, a deflagration is formed that propagates at about half the CJ detonation velocity of the mixture. This “critical” deflagration (which propagates at a relatively constant velocity for a certain duration prior to the onset of detonation) is comprised of a leading shock wave followed by an extended turbulent reaction zone. The critical deflagration speed is not dependent on the turbulence characteristics of the perforated plate but rather on the energetics of the mixture like a CJ detonation (i.e., the deflagration front is driven by the expansion of the combustion products). Hence, the critical deflagration is identified as a CJ deflagration. The high intensity turbulence that is required to sustain its propagation is maintained via chemical instabilities in the reaction zone due to the coupling of pressure fluctuations with the energy release. Therefore, in “unstable” mixtures, critical deflagrations can be supported for long durations, whereas in “stable” mixtures, deflagrations decay as the initial plate generated turbulence decays. The eventual onset of detonation is postulated to be a result of the amplification of pressure waves (i.e., turbulence) that leads to the formation of local explosion centers via the SWACER mechanism during the pre-detonation period.     

Asymptotic description of solitary wave trains in fully nonlinear shallow-water theory

We derive an asymptotic formula for the amplitude distribution in a fully nonlinear shallow-water solitary wave train which is formed as the long-time outcome of the initial-value problem for the Su–Gardner (or one-dimensional Green–Naghdi) system. Our analysis is based on the properties of the characteristics of the associated Whitham modulation system which describes an intermediate “undular bore” stage of the evolution. The resulting formula represents a “non-integrable” analogue of the well-known semi-classical distribution for the Korteweg–de Vries equation, which is usually obtained through the inverse scattering transform. Our analytical results are shown to agree with the results of direct numerical simulations of the Su–Gardner system. Our analysis can be generalised to other weakly dispersive, fully nonlinear systems which are not necessarily completely integrable.