Colliding Blast Waves Driven by the Interaction of a Short-Pulse Laser with a Gas of Atomic Clusters

Collisions between shocks are commonly found in many astrophysical objects, however robust numerical models or laboratory analogues of these complex systems remain challenging to implement. We report on the development of scaled laboratory experiments which employ new techniques for launching and diagnosing colliding shocks and high Mach number blast waves, scalable to a limited subset of astrophysically-relevant regimes. Use of an extended medium of atomic clusters enables efficient (>80%) coupling of 700 fs, 1 J, 1054 nm laser pulses to a “cluster” gas with an average density of ≈10¹⁹ particles cm⁻³, producing an initial energy density >10⁵ J cm⁻³, equivalent to ≈5×10⁹ J/g. Multiple laser foci are used to tailor the spatial profile of energy deposition, or to launch pairs of counter-propagating cylindrical shocks which then collide. By probing the collision interferometrically at multiple view angles in 5^{^} increments and applying an inverse Radon transform to the resulting phase projections we have been able to tomographicall reconstruct the full three-dimensional, time-framed electron density profile of the system.     

Directed blasts and blast-generated pyroclastic density currents: a comparison of the Bezymianny 1956, Mount St Helens 1980, and Soufrière Hills,Montserrat 1997 eruptions and deposits

We compare eruptive dynamics, effects and deposits of the Bezymianny 1956 (BZ), Mount St Helens 1980 (MSH), and Soufrière Hills volcano, Montserrat 1997 (SHV) eruptions, the key events of which included powerful directed blasts. Each blast subsequently generated a high-energy stratified pyroclastic density current (PDC) with a high speed at onset. The blasts were triggered by rapid unloading of an extruding or intruding shallow magma body (lava dome and/or cryptodome) of andesitic or dacitic composition. The unloading was caused by sector failures of the volcanic edifices, with respective volumes for BZ, MSH, and SHV c. 0.5, 2.5, and 0.05 km³. The blasts devastated approximately elliptical areas, axial directions of which coincided with the directions of sector failures. We separate the transient directed blast phenomenon into three main parts, the burst phase, the collapse phase, and the PDC phase. In the burst phase the pressurized mixture is driven by initial kinetic energy and expands rapidly into the atmosphere, with much of the expansion having an initially lateral component. The erupted material fails to mix with sufficient air to form a buoyant column, but in the collapse phase, falls beyond the source as an inclined fountain, and thereafter generates a PDC moving parallel to the ground surface. It is possible for the burst phase to comprise an overpressured jet, which requires injection of momentum from an orifice; however some exploding sources may have different geometry and a jet is not necessarily formed. A major unresolved question is whether the preponderance of strong damage observed in the volcanic blasts should be attributed to shock waves within an overpressured jet, or alternatively to dynamic pressures and shocks within the energetic collapse and PDC phases. Internal shock structures related to unsteady flow and compressibility effects can occur in each phase. We withhold judgment about published shock models as a primary explanation for the damage sustained at MSH until modern 3D numerical modeling is accomplished, but argue that much of the damage observed in directed blasts can be reasonably interpreted to have been caused by high dynamic pressures and clast impact loading by an inclined collapsing fountain and stratified PDC. This view is reinforced by recent modeling cited for SHV. In distal and peripheral regions, solids concentration, maximum particle size, current speed, and dynamic pressure are diminished, resulting in lesser damage and enhanced influence by local topography on the PDC. Despite the different scales of the blasts (devastated areas were respectively 500, 600, and >10 km² for BZ, MSH, and SHV), and some complexity involving retrogressive slide blocks and clusters of explosions, their pyroclastic deposits demonstrate strong similarity. Juvenile material composes >50% of the deposits, implying for the blasts a dominantly magmatic mechanism although hydrothermal explosions also occurred. The character of the magma fragmented by explosions (highly viscous, phenocryst-rich, variable microlite content) determined the bimodal distributions of juvenile clast density and vesicularity. Thickness of the deposits fluctuates in proximal areas but in general decreases with distance from the crater, and laterally from the axial region. The proximal stratigraphy of the blast deposits comprises four layers named A, B, C, D from bottom to top. Layer A is represented by very poorly sorted debris with admixtures of vegetation and soil, with a strongly erosive ground contact; its appearance varies at different sites due to different ground conditions at the time of the blasts. The layer reflects intense turbulent boundary shear between the basal part of the energetic head of the PDC and the substrate. Layer B exhibits relatively well-sorted fines-depleted debris with some charred plant fragments; its deposition occurred by rapid suspension sedimentation in rapidly waning, high-concentration conditions. Layer C is mainly a poorly sorted massive layer enriched by fines with its uppermost part laminated, created by rapid sedimentation under moderate-concentration, weakly tractive conditions, with the uppermost laminated part reflecting a dilute depositional regime with grain-by-grain traction deposition. By analogy to laboratory experiments, mixing at the flow head of the PDC created a turbulent dilute wake above the body of a gravity current, with layer B deposited by the flow body and layer C by the wake. The uppermost layer D of fines and accretionary lapilli is an ash fallout deposit of the finest particles from the high-rising buoyant thermal plume derived from the sediment-depleted pyroclastic density current. The strong similarity among these eruptions and their deposits suggests that these cases represent similar source, transport and depositional phenomena.     

Sensitivity equations for hyperbolic conservation law-based flow models

The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.     

辉长岩在冲击加载下的相变及物态方程

岩石相变关系对限定地球内部物质成分和动力学关系有着重要意义.利用光探针在14~39 GPa压力区间对热压烧结的辉长岩样品进行了冲击波加载实验.从16 GPa开始,D-u线斜率发生明显变化,说明样品内发生了相变.利用原始的Hugoniot曲线结合地震波状态方程（Seismic Equation of State）得到了Birch-Murnaghan型的等熵压缩线和相变能.结合前人的实验结果,得到了辉长岩的Grü neisen参数γ0,约为2.3.计算得到相变后的矿物集合体在零压下的密度为3.41 g·cm-3.根据以上参数分析得出,在此压力段辉长岩相变主要由长石的分解反应和石英的高压相变控制.     

哈达山水利枢纽坝基砂层震动液化研究

砂层震动液化问题是哈达山水利枢纽工程首要的工程地质问题。根据工程勘察取得的大量试验资料,通过标贯、GJS原状取砂、OYO 3000系列密度测井及瑞利波测试等综合方法,对砂层震动液化问题进行了深入的研究,为坝基加固处理提供了重要依据。其中瑞利波测试技术用于砂土震动液化是一个新的尝试。     

大跨度悬索桥地震响应控制分析

以主跨为1 490 m的润扬长江大桥为背景,采用ANSYS软件建立结构有限元模型,计算大跨度悬索桥动力特性。对黏滞阻尼器和软钢阻尼器进行参数敏感性分析,得出控制大跨度悬索桥地震响应的最优参数值,并分析一致激励和行波激励作用下黏滞阻尼器和软钢阻尼器减震效果。研究结果表明,黏滞阻尼器和软钢阻尼器对塔梁相对位移有较好的控制效果,但会使塔底内力有所增加;就位移控制而言,软钢阻尼器的效果更好;在低视波速区间内,黏滞阻尼器和软钢阻尼器减震效果明显存在波动特征;随着视波速的逐渐增大,黏滞阻尼器和软钢阻尼器减震效果受视波速的影响逐渐减小,悬索桥地震响应逐渐平缓,并趋于一致激励作用下的对应值。     

A Raman spectroscopic study of shock-wave densification of vitreous silica

The densification processes in SiO₂ glass induced by shock-wave compression up to 43.4 GPa are investigated by Raman spectroscopy. At first, densification increases with increasing shock pressure. A maximum densification of 11% is obtained for a shock pressure of 26.3 GPa. This densification is attributed to the reduction of the average Si−O−Si angle, which occurs first by the collapse of the largest ring cavities, then by further reduction of the average ring size. For higher shock pressures, a different structural modification is observed, resulting in decreasing densification with increasing shock pressure. Indeed, the recovered densification becomes very small, with values of 1.8 and 0.5% at 32 and 43.4 GPa, respectively. This is attributed to partial annealing of the samples due to high after shock residual temperatures. The study of the annealing process of the most densified glass by in situ high temperature Raman spectroscopy confirms that relaxation of the Si−O−Si angle starts at a lower temperature (about 800 K) than that of the siloxane rings (about 1000 K), thus explaining the high intensity of the siloxane defect bands in the samples schocked at compressions of 32 and 43.4 GPa. The large intensity of the siloxane bands in the nearly undensified samples shocked by compressions above 30 GPa may be explained by the relaxation during decompression of five- and six-fold coordinated silicon species formed at high pressure and high temperature during the shock event. Received: March 30, 1998 / Revised, accepted: August 21, 1998     

Diagnostics of supernova remnant shock waves

Fairly simple models can explain the emission from non-radiative shock waves in supernova remnants. This talk reviews some of the more robust diagnostics of shock parameters and some of the implications for the physics of collisionless shock waves in interstellar gas. The Hopkins Ultraviolet Telescope observed several non-radiative shocks during the ASTRO-2 mission, and the spectra have interesting implications for the physics of collisionless shocks.     

The deflection of jets by clouds

The role of collisions between extragalactic jets and dense clouds in determining the appearance of high-redshift radio galaxies is discussed and investigated through numerical hydrodynamic simulations in three dimensions. The code has the facility to track jet material separately from ambient material. This allows us to use simplifying assumptions to calculate synthetic radio images. The results indicate that the most powerful radio sources are likely to be observed during or shortly after an interaction, and that such interactions can explain both the radio structures and the spatial association between optical and radio light found in powerful radio galaxies. In some cases such a scenario may provide an alternative explanation of jet properties to mechanisms based on variations in the source or fluid-dynamical instabilities.This author is supported by a PPARC research studentship