The numerical simulations of explosion and implosion in air: use of a modified Harten's TVD scheme

Numerical simulations of explosion and implosion in air are carried out with a modified Harten's TVD scheme. The new scheme has a high resolution for contact discontinuities in addition to maintaining the good features of Harten's TVD scheme. In the numerical experiment of spherical explosion in air, the second shock wave (which does not exist in the one‐dimensional shock tube problem) and its subsequent implosion on the origin have been successfully captured. The positions of the main shock wave, the contact discontinuity and the second shock wave have shown satisfactory agreement with those predicted from previous analysis. The numerical results are also compared with those obtained experimentally. Finally, simulations of a cylindrical explosion and implosion in air are carried out. Results of the cylindrical implosion in air are compared with those of previous work, including the interaction of the reflected main shock wave with the contact discontinuity and the formation of a second shock wave. All these attest to the successful use of the modified Harten's TVD scheme for the simulations of shock waves arising from explosion and implosion. Copyright © 1999 John Wiley & Sons, Ltd.     

Reflection of a ring shock wave from a rigid wall

A new type of an implosion has been observed experimentally and simulated in a numerical calculation: the formation of a quasi-spherical converging shock wave after the reflection of a ring shock wave from a solid wall. The conversion of the ring shock wave into the quasi-spherical converging shock wave intensifies the local implosion properties.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.     

基于面积折减等效模型的光电倍增管水下内爆机理研究

光电倍增管(photomultiplier tube，PMT)是中微子探测器的核心部件，是由玻璃材料制成的内部真空的薄壳结构，排列在深水中工作，若一个PMT被压溃会产生内爆冲击波，会引起周围PMT发生殉爆。针对PMT内爆，建立了PMT内爆数值计算简化模型，并将计算与试验结果进行对比，验证简化模型的合理性。在此基础上，提出了基于面积折减等效模型的PMT内爆计算方法，通过等效模型分析了防护装置破口面积对PMT内爆的影响，得出随着防护装置破口面积的减小，水流碰撞发生PMT内爆的时刻相应提前，内爆产生的冲击波脉宽基本保持不变，冲击波压力峰值明显减小。该研究有利于找到有效的PMT内爆防护方法。     

Formation of blast waves as a result of a low pressure domain decomposition

In this article the formation and propagation of blast wave as a result of a focusing of shock wave in a domain with low pressure and density are examined in the frame of an ideal gas model. We consider the decomposition of a pressure-density discontinuity on the boundary of the spherical (or cylindrical) domain which is filled by a gas whose pressure and density are lower than the pressure and density of gas that filling the external space. At initial moments of this decomposition the rarefaction wave propagates in the external space and the converging shock wave is formed in the domain of low pressure and density. The intensity of the converging shock wave gradually increases, and the wave transitions to the self-similar regime. After implosion, a diverging shock wave is formed which propagates through the disturbed gas. The values of wave intensity, wave impulse and other parameters for some magnitudes of initial parameters have been determined by means of numerical calculations. Received 10 August 1997 / Accepted 13 July 1998     

Super-spherical cumulation: shock collapse which is more intense than the spherical one

A shock wave implosion in an axisymmetric chamber with a convex bounding wall is studied experimentally, analytically and numerically. The converging shock front area in this geometry shrinks quickly as the shock wave approaches the center point. The analytical theory predicts that the corresponding rate of post-shock pressure and density increase in this case exceeds essentially that achieved in the classical cylindrical or spherical shock implosions, hence, the phenomenon is referred to as “super-spherical cumulation”. The experiments confirm higher intensity of the super-spherical implosion compared with the cylindrical one, both driven by identical high-current pulsed electric discharges. The converging shock stability is analyzed in the framework of the CCW theory. The numerical results obtained using a locally-adaptive unstructured grid technique agree well with the theoretical predictions of the converging shock wave intensity. Received 28 January 1998 / Accepted 6 November 1998     

聚心火焰与激波相互作用的数值研究

基于带化学反应的2维轴对称Euler方程,利用带有monotonized centered(MC)限制器的波传播算法,在两端敞开的圆桶中对惰性介质的聚心激波和氢气-空气混合物的聚心火焰与激波的相互作用进行了数值模拟。数值结果表明,在惰性介质中激波在轴心的每次汇聚均可成长为马赫干,马赫干的追赶使激波得到一定程度的增强,但整体还呈下降趋势。在氢气-空气混合物中,燃烧诱导的激波,由于与火焰的反复作用,使激波在轴心处产生马赫干的频率和强度皆高于惰性介质中的情形。同时,火焰在与激波的相互作用过程中发生失稳变形,使其形状呈扁平头部的蘑菇云。     

Simulation of converging cylindrical GPa-range shock waves generated by wire array underwater electrical explosions

Results of one-dimensional numerical simulations of the parameters of the converging strong shock wave generated by electrical underwater explosions of a cylindrical wire array with different array radii and different deposited energies are presented. It was shown that for each wire array radius there exists an optimal duration of the energy deposition into the exploding array, which allows one to maximize the shock wave pressure and temperature in the vicinity of the implosion axis. The simulation results agree well with the 130-GPa pressure in the vicinity of the implosion axis that was recently obtained, which strongly indicates the azimuthal symmetry of the converging shock wave at these extreme conditions. Also, simulations showed that using a pulsed power generator with a stored energy of ~200 kJ, the pressure and temperature at the shock wave front reaches ~220 GPa and 1.7 eV at 0.1 mm from the axis of implosion in the case of a 2.5 mm radius wire array explosion. It was found that, in spite of the complicated equation of state of water, the maximum pressure at the shock wave front at radius r can be estimated as P ≈ (P*(r*/r) ^α , where P* is the known value of pressure at the shock wave front at radius r* ≥ r and α is a parameter that equals 0.62±0.02. A rough estimate of the implosion parameters of the hydrogen target after the interaction with the converging strong shock wave is presented as well.     

Numerical investigation of axisymmetric shock wave focusing over paraboloidal reflectors

The problem of a plane shock wave incident to a paraboloidal reflector is numerically investigated. The numerical solver used is developed by an improved, implicit, upwind total variation diminishing scheme in a finite-volume approach. The real-gas effect is taken into account if high temperature occurs. The solver is validated on four test problems. The complicated flow fields of axisymmetric shock wave focusing for different-depth reflectors at various incident shock Mach numbers are studied. An interesting result of a maximum pressure happening at the reflector center is found. This is due to the occurrence of an implosion phenomenon. A maximum temperature might occur at the reflector center or at other locations, depending on the incident shock Mach number and the reflector depth. Moreover, vortical flows induced by shock wave focusing and their formation mechanism are explored. It was found that the vortices near the reflector are caused by a ring-shaped shock/slipline interaction. Owing to the slipline on the symmetry axis, a jet flow is induced, resulting in the formation of vortices near the symmetry axis. Received 13 January 1998 / Accepted 10 November 1998     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序,对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为,进行了数值模拟研究。结果表明,首次冲击后至约12 μs,界面发展以RM（Richtmyer-Meshkov）不稳定性为主;12 μs后至冲击波聚心反弹加载前,界面聚心运动处于加速减速状态,界面发展由RT (Rayleigh-Taylor)不稳定性主导;冲击波聚心反弹加载后,界面发展又由RM不稳定性主导。另外,还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大;初始波长较小（模数较大）时振幅增长较小,而且存在一个截止波长;钢壳厚度会抑制扰动增长,也存在一个截止厚度;几何汇聚效应会使扰动增长速度更快。     

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.     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Shock tube investigation of hydrodynamic issues related to inertial confinement fusion

A shock tube investigation of two hydrodynamic issues related to inertial confinement fusion (ICF) is undertaken. ICF is a promising source of energy for the future. There has been a considerable increase in the interest in ICF with the development of the National Ignition Facility (NIF). However, much remains to be investigated before a useful yield is obtained from a fusion reaction for power generation. The physics involved in carrying out a fusion reaction combines hydrodynamics, plasma physics and radiation effects superimposed on each other, at extremely small scales, making the problem very complex. One such phenomenon occurring in the deuterium-tritium pellet implosion is the Richtmyer-Meshkov instability occuring at each layer of the fuel which results in the mixing of the ablator with the fuel. This causes dilution of the fuel and reduces the yield of the reaction. Another issue is the impulsive loading of ICF reactor cooling tubes due to the shock wave produced as a result of the fusion reaction. These tubes must withstand the impulse of the shock wave. A shock tube provides an ideal environment to study these issues at large geometric scales with the isolation of hydrodynamics from other effects. A new vertical, square shock tube has been designed specifically for the purpose of studying these fluid flow phenomena from a fundamental point of view. The shock tube is vertical, with a large square inner cross-section and is designed to allow for the release of a shock into air at atmospheric pressure. In this paper, we describe the new shock tube and related instrumentation in detail and present a few preliminary results on the Richtmyer-Meshkov instability and shock-cylinder interactions. Received 5 January 1999 / Accepted 10 July 2000     

A STUDY ON JET PHENOMENON OF R22 GAS CYLINDER UNDER THE IMPACT OF SHOCK WAVE

重物质气泡在平面激波作用下出现射流是一个非常有趣的物理现象, 它与气体属性、激波强度密切相关, 气泡内外波系结构复杂, 值得进行深入探讨. 论文采用水平集方法结合真实虚拟流体方法捕捉物质界面, 对R₂₂圆柱气泡在平面激波作用下的射流现象进行了数值模拟, 详细阐述了波系结构在气泡内外的演变过程, 总结了导致气泡射流现象形成的典型波系结构特征, 从波系结构角度解释了射流形成的原因.论文也对不同强度的激波冲击R₂₂圆柱气泡的情况进行了归纳, 发现激波强度对射流现象的出现以及射流结构的尺寸特征有重要影响, 得出了 R₂₂圆柱气泡能否出现射流结构的临界激波强度, 以期对实验结果起到指导与验证作用.     

Interaction of a shock wave with a particle cloud of finite size

We present the results of the numerical modelling of the interaction of a shock wave with a cloud of finite size particles. The computations were carried out within the framework of continuum/discrete model with the use of the techniques of digital diagnostics and pattern recognition. The shock wave and vortex formation behind the cloud of particles as well as the formation of a dense layer in the cloud have been revealed. For this reason, the use of a cloud of particles for relaxing the shock wave may prove to be inefficient.     

Collective Bow Shock Ahead of a Transverse System of Spheres in a Supersonic Flow Behind a Moving Shock Wave

The results of investigating the dynamics and physical conditions of formation of a collective bow shock ahead of a system of spheres with the line of centers normal to the supersonic flow behind a traveling shock wave are presented. Two types of shock-wave patterns that necessarily precede the formation of the collective shock wave and correspond to regular and Mach interaction of the bow shocks were detected experimentally. On the basis of a local gasdynamic-discontinuity interference theory, quantitative criteria of the existence of these regimes and of the formation of a common shock wave are determined. These criteria are confirmed in a series of experiments for the transitional regimes.     

激波喷出过程和在林带中的传播

运用激波管技术研究了激波从管口喷出和在林带中传播这两个基本问题。给出了三方面的结果 :揭示了激波从管口喷出时的复杂流场 ,显示了主激波弯曲和衰减 ,二次激波的形成和合拢 ,涡环的长大和发展 ,以及涡环前缘出现正激波等基本物理现象 ;显示了激波在林带中的传播过程 ,取得了激波遇林地的反射 ,马赫杆形成 ,激波与林冠的相遇 ,林冠对涡环的阻滞效应等纹影照片 ;沿激波在林地的传播方向测量了地面压力 ,证实林带有明显的消波效应。     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成; (2) U型空泡团头部溃灭; (3) U型空泡团腿部溃灭.在U 型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象, 抑制了下一周期的空化发展.      

Study of unsteady flow past bodies behind shock waves

Several theoretical and experimental studies have been devoted to the problem of the nonstationary action of the stream behind a shock wave on bodies of varied shape. In particular, in [1], the pressure and density are calculated for flow about bodies of the more typical shapes in the initial stage of the process. The basic relations which accompany the interaction of shock waves are considered in [2, 3]. The analysis of the phenomena of diffraction of shock waves on the sphere, cylinder, and cone is presented in [4]. Problems of unsteady flow about a wing are examined in [5, 6]. A detailed review of the foreign studies on unsteady flow is given in [7]. Of great practical interest is the question of the time for flow formation and the magnitudes of the unsteady loads during this period. Experimental investigations have been made recently [8, 9] in which some criteria are presented for estimating the bow shock formation time for supersonic flow about the sphere and cylinder with flat blunting. However the question of the formation time of the stationary pressure on the body surface is not referred to in these studies and no relationship is shown between the transient position of the reflected wave and the corresponding unsteady pressure on the surface. Moreover, in [8] the dimensionless time criterion is determined very approximately, independently of the Mach number of the shock wave. The present study was undertaken with the object of determining the basic criteria which characterize unsteady flow about bodies behind a plane shock wave which has time-independent parameters, and clarification of the shock wave reflected from the body and the pressure on the surface of the body during the transient period. The most typical body shapes were studied: 1) a cylinder with flat face aligned with the stream; 2) a spherically-blunted cylinder; and 3) a cylinder transverse to the stream. The experiments were conducted in a conventional shock tube using the single-diaphragm scheme. The measurements of the pressure on the models and the velocity of the incident shock wave were made using the technique analogous to that of [10, 11]. A highspeed movie camera was used to record the pattern of the wave diffraction on the body. The Mach number of the incident shock wave varied in the range from M=1.5 to M≈6.0, which corresponded to a range of Mach numbers M_∞ of the stream behind the shock wave from 0.6 to 2.1. The calculations of the required gas dynamic parameters for high temperatures were made with account for equilibrium dissociation of the air on the basis of the data of [10, 12, 13]. The magnitude of the relative maximal shock wave standoff Δ at the stagnation point obtained in the present experiments was compared with the values of Δ from other studies. In the case of the flat-blunted cylinder it was in good agreement with the results of [8–14], and in the case of the spherically-blunted cylinder and the transverse cylinder it was in agreement with the results of [15].     

Effects of a Single-pulse Energy Deposition on Steady Shock Wave Reflection

The effects of energy deposition in the free stream on steady regular and Mach shock wave reflections are studied numerically. A short-duration laser pulse is focused upstream of the incident shock waves. It causes formation of the expanding blast wave and the residual hot-spot interacting in a complex way with the steady shock wave reflection. It was found that the laser energy addition in the free stream may force the transition from regular to Mach reflection in the dual solution domain. In contrast to previously reported numerical results, the transition from Mach to regular reflection has not been reproduced in our refined computations since the Mach reflection is restored after the flow perturbation.     

Half a century of continuous shock interaction investigations in the Joint Institute for High Temperatures of Russian Academy of Sciences

This article describes the history of the investigations of shock wave interactions at the Physical Gasdynamic Department, starting from the early 50s of the last century, when the first research related to missile reentry was made. The review focuses on a number of topics studied over more than 50 years and includes the study of strong shock waves, where it is necessary to take into account the physicochemical transformations in gases, shock wave reflection, diffraction, interaction with the boundary layer and with the nozzle, as well as detonation wave formation and interactions. The investigation of shock wave interactions is a current topic at the Joint Institute for High Temperatures of the Russian Academy of Sciences. Some new results are observed: the formation of impulse jets and the self-ignition of a cold hydrogen jet, diffraction of 3D shock waves, and the effect of an impulse jet and diffracted shock wave on an obstacle.