Numerical investigations of transition between regular and Mach reflections caused by free-stream disturbances

Numerical simulations have been performed to study the influence of the free-stream disturbances on the alternation of the steady shock wave reflection configurations in the dual solution domain. Different types of disturbances have been considered. The analysis of interaction between disturbances and the incident shock wave can be substantially simplified for the localized density disturbances. It is shown that such disturbances can indeed cause the transition from regular reflection to Mach reflection and back, so that within a certain range of angles of incidence the shock wave reflection configuration can be considered as a bi-stable system. The threshold amplitude of the localized density disturbance, able to induce the transition, has been estimated theoretically. The results of numerical computations convince of higher stability of the Mach reflection in the dual solution domain compared to the regular reflection, which is in accordance with available experimental data. Received 10 May 2001 / Accepted 15 November 2001 Published online 8 July 2002     

Stability and analogy of shock wave reflection in steady flow

Numerical simulations have been performed to investigate the stability of shock wave reflection in supersonic steady flow. Wall deflection control has been applied just downstream of the reflection point in the regular reflection configuration. The results provide the magnitude of the disturbance required to cause transition from one configuration to the other throughout the range of incident shock angle. An argument focusing on the subsonic region generated behind the Mach stem in the Mach reflection configuration explains the mechanism of the transition. Numerical results show that both regular and Mach reflections are possible in the dual-solution domain, and also indicate the presence of the hysteresis effect. The transition processes and the stability of the possible states are shown to be described consistently by an analogy based on the potential energy of a particle on a surface. The necessity of more sophisticated experimental investigations is emphasized to verify the argument about the stability of shock reflections and proposed analogy. Received 17 March 1997 / Accepted 26 February 1998     

Effects of a Single-pulse Energy Deposition on Steady Shock Wave Reflection

The effects of energy deposition in the free stream on steady regular and Mach shock wave reflections are studied numerically. A short-duration laser pulse is focused upstream of the incident shock waves. It causes formation of the expanding blast wave and the residual hot-spot interacting in a complex way with the steady shock wave reflection. It was found that the laser energy addition in the free stream may force the transition from regular to Mach reflection in the dual solution domain. In contrast to previously reported numerical results, the transition from Mach to regular reflection has not been reproduced in our refined computations since the Mach reflection is restored after the flow perturbation.     

Numerical confirmation of the hysteresis phenomenon in the regular to the Mach reflection transition in steady flows

Numerical calculations based on the Navier-Stokes equations are carried out to investigate the reflection of shock waves over straight reflecting surfaces in steady flows. The results for a flow Mach number of M₀=4.96 confirm the recent experimental findings of Chpoun et al. (1995) concerning the transition from regular to Mach reflection. Numerical calculations as well as experimental results show a hysteresis phenomenon during this transition and the regular reflection is found to be stable in the dual-solution domain in which theoretically both regular and Mach reflection wave configurations are possible.     

Least energy as a criterion for transition between regular and Mach reflection

A new criterion is suggested to define the point of transition between regular and Mach reflection. The suggested criterion is based on the natural tendency of a physical system to minimize its energy. The increases of the specific energy behind the reflected shock of a regular reflection and behind the Mach stem of a Mach reflection are calculated. It is hypothesized that the type of reflection that will occur is that which produces the smaller change of specific energy. The transition angles predicted using this criterion show better agreement with experimental results than those predicted using the detachment criterion for incident shock waves with Mach numbers between 1.1 and 2.0.This article was processed using Springer-Verlag TEX Shock Waves macro package     

Steady shock wave reflections in thermochemical nonequilibrium flows

The shock wave reflection phenomena in hypersonic steady air flows, including thermochemical nonequilibrium effects, are investigated. The main objectives are to study the influence of these effects on the two shock wave reflections (regular and Mach reflections), on the Mach stem height and on the hysteresis behavior. The air computations are performed using a multi-block MUSCL-TVD finite-volume scheme. The computational results with and without thermochemical effects in the air mixture flow at an upstream Mach number equal to 7 are compared. The comparison reveals a strong dependence of the transition angles, of the height and location of the Mach stem on the physical modeling of the gas flow. Received 17 February 2000 / Accepted 30 August 2000     

带隔板装药爆轰波马赫反射理论研究和数值模拟

基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。     

Investigation of the hysteresis phenomena in steady shock reflection using kinetic and continuum methods

The problem of transition of planar shock waves over straight wedges in steady flows from regular to Mach reflection and back was numerically studied by the DSMC method for solving the Boltzmann equation and finite difference method with FCT algorithm for solving the Euler equations. It is shown that the transition from regular to Mach reflection takes place in accordance with detachment criterion while the opposite transition occurs at smaller angles. The hysteresis effect was observed at increasing and decreasing shock wave angle. Received September 1, 1995 / Accepted November 20, 1995     

Analysis of double-Mach-reflection wave configurations with convexly curved Mach stems

The wave configuration of a double Mach reflection (DMR) with a convexly curved Mach stem and the resulted flow fields are analyzed. An analytical model describing the formation of a with a curved Mach stem and predicting its wave configuration is proposed. The transition criterion from a with a straight Mach stem to a with a curved Mach stem is also suggested. Predictions based on the analytical model are compared to experimental results. The agreement is found to be good to excellent. Received 5 July 1996 / Accepted 18 March 1998     

Diffuse holographic interferometric observation of shock wave reflection from a skewed wedge

The pattern of shock wave reflection over a wedge is, in general, either a regular reflection or a Mach reflection, depending on wedge angles, shock wave Mach numbers, and specific heat ratios of gases. However, regular and Mach reflections can coexist, in particular, over a three-dimensional wedge surface, whose inclination angles locally vary normal to the direction of shock propagation. This paper reports a result of diffuse double exposure holographic interferometric observations of shock wave reflections over a skewed wedge surface placed in a 100 × 180 mm shock tube. The wedge consists of a straight generating line whose local inclination angle varies continuously from 30° to 60°. Painting its surface with fluorescent spray paint and irradiating its surface with a collimated object beam at a time interval of a few microseconds, we succeeded in visualizing three-dimensional shock reflection over the skewed wedge surface. Experiments were performed at shock Mach numbers, 1.55, 2.02, and 2.53 in air. From reconstructed holographic images, we estimated critical transition angles at these shock wave Mach numbers and found that these were very close to those over straight wedges. This is attributable to the flow three-dimensionality.      

Three-dimensional effects on regular reflection in steady supersonic flows

The reflection of shock waves between two symmetrical wedges is investigated for the case of three-dimensional flows. Oblique shadowgraphs at various optical angles of yaw and pitch were used to examine the nature of fully three-dimensional flows, with wedge aspect ratios as low as 0.25 being considered. These images were used to construct surface models of the overall flow field for various reflection patterns and aspect ratios, which provides a visual indication of the flow field shape. For a Mach number of 3.1, and suitable wedge angles, the flow field with regular reflection on the tunnel centreline and Mach reflection further out is examined. The point of transition from regular reflection to the peripheral Mach surfaces is identified for various wedge angles and aspect ratios. It is shown that the transition points move outwards from the central plane as the aspect ratio decreases. This shows that three-dimensional flows favor regular reflection, because of the increasing curvature of the incident shock as the wedge becomes narrower, causing a decrease in the local angle of incidence. The height of the Mach stem is shown to be highly dependent on the geometry of the test wedge models. The Mach stem height decreases with aspect ratio due to the three-dimensional relieving effect, where the increase in lateral flow relieves the pressure over the surfaces of the wedges. Experimental evidence of the existence of the strong oblique shock solution in steady flows is presented.Received: 7 July 2003, Revised: 20 October 2003, Accepted: 6 November 2003, Published online: 10 February 2004PACS: 47.40.Nm Correspondence to: B.W. Skews     

Numerical simulation of Mach reflections

In this paper, the “FLIC” difference method with triangular mesh is adopted to numerically simulate the regular and Mach reflections that occur when a shock wave pass around a wedge. The compuational result is compared with the shock tube experimental results of G. Ben-Dor and I. I. Glass. The comparison shows that the position, shape of shock wave and height of Mach stem all show a good agreement. Consequently, the “FLIC” difference method with triangular mesh is quite satisfactory in numerical simulation of the regular and Mach reflections.     

球面激波在固壁的马赫反射

本文对金属丝电爆炸产生球面激波在固壁上的马赫反射作了研究。由高速纹影摄影系统将反射图像拍摄下来,从而测量出马赫杆在壁面的马赫数,进而得到超压以及三波点迹线,实验结果与我们的数值计算结果作了比较,二者比较吻合。     

Simulation of the Mach reflection in supersonic flows by the CE/SE method

This study employs the Space-Time Conservation Element and Solution Element (CE/SE) method to determine the influence of downstream flow conditions on Mach stem height. The results indicate that the Mach stem height depends on the incident shock wave angle and the distance between the trailing edge and the symmetry plane. Furthermore, it is shown that the downstream length ratio and the trailing edge angle do not affect the Mach stem height nor the Mach reflection (MR) configuration, and the Space-Time Conservation Element and Solution Element method is able to simulate the MR as well as many other numerical schemes. Communicated by K. Takayama PACS 47.40.Nm     

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.     

Analytical and numerical investigations of the reflection of asymmetric nonstationary shock waves

The reflection of asymmetric nonstationary shock waves is analytically and numerically studied in this paper. An analytical approach, which is a combination of the shock dynamic and shock polar methods, is advanced to predict the reflection wave configurations. The numerical simulations are performed by the finite volume method based on the second-order MUSCL-Hancock scheme and the HLLC approximate Riemann solver, with the self-adaptive unstructured mesh. It is found that the transition between the overall regular reflection and overall Mach reflection in the asymmetric nonstationary reflection agrees with the detachment criterion, which is analogous to the reflection in pseudo-steady flows (i.e. shock reflection over a wedge). Some special reflection wave configurations, which have never been observed in steady or nonstationary shock reflections so far, are found to exist in this asymmetric reflection. Furthermore, the domains and boundaries of various overall reflection wave configurations are analytically predicted, and the effect of mis-synchronization is also discussed.     

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

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

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.

---

     

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 study of shock wave interaction with a rotating cylinder

Shock wave reflection over a rotating circular cylinder is numerically and experimentally investigated. It is shown that the transition from the regular reflection to the Mach reflection is promoted on the cylinder surface which rotates in the same direction of the incident shock motion, whereas it is retarded on the surface that rotates to the reverse direction. Numerical calculations solving the Navier-Stokes equations using extremely fine grids also reveal that the reflected shock transition from RR MR is either advanced or retarded depending on whether or not the surface motion favors the incident shock wave. The interpretation of viscous effects on the reflected shock transition is given by the dimensional analysis and from the viewpoint of signal propagation.Received: 24 April 2002, Accepted: 16 August 2002, Published online: 25 March 2003