Modification of Whitham's ray shock theory for analyzing the reflection of weak shock waves over small wedge angles

The weak Mach reflection phenomenon has been analyzed by applying both the shock dynamics approach and the disturbance propagation concept. The analysis which is based on modified Whitham's ray shock theory results in analytical expressions for the triple point trajectory angle,, and the shape of the curved Mach stem, which are functions of the incident shock wave Mach number,M _i, and the reflecting wedge angle, _w. The analytical results were found to be in good agreement with experimental results.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.     

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.     

Mach reflection of a plane shock wave passing a mountain of 45° inclination

The transition from regular reflection (RR) to Mach reflection (MR) as a plane shock wave diffracts around a triangular mountain of 45° inclination is analysed in this paper, both by optical measurement in a shock tube and by numerical simulation the numerical method developed by Li Yingfan^[1] is of the FLIC type with triangular mesh. The dependence of the critical transition point Lk ofRR→MR on shock Mach numberM _i is analyzed and the variations of the incidence angle ω _i of the impinging shock and the reflection angle ω _r with the distanceL ^* are investigated. Our experimental and numerical results agree well with the theoretical results of Iton and Italya.     

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.     

On the influence of low initial pressure and detonation stochastic nature on Mach reflection of gaseous detonation waves

The two-dimensional, time-dependent and reactive Navier–Stokes equations were solved to obtain an insight into Mach reflection of gaseous detonation in a stoichiometric hydrogen-oxygen mixture diluted by 25 % argon. This mixture generates a mode-7 detonation wave under an initial pressure of 8.00 kPa. Chemical kinetics was simulated by an eight-species, forty-eight-reaction mechanism. It was found that a Mach reflection mode always occurs for a planar detonation wave or planar air shock wave sweeping over wedges with apex angles ranging from \(5^\circ \) to \(50^\circ \) . However, for cellular detonation waves, regular reflection always occurs first, which then transforms into Mach reflection. This phenomenon is more evident for detonations ignited under low initial pressure. Low initial pressure may lead to a curved wave front, that determines the reflection mode. The stochastic nature of boundary shape and transition distance, during deflagration-to-detonation transition, leads to relative disorder of detonation cell location and cell shape. Consequently, when a detonation wave hits the wedge apex, there appears a stochastic variation of triple point origin and variation of the angle between the triple point trajectory and the wedge surface. As the wedge apex angle increases, the distance between the triple point trajectory origin and the wedge apex increases, and the angle between the triple point trajectory and the wedge surface decreases exponentially.     

An experimental investigation of the sonic criterion for transition from regular to Mach reflection of weak shock waves

The signal speed, namely the local sound speed plus the flow velocity, behind the reflected shocks produced by the interaction of weak shock waves (M _i < 1.4) with rigid inclined surfaces has been measured for several shock strengths close to the point of transition from regular to Mach reflection. The signal speed was measured using piezo-electric transducers, and with a multiple schlieren system to photograph acoustic signals created by a spark discharge behind a small aperture in the reflecting surfaces. Both methods yielded results with equal values within experimental error. The theoretical signal speeds behind regularly reflected shocks were calculated using a non-stationary model, and these agreed with the measured results at large angles of incidence. As the angle of incidence was reduced, for the same incident shock Mach number, so as to approach the point of transition from regular to Mach reflection, the measured values of the signal speed deviated significantly from the theoretical predictions. It was found, within experimental uncertainty, that transition from regular to Mach reflection occurred at the experimentally observed sonic point, namely, when the signal speed was equal to the speed of the reflection point along the reflecting surface. This sonic condition did not coincide with the theoretical value.     

Shock wave triple-point morphology

The existence and characteristics of shock wave triple points are examined. The analysis utilizes a single flow plane for the three shocks and is local to the triple point. It applies when the flow is unsteady, three-dimensional, and the upstream flow is nonuniform. Under more restrictive conditions, a relation is also derived for the ratio of the curvature of the Mach stem to that of the reflected shock. For given values of the ratio of specific heats, γ, and the upstream Mach number, M ₁, a solution window is established. A parametric set of solutions is generated within the window for γ = 1, 1.4, and 5/3 and for 16 values of M ₁ ranging from solution onset to M ₁ = 6.A solution can be one of three types, these stem from the velocity tangency condition along the slip stream. Topics are addressed such as solution multiplicity, shock wave and slip stream orientation, the nature of the reflected wave (weak, strong, inverted, normal), the nature of the Mach stem (weak, strong, normal), and differences due to changes in γ and M ₁.     

Normal shock interaction with a yawed wedge moving at supersonic speed

In this article, the interaction of a normal shock with a yawed wedge moving at supersonic speed has been considered. The vorticity distribution of a particle over the diffracted shock wave for various combinations of yawed angles, Mach number of the shock wave and Mach number of the moving wedge have been obtained. Further triple point angle χ in Mach reflection has been calculated for the various parameters.      

Effect of flow regime change from subsonic to transonic on the air loads of an oscillating airfoil

In this research, the effect of flow regime change from subsonic to transonic on the air loads of a pitching NACA0012 airfoil is investigated. To do this, the effect of change in flow regime on the lift and pitching moment coefficients hysteresis cycles is studied. The harmonic balance approach is utilized for numerical calculation due to its low computational time. Verifications are also made with previous works and good agreements are observed. The assessment of flow regime change on the aforementioned hysteresis cycles is accomplished in the Mach number range of M=0.65–0.755. The reduced frequency and pitch amplitude also vary from k=0.03 to 0.1 and α₀=1–2.51°, respectively. Results show that the effect of increase in Mach number is to increase and decrease the lift coefficient during downstroke and upstroke, respectively, whereas at low reduced frequencies, the effect of increase in Mach number may lead to a reverse manner when airfoil moves toward its extremum angle of attack. Results also reveal that as the pitch amplitude varies, the shape of lift coefficient hysteresis cycle depends more on the pitch amplitude than on the appearance of shock. It is shown that as the Mach number increases, the incidence angles correspond to the extremum pitching moment, and depending on the reduced frequency, lie between zero and extremum angle of attack. These incidence angles shift toward the extremum angle of attack as the reduced frequency decreases. Results also show that the increase in pitch amplitude at low Mach number, in such a way that leads to the formation of shock around the extremum angle of attack, causes the extremum pitching moment to appear around these angles and at high Mach number, depending on the reduced frequency, the extremum pitching moment incidence angles would be between zero and extremum incidence angle.     

Wave dynamic processes in cellular detonation reflection from wedges

When the cell width of the incident detonation wave (IDW) is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H ₂/O ₂ diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes. The English text was polished by Yumming Chen.     

Incidence of plane supersonic jet on a plane at an arbitrary angle

A solution method is described and results are presented of numerical calculations for the problem of determining the subsonic part of the flow with incidence of a plane uniform supersonic jet on a plane at an arbitrary angle, which corresponds to the flow regime with a a detached shock wave. For the problem solution we use the method of integral relations in the first approximation in a polar coordinate system. The calculation results (pressure distribution on the plate, shock wave shape, and velocity gradient magnitude at the stagnation point) are presented for Mach numbers of 5 and 20, in a range of incidence angles from 0° to 35°, and also for M=3 for an incidence angle 0° (angles measured from normal to the plate).In conclusion the authors wish to thank G. I. Taganov for guidance in this work.     

Oblique reflection of shock waves from rigid porous materials

A study to determine the general gas dynamic behaviour associated with the impact of a shock wave on a porous wedge has been undertaken. A number of interesting features are noted. The pattern of wave reflection is shown to be significantly affected by the inflow of gas into the wedge. This has the effect of reducing the triple point trajectory angle for cases of Mach reflection and for strongly reducing the reflection angle in regular reflection. The permeability of the wedge has a significant effect on the strength of the reflected wave and in some cases this wave can be attenuated to the extent that it is almost eradicated. Pressure measurements taken under the wedge are characterized by oscillations which are of similar shape, for a given wedge, over a range of shock wave Mach numbers. It is shown that the wave transmitted into the wedge is attenuated to varying degrees depending on the material properties, and that for weak incident waves the mean propagation velocity can be less than the sound speed in the pore fluid. Photographs taken using a specially constructed wedge which allows the transmitted wave to be visualised, show that the transmitted wave is nearly plane.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.     

Vorticity and its rate of change just downstream of a curved shock

A theory is developed for the vorticity and its substantial derivative just downstream of a curved shock wave, the resulting formulas are exact, algebraic, and explicit. Analysis is for a cylinder-wedge or sphere-cone body, at zero incidence, whose downstream half-angle is θ_b. Derived formulas directly depend only on the ratio of specific heats, γ, the freestream Mach number, M ₁, the local slope and curvature of the shock, and the dimensionality parameter, σ, which is zero for a two-dimensional shock and unity for an axisymmetric shock. In turn, the slope and curvature depend on γ, M ₁, and θ_b. Numerical results are provided for a bow shock in which θ_b is 5°, 10°, or 15°, M ₁ is 2, 4, or 6, and γ = 1.4. There is little dependence on the half angle but a strong dependence on the freestream Mach number and on dimensionality. For vorticity and its substantial derivative, the dimensionality dependence gradually decreases with increasing Mach number. In comparison to the two-dimensional case, an axisymmetric shock generates considerable vorticity in a region relatively close to the symmetry axis. Moreover, the magnitude of the vorticity, in this region, is further enhanced in the flow downstream of the shock. This dimensionality difference in vorticity and its substantial derivative is attributed to the three-dimensional relief effect in an axisymmetric flow.

---

     

Experimental study of the shock–bubble interaction with reshock

The interaction of a planar shock wave with a spherical density inhomogeneity is studied experimentally under reshock conditions. Reshock occurs when the incident shock wave, which has already accelerated the spherical bubble, reflects off the tube end wall and reaccelerates the inhomogeneity for a second time. These experiments are performed at the Wisconsin Shock Tube Laboratory, in a 9m-long vertical shock tube with a large square cross section (25.4×25.4 cm²). The bubble is prepared on a pneumatically retracted injector and released into a state of free fall. Planar diagnostic methods are used to study the bubble morphology after reshock. Data are presented for experiments involving two Atwood numbers (A = 0.17 and 0.68) and three Mach numbers (1.35 < M < 2.33). For the low Atwood number case, a secondary vortex ring appears immediately after reshock which is not observed for the larger Atwood number. The post-reshock vortex velocity is shown to be proportional to the incident Mach number, M, the initial Atwood number, A, and the incident shock wave speed, W _i.     

Visualization of shock-wave formation processes during shock reflection at obstacles with multiple steps

According to standard textbooks on compressible fluid dynamics, a shock wave is formed by an accumulation of compression waves. However, the process by which an accumulated compression wave grows into a shock wave has never been visualized. In the present paper, the authors tried to visualize this process using a model wedge with multiple steps. This model is useful for generating a series of compression waves and can simulate a compression process that occurs in a shock tube. By estimating the triple-point trajectory angle, we demonstrated visually that an accumulated compression wave grows into a shock wave. Further reflection experiments over a rough-surface wedge confirmed the tendency for the triple point trajectory angle to reach the asymptotic value _s in the end.This work was first presented at the Symposium on Shock Waves, Japan 2002     

Effect of Partial Closure of the Channel on the Pressure Pulse in a Shock Wave Emerging from a Duct

The temporal-spatial development of the flow pattern behind a shock wave emerging from open-ended and half-closed ducts at M₀=1.15 to 3.0 Mach number is investigated both experimentally and numerically and the action of the diffracted wave on a barrier placed at different distances from the duct exit is studied. Flow toeplerograms are obtained and the pressure on the barrier is measured. The distinctive features of the interaction between the shock wave and a barrier mounted perpendicular to the duct axis are established. It is found that partial closing of the duct results in a decrease in the barrier pressure in the case of a strong shock (M₀ > 2.2) and leads to a pressure decrease when a weak shock is diffracted (M₀=1.1 to 1.7). A dependence characterizing the dynamic action of the shock wave on the barrier and specifying the threshold value of a combination of the shock Mach number and the distance from the barrier determining whether the pressure pulse on the barrier increases or decreases, is obtained.     

Nonlinear Strong Shock Interactions: A Shock-Fitted Approach

This paper addresses nonlinear effects which result from the interaction of shock waves with vortices. A series of experiments are carried out, which involve the interaction of a strong shock wave with a single plane vorticity wave and a randomly distributed wave system. These experiments are first conducted in the linear regime to obtain a mutual verification of theory and computation. They are subsequently extended into the nonlinear regime. A systematic study of the interaction of a plane shock wave and a single vortex is then conducted. Specifically, we investigate the conditions under which nonlinear effects become important, both as a function of shock Mach number, M ₁, and incident vortex strength (characterized by its circulation Γ). The shock Mach number is varied from 2 to 8, while the circulation of the vortex is varied from infinitesimally small values (linear theory) to unity. Budgets of vorticity, dilatation, and pressure are obtained. They indicate that nonlinear effects become more significant as both the shock Mach number and the circulation increase. For Mach numbers equal to 5 and above, the dilatation in the vortex core grows quadratically with circulation. An acoustic wave propagates radially outward from the vortex center. As circulation increases, its upstream-facing front steepens at low Mach numbers, and its downstream-facing front steepens at high Mach numbers. A high Mach number asymptotic expansion of the Rankine--Hugoniot conditions reveals that nonlinear effects dominate both the shock motion and the downstream flow for ΓM ₁ > 1. Received 28 June 1997 and accepted 25 November 1997     

Effect of incident shock wave strength on the decay of Richtmyer–Meshkov instability-introduced perturbations in the refracted shock wave

The effect of incident shock wave strength on the decay of interface introduced perturbations in the refracted shock wave was studied by performing 20 different simulations with varying incident shock wave Mach numbers (M ~ 1.1? 3.5). The analysis showed that the amplitude decay can be represented as a power law model shown in Eq.7, where A is the average amplitude of perturbations (cm), B is the base constant (cm^?(E?1), S is the distance travelled by the refracted shockwave (cm), and E is the power constant. The proposed model fits the data well for low incident Mach numbers, while at higher mach numbers the presence of large and irregular late time oscillations of the perturbation amplitude makes it hard for the power law to fit as effectively. When the coefficients from the power law decay model are plotted versus Mach number, a distinct transition region can be seen. This region is likely to result from the transition of the post-shock heavy gas velocity from subsonic to supersonic range in the lab frame. This region separates the data into a high and low Mach number region. Correlations for the power law coefficients to the incident shock Mach number are reported for the high and low Mach number regions. It is shown that perturbations in the refracted shock wave persist even at late times for high incident Mach numbers.     

气相爆轰波在障碍物上Mach反射的实验验证

本文公布了气相爆轰波沿收缩管道传播时发生Mach反射的实验证据。在爆轰波通过的管道中安装不同楔角的楔块,形成管道的收缩。爆轰波在通过楔块时会发生Mach反射。利用烟熏玻璃片记录到了爆轰波Mach反射时形成的三波点迹线及其两侧胞格尺寸和密度的变化。据我们掌握的资料,这是首次用胞格结构变化的记录证实,气相爆轰波与无化学反应的空气中的冲击波一样,在一定的入射条件下会发生Mach反射。这一实验结果可使我们更深入了解爆轰波的本质,也为数值模拟气相爆轰波在障碍物上Mach反射现象提供了可对比的依据。     

Interaction of a shock with a sphere suspended in a vertical shock tube

Shock wave interaction with a sphere is one of the benchmark tests in shock dynamics. However, unlike wind tunnel experiments, unsteady drag force on a sphere installed in a shock tube have not been measured quantitatively. This paper presents an experimental and numerical study of the unsteady drag force acting on a 80 mm diameter sphere which was vertically suspended in a 300 mm x 300 mm vertical shock tube and loaded with a planar shock wave of M _s = 1.22 in air. The drag force history on the sphere was measured by an accelerometer installed in it. Accelerometer output signals were subjected to deconvolution data processing, producing a drag history comparable to that obtained by solving numerically the Navier-Stokes equations. A good agreement was obtained between the measured and computed drag force histories. In order to interpret the interaction of shock wave over the sphere, high speed video recordings and double exposure holographic interferometric observations were also conducted. It was found that the maximum drag force appeared not at the time instant when the shock arrived at the equator of the sphere, but at some earlier time before the transition of the reflected shock wave from regular to Mach reflection took place. A negative value of the drag force was observed, even though for a very short duration of time, when the Mach stem of the transmitted shock wave relfected and focused at the rear stagnation point of the sphere.Received: 31 March 2003, Accepted: 7 July 2003, Published online: 2 September 2003