Forward-Running Detonation Drivers for High-Enthalpy Shock Tunnels

To improve the quality of driving flows generated with detonation-driven shock tunnels operated in the forward-running mode, various detonation drivers with specially designed sections were examined. Four configurations of the specially designed section, three with different converging angles and one with a cavity ring, were simulated by solving the Euler equations implemented with a pseudo kinetic reaction model. From the first three cases, it is observed that the reflection of detonation fronts at the converging wall results in an upstream-traveling shock wave that can increase the flow pressure that has decreased due to expansion waves, which leads to improvement of the driving flow. The configuration with a cavity ring is found to be more promising because the upstream-traveling shock wave appears stronger and the detonation front is less overdriven. Although pressure fluctuations due to shock wave focusing and shock wave reflection are observable in these detonation-drivers, they attenuate very rapidly to an acceptable level as the detonation wave propagates downstream. Based on the numerical observations, a new detonation-driven shock tunnel with a cavity ring is designed and installed for experimental investigation. Experimental results confirm the conclusion drawn from numerical simulations. The generated driving flow in this shock tunnel could maintain uniformity for as long as 4 ms. Feasibility of the proposed detonation driver for high-enthalpy shock tunnels is well demonstrated.

Oxyhydrogen combustion and detonation driven shock tube

The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incident shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation.

Analysis of residual stress gradient in MEMS multi-layer structure

Residual stress and its gradient through the thickness are among the most important properties of as-deposited films. Recently, a new mechanism based on a revised Thomas-Fermi-Dirac (TFD) model was proposed for the origin of intrinsic stress in solid film

Abel Inversion of Axially-Symmetric Shock Wave Flows

Finite-fringe interferograms produced for axisymmetric shock wave flows are analyzed by Fourier transform fringe analysis and an Abel inversion method to produce density field data for the validation of numerical models. For the Abel inversion process, we use basis functions to model phase data from axially-symmetric shock wave structure. Steady and unsteady flow problems are studied, and compared with numerical simulations. Good agreement between theoretical and experimental results is obtained when one set of basis functions is used during the inversion process, but the shock front is smeared when another is used. This is because each function in the second set of basis functions is infinitely differentiable, making them poorly-suited to the modelling of a step function as is required in the representation of a shock wave.

Simulating shock to detonation transition: Algorithm and results

An algorithm based on flux-corrected transport and the Lagrangian finite element method is presented for solving the problem of shock dynamics. It is verified through the model problem of one-dimensional strain elastoplastic shock wave propagation that the algorithm leads to stable, non-oscillatory results. Shock initiation and detonation wave propagation is simulated using the algorithm, and some interesting results are obtained. (C) 1999 Academic Press.

Shear Strength Measurements in LY-12 Aluminium Alloy During Shock Loading

Lateral stress of LY-12 alummium alloy under plate impact shock loading was measured. Based on the measured data, the Hugoniot relation and shear strength were obtained. The result has demonstrated that the shear strenath of the tested material increases remarkably with the increasing longitudinal stress. This means that the assumption of constant shear strength usually adopted in shock stress calculation is not suitable for the present material.

Shock wave diffraction by a square cavity filled with dusty gas

A simple two-dimensional square cavity model is used to study shock attenuating effects of dust suspension in air. The GRP scheme for compressible flows was extended to simulate the fluid dynamics of dilute dust suspensions, employing the conventional two-phase approximation. A planar shock of constant intensity propagated in pure air over Aat ground and diffracted into a square cavity filled with a dusty quiescent suspension. Shock intensities were M-s = 1.30 and M-s = 2.032, dust loading ratios were alpha = 1 and alpha = 5, and particle diameters were d = 1, 10 and 50 mum. It was found that the diffraction patterns in the cavity were decisively attenuated by the dust suspension, particularly for the higher loading ratio. The particle size has a pronounced effect on the flow and wave pattern developed inside the cavity. Wall pressure historics were recorded for each of the three cavity walls, showing a clear attenuating effect of the dust suspension.

Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Flow Characteristics and Its Friction Coefficient of Laminar Near-Wall Boundary Layer in Dusty Atmosphere

采用双向耦合的双流体模型,研究了大气悬浮沙尘的存在对大气边界层中层流底层流动特性及摩阻系数的影响,计算并讨论了不同沙尘含量下含尘大气相对于无尘大气摩阻系数的变化。结果表明：摩阻系数的变化取决于悬浮沙尘的初始运动状态和质量载荷率。

Wave dynamic processes induced by a supersonic projectile discharging from a shock tube

A numerical study on wave dynamic processes occurring in muzzle blast flows, which are created by a supersonic projectile released from the open-end of a shock tube into ambient air, is described in this paper. The Euler equations, assuming axisymmetric flows, are solved by using a dispersion-controlled scheme implemented with moving boundary conditions. Three test cases are simulated for examining friction effects on the muzzle flow. From numerical simulations, the wave dynamic processes, including two blast waves, two jet flows, the bow shock wave and their interactions in the muzzle blasts, are demonstrated and discussed in detail. The study shows that the major wave dynamic processes developing in the muzzle flow remain similar when the friction varies, but some wave processes, such as shock-shock interactions, shock-jet interactions and the contact surface instability, get more intensive, which result in more complex muzzle blast flows.

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.

Influence of the thickness ratio of air layer to liquid layer on the onset of Benard-Marangoni convection

An experimental investigation of the onset of Benard-Marangoni convection has been performed in a liquid layer of rectangular configuration. The critical temperature difference was measured via the detections of both temperature field pattern (IR-imaging) on the free surface and fluid convection (PIV) in the liquid layer. The critical temperature difference or the critical Marangoni number was given. The experiments were performed for a fixed depth of air layer and a changeable depth of the liquid layer, and then the influence of the thickness ratio of the air layer to liquid layer on the Marangoni instability was studied.

Coherent Structures in Transition of a Flat-Plate Boundary Layer at Ma=0.7

Direct numerical simulation (DNS) of a spatially evolving flat-plate boundary layer transition process at free stream Mach number 0.7 is performed. Tollmien-Schlichting (T-S) waves are added on the inlet boundary as the disturbances before transition. Typical coherent structures in the transition process are investigated based on the second invariant of velocity gradient tensor. The instantaneous shear stress and the mean velocity profile in the transition region are studied. In our view, the fact that the peak value of shear stress in the stress concentration area increases and exceeds a threshold value during the later stage of the transition process plays an important role in the laminar breakdown process.

Investigation of High-Speed Combustible Gas Ignited by a Hot Gas Jet Produced in the Shock Tube

The high-speed combustible gas ignited by a hot gas jet, which is induced by shock focusing, was experimentally investigated. By use of the separation mode of shock tube, the test section of a single shock tube is split into two parts, which provide the high-speed flow of combustible gas and pilot flame of hot gas jet, respectively. In the interface of two parts of test sections the flame of jet was formed and spread to the high-speed combustible gas. Two kinds of the ignitions, 3-D “line-flame ignition” and 2-D “plane-flame ignition”, were investigated. In the condition of 3-D “lineflame ignition” of combustion, thicker hot gas jet than pure air jet, was observed in schlieren photos. In the condition of 2-D “plane-flame ignition” of combustion, the delay time of ignition and the angle of flame front in schlieren photos were measured, from which the velocity of flame propagation in the high-speed combustible gas is estimated in the range of 30–90 m/s and the delay time of ignition is estimated in the range of 0.12–0.29 ms.

Hybrid Finite Compact-WENO Schemes for Shock Calculation

Hybrid finite compact (FC)-WENO schemes are proposed for shock calculations. The two sub-schemes (finite compact difference scheme and WENO scheme) are hybridized by means of the similar treatment as in ENO schemes. The hybrid schemes have the advantages of FC and WENO schemes. One is that they possess the merit of the finite compact difference scheme, which requires only bi-diagonal matrix inversion and can apply the known high-resolution flux to obtain high-performance numerical flux function; another is that they have the high-resolution property of WENO scheme for shock capturing. The numerical results show that FC-WENO schemes have better resolution properties than both FC-ENO schemes and WENO schemes. In addition, some comparisons of FC-ENO and artificial compression method (ACM) filter scheme of Yee et al. are also given.