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

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

Structure of an oblique detonation wave induced by a wedge

The structure of an oblique detonation wave (ODW) induced by a wedge is investigated via numerical simulations and Rankine–Hugoniot analysis. The two-dimensional Euler equations coupled with a two-step chemical reaction model are solved. In the numerical results, four configurations of the Chapman–Jouguet (CJ) ODW reflection (overall Mach reflection, Mach reflection, regular reflection, and non-reflection) are observed to take place sequentially as the inflow Mach number increases. According to the numerical and analytical results, the change of the CJ ODW reflection configuration results from the interaction among the ODW, the CJ ODW, and the centered expansion wave.     

Theoretical considerations on shock reflections and their implications on the evaluation of air intake performance

Many theoretical studies have shown the existence of a hysteresis effect in the solution of oblique shock reflections. In fact, a wide domain of free-stream Mach number and shock angle values exists where regular reflection and Mach reflection are both possible solutions for the same flow conditions. Part of this domain overlaps the typical operating conditions of supersonic air intakes, and therefore it is of practical interest to obtain a deeper understanding of the theoretical problem. Indeed, although both solutions are theoretically possible, they yield very different flowfields and consequently large discrepancies in the evaluation of the air intake performance. Numerical solutions for steady configurations have been carried out and compared with the flow evolution obtained for time-dependent cases. The results have confirmed numerically the existence of the multiple solution domain where hysteresis takes place in time-dependent simulations. The analysis of the physical and numerical problems encountered has provided indications for a correct simulation in practical applications. Received 10 August 1999 / Accepted 6 October 2000     

不同状态方程对固体介质中斜激波反射的影响

相比气体,固体介质在高压下的状态方程更为复杂,形式也多种多样.现有关于固体介质中激波反射的理论研究,一般直接采用某种状态方程,缺乏对采用不同状态方程得到的结果的对比.本项工作采用激波极曲线的理论分析方法,选择4种不同组合形式的状态方程(一次冲击激波采用线性的冲击波速度与粒子速度关系式,二次冲击激波采用Gr(u|¨)neisen状态方程;一次冲击和二次冲击激波均采用冲击波速度与粒子速度关系式:一次冲击激波采用线性冲击波速度与粒子速度关系式,二次冲击激波采用刚性气体状态方程;以及一次冲击激波和二次冲击激波均采用刚性气体状态方程),研究固体介质中的斜激波反射,比较了采用不同组合形式的状态方程对反射激波波后压力的影响.利用量纲分析方法讨论了简化状态方程达到较高精度的条件.此外,用ANSYS/LS-DYNA软件,对激波极曲线理论给出的结果进行了验证.本项工作可为固体介质中激波反射问题状态方程的选取提供一定的指导.     

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.     

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.

---

     

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.      

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     

An experiment on imploding conical shock waves

Over the past decade and a half there have been a number of numerical studies of the reflection of oblique axisymmetric shock waves from the axis of symmetry. Many of these have shown a complex Mach reflection pattern together with a strong toroidal vortex which significantly distorts the Mach disk. This geometry has never been captured experimentally. This note describes a special rig for the generation of strong imploding conical shocks, and presents photographs of the predicted reflection pattern. Received 4 October 2001 / Accepted 13 November 2001     

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.     

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.     

STUDY ON THE FUEL AIR MIXING INDUCED BY A SHOCK WAVE PROPAGATING INTO A H_2-AIR INTERFACE

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｈｅａｒｌａｙｅｒｆｌｏｗｉｓａｎｉｍｐｏｒｔａｎｔｐｒｏｂｌｅｍｉｎｓｔｕｄｙｉｎｇｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗａｎｄｓｕｐｅｒｓｏｎｉｃｃｏｍｂｕｓｔｉｏｎ ,ｓｕｃｈａｓｆｕｅｌｄｉｆｆｕｓｉｏｎａｎｄｍｉｘｉｎｇｉｎａＳｃｒａｍｊｅｔ[1].Ｒｏｓｈｋｏ[2 ]ｓｔｕｄｉｅｄｅｘｐｅｒｉｍｅｎｔａｌｌｙｔｈｅｐｈｅｎｏｍｅｎｏｎｏｆｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｕｂｓｏｎｉｃｓｈｅａｒｆｌｏｗａｎｄｆｏｕｎｄｌａｒｇｅ ,ｃｏｈ…     

Estimation of shock induced vorticity on irregular gaseous interfaces: a wavelet-based approach

We study the interaction of a shock with a density-stratified gaseous interface (Richtmyer–Meshkov instability) with localized jagged and irregular perturbations, with the aim of developing an analytical model of the vorticity deposition on the interface immediately after the passage of the shock. The jagged perturbations, meant to simulate machining errors on the surface of a laser fusion target, are characterized using Haar wavelets. Numerical solutions of the Euler equations show that the vortex sheet deposited on the jagged interface rolls into multiple mushroom-shaped dipolar structures which begin to merge before the interface evolves into a bubble-spike structure. The peaks in the distribution of x-integrated vorticity (vorticity integrated in the direction of the shock motion) decay in time as their bases widen, corresponding to the growth and merger of the mushrooms. However, these peaks were not seen to move significantly along the interface at early times i.e. t < 10 τ, where τ is the interface traversal time of the shock. We tested our analytical model against inviscid simulations for two test cases – a Mach 1.5 shock interacting with an interface with a density ratio of 3 and a Mach 10 shock interacting with a density ratio of 10. We find that this model captures the early time (t/τ ∼ 1) vorticity deposition (as characterized by the first and second moments of vorticity distributions) to within 5% of the numerical results. PACS 47.40.Nm; 47.20.Ma     

Linear relationship between Hugoniot temperature and shock Mach number

It is the purpose of this publication to discuss further the apparent validity of a linear relationship between the Hugoniot temperature and the shock Mach number, when used as an independent variable in the thermodynamics of very high pressures. Additional evidence for seventeen different materials is presented. Some of the materials discussed might present phase transitions within the ranges of pressure and temperature here studied. The case of molybdenum is discussed in particular because experimental data on phase transitions are available within the ranges of pressure and temperature considered. Equation of state results for a few materials, obtained using an exact analytical equation of state, are compared with those computed employing an approximate form of the equation, consequence of the linear relationship between the Hugoniot temperature and shock Mach number. The excellent agreement shows that this approximate and very simple equation of state can be very reliable and useful. Received 17 June 1997 / Accepted 4 November 1997     

Pseudo-steady oblique shock wave reflections over water wedges

A preliminary experimental investigation of the reflection phenomenon of a planar shock wave over a water wedge is reported. The experimental investigation was conducted using a specially designed shock tube which was capable of being tilted from a horizontal situation to a vertical one in an accuracy of 0.1°. The four well-known types of pseudo-steady oblique shock wave reflections were observed. In addition, it has been found that the actual wedge angles at which the regular reflection transitions to a Mach reflection differ from those measured over a solid wedge.     

铅飞层中斜冲击波对碰马赫反射行为实验研究

采用线阵多普勒光纤探针测速技术(Doppler pins system,DPS)和高速光电分幅相机照相两种精密诊断技术,对铅飞层中斜冲击波对碰后的反射行为进行了观测。获得了飞层对碰部位速度-时间历史曲线和凸起形貌演化图像,给出了凸起轮廓发展演化过程、压力分布等实验数据和信息。结合冲击波反射理论,对铅飞层对碰区动力学现象进行了分析和解释,证实铅飞层中斜冲击波对碰后发生了马赫反射。     

平面激波加载下砂墙结构的冲击响应特性

砂墙结构在爆炸安全防护领域具有广泛应用,为了研究激波加载下砂墙结构的冲击响应特性,基于水平激波管实验装置,开展平面激波冲击砂墙结构系列实验,采用高速纹影摄像系统捕捉流场中激波波系的演化过程和砂墙结构的运动过程。入射激波马赫数为1.827~2.413,相应入射激波载荷强度为0.378~0.724 MPa。砂墙结构利用铁砂、矾土、石英砂3种实验用砂制备,所制备砂墙结构孔隙度分别为56.6%、69.3%、56.6%。高速纹影照片显示:平面激波冲击砂墙结构发生反射和透射,伴随入射激波和透射激波的传播,在百微秒内,砂墙未产生显著运动,表现出显著的类固体动力学响应特性。基于冲击理论,确定了铁砂墙、矾土砂墙、石英砂墙的线性冲击关系,冲击关系中线性常数λ值量级为100,根据凝聚介质实用状态方程推断:较低强度载荷冲击作用下,砂墙主要产生体积变形,而由冲击引起的热能效应则可以忽略。     

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     

Numerical simulation for nonstationary Mach reflection of a shock wave: A kinetic-model approach

A numerical simulation was performed for the process of formation of single Mach reflection on a wedge by solving a BGK type kinetic equation for the reduced distribution function with a finite difference scheme. The calculations were carried out for a shock Mach number 2.75 and wedge angle 25° in a monatomic gas, which corresponds to the conditions of single Mach reflection in the classical von Neumann theory. The calculations were performed for both diffuse and specular reflection of molecules at the wall surface. It is concluded that the diffuse reflection of molecules at the wall surface or the existence of the viscous or thermal layer is an essential factor for a nonstationary process at the initial stage of Mach reflection. Furthermore, the numerical results for diffuse reflection are found to simulate the experimental results very well, such as a transient process from regular reflection to Mach reflection along with shock propagation.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

HNS-Ⅳ炸药的点火增长模型

进行直径?10mm 的圆筒实验,根据冲击Hugoniot关系,拟合了HNS-Ⅳ炸药的未反应时的JWL 状态方程,并通过数值模拟确定了HNS-Ⅳ炸药爆轰产物JWL状态方程参数;测量炸药/窗口界面粒子速度 历程,结合数值模拟获得HNS-Ⅳ炸药的点火增长模型反应速率方程参数。研究表明,HNS-Ⅳ炸药点火增长 模型能够描述反应过程并与实验结果吻合较好。