Study on the Fuel Air Mixing Induced by a Shock Wave Propagating into a H2-Air Interface

2nd-order upwind TVD scheme was used to solve the laminar, fully Navier-Stokes equations. The numerical simulations were done on the propagation of a shock wave with Ma _S = 2 and 4 into a hydrogen and air mixture in a duct and a duct with a rearward step. The results indicate that a swirling vortex may be generated in the lopsided interface behind the moving shock. Meanwhile, the complex shock system is also formed in this shear flow region. A large swirling vortex is produced and the fuel mixing can be enhanced by a shock wave at low Mach number. But in a duct with a rearward step, the shock almost disappears in hydrogen for Ma _S = 2. The shock in hydrogen will become strong if Ma _S is large. Similar to the condition of a shock moving in a duct full of hydrogen and air, a large vortex can be formed in the shear flow region. The large swirling vortex even gets through the reflected shock and impacts on the lower wall. Then, the distribution of hydrogen behind the rearward step is divided into two regions. The transition from regular reflection to Mach reflection was observed as well in case Ma _S = 4.     

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     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象．数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较．结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征．提出了判断流向涡与运动激波相互作用中旋涡破碎的准则     

Shock wave deformation in shock-vortex interactions

An initially planar shock wave can undergo significant distortion to its shape along with changes in its strength during the period of its interaction with a compressible vortex. This phenomenon is studied by numerically simulating the shock wave-vortex interaction with a high resolution shock-capturing scheme. Incident shock waves of various Mach numbers are made to interact with a compressible vortex and the dependence of the shock wave distortion on the strength of the incident shock wave is studied in detail. It is known that the type of complex shock structure formed in the later stages of a compressible vortex-shock wave interaction is dependent on the Mach number of the incident shock wave. A simple physical model based on the principle of shock wave reflection is proposed to explain this complex shock structure formation and its dependence on the relative strengths of the interacting vortex and shock wave. Received July 28, 1997 / Accepted November 17, 1997     

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     

爆轰波在突扩通道中传播的数值模拟研究

建立了描述甲烷 空气混合物爆轰波传播的单步化学反应爆轰模型 ,通过数值模拟研究了在二维突扩通道中爆轰波的强度变化和各种波行为。结果表明 :爆轰波在进入突扩通道初始阶段的衍射使爆轰波局部向爆燃转变 ;爆炸波在壁面发生马赫反射形成的高温高压区域将直接诱导自持爆轰波的重新形成。     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

Shear layer behavior resulting from shock wave diffraction

All previous studies on shock wave diffraction in shock tubes have spatial and temporal limitations due to the size of the test sections. These limitations result from either the reflection of the expansion wave, generated at the corner, from the top wall and/or of the reflection of the incident diffracted shock from the bottom wall of the test section passing back through the region of interest. This has limited the study of the evolution of the shear layer and its associated vortex, which forms a relatively small region of the flow behind the shock with an extent of only a few centimeters, and yet is a region of significant interest. A special shock tube is used in the current tests which allow evolution of the flow to be examined at a scale about an order of magnitude larger than in previously published results, with shear layer lengths of up to 250 mm being achieved without interference from adjacent walls. Tests are presented for incident shock wave Mach numbers of nominally 1.3–1.5. Studies have been undertaken with wall angles of 10, 20, 30 and 90°. Significant changes are noted as the spatial and temporal scale of the experiment increases. For a given wall angle, the flow behind the incident shock is not self-similar as is usually assumed. Both shear layer instability and the development of turbulent patches become evident, neither of which have been noted in previous tests.     

Characteristics of Oscillations in Supersonic Open Cavity Flows

Characteristics of oscillations in supersonic open cavity flows are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier-Stokes/Large Eddy Simulation) method. The oscillation regimes and feedback mechanisms for the supersonic cavity flows are identified and analyzed. The calculation captures a mixed shear-layer/wake oscillation mode in the flow of Ma = 1.75, where these two modes occur alternately. The shear-layer mode and wake mode are driven by vortex convection-acoustic feedback and absolute instability, respectively. In particular, the results indicate that the feedback-acoustic-wave in the shear-layer mode is probably generated by the reflection of the downstream-traveling pressure wave, associated with the shed vortex in the shear layer, on the aft wall. The cavity flow of Ma = 2.52 is then simulated to see the influence of Mach number. It is found that the increase of Mach number may decrease the amplitude of the fluctuations in the shear layer, inhibiting the transition to wake mode. Furthermore, the influence of upstream injection is also studied, where the results show that the injection only weakens the oscillations and faintly shifts the resonant frequencies.     

A numerical study of oblique shock wave reflections over wedges in an ideal quantum gas

The various oblique shock wave reflection patterns generated by a moving incident shock on a planar wedge using an ideal quantum gas model are numerically studied using a novel high resolution quantum kinetic flux splitting scheme. With different incident shock Mach numbers and wedge angles as flow parameters, four different types of reflection patterns, namely, the regular reflection, simple Mach reflection, complex Mach reflection and the double Mach reflection as in the classical gas can be classified and observed. Both Bose–Einstein and Fermi–Dirac gases are considered.      

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程,研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体,在初始压强20kPa的条件下, 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的,第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰,起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构,同时影响了爆轰波传播的波动 力学过程.     

A note on numerical simulation of vortical structures in shock diffraction

In numerical simulation of the Euler equations, the slipstream or shear layer that appears behind a diffracted shock wave may develop small discrete vortices using fine computational meshes. Similar phenomena were also observed in the simulation of a Mach reflection that is accompanied by a shear layer. However, these small vortices have never been observed in any shock-tube experiment, although the wave pattern and the shape of the main vortex agree very well with visualization results. Numerical solutions obtained with coarse grids may agree better with experimental photos than those with very fine grids because of the pollution of the small vortices. This note tries to investigate the effect of viscosity on the small vortices by comparing the solutions of the laminar Navier-Stokes equations and the turbulence model. It is found that the small vortices are still observed in the solution of the laminar Navier-Stokes equations, although they can be suppressed by using the turbulence model. Numerical and experimental factors that are responsible for the deviation of the laminar solutions from experimental results are discussed. The secondary vortex in shock diffraction is successfully simulated by solving the Navier-Stokes equations.Received: 28 March 2003, Accepted: 6 May 2003, Published online: 11 June 2003     

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.      

Self-intensification in shock wave and vortex interaction

Computer simulation has been performed for the interaction between a shock wave and a vortex ring moving toward the wave. The computed density contours are compared with the pattern of shadowgraphs. A remarkable property found in the simulation is that, during the passage of the shock wave over the vortex ring, the part of the wave propagating through the inside of the ring-vortex is intensified spontaneously at a localized region. Maximum pressure occurs inside the vortex and attains a high value, about several times that of the impinging shock for incident Mach numbers of around 1.2 with the vortex translation Mach number 0.60. This is due to a double-step mechanism of intensification within the flow field by the shock-vortex interaction.     

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.      

可压缩流向涡与反向运动激波相互作用的实验

对可压缩流向涡与反向运动激波相互作用的现象进行了实验研究．实验在９４ｍｍ×９４ｍｍ的方截面激波管中进行．在实验段上游安装了一个有限翼展平直机翼．当入射激波通过机翼后，波后２区气流在模型翼尖诱导出一条流向涡．入射激波在激波管端壁反射后，形成的反射激波在观察窗处和流向涡发生作用．实验中拍摄了激波与流向涡作用全过程的纹影照片，观察到了一些和定常激波与旋涡相互作用不同的现象，并与数值计算结果进行了初步比较     

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.     

Imploding conical shock waves

The behaviour of conical shock waves imploding axisymmetrically was first studied numerically by Hornung (J Fluid Mech 409:1–12, 2000) and this prompted a limited experimental investigation into these complex flow patterns by Skews et al. (Shock Waves 11:323–326, 2002). Modification of the simulation boundary conditions, resulting in the loss of self-similarity, was necessary to image the flow experimentally. The current tests examine the temporal evolution of these flows utilising a converging conical gap of fixed width fed by a shock wave impinging at its entrance, supported by CFD simulations. The effects of gap thickness, angle and incident shock strength were investigated. The wave initially diffracts around the outer lip of the gap shedding a vortex which, for strong incident shock cases, can contain embedded shocks. The converging shock at exit reflects on the axis of symmetry with the reflected wave propagating outwards resulting in a triple point developing on the incident wave together with the associated shear layer. This axisymmetric shear layer rolls up into a mushroom-shaped toroidal vortex ring and forward-facing jet. For strong shocks, this deforms the Mach disk to the extent of forming a second triple point with the primary shock exhibiting a double bulge. Separate features resembling the Richtmeyer–Meshkov and Kelvin–Helmholtz instabilities were noted in some tests. Aside from the incident wave curvature, the reflection patterns demonstrated correspond well with the V- and DV-types identified by Hornung although type S was not clearly seen, possibly due to the occlusion of the reflection region by the outer diffraction vortex at these early times. Some additional computational work explicitly exploring the limits of the parameter space for such systems has demonstrated the existence of a possible further reflection type, called vN-type, which is similar to the von Neumann reflection for plane waves. It is recommended that the parameter space be more thoroughly explored experimentally.     

Three-dimensional supersonic internal flows

In order to examine the transition between regular and Mach reflection in a three-dimensional flow, a range of special geometry test pieces, and inlets, were designed. The concept is to have a geometry consisting of two plane wedges which results in regular reflection between the incident waves off the top and bottom of the inlet capped by two curved end sections causing Mach reflection. The merging of these two reflection patterns and the resulting downstream flow are studied using laser vapor screen and shadowgraph imaging supported by numerical simulation. An angled Mach disc is formed which merges with the line of regular reflection. A complex wave pattern results with the generation of a bridging shock connecting the reflected wave from the Mach reflection with the reflected waves from the regular reflection. In order to experimentally access the flow within the duct, a number of tests were conducted with one end cap removed. This resulted in a modified flow due to the expansive flow at the open end the influence of which was also studied in more detail.     

激波喷出过程和在林带中的传播

运用激波管技术研究了激波从管口喷出和在林带中传播这两个基本问题。给出了三方面的结果 :揭示了激波从管口喷出时的复杂流场 ,显示了主激波弯曲和衰减 ,二次激波的形成和合拢 ,涡环的长大和发展 ,以及涡环前缘出现正激波等基本物理现象 ;显示了激波在林带中的传播过程 ,取得了激波遇林地的反射 ,马赫杆形成 ,激波与林冠的相遇 ,林冠对涡环的阻滞效应等纹影照片 ;沿激波在林地的传播方向测量了地面压力 ,证实林带有明显的消波效应。