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     

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.     

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.      

Experiments on Guderley Mach reflection

Experiments have been conducted in a large shock tube to examine the four-wave shock reflection pattern, now known as Guderley reflection (GR). The fourth wave, an expansion, is clearly identified, as is the supersonic patch behind the reflected wave. A shocklet terminating the supersonic patch behind the reflected wave is identified, which forms a second triple point further down the Mach stem. Evidence is presented showing the presence of more than one expansion wave and more than one shocklet, thus indicating the existence of more than one supersonic patch. In order to distinguish between cases with a single patch without the shocklet as originally proposed by Guderley and found in some computations, and the indications of a multi-patch geometry found here, and also in other computations, this latter case is designated Guderley Mach reflection (GMR). Multi-exposure images of the shock propagation superimposed on a single image frame enable estimates to be made of the strength of the major waves, and it is shown that the reflected wave is very weak.      

Three-dimensional steady supersonic duct flow using Lagrangian formulation

In this paper, numerical simulation of three-dimensional supersonic flow in a duct is presented. The flow field in the duct is complex and can find its applications in the inlet of air-breathing engines. A unique streamwise marching Lagrangian method is employed for solving the steady Euler equations. The method was first initiated by Loh and Hui (1990) for 2-D steady supersonic flow computations and then extended to 3-D computation by the present authors Loh and Liou (1992). The new scheme is shown to be capable of accurately resolving complicated shock or contact discontinuities and their interactions. In all the computations, a free stream of Mach numberM=4 is considered.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.     

Investigation on onset of shock-induced separation

A great number of experimental data indicating shock wave/boundary layer interactions in internal or external supersonic flows were reviewed to make clear the mechanism of the interaction and to decide the onset of shock-induced separation. The interesting conclusions were obtained for the considerably wide range of flow geometries that the onset of separation is independent of the flow geometries and the boundary layer Reynolds number. It is found that the pressure rise necessary to separate the boundary layer in supersonic external flows could be applied to such internal flows as overexpanded nozzles or diffusers. This is due to the fact that the separation phenomenon caused by shock wave/boundary layer interactions is processed through a supersonic deceleration. The shock-induced separation in almost all of interacting flow fields is governed by the concept of free interaction, and the onset of shock-induced separation is only a function of the Mach number just upstream of shock wave. However, physical scales of the produced separation are not independent of the downstream flow fields.     

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     

Collective Bow Shock Ahead of a Transverse System of Spheres in a Supersonic Flow Behind a Moving Shock Wave

The results of investigating the dynamics and physical conditions of formation of a collective bow shock ahead of a system of spheres with the line of centers normal to the supersonic flow behind a traveling shock wave are presented. Two types of shock-wave patterns that necessarily precede the formation of the collective shock wave and correspond to regular and Mach interaction of the bow shocks were detected experimentally. On the basis of a local gasdynamic-discontinuity interference theory, quantitative criteria of the existence of these regimes and of the formation of a common shock wave are determined. These criteria are confirmed in a series of experiments for the transitional regimes.     

Shock wave diffraction in the presence of a supersonic co-flow jet

The interaction between a diffracting shock wave and a uniform jet is a case that so far has only been partially investigated. This interaction is extremely important for the control of noise generation and improvement of combustor performance. To fill this knowledge gap, three geometries of the diffracting corner, namely a straight ramp, a serrated ramp, and a rounded corner, have been tested experimentally to study the interaction of shock diffraction with a supersonic co-flow jet at incident Mach numbers of 1.31 and 1.59, with Reynolds numbers of \(1.08\times 10^{6}\) and \(1.68\times 10^{6}\), respectively. Schlieren photography was employed to analyse the evolution of the flow phenomena. The aim is to provide a qualitative understanding of the interaction between the diffracting shock wave and the uniform jet relevant to future high-speed transport. The results show that the flow field evolves more rapidly and develops stronger structures for a higher shock Mach number. The diffraction around a rounded splitter develops a periodical vortical structure which continues after the disturbance introduced by the passage of the shock wave is removed.     

Experimental and numerical investigation of inclined air/SF6 interface instability under shock wave

The shock tube experiments of inclined air/SF6 interface instability under the shock wave with the Mach numbers 1.23 and 1.41 are conducted. The numerical simulation is done with the parallel algorithm and the multi-viscous-fluid and turbulence (MVFT) code of the large-eddy simulation (LES). The developing process of the interface accelerated by the shock wave is reproduced by the simulations. The complex wave structures, e.g., the propagation, refraction, and reflection of the shock wave, are clearly revealed in the flows. The simulated evolving images of the interface are consistent with the experimental ones. The simulated width of the turbulent mixing zone (TMZ) and the displacements of the bubble and the spike also agree well with the experimental data. Also, the reliability and effectiveness of the MVFT in simulating the problem of interface instability are validated. The more energies are injected into the TMZ when the shock wave has a larger Mach number. Therefore, the perturbed interface develops faster.     

圆球诱发斜爆轰波的数值研究

斜爆轰发动机是飞行器在高马赫数飞行条件下的一种新型发动机,具有结构简单、成本低和比冲高等优点.但是斜爆轰发动机的来流马赫数范围广,来流条件复杂,为实现斜爆轰波的迅速、可靠引发,采用钝头体来诱发.利用Euler方程和氢氧基元反应模型,对超声速氢气/空气混合气体中圆球诱导的斜爆轰流场进行了数值研究.不同于楔面诱发的斜爆轰波,球体首先会在驻点附近诱发正激波/爆轰波,然后在稀疏波作用下发展为斜激波/爆轰波.模拟结果显示,经过钝头体压缩的预混气体达到自燃温度后,会出现两种流场:当马赫数较低时,由于稀疏波的影响,燃烧熄灭,钝头体下游不会出现燃烧情况;而当马赫数较高时,燃烧阵面能传到下游.分析表明,当钝头体的尺度较小时,驻点附近的能量不足以诱发爆轰波,只会形成明显的燃烧带与激波非耦合结构;当钝头体的尺度较大时,流场中不会出现燃烧带与激波的非耦合现象,且这一特征与马赫数无关.通过调整球体直径,获得了激波和燃烧带部分耦合的燃烧流场结构,这一流场结构在楔面诱发的斜爆轰波中并不存在,说明稀疏波与爆轰波面的相互作用是决定圆球诱发斜爆轰波的关键.     

Shock wave standoff distance for a sphere slightly above Mach one

This paper describes a numerical simulation of bow shock formation ahead of a sphere at steady supersonic flow in the Mach number range of 1.025–1.20. Turbulent viscous flow results are presented using the Spalart–Allmaras turbulence model. The purpose of this study is to determine the shock standoff distance for a spherical projectile at slightly supersonic free flight speeds. Results are compared to experimental data, including double exposure holographic interferograms obtained from a 40 mm polycarbonate sphere launched by a light gas gun. The shock standoff distance was determined from the interferograms. The present numerical simulations were found to agree with previously published data, and reached down to M = 1.025—a range where almost no previously published data exists. The computed flow structure and shock wave locations agree well with recently obtained free-flight interferograms.     

Mach number effects on shock-bubble interaction

An investigation of Mach number effects on the interaction of a shock wave with a cylindrical bubble, is presented. We have conducted simulations with the Euler equations for various incident shock Mach numbers () in the range of , using high-resolution Godunov-type methods and an implicit solver. Our results are found in a very good agreement with previous investigations and further reveal additional gasdynamic features with increasing the Mach number. At higher Mach numbers larger deformations of the bubble occur and a secondary-reflected shock wave arises upstream of the bubble. Negative vorticity forms at all Mach numbers, but the “c-shaped” vortical structure appeared at gives its place to a circular-shaped structure at higher Mach numbers. The computations reveal that the (instantaneous) displacements of the upstream, downstream and jet interfaces are not significantly affected by the incident Mach number for values (approximately) greater than . With increasing the incident Mach number, the speed of the jet (arising from the centre of the bubble during the interaction) also increases. Received 21 December 2000 / Accepted 23 April 2001     

Numerical investigation of inviscid shock wave dynamics in an expansion tube

Highly complicated shock wave dynamics has been numerically calculated by solving the Euler equations for a circular shock tube suddenly expanded three times of the original tube diameter atx=0. Shock waves of different shock Mach number,M _s =1.5 and 2.0, have produced remarkably distinct blast jet structures. A planar shock wave took its final form after the blast by repeated Mach reflections of the blast wave: the first one at the wall and the second one at the central axis. The central Mach disc overtook and merged with the annular Mach stem before the planar shock wave was formed. In contrast to the blast wave which would propagate spherically in an open space, the present blast wave undergoes complex morphological transformation in the restricted flow passage, resulting in an unstable and oscillatory blast jet structure of highly rotational nature. The slipstream tube emanating from the shock tube exit corner decomposed into a chain of small vortex rings that interacted with the barrel shock of the jet, which caused periodic collapse of the jet structure. The finite volume-FCT formulation equipped with the time-dependenth-refinement adaptive unstructured triangular mesh technique in the present paper has contributed to resolution of the intricate physical discontinuities developing in the blast flow fields.     

Analogy between soap film and gas dynamics. II. Experiments on one-dimensional motion of shock waves in soap films

This paper presents an experimental investigation of one-dimensional moving shock waves in vertical soap films. The shock waves were generated by bursting the films with a perforating spark. Images of propagating shock waves and small disturbances were recorded using a fast line scan CCD camera. An aureole and a shock wave preceding the rim of the expanding hole were clearly observed. These images are similar to the x-t diagrams in gas dynamics and give the velocities of shock and sound waves. The moving shock waves cause jumps in thickness. The variations of the induced Mach number, M₂ and the ratio of film thickness across the shock wave, δ ₂/δ ₁, are plotted versus the shock Mach number, M _s. Both results suggest that soap films are analogous to compressible gases with a specific heat ratio of γ≅1.0. Published online: 15 October 2002     

Effects of Compressibility and Shock-Wave Interactions on Turbulent Shear Flows

Compressibility effects are present in many practical turbulent flows, ranging from shock-wave/boundary-layer interactions on the wings of aircraft operating in the transonic flight regime to supersonic and hypersonic engine intake flows. Besides shock wave interactions, compressible flows have additional dilatational effects and, due to the finite sound speed, pressure fluctuations are localized and modified relative to incompressible turbulent flows. Such changes can be highly significant, for example the growth rates of mixing layers and turbulent spots are reduced by factors of more than three at high Mach number. The present contribution contains a combination of review and original material. We first review some of the basic effects of compressibility on canonical turbulent flows and attempt to rationalise the differing effects of Mach number in different flows using a flow instability concept. We then turn our attention to shock-wave/boundary-layer interactions, reviewing recent progress for cases where strong interactions lead to separated flow zones and where a simplified spanwise-homogeneous problem is amenable to numerical simulation. This has led to improved understanding, in particular of the origin of low-frequency behaviour of the shock wave and shown how this is coupled to the separation bubble. Finally, we consider a class of problems including side walls that is becoming amenable to simulation. Direct effects of shock waves, due to their penetration into the outer part of the boundary layer, are observed, as well as indirect effects due to the high convective Mach number of the shock-induced separation zone. It is noted in particular how shock-induced turning of the detached shear layer results in strong localized damping of turbulence kinetic energy.     

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 analysis on flutter of Busemann-type supersonic biplane airfoil

The Busemann-type supersonic biplane can effectively reduce the wave drag through shock interference effect between airfoils. However, considering the elastic property of the wing structure, the vibration of the wings can cause the shock oscillation between the biplane, which may result in relative aeroelastic problems of the wing. In this research, fluid–structure interaction characteristics of the Busemann-type supersonic biplane at its design condition have been studied. A theoretical two-dimensional structure model has been established to consider the main elastic characteristics of the wing structure. Coupled with unsteady Navier–Stokes equations, the fluid–structure dynamic system of the supersonic biplane is studied through the two-way computational fluid dynamics/computational structural dynamics (CFD/CSD) coupling method. The biplane system has been simulated at its design Mach number with different nondimensional velocities. Different initial disturbance has been applied to excite the system and the effects of the position of the mass center on the system’s aeroelastic stability is also discussed. The results reveal that the stability of the airfoil in supersonic biplane system is decreased compared with that of the airfoil isolated in supersonic flow and such stability reduction effect should be given due attention in practical design.     

Experimental investigation of the interaction of a shock wave with a magnetic field

We present certain results of an experimental investigation of the propagation of a shock wave S through a magnetohydrodynamic channel of the Faraday type. Under conditions of short circuiting of the induced currents in the supersonic stream that follows the shock wave, we registered the occurrence of a shock front T. The x-t diagrams of the motion of the S and T shock waves in the channel and behind it are determined. For a number of fixed sections in the channel, we have measured the density and degree of ionization of the gas and determined their time dependence. The investigations were performed in argon, and the ionizing shock wave propagated with Mach numbers 12–13. The magnetic field intensity was 1.5 T.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 153–183, May–June, 1971.     

等压比热在基于壁面放电的激波控制中的影响

放电等离子体对流动的控制机理可按热效应和非热效应分为两大类,其中放电等离子体的热效应对流场中激波结构有着明显的控制作用. 目前在放电等离子体热激励对激波控制的数值模拟过程中,通常采用等效热源的方式来实现放电的热效应,数值模拟和实验的结果显示放电产生的局部温度可达到上万度. 如果数值模拟的过程中没有考虑到气体等压比热随温度的非线性变化,计算得到的结果是有失真实性的. 本文以5 马赫的超音速进气道为平台,对基于壁面放电的激波控制过程进行了数值模拟. 选取了随温度非线性变化的等压比热,并且将其结果与定等压比热的计算结果进行了对比. 结果发现:(1) 两种等压比热下,计算结果显示放电热激励在激波控制上都有着显著的效果;(2) 两种计算结果在模拟与温度相关的参数(温度、马赫数和总压恢复系数)上的差别非常明显. 因此,为了获得壁面放电对激波控制更真实的计算结果,必须考虑到等压比热随温度发生非线性变化效应的影响.