Unsteady shock generated vorticity

Abstract. A two-dimensional theory is developed for the vorticity just downstream of a curved, unsteady shock wave. By utilizing Crocco's equation, an explicit formula is obtained for the vorticity that does not require a perfect gas and that holds for arbitrary conditions upstream of the shock wave. The analysis is applied to the flow just downstream of the reflected shock that occurs in a single-Mach reflection pattern. Flow conditions are based on an interferometric photograph of Ben-Dor and Glass (1978). In this case, the reflected shock is weak everywhere from its upstream intersection with the wall to the triple point. The vorticity has a singularity and a change of sign near the triple point that indicates the presence of a weak shear layer downstream of this location. Received 20 December 1999 / Accepted 19 March 2000     

A shock dynamics theory based analytical solution of double Mach reflections

The shock-dynamics based governing equations of double Mach reflection (DMR) wave configurations were developed and solved. The analytical results were compared to experimental results from various sources. Very good to excellent agreement was evident.     

Shock wave calibration of under-expanded natural gas fuel jets

Natural gas, a fuel abundant in nature, cannot be used by itself in conventional diesel engines because of its low cetane number. However, it can be used as the primary fuel with ignition by a pilot diesel spray. This is called dual-fuelling. The gas may be introduced either into the inlet manifold or, preferably, directly into the cylinder where it is injected as a short duration, intermittent, sonic jet. For accurate delivery in the latter case, a constant flow-rate from the injector is required into the constantly varying pressure in the cylinder. Thus, a sonic (choked) jet is required which is generally highly under-expanded. Immediately at the nozzle exit, a shock structure develops which can provide essential information about the downstream flow. This shock structure, generally referred to as a “barrel” shock, provides a key to understanding the full injection process. It is examined both experimentally and numerically in this paper.

---

     

Investigation of shock waves interference and associated hysteresis effect at variable-Mach-number upstream flow

Experimental investigations and numerical simulations have been performed to study the transition between steady regular and Mach reflections induced by flow Mach number variation. The experiments have been carried out in the supersonic wind tunnel SIGMA 4B, at the Institut Aérotechnique (IAT) in St Cyr L'Ecole, France. Symmetric and asymmetric arrangements of wedges have been tested. No significant hysteresis phenomenon has been detected experimentally. However, this phenomenon has been revealed by numerical computations obtained by solving the Navier-Stokes equations.Received: 9 October 2001, Accepted: 11 October 2002, Published online: 21 February 2003     

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     

A weak shock wave reflection over wedges

When a weak shock wave reflects from wedges its reflection pattern does not appear to be a simple Mach reflection. This reflection pattern is known to be von Neumann Mach reflection in which a Mach stem can not necessarily be straight. In this paper the local change of the Mach stem curvature was experimentally and numerically investigated. A distinct triple point, at which the curvature becomes infinite as appears in a simple Mach reflection, was not observed but the Mach stem curvature became a maximum between foot of the Mach stem and a point, P₁, at which an incident shock met with a reflected shock. Maximum curvature point P₂ and P₁ do not coincide for small wedge angles and tend to merge over a certain wedge angle. Experimental results agreed with numerical results. The trajectory angle of P₂ was found to be expressed well by Whitham's shock-shock angle.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.     

Numerical simulations of shock wave reflection phenomena in non-stationary flows using regularized smoothed particle hydrodynamics (RSPH)

In this paper we wish to demonstrate to what extent the numerical method regularized smoothed particle hydrodynamics (RSPH) is capable of modelling shocks and shock reflection patterns in a satisfactory manner. The use of SPH based methods to model shock wave problems has been relatively sparse, both due to historical reasons, as the method was originally developed for studies of astrophysical gas dynamics, but also due to the fact that boundary treatment in Lagrangian methods may be a difficult task. The boundary conditions have therefore been given special attention in this paper. Results presented for one quasi-stationary and three non-stationary flow tests reveal a high degree of similarity, when compared to published numerical and experimental data. The difference is found to be below 5, in the case where experimental data was found tabulated. The transition from regular reflection (RR) to Mach reflection (MR) and the opposite transition from MR to RR are studied. The results are found to be in close agreement with the results obtained from various empirical and semi-empirical formulas published in the literature. A convergence test shows a convergence rate slightly steeper than linear, comparable to what is found for other numerical methods when shocks are involved.     

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.     

Theoretical considerations on shock reflections and their implications on the evaluation of air intake performance

Many theoretical studies have shown the existence of a hysteresis effect in the solution of oblique shock reflections. In fact, a wide domain of free-stream Mach number and shock angle values exists where regular reflection and Mach reflection are both possible solutions for the same flow conditions. Part of this domain overlaps the typical operating conditions of supersonic air intakes, and therefore it is of practical interest to obtain a deeper understanding of the theoretical problem. Indeed, although both solutions are theoretically possible, they yield very different flowfields and consequently large discrepancies in the evaluation of the air intake performance. Numerical solutions for steady configurations have been carried out and compared with the flow evolution obtained for time-dependent cases. The results have confirmed numerically the existence of the multiple solution domain where hysteresis takes place in time-dependent simulations. The analysis of the physical and numerical problems encountered has provided indications for a correct simulation in practical applications. Received 10 August 1999 / Accepted 6 October 2000     

Supersonic liquid jets: Their generation and shock wave characteristics

The generation of high-speed liquid (water and diesel fuel) jets in the supersonic range using a vertical single-stage powder gun is described. The effect of projectile velocity and mass on the jet velocity is investigated experimentally. Jet exit velocities for a set of nozzle inner profiles (e.g. straight cone with different cone angles, exponential, hyperbolic etc.) are compared. The optimum condition to achieve the maximum jet velocity and hence better atomization and mixing is then determined. The visual images of supersonic diesel fuel jets (velocity about 2000 m/s) were obtained by the shadowgraph method. This provides better understanding of each stage of the generation of the jets and makes the study of their characteristics and the potential for auto-ignition possible. In the experiments, a pressure relief section has been used to minimize the compressed air wave ahead of the projectile. To clarify the processes inside the section, additional experiments have been performed with the use of the shadowgraph method, showing the projectile travelling inside and leaving the pressure relief section at a velocity of about 1100 m/s. Received 23 January 2001 / Accepted 2 July 2001     

On the RR→MR transition of asymmetric shock waves in steady flows

In this paper, RR→MR transition of asymmetric shock waves has been theoretically studied. The transition can occur between the sonic-point and maximum-deflection criteria due to the the effects of expansion fans which are inherent flow structures. Comparison shows a better agreement among experiments and the analytical results. Some discrepancies reported in previous studies among experiments and theory have also been explained based on the threshold for RR→MR transition.      

Analysis of the shock structures in a regular detonation

Time-dependent two-dimensional numerical simulations have been used to investigate the detailed shock structures and patterns of energy release in the regions of the triple points and transverse waves in a planar detonation. As the system of shock triple points evolves between collisions, they trace a well shaped cellular pattern characteristic of detonations in argon-diluted, low-pressure mixtures of hydrogen and oxygen. In the region of the triple points, the shock structure evolves continuously from a single Mach structure to a double Mach structure and finally to a complex Mach structure characteristic of spinning detonations. Most of the energy released in the region of the triple points. The amount of energy release increases as the triple point comes closer to a collision with a wall or another triple point. Just before the collision, there is a large region of energy release that covers the length of the interacting transverse waves. The result is a rectangular high-energy region which boosts the propagation of the new detonation cell.     

Study on Mach stems induced by interaction of planar shock waves on two intersecting wedges

The properties of Mach stems in hypersonic corner flow induced by Mach interaction over 3D inter-secting wedges were studied theoretically and numerically. A new method called “spatial dimension reduction” was used to analyze theoretically the location and Mach num-ber behind Mach stems. By using this approach, the problem of 3D steady shock/shock interaction over 3D intersecting wedges was transformed into a 2D moving one on cross sec-tions, which can be solved by shock-polar theory and shock dynamics theory. The properties of Mach interaction over 3D intersecting wedges can be analyzed with the new method, including pressure, temperature, density in the vicinity of triple points, location, and Mach number behind Mach stems. Theoretical results were compared with numerical results, and good agreement was obtained. Also, the influence of Mach number and wedge angle on the properties of a 3D Mach stem was studied.     

浮动冲击平台横向冲击谱受外形结构的影响

基于ABAQUS软件中的声固耦合法,采用设备、平台一体化分析方法,展开对中型浮动冲击平台的研究,探讨平台外形结构形式对平台冲击环境的影响,并提出了提高平台横向冲击谱谱值的设计方案。首先对外形结构进行初步设计,分析其对平台冲击谱的影响,并找到决定影响程度的关键因素。然后针对结构进行优化,使其更大幅度的提高平台冲击谱谱值。计算表明：在平台外部下方加装挡板结构不会明显影响平台垂向冲击谱,但可以提高平台横向冲击谱;舷侧底端向下延伸加装竖直挡板由于受到冲击波绕射和阻力的影响,增加横向谱值的效果不是很明显;平台底部流线型挡板可以有效增加平台对爆炸载荷的接收效果,同时可以尽量减小阻力影响,从而明显提高平台横向冲击谱谱值。

     

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.     

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     

Measurement of laser driven shock wave transit time through thin aluminium targets by optical shadowgraphy

Laser driven shock wave transit time in thin aluminium targets was experimentally estimated by determining the shock emergence time at the rear of thin aluminium foils of varying thickness from 5 to 35 μm. A 20 J, 5 ns Nd:glass laser was focused to produce laser intensity of 10¹² to 5 × 10¹³ W/cm² on the targets which were placed in vacuum. Target foil movement was measured to an accuracy of 10 μm using optical shadowgraphy technique. This technique was used to accurately measure the shock transit time by recording the optical shadowgrams at various instants of time and thus identify the instant at which the foil is just set into motion. Shock transit time measured in foils of different thickness can give the value of shock velocity at a given laser intensity. Target motion recorded by shadowgraphy can also give the target foil velocity from which shock pressure can be estimated. Experimental values of shock transit time, shock velocity and shock pressure were observed to agree well with the values using one-dimensional multi-group radiation hydrodynamic simulations. PACS 52.50Jm; 52.50Lp; 52.25 Communicated by K. Takayama     

Numerical investigation of toroidal shock wave focusing in a cylindrical chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics. PACS 47.40.Ki; 47.40.Nm; 52.35.Tc     

An experiment on imploding conical shock waves

Over the past decade and a half there have been a number of numerical studies of the reflection of oblique axisymmetric shock waves from the axis of symmetry. Many of these have shown a complex Mach reflection pattern together with a strong toroidal vortex which significantly distorts the Mach disk. This geometry has never been captured experimentally. This note describes a special rig for the generation of strong imploding conical shocks, and presents photographs of the predicted reflection pattern. Received 4 October 2001 / Accepted 13 November 2001     

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.