Shock-tube experiments on the stability of regular reflection in the dual-solution domain

Regular reflection (RR) and Mach reflection (MR) are theoretically both possible in the dual-solution domain of oblique shock reflection. The physical difference between the two types of reflection is the pressure behind the reflected shock wave: that of MR is lower than that of RR for strong shock reflection. The magnitude relation of these pressures is inverted for weak shock reflection. In the present paper, we performed two kinds of experiment, depending on whether the oblique shock reflection is weak or strong. For strong shock reflection, we decreased the pressure behind the reflected wave of RR using a convex double-wedge. For weak shock reflection, we increased this pressure using a concave double-wedge. Thus, we investigated the stability of RR against pressure disturbances. The results indicate that RR in a shock tube is stable, in the dual-solution domain, for both weak and strong incident shocks.     

Existence criteria for reflection configurations in shock–vortex interactions

The length scale criteria is widely accepted as an explanation for transition and hence existence of different shock wave reflection configurations in pseudo-steady flows. However, there has not been any attempt to validate this criteria using information obtained from a time-dependent numerical simulation. A high resolution time-dependent numerical simulation in pseudo-steady flow is carried out in the present work. Time-dependent numerical data is used to calculate flow features in a laboratory frame of reference to verify validity of the length scale criteria for existence of different shock wave reflection configurations in pseudo-steady flow. This analysis is then extended to the study of unsteady shock wave reflection configurations in shock–vortex interactions. It is shown that the existence of regular reflection (RR) and Mach reflection (MR) configurations in an unsteady flowfield resulting from shock–vortex interactions can also be explained locally based on limiting conditions similar to that prescribed by the length scale criteria for pseudo-steady flow.

---

     

Head-on Collision of a Detonation with a Planar Shock Wave

The phenomenon that occurs when a Chapman–Jouguet (CJ) detonation collides with a shock wave is discussed. Assuming a one-dimensional steady wave configuration analogous to a planar shock–shock frontal interaction, analytical solutions of the Rankine–Hugoniot relationships for the transmitted detonation and the transmitted shock are obtained by matching the pressure and particle velocity at the contact surface. The analytical results indicate that there exist three possible regions of solutions, i.e. the transmitted detonation can have either strong, weak or CJ solution, depending on the incident detonation and shock strengths. On the other hand, if we impose the transmitted detonation to have a CJ solution followed by a rarefaction fan, the boundary conditions are also satisfied at the contact surface. The existence of these multiple solutions is verified by an experimental investigation. It is found that the experimental results agree well with those predicted by the second wave interaction model and that the transmitted detonation is a CJ detonation. Unsteady numerical simulations of the reactive Euler equations with both simple one-step Arrhenius kinetic and chain-branching kinetic models are also carried out to look at the transient phenomena and at the influence of a finite reaction thickness of a detonation wave on the problem of head-on collision with a shock. From all the computational results, a relaxation process consisting of a quasi-steady period and an overshoot for the transmitted detonation subsequent to the head-on collisions can be observed, followed by the asymptotic decay to a CJ detonation as predicted theoretically. For unstable pulsating detonations, it is found that, due to the increase in the thermodynamic state of the reactive mixture caused by the shock, the transmitted pulsating detonation can become more stable with smaller amplitude and period oscillation. These observations are in good agreement with experimental evidence obtained from smoked foils where there is a significant decrease in the detonation cell size after a region of relaxation when the detonation collides head-on with a shock wave.     

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.     

Dynamics of stress wave propagation in a chain of photoelastic discs impacted by a planar shock wave; Part I, experimental investigation

The propagation of stress waves through a chain of discs has been studied experimentally. Optically transparent 20-mm diameter discs, made of epoxy, were loaded dynamically by head-on collision with an incident planar shock wave. The loading was done in a vertical shock tube. The head-on collision between the punch-plate, placed on top of the chain of discs, and the incident shock wave resulted in a head-on reflected shock wave inducing behind it a fairly uniform step-wise pressure pulse having duration of about 6 ms. The recorded fringe patterns of the stress field, in the discs-chain, show that the input pressure pulse was broken into several oscillating cycles. The back and forth bouncing of stress waves gave rise to two different modes of the contact stress oscillations, which continued until the overall stress reaches equilibrium with the input conditions. The registered propagation velocity of the stress wave was significantly lower than the appropriate speed of sound in the material from which the discs were made.      

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.      

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.     

Transition between regular and Mach reflection of shock waves: new numerical and experimental results

New numerical and experimental results on the transition between regular and Mach reflections of steady shock waves are presented. The influence of flow three-dimensionality on transition between steady regular and Mach reflection has been studied in detail both numerically and experimentally. Characteristic features of 3D shock wave configuration, such as peripheral Mach reflection, non-monotonous Mach stem variation in transverse direction, the existence of combined Mach-regular-peripheral Mach shock wave configuration, have been found in the numerical simulations. The application of laser sheet imaging technique in streamwise direction allowed us to confirm all the details of shock wave configuration in the experiments. Close agreement of the numerical and experimental data on Mach stem heights is shown. Received 23 November 2000 / Accepted 25 April 2001     

Numerical simulation of the reflection of a planar shock wave over a double wedge

A numerical code based on the upwind TVD scheme for simulating the various reflection processes of a planar shock wave over a concave or convex double wedge has been developed. The numerical results were compared with actual experiments and excellent agreement was obtained. The excellent agreement serves also as a validation of the shock-capturing performance of the numerical scheme.     

Three-dimensional shock wave reflection over a corner of two intersecting wedges

This paper presents an experimental and numerical investigation of three-dimensional shock wave reflections over a corner of two wedges intersecting perpendicularly in a shock tube. Experiments were conducted in a diaphragmless shock tube equipped with double-exposure diffuse holographic interferometry in which the time interval between the first and second exposure was set to be . This arrangement clearly visualized complex configurations of three-dimensional shock wave reflections. A numerical study was also carried out for interpreting these holographic interferometric observations by using the Weighted Average Flux (WAF) method to solve the three-dimensional unsteady compressible Euler equations. It was found that along the line of the intersection of these two wedges, two Mach stems intersected each other resulting in the formation of a Mach stem which leaned forward. Received 30 June 1996 / Accepted 6 October 1996     

Formation of a jet and vortices behind a shock wave reflected from a concave body

A jet and vortices have been observed when a plane shock wave reflects from a concave body in a shock tube. If the cavity is deep enough then two reflected shocks appear near its edges. Air, carbon tetrafluoride (CF) and dichlorodifluoromethane (CClF) were chosen as test gases. The flow was visualized with the aid of a conventional shadow technique. Pressure measurements at the body surface were also obtained. Numerical studies have been conducted using a two-dimensional inviscid model. There is a good qualitative agreement between the experimental and numerical results. Received 8 February 1996 / Accepted 30 June 1997     

LES of spatially developing 3D compressible mixing layer

The spatial development of a turbulent compressible mixing layer is investigated by means of large eddy simulation (LES). The subgrid viscosity is represented by the so-called mixed-scale model, adapted to compressible flows. Two different shock capturing schemes and three sets of inlet white-noise perturbations are investigated. The comparison between numerical and experimental results gives an overall good agreement.     

Numerical investigation on three-dimensional shock wave reflection over two perpendicularly intersecting wedges

The three-dimensional (3D) shock wave reflections over two perpendicularly intersecting wedges are numerically studied in this paper, using the finite volume method which is based on the MUSCL-Hancock interpolation technique and self-adaptive unstructured mesh. Two kinds of 3D Mach stem structures are demonstrated by the numerical simulations for different shock Mach numbers and wedge angles. A four-shock or three-shock wave configuration appears in the vicinity of the corner of the wedges.     

非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

Height of burst explosions: a comparative study of numerical and experimental results

In the current work, we use the Constant Volume model and the numerical method, Regularized Smoothed Particle Hydrodynamics (RSPH) to study propagation and reflection of blast waves from detonations of the high explosives C-4 and TNT. The results from simulations of free-field TNT explosions are compared to previously published data, and good agreement is found. Measurements from height of burst tests performed by the Norwegian Defence Estates Agency are used to compare against numerical simulations. The results for shock time of arrival and the pressure levels are well represented by the numerical results. The results are also found to be in good agreement with results from a commercially available code. The effect of allowing different ratios of specific heat capacities in the explosive products are studied. We also evaluate the effect of changing the charge shape and height of burst on the triple point trajectory.      

The Mach reflection triple-point locus for internal and external conical diffraction of a moving shock wave

This paper examines the different behavior that occurs for the Mach reflection triple-point loci between the two fundamental axisymmetric cases, these being the external diffraction by a cone and the internal diffraction within a conically contracting channel. From equations derived in this paper using a shock dynamics approach, it has been shown that, for external diffraction over a cone, a possible solution is that the triple-point locus is a straight line which corresponds to the experimental results available, while for internal diffraction along a conically converging channel, it cannot be straight and is, in fact, a convex curve. In the latter case, a transition point is noted on the triple-point locus before which the locus is nearly straight but after which the curvature becomes marked. The second region diminishes as a proportion of the total locus with decreasing half cone angle.For the external case, a set of simple, axisymmetric equations are derived which allow a rapid estimation of the triple point locus angle and the Mach stem strength for any incident shock Mach number and cone angle combination. The equations for internal diffraction are similar and allow a quick computation of both the curved triple-point locus and the strength of the diffracting front of the shock wave. A comparison with experiment has been carried out and agreement is good.     

金属铅在斜冲击波对碰加载下动态行为的数值模拟

高强度冲击加载作用下金属材料的动态物理行为是当前冲击波领域基础研究和工程应用最为关注的焦点。采用光滑粒子法(SPH)开展不同位置起爆诱发的斜冲击波对碰加载金属铅的二维数值模拟研究,得到了金属铅内入射斜冲击波的角度和强度,并利用极曲线方法理论上导出发生马赫反射时的临界入射角和入射马赫数关系。根据计算结果可知,金属铅内入射斜波对碰后将发生马赫反射。随着起爆位置与金属铅表面距离的增加,不仅金属铅内入射冲击波强度和入射角增加,而且形成的马赫杆宽度也在增加。由自由面速度剖面给出了马赫杆宽度及张角,结果与理论预测的结果吻合较好。     

Numerical solutions of the Boltzmann equation: comparison of different algorithms

In the paper we compare different algorithms for numerical solutions of the Boltzmann equation. For this comparison we have taken the standard problem of the shock wave structure in a mono-atomic rarefied gas. Different parameters characterizing the shock structure have been calculated by a Monte Carlo simulation (DSMC), a second order time-relaxed Monte Carlo method (TRMC2), a fully deterministic discrete velocity method (DV), a discrete velocity method with Monte Carlo calculations of collision integral (DVMC) and a molecular dynamics method (MD). Results of these calculations have been compared with the shock wave structure obtained in experiments in a shock tube. The results of the comparison are not conclusive. We have observed general agreement between numerical and experimental results but there is no single numerical method which fits best to the experimental measurements.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison with experiment for TVD calculations of blast waves from a shock tube

Harten's second-order-accurate total-variation-diminishing (TVD) scheme is applied to calculation of flow from the open end of a shock tube. Comparison of numerical results with available experimental data for overpressure at selected points around the shock tube exit shows good agreement. Numerically indicated positions of the moving shock front and Mach stem also compare well with flow shadowgraph data. Where the problem geometry is sufficiently simple and rectangular gridding can be used, Harten's method affords a good choice for blast wave calculations.