Head on collision of a planar shock wave with a dusty gas layer

The head-on collision of a planar shock wave with a dust-air suspension is studied numerically. In this study the suspension is placed inside a conduit adjacent to its rigid end-wall. It is shown that as a result of this collision two different types of transmitted shock waves are possible, depending on the strength of the incident shock wave and the dust loading ratio in the suspension. One possibility is a partially dispersed shock wave, the other is a compression wave. The flow fields resulting in these two options are investigated. It is shown that in both cases, at late times after the head-on reflection of the transmitted shock wave from the conduit end-wall a negative flow (away from the end-wall) is evident. The observed flow behavior may suggest a kind of dust particle lifting mechanism that could shed new light on the complex phenomenon of dust entrainment behind sliding shock waves.      

Shock wave attenuation by grids and orifice plates

The interaction of weak shock waves with porous barriers of different geometries and porosities is examined. Installing a barrier inside the shock tube test section will cause the development of the following wave pattern upon a head-on collision between the incident shock wave and the barrier: a reflected shock from the barrier and a transmitted shock propagating towards the shock tube end wall. Once the transmitted shock wave reaches the end wall it is reflected back towards the barrier. This is the beginning of multiple reflections between the barrier and the end wall. This full cycle of shock reflections/interactions resulting from the incident shock wave collision with the barrier can be studied in a single shock tube test. A one-dimensional (1D), inviscid flow model was proposed for simulating the flow resulting from the initial collision of the incident shock wave with the barrier. Fairly good agreement is found between experimental findings and simulations based on a 1D flow model. Based on obtained numerical and experimental findings an optimal design procedure for shock wave attenuator is suggested. The suggested attenuator may ensure the safety of the shelter’s ventilation systems.     

A holographic interferometric study of shock wave focusing in a circular reflector

A holographic interferometric study was made of the focusing of reflected shock waves from a circular reflector. A diaphragmless shock tube was used for incident shock Mach numbers ranging from 1.03 to 1.74. Hence, the process of reflected shock wave focusing was quantitatively observed. It is found that a converging shock wave along the curved wall undergoes an unsteady evolution of mach reflection and its focusing is, therefore, subject to the evolution of the process of shock wave reflections. The collision of triple points terminates the focusing process at the geometrical focus. In order to interprete quantitatively these interferograms, a numerical simulation using an Eulerian solver combined with adaptive unstructured grids was carried out. It is found numerically that the highest density appears immediately after the triple point collision. This implies that the final stage of focusing is mainly determined by the interaction between shock waves and vortices. The interaction of finite strength shock waves, hence, prevents a curved shock wave from creating the infinite increase of density or pressure at a focal point which is otherwise predicted by the linear acoustic theory.     

Dynamics of stress wave propagation in a chain of photoelastic discs impacted by a planar shock wave: Part II, numerical investigation

The propagation of stress waves through a chain of discs has been studied experimentally in Part I (Glam et al. [1]) and is completed here with numerical investigation using the standard package ABAQUS. A fair agreement is found between experimental findings and their simulations. Based on this agreement, parametric study of wave propagation through disc-chains was conducted. Specifically, effects associated with changes in the disc diameter, material density, stiffness/rigidity and the number of discs in the chain on the stressed chain have been studied. It was found that the propagation velocity of the evolved waves increases with improving contacts between the chain’s discs by exposing the chain to a static load before its dynamic loading. The wave- propagation velocity decreases with increase in the discs material density and it increases when its diameter increases. In case of a chain composed of small diameter discs and/or small material density, the transmitted stress wave is first strengthened and only at discs further down the chain it starts decaying. When checking the influence of the dynamic-loading duration it was found that long dynamic-load duration dissolves quickly into short pulses. It was also found that there is a ‘characteristic’ wave for a given chain. This wave propagates with minimal dispersion. Dynamic loads having shorter time duration than the ‘characteristic’ one experiences significant attenuation.     

Shock wave interaction with cellular materials

The equations governing the head-on collision of a planar shock wave with a cellular material and a numerical scheme for solving the set of the governing equations were outlined. In addition, the condition for the transmitted compression waves to transform into a shock wave, inside the cellular material was introduced. It was proved analytically that a cellular material cannot be used as a means of reducing the pressure load acting on the end-wall of the shock tube. In subsequent papers, the interaction of planar shock waves with specific cellular materials, e.g., foams and honeycombs will be described in detail.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.     

矩形凹槽管道中激波传播的数值研究

基于二维Euler方程,结合五阶加权基本无振荡(weighted essentially nonoscillatory,WENO)格式以及自适应网格加密(adaptive mesh refinement,AMR)技术对入射激波在矩形凹槽管道内传播过程进行了数值模拟。数值结果清晰地显示了入射激波传播过程中与多个矩形凹槽作用以及在凹槽内变化的整个过程,且与已有的实验结果吻合较好。另外,结果还揭示了入射激波与单个凹槽作用时,会发生绕射产生膨胀波,还会发生碰撞从而诱导反射激波。膨胀波会导致入射激波压力降低,而反射激波则导致其升高,但膨胀波的影响占主导作用,因而入射激波波阵面强度出现振荡下降。     

Head-on collision between normal shock waves and a rubber-supported plate,a parametric study

The equations governing the flow resulting from a head-on collision between a normal shock wave and a rubber-supported plate are listed. The non-dimensional parameters that may affect the resulting flow are specified and their influence on the post-collision flow and waves is studied numerically. It is shown that changes in: the area-ratio between the gas and the rubber cross-sections, the incident shock wave Mach number and the mass ratio between the rubber and the plate it supports, all have significant effects on the post-collision gas and rubber responses. Changes in the rubber elasticity constant also affect the post-collision flow. The extent of the effect that changes in the above mentioned parameters have on the post-collision flow responses depends upon the loading mode used. Three different modes were studied; uni-axial stress loading, bi-axial stress loading and uni-axial strain loading.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.     

The formation of a secondary shock wave behind a shock wave diffracting at a convex corner

This paper deals with the formation of a secondary shock wave behind the shock wave diffracting at a two-dimensional convex corner for incident shock Mach numbers ranging from 1.03 to 1.74 in air. Experiments were carried out using a 60 mm 150 mm shock tube equipped with holographic interferometry. The threshold incident shock wave Mach number () at which a secondary shock wave appeared was found to be = 1.32 at an 81° corner and = 1.33 at a 120° corner. These secondary shock waves are formed due to the existence of a locally supersonic flow behind the diffracting shock wave. Behind the diffracting shock wave, the subsonic flow is accelerated and eventually becomes locally supersonic. A simple unsteady flow analysis revealed that for gases with specific heats ratio the threshold shock wave Mach number was = 1.346. When the value of is less than this, the vortex is formed at the corner without any discontinuous waves accompanying above the slip line. The viscosity was found to be less effective on the threshold of the secondary shock wave, although it attenuated the pressure jump at the secondary shock wave. This is well understood by the consideration of the effect of the wall friction in one-dimensional duct flows. In order to interpret the experimental results a numerical simulation using a shock adaptive unstructured grid Eulerian solver was also carried out. Received 1 May 1996 / Accepted 12 September 1996     

水下连续脉冲冲击波的声学特性

为探求金属爆炸索在水下爆炸声源研究领域的应用前景,设计了一种可以连续产生若干个脉冲冲击波的装置,称之为水下连续脉冲冲击波发生装置。利用小波分析对该装置产生的连续脉冲冲击波信号进行分解与重构,考察其频谱特性,并进一步分析了信号的声压级特性。结果表明:该装置产生的信号声压级较高,具有很强的声功率;信号包含频率十分丰富,雷管和金属爆炸索由于装药结构及传爆方式的不同,爆炸所产生的冲击波频谱特性也有所差异。雷管爆炸产生的冲击波主要分布在15.6 kHz以下的频带内,金属爆炸索爆炸产生的脉冲冲击波信号则主要分布在62.5 kHz以下的频带内;脉冲冲击波的个数和声持续时间可由爆炸索的排列方式和长度控制,脉冲冲击波间的时间间隔可调,发生装置稳定易控     

Numerical simulation of the head-on reflection of a regular reflection

The head-on collision and subsequent reflection of a Regular Reflection (RR) from the end-wall of a shock tube has been investigated both experimentally and numerically for two different incident shock wave Mach numbers and two different reflecting wedge angles. The agreement between the double-exposure holographic interferograms and the numerical simulations which were obtained using a GRP based numerical code, was found to be excellent in the RR region and very good behind the head-on reflected RR. The overall good agreement between the computed and experimental constant-density contours (isopycnics) constitutes a validation of the computational method, including the oblique-wall boundary condition.     

The reflection of a shock wave over a cone

An investigation was made of the reflection of planar shock waves from cones. 86 cones, the half apex angle of which varied from 10° to 52° at every 0.5°, were installed in a 60 mm×150 mm diaphragmless shock tube equipped with holographic interferometry. The diaphragmless shock tube had a high degree of reproducibility with which the scatter of shock wave Mach number was within ±0.25% for shock wave Mach number ranging from 1.16 to approximately 2.0. The reflection of shock waves over cones was visualized using double exposure holographic interferometry. Whitham's geometrical shock wave dynamics was used to analyse the motion of Mach stems over cones. It is found that for relatively smaller apex angles of cones trajectory angles of resulting irregular reflections coincide with the so-called glancing incidence angles and their Mach stems appear to be continuously curved from its intersection point with the incident shock wave, which shows the chractericstic of von Neumann reflection. The domain of the existence of the von Neumann reflection was analytically obtained and was found to be broadened much more widely than that of two-dimensional reflections of shock waves over wedges.     

Shock wave propagation in a branched duct

The propagation of a planar shock wave in a 90° branched duct is studied experimentally and numerically. It is shown that the interaction of the transmitted shock wave with the branching segment results in a complex, two-dimensional unsteady flow. Multiple shock wave reflections from the duct's walls cause weakening of transmitted waves and, at late times, an approach to an equilibrium, one-dimensional flow. While at most places along the branched duct walls calculated pressures are lower than that existing behind the original incident shock wave, at the branching segment's right corner, where a head on-collision between the transmitted wave and the corner is experienced, pressures that are significantly higher than those existing behind the original incident shock wave are encountered. The numerically evaluated pressures can be accepted with confidence, due to the very good agreement found between experimental and numerical results with respect to the geometry of the complex wave pattern observed inside the branched duct. Received 15 July 1996 / Accepted 20 February 1997     

爆轰波与激波对撞的实验研究

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波,以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.      

Evolution of a laser-generated shock wave in iron and its interaction with martensitic transformation and twinning

Effects of shock waves (generated by a nanosecond laser pulse in plates of Armco-iron) on structural changes are analysed. Localisation of processes of martensitic transformation and twinning – for various values of laser pulse duration – is studied both experimentally and numerically. A proposed model accounts for interaction of shock wave propagation and structure changes. Realisation of martensitic transformation and twin formation influences wave front modification. A stress amplitude decrease with increasing distance from a microcrater determines, together with the pulse duration, a character of spatial localisation of structural changes. Numerical results are compared with experimental data and serve as a basis for additional interpretation of phenomena. Received 9 August 1994 / Accepted 30 June 1997     

Interaction of a shock wave with a loose dusty bulk layer

The interaction of a planar shock wave with a loose dusty bulk layer has been investigated both experimentally and numerically. Experiments were conducted in a shock tube. The incident shock wave velocity and particle diameters were measured with the use of pressure transducers and a Malvern particle sizer, respectively. The flow fields, induced by shock waves, of both gas and granular phase were visualized by means of shadowgraphs and pulsed X-ray radiography with trace particles added. In addition, a two-phase model for granular flow presented by Gidaspow is introduced and is extended to describe such a complex phenomenon. Based on the kinetic theory, such a two-phase model has the advantage of being able to clarify many physical concepts, like particulate viscosity, granular conductivity and solid pressure, and deduce the correlative constitutive equations of the solid phase. The AUSM scheme was employed for the numerical calculation. The flow field behind the shock wave was displayed numerically and agrees well with our corresponding experimental results.      

爆轰波与激波对撞爆轰状态演化机理研究

对当量比氢氧混合气体中爆轰波与激波的正面对撞过程进行了二维数值研究. 采用了二阶精度NND差分格式与改进的两阶段化学反应模型,并以数值x-t纹影图以及烟迹图记录了对撞过程. 数值研究表明,透射爆轰波受到膨胀影响首先会衰减,甚至产生局部解耦现象;然后由于三波点的碰撞又能再次重新耦合. 在爆轰波对撞过程中,由于燃烧不均匀性而产生的弱横波对爆轰胞格的形成起着重要作用.      

Shock wave in a gas-liquid medium

The propagation of weak shock waves and the conditions for their existence in a gas-liquid medium are studied in [1]. The article [2] is devoted to an examination of powerful shock waves in liquids containing gas bubbles. The possibility of the existence in such a medium of a shock wave having an oscillatory pressure profile at the front is demonstrated in [3] based on the general results of nonlinear wave dynamics. It is shown in [4, 5] that a shock wave in a gas-liquid mixture actually has a profile having an oscillating pressure. The drawback of [3–5] is the necessity of postulating the existence of the shock waves. This is connected with the absence of a direct calculation of the dissipative effects in the fundamental equations. The present article is devoted to the theoretical and experimental study of the structure of a shock wave in a gas-liquid medium. It is shown, within the framework of a homogeneous biphasic model, that the structure of the shock wave can be studied on the basis of the Burgers-Korteweg-de Vries equation. The results of piezoelectric measurements of the pressure profile along the shock wave front agree qualitatively with the theoretical representations of the structure of the shock wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 65–69, May–June, 1973.     

The interaction between shock waves and solid spheres arrays in a shock tube

When a shock wave interacts with a group of solid spheres, non-linear aerodynamic behaviors come into effect. The complicated wave reflections such as the Mach reflection occur in the wave propagation process. The wave interactions with vortices behind each sphere‘s wake cause fluctuation in the pressure profiles of shock waves. This paper reports an experimental study for the aerodynamic processes involved in the interaction between shock waves and solid spheres. A schlieren photography was applied to visualize the various shock waves passing through solid spheres. Pressure measurements were performed along different downstream positions. The experiments were conducted in both rectangular and circular shock tubes. The data with respect to the effect of the sphere array, size, interval distance, incident Mach number, etc., on the shock wave attenuation were obtained.     

SHPB试验中试件的轴向应力均匀性

针对SHPB试验中试件的轴向应力均匀性问题,采用一维弹性波理论,推导了具有任意形状前沿的入射波加载下,试件内应力的时空分布计算公式。以脉冲前沿的上升时间为参数,将矩形、梯形和坡形3种典型的输入脉冲统一表示为梯形波的形式,计算了不同入射波上升时间和不同试件-压杆波阻抗比情况下试件中的应力传播过程,得到了相应的应力均匀度时程曲线以及应力平衡时间。分析了入射波上升时间和试件-压杆波阻抗比对应力平衡时间的影响,得到了一些有意义的认识,为SHPB试验的设计与分析提供了一定的理论依据。     

Head-on collision of normal shock waves with rigid porous materials

The head-on collision of a planar shock wave with a rigid porous material has been investigated experimentally in a 75 mm × 75 mm shock tube. The experimental study indicated that unlike the reflection from a flexible porous material (e.g., polyurethane foam) where the transmitted compression waves do not converge to a sharp shock wave, in the case of a rigid porous material (e.g., alumina) the transmitted compression waves do converge to a sharp shock wave, which decays as it propagates along the porous material. In addition to this major difference, many other differences were observed. They are outlined in the following sections.Based on these observations a suggestion modifying the phenomenology of the reflection/interaction process in the case a porous material with large permeability is proposed.