Transition Between Detonation and Deflagration in a Tube with a Cavity

In this paper, the transition of a detonation from deflagration was investigated numerically while a detonation wave propagates in a tube with a sudden change in cross section, referred to as the expansion cavity. The dispersion-controlled scheme was adopted to solve Euler equations of axis-symmetric flows implemented with detailed chemical reaction kinetics of hydrogen-oxygen (or hydrogen-air) mixture. The fractional step method was applied to treat the stiff problems of chemical reaction flow. It is observed that phenomena of detonation quenching and reigniting appear when the planar detonation front diffracts at the vertex of the expansion cavity entrance. Numerical results show that detonation front in mixture of higher sensitivity keeps its substantial coupled structure when it propagates into the expansion cavity. However, the leading shock wave decouples with the combustion zone if mixture of lower sensitivity was set as the initial gas.

郭永怀文集

冲击压缩下玻璃等脆性材料中失效波的实验和理论研究

本文首先回顾了失效波的研究进展和存在的问题，然后在现有实验手段改进的基础上对目前失效波研究中尚存的问题进行了较系统的实验研究，并建立了理论模型，进行了失效波传播的数值模拟。本文实验工作给出了K9玻璃中失效波速度和冲击载荷之间的定量关系和产生失效波的最小载荷阈值；通过：（1）玻璃样品和飞片表面状况对萌生失效波的载荷阈值的影响；（2）高速摄影观测冲击压缩下玻璃样品表面和内置界面处失效波的萌生和冲击波后压缩区内损伤演化情况；（3）回收样品的显微分析和Ｘ光衍射分析等实验，证实失效波的萌生基本与冲击相变无关，主要与玻璃样品表面状况有关，失效波的本质是玻璃样品表面固有微裂纹和冲击瞬间在此处萌生的微裂纹系统向玻璃样品中扩展的宏观统计表现。高速摄影观测到冲击波后有破碎界面在移动，其速度明显高于同等加载条件下VISAR测试的失效波速度，据此推测失效波是由大量裂纹扩展的宏观表现一破碎界面和其后方声阻抗明显降低的移动界面组成。实验同时研究了微晶玻璃、高纯度石英玻璃以及碱石灰玻璃在冲击压缩下的动力学响应特性。本文关于失效波方面的研究工作结果大部分未见相关文献报道。实验工作对于深入了解玻璃中失效波的萌生机制和失效波的力学性质等是非常重要的，并有助于建立与失效波萌生、传播有关的理论模型。改进有关实验技术，是完成本文实验研究的必要条件，也是本文工作的重要组成部分。本文设计了一种高接收效率、景深可调的新型VISAR探头，其技术指标接近国外同类产品，成本不到国外同类产品的二十分之一，且结构是国内外同类产品中最简单的。本文设计高速摄影的阴影和纹影光路，用于观察冲击压缩下玻璃样品中冲击损伤和失效波的演化及发展规律，其技术优于国外同类实验。本文设计了两种冲击压力低于玻璃样品HEL值的爆轰加载装置，用于配合高速摄影诊断实验。爆轰驱动厚飞片装置的设计，在飞片的炸药透镜之间增加了一个空腔，既降低了飞片的速度，又有效避免了以往同类装置驱动厚飞片时经常遇到的层裂问题，并且成功地进行了与爆轰驱动有关的数值模拟及设计工作。本文在理想微裂系统演化理论的基础上，建立了描述失效波的理论模型，分别讨论了表面损伤、微裂纹扩展和微裂纹形核、长大对失效波萌生、传播的影响，给出了描述失效波扩展的损伤演化方程，结合冲击压缩下材料的Resende压缩损伤本构方程，进行了各种加载条件下玻璃样品后自由表面速度和玻璃体内应力分存的数值模拟，计算结果与实验结果相符，表明本文建立的理论模型是可靠的，能够反映冲击压缩下失效波的传播特性的基本力学特性。用飞片碰玻璃样品时，实验测量的表面速度时程曲线出现“过冲现象”。本文通过高速摄影和VISAR等测试手段从实验上对其进行了较系统的研究，发现这种现象与玻璃样品后自由表面的破碎有关，并且其萌生所需的载荷阈值与萌生失效波的载荷阈值接近。本文建立了描述这一破坏现象的损伤演化方程，对其进行了数值模拟。计算结果和实验结果吻合，表明本文对这一现象的解释是合理的。

Development And Application Of Flexible Substrate Sensors In Instantaneous Heat Flux Measurement

A new type of sensor with the flexible substrate is introduced. It is applicable in measuring instantaneous heat flux on the model surface in a hypersonic shock tunnel. The working principle, structure and manufacture process of the sensor are presented. The substrate thickness and the dynamic response parameter of the sensor are calculated. Because this sensor was successfully used in measuring the instantaneous heat flux on the surface of a flat plate in a detonation-driven shock tunnel, it may be effective in measuring instantaneous heat flux on the model surface.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

MHD控制激波诱导边界层分离的机理

激波边界层相互作用是高超声速飞行器面临的重大问题,激波入射到平板引起的边界层分离是其中最具代表性的一种.用加权的3阶ENN格式计算了小磁雷诺数近似的MHD方程,研究了 MHD控制层流边界层分离的机理.数值结果显示,通过局部电离空气并施加洛伦兹力,能使分离点向下游移动,分离区尺寸减小,从而抑制和缓解由于激波-边界层相互作用而引起的分离.

Planar reflection of gaseous detonations

Pipes containing flammable gaseous mixtures may be subjected to internal detonation. When the detonation normally impinges on a closed end, a reflected shock wave is created to bring the flow back to rest. This study built on the work of Karnesky (2010) and examined deformation of thin-walled stainless steel tubes subjected to internal reflected gaseous detonations. A ripple pattern was observed in the tube wall for certain fill pressures, and a criterion was developed that predicted when the ripple pattern would form. A two-dimensional finite element analysis was performed using Johnson-Cook material properties; the pressure loading created by reflected gaseous detonations was accounted for with a previously developed pressure model. The residual plastic strain between experiments and computations was in good agreement.

During the examination of detonation-driven deformation, discrepancies were discovered in our understanding of reflected gaseous detonation behavior. Previous models did not accurately describe the nature of the reflected shock wave, which motivated further experiments in a detonation tube with optical access. Pressure sensors and schlieren images were used to examine reflected shock behavior, and it was determined that the discrepancies were related to the reaction zone thickness extant behind the detonation front. During these experiments reflected shock bifurcation did not appear to occur, but the unfocused visualization system made certainty impossible. This prompted construction of a focused schlieren system that investigated possible shock wave-boundary layer interaction, and heat-flux gauges analyzed the boundary layer behind the detonation front. Using these data with an analytical boundary layer solution, it was determined that the strong thermal boundary layer present behind the detonation front inhibits the development of reflected shock wave bifurcation.

Techniques for strength measurement at high pressures and strain-rates using transverse waves

The study of the strength of a material is relevant to a variety of applications including automobile collisions, armor penetration and inertial confinement fusion. Although dynamic behavior of materials at high pressures and strain-rates has been studied extensively using plate impact experiments, the results provide measurements in one direction only. Material behavior that is dependent on strength is unaccounted for. The research in this study proposes two novel configurations to mitigate this problem.

The first configuration introduced is the oblique wedge experiment, which is comprised of a driver material, an angled target of interest and a backing material used to measure in-situ velocities. Upon impact, a shock wave is generated in the driver material. As the shock encounters the angled target, it is reflected back into the driver and transmitted into the target. Due to the angle of obliquity of the incident wave, a transverse wave is generated that allows the target to be subjected to shear while being compressed by the initial longitudinal shock such that the material does not slip. Using numerical simulations, this study shows that a variety of oblique wedge configurations can be used to study the shear response of materials and this can be extended to strength measurement as well. Experiments were performed on an oblique wedge setup with a copper impactor, polymethylmethacrylate driver, aluminum 6061-t6 target, and a lithium fluoride window. Particle velocities were measured using laser interferometry and results agree well with the simulations.

The second novel configuration is the y-cut quartz sandwich design, which uses the anisotropic properties of y-cut quartz to generate a shear wave that is transmitted into a thin sample. By using an anvil material to back the thin sample, particle velocities measured at the rear surface of the backing plate can be implemented to calculate the shear stress in the material and subsequently the strength. Numerical simulations were conducted to show that this configuration has the ability to measure the strength for a variety of materials.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Metaconcrete: Engineered aggregates for enhanced dynamic performance

This work presents the development and investigation of a new type of concrete for the attenuation of waves induced by dynamic excitation. Recent progress in the field of metamaterials science has led to a range of novel composites which display unusual properties when interacting with electromagnetic, acoustic, and elastic waves. A new structural metamaterial with enhanced properties for dynamic loading applications is presented, which is named metaconcrete. In this new composite material the standard stone and gravel aggregates of regular concrete are replaced with spherical engineered inclusions. Each metaconcrete aggregate has a layered structure, consisting of a heavy core and a thin compliant outer coating. This structure allows for resonance at or near the eigenfrequencies of the inclusions, and the aggregates can be tuned so that resonant oscillations will be activated by particular frequencies of an applied dynamic loading. The activation of resonance within the aggregates causes the overall system to exhibit negative effective mass, which leads to attenuation of the applied wave motion. To investigate the behavior of metaconcrete slabs under a variety of different loading conditions a finite element slab model containing a periodic array of aggregates is utilized. The frequency dependent nature of metaconcrete is investigated by considering the transmission of wave energy through a slab, which indicates the presence of large attenuation bands near the resonant frequencies of the aggregates. Applying a blast wave loading to both an elastic slab and a slab model that incorporates the fracture characteristics of the mortar matrix reveals that a significant portion of the supplied energy can be absorbed by aggregates which are activated by the chosen blast wave profile. The transfer of energy from the mortar matrix to the metaconcrete aggregates leads to a significant reduction in the maximum longitudinal stress, greatly improving the ability of the material to resist damage induced by a propagating shock wave. The various analyses presented in this work provide the theoretical and numerical background necessary for the informed design and development of metaconcrete aggregates for dynamic loading applications, such as blast shielding, impact protection, and seismic mitigation.

Three-dimensional SBLI control for transonic airfoils

The application of shock control to transonic airfoils and wings has been demonstrated widely to have the potential to reduce wave drag. Most of the suggested control devices are two-dimensional, that is they are of uniform geometry in spanwise direction. Examples of such techniques include contour bumps and passive control. Recently it has been observed that a spanwise array of discrete three-dimensional controls can have similar benefits but also offer advantages in terms of installation complexity and drag. This paper describes research carried out in Cambridge into various three-dimensional devices, such as slots, grooves and bumps. In all cases the control device is applied to the interaction of a normal shock wave (M=1.3) with a turbulent boundary layer. Theoretical considerations are proposed to determine how such fundamental experiments can provide estimates of control performance on a transonic wing. The potential of each class of three-dimensional device for wave drag reduction on airfoils is discussed and surface bumps in particular are identified as offering potential drag savings for typical transonic wing applications under cruise conditions.

Micro-vortex generator flow control for supersonic engine inlets

To investigate the flow control potential of micro-vortex generators for supersonic mixed-compression inlets, a basic model experiment has been designed which combines a normal shock wave with a subsonic diffuser. The diffuser is formed by a simple expansion corner, with a divergence angle of 6 degrees. The diffuser entry Mach numbers were M=1.3 and M=1.5 and a number of shock locations relative to the corner position were tested. Flow control was applied in the form of counter-rotating micro-vanes with heights of approximately 20% of boundary layer thickness. Furthermore, corner fences where employed to reduce sidewall effects. It was found that micro-vortex generators were able to significantly reduce the extent of flow separation under all conditions, but could not eliminate it altogether. Corner fences also demonstrated potential for improving the flow in rectangular cross section channels and the combination of corner fences with micro-vortex generators was found to give the greatest benefits. At M=1.3 the combination of corner fences and micro-vanes placed close to the diffuser entry could prevent separation for a wide range of conditions. At the higher diffuser entry Mach number the benefits of flow control were less significant although a reduction of separation size and an improved pressure recovery was observed. It is thought that micro-vortex generators can have significant flow control potential if they are placed close to the expected separation onset and when the adverse pressure gradient is not too far above the incipient separation level. The significant beneficial effects of corner fences warrant a more comprehensive further investigation. It is thought that the control methods suggested here are capable of reducing the bleed requirement on an inlet, which could provide significant performance advantages.

Dynamic compressive response of composite square honeycombs

Carbon fibre-epoxy composite square honeycombs, and the parent composite material, were tested in quasi-static compression at a strain rate of 10 -3 s -1 and in dynamic compression at strain rates of 10 3-10 4 s -1 using an instrumented Kolsky bar arrangement. Taken together, these tests provide an assessment of the potential of this composite topology for use as a lightweight sandwich core. The honeycombs had two relative densities, 0.12 and 0.24, and two material orientations, ±45° and 0/90° with respect to the prismatic, loading direction of the honeycomb. Honeycomb manufacture was by slotting, assembling and bonding together carbon fibre/epoxy woven plies of composite sheets of 2 × 2 twill weave construction. The peak value of wall stress in the honeycombs was about one third that of the parent material, for all strain rates. An elastic finite element analysis was used to trace the source of this knock-down in strength: a stress concentration exists at the root of the slots and leads to premature failure by microbuckling. Shock-wave effects were evident at impact velocities exceeding 50 ms -1 for the honeycomb of relative density 0.12. This was traced to stubbing of the buckled cell walls against the face of the Kolsky bar. © 2011 Elsevier Ltd. All rights reserved.

A canonical normal SBLI flow relevant to external compression inlets

The normal shock wave / boundary layer interaction (normal SBLI) is important to the operation and performance of a supersonic inlet, and the normal SBLI is particularly prominent in external compression inlets. To improve our understanding of such interactions, it is helpful to make use of fundamental flows which capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental fiow-fleld configurations have been considered as possible test cases to represent the normal SBLI aspects found in typical external compression inlets, and it was found that the spillage-diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flow-fleld allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to be all held approximately constant mid independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage-diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, the width around twice or three times the height, and with the area expansion just downstream of the shock on the conservative side of the stall limit for incompressible diffusers. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Research on the mechanics of underwater supersonic gas jets

An experimental research was carried out to study the fluid mechanics of underwater supersonic gas jets. High pressure air was injected into a water tank through converging-diverging nozzles (Laval nozzles). The jets were operated at different conditions of over-, full-and under-expansions. The jet sequences were visualized using a CCD camera. It was found that the injection of supersonic air jets into water is always accompanied by strong flow oscillation, which is related to the phenomenon of shock waves feedback in the gas phase. The shock wave feedback is different from the acoustic feedback when a supersonic gas jet discharges into open air, which causes screech tone. It is a process that the shock waves enclosed in the gas pocket induce a periodic pressure with large amplitude variation in the gas jet. Consequently, the periodic pressure causes the jet oscillation including the large amplitude expansion. Detailed pressure measurements were also conducted to verify the shock wave feedback phenomenon. Three kinds of measuring methods were used, i.e., pressure probe submerged in water, pressure measurements from the side and front walls of the nozzle devices respectively. The results measured by these methods are in a good agreement. They show that every oscillation of the jets causes a sudden increase of pressure and the average frequency of the shock wave feedback is about 5-10 Hz.