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     

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.     

Interaction of a shock wave with a two-phase interface

The acceleration by an incident shock of a planar interface between a gas and a particle-gas mixture has been investigated experimentally and numerically. The experiments were conducted in a newly developed vertical shock tube in which the planar interface of the particle-gas mixture was generated and its particle concentration history was measured. Polydisperse corn starch particles with a mean diameter of 10m were used. We recorded the motion of the interface, as well as of the incident and reflected shock by using a 4 channel spark shadowgraph. The experimental conditions were Mach numberM _s=5.15 and initial pressurep ₁=50kPa for various particle concentrations in nitrogen. The reflected shock appears with a delay after the incident shock enters the particle-gas mixture. Numerical methods were employed to solve the two-phase governing equations. Experiments and numerical solutions are in good agreement.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.     

Three-dimensional propagation of the transmitted shock wave in a square cross-sectional chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.PACS: 43.40.Nm     

Shock wave interaction with a cloud of particles

The present paper is devoted to experimental and theoretical investigation of the shock wave (SW) propagation in a mixture of gas and solid particles in the presence of explicit boundaries of the two-phase region (cloud of particles). The effect of the qualitative change in the supersonic flow behind the SW in a cloud of particles within the range of the volume concentration of the disperse phase 0.1-3% is experimentally shown and theoretically grounded. Received 15 April 1996 / Accepted 3 June 1996     

A non-equilibrium multiscale simulation of shock wave propagation

A novel non-equilibrium multiscale dynamics (NEMSD) is proposed to simulate non-equilibrium thermal–mechanical processes. The model couples coarse-grain thermodynamics with a fine scale molecular dynamics. A Distributed Nośe-Hoover Thermostat Network is used, which regulates the temperature in each coarse scale Voronoi cell according to the finite element (FE) nodal temperature. The atoms in each element-cell, namely Voronoi cell-ensemble, are assumed to be in a local equilibrium state within one coarse scale time step. The change of FE nodal temperature provides a source of random forces, which drive the system out of equilibrium. The proposed NEMSD can successfully simulate shock wave propagation in a cubic lattice.     

Numerical simulations of shock wave propagation in condensed multiphase materials

We propose to find out numerical solutions of a travelling shock wave in condensed mixtures by using a direct numerical simulation. Condensed multiphase materials under shock wave conditions are mechanically characterized by a unique pressure and a unique velocity. In this study, the mixture is considered as a collection of grains separated by interface between each material: this problem of interfaces is solved by a diffuse interface method. The results will be compared to existing one-dimensional numerical models, analytical solutions and also experimental data. The volume fraction (or the phase temperature) is not measured in experiments and it is then important to verify the behaviour of a phase quantity through various methods. A non-monotonous evolution of the volume fraction is obtained with analytical solution as well as numerical simulation.      

Gas combination for shock wave enhancement

It was recently demonstrated that shock wave enhancement could be achieved when a shock propagates in a constant cross-section duct through pairs of air–helium layers having a continually decreasing width (Igra and Igra in Shock Waves 16(3):199–207). A parametric study was conducted aimed at finding a two-layered, light–heavy gas arrangement that yields maximal shock enhancement; the heavy and the light gases used were air and helium, respectively. Effects associated with changes in following parameters were investigated: the number of alternating heavy/light gas layers, the applied reduction ratio between successive layers thickness, and the initial shock wave Mach number.      

Simulation of blast wave propagation over a cylinder

A numerical study of the interaction of plane blast waves with a cylinder is presented. Computations are carried out for various blast-wave durations and comparisons are obtained with the corresponding results of planar shock-wave. Both inviscid and viscous results based on the solution of the Euler and Navier-Stokes equations are presented. The equations are solved by an adaptive-grid method and a second-order Godunov scheme. The shock wave diffraction over the cylinder is investigated by means of various contour plots, as well as, pressure and skin-friction histories. The study reveals that the blast-wave duration significantly influences the unsteady flow over the cylinder. The differences between the viscous and inviscid results are also discussed. Received 2 March 1996 / Accepted 28 February 1997     

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.     

Experimental and numerical studies of underwater shock wave attenuation

The attenuation of an underwater shock wave by a thin porous layer is studied both experimentally and numerically. The shock waves are generated by exploding 10 mg silver azide pellets and the pressures at different distances from the explosion center are measured. Measurements are also carried out with a gauze layer placed between the explosion source and the pressure gauge. The results with and without the gauze layer are compared evaluating the shock wave attenuation. Numerical simulations of the phenomenon are also carried out for a simple wave attenuation model. The results are compared with the experimental data. Despite the simple mathematical model of wave attenuation, the agreement between the experimental and numerical results is reasonable.Received: 22 October 2002, Accepted: 17 June 2003, Published online: 5 August 2003PACS: 47.11.+j, 47.40.Nm, 47.55.Mh     

激波在带扩容室的变截面管中传播和消波效果的数值分析

在利用实验结果验证程序的基础上,对激波在具有复杂结构管道内的传播特性进行了数值模拟,比较了不同入射激波条件和特征尺寸对扩容室消波效应的影响,并考虑了二维(轴对称、平面)模型和三维模型的影响。计算结果表明,同一尺寸的扩容室对不同马赫数入射激波的消波效果是相近的,扩容室特征尺寸不同则消波效果差别很大。     

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.     

Color schlieren methods in shock wave research

Schlieren methods are widely known and well established to visualize refractive index variations in transparent media. The use of color allows one to obtain more data and previously inaccessible information from a picture taken with this technique. In general, a hue can be related to a certain strength or a certain direction of a refractive index gradient. While the first case essentially corresponds to the usual black- and-white system the latter correlation cannot be made adequately evident without the use of color. Two color schlieren techniques are presented here, which reach or even exceed the quality and sensitivity range of conventional black- and-white methods. Using a powerful short duration light source these methods are applied to visualize transient flow phenomena in a shock tube.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Shock wave diffraction by a square cavity filled with dusty gas

A simple two-dimensional square cavity model is used to study shock attenuating effects of dust suspension in air. The GRP scheme for compressible flows was extended to simulate the fluid dynamics of dilute dust suspensions, employing the conventional two-phase approximation. A planar shock of constant intensity propagated in pure air over flat ground and diffracted into a square cavity filled with a dusty quiescent suspension. Shock intensities were and , dust loading ratios were and , and particle diameters were and {\rm \mu}$m. It was found that the diffraction patterns in the cavity were decisively attenuated by the dust suspension, particularly for the higher loading ratio. The particle size has a pronounced effect on the flow and wave pattern developed inside the cavity. Wall pressure histories were recorded for each of the three cavity walls, showing a clear attenuating effect of the dust suspension. Received 15 November 1999 / Accepted 25 October 2000     

Discrete element simulation of shock wave propagation in polycrystalline copper

The implementation of the characteristic of compressive plasticity into the Discrete Element Code, DM2, while maintaining its quasi-molecular scheme, is described. The code is used to simulate the shock compression of polycrystalline copper at 3.35 and 11.0 GPa. The model polycrystal has a normal distribution of grain sizes, with mean diameter 14 μm, and three distinct grain orientations are permitted with respect to the shock direction; 〈1 0 0〉, 〈1 1 0〉, and 〈1 1 1〉. Particle velocity dispersion (PVD) is present in the shock-induced flow, attaining its maximum magnitude at the plastic wave rise. PVD normalised to the average particle velocity of and are yielded for the 3.35 and 11.0 GPa shocks, respectively, and are of the same order as those seen in the experiment. Non-planar elastic and plastic wave fronts are present, the distribution in shock front position increasing with propagation distance. The rate of increase of the spread in shock front positions is found to be significantly smaller than that seen in probabilistic calculations on nickel polycrystals, and this difference is attributed, in the main, to grain interaction. Reflections at free surfaces yield a region of tension near to the target free surface. Due to the dispersive nature of the shock particle velocity and the non-planarity of the shock front, the tensile pressure is distributed. This may have implications for the spall strength, which are discussed. Simulations reveal a transient shear stress distribution behind the shock front. Such a distribution agrees with that put forward by Lipkin and Asay to explain the quasi-elastic reloading phenomenon. Simulation of reloading shocks show that the shear stress distribution can give rise to quasi-elastic reloading on the grain scale.     

Attenuation of shock waves propagating in polyurethane foams

Shock wave attenuation in polyurethane foams is investigated experimentally and numerically. This study is a part of research project regarding shock propagation in polyurethane foams with high-porosities = 0.951 ~ 0.977 and low densities of ρc = 27.6 ~55.8 kg/m³. Sixty Millimeter long cylindrical foams with various cell numbers and foam insertion condition were installed in a horizontal shock tube of 50 mm i.d. and 5.4 mm in length. Results of pressure measurements in air/foam combination are compared with CFD simulation solving the one-dimensional Euler equations. In the case of a foam B fixed on shock tube wall, pressures at the shock tube end wall increases relatively slowly comparing to non-fixed foam, free to move and a foam A fixed on shock tube wall. This implies that elastic inertia hardly contributes to pressure build up. Pressures behind a foam C fixed on shock tube wall decrease indicating that shock wave is degenerated into compression wave. Dimensionless impulse and attenuation factor decrease as the initial cell number increases. The momentum loss varies depending on cell structure and cell number.     

Oblique shock wave with sweep

An attached planar, oblique shock with sweep is investigated for the inviscid flow of a perfect gas. The ratio of specific heats, freestream Mach number, and wedge angle in the plane of the freestream velocity are prescribed, with the sweep angle as a free parameter. Explicit relations are provided for jump and detachment conditions. A number of trends, some non-intuitive, are discussed, e.g., the downstream Mach number may increase with sweep.     

Experimental and numerical investigation of the shock-induced fluidization of a particles bed

An experimental study on unsteady two phase flow is conducted in a vertical shock tube. Shock Mach numbers range from 1.3 to 1.5 in 1 bar. The particles are initially positioned in horizontal beds of various thicknesses. Our research covers a large domain of void fraction from 1 (single particles) to 0.35 (compact beds). The experiments provide shadowgraph images for the recording of particle trajectories (effect of the gas on the particles) and side-wall pressures (action of the particles on the gas). A dense two phase flow model has been elaborated and numerically solved using a finite difference scheme with pseudoviscosity. The simulated shock-induced fluidization of a 2 cm thick bed of 1.5 mm diameter glass particles is compared to the experiment. Received 10 September 1996 / Accepted 4 January 1997