Shock-wave propagation and cavitation bubble oscillation by Nd:YAG laser ablation of a metal in water.

PubMed

Chen, Xiao; Xu, Rong-Qing; Chen, Jian-Ping; Shen, Zhong-Hua; Jian, Lu; Ni, Xiao-Wu

2004-06-01

A highly sensitive fiber-optic sensor based on optical beam deflection is applied for investigating the propagation of a laser-induced plasma shock wave, the oscillation of a cavitation bubble diameter, and the development of a bubble-collapse-induced shock wave when a Nd:YAG laser pulse is focused upon an aluminum surface in water. By the sequence of experimental waveforms detected at different distances, the attenuation properties of the plasma shock wave and of the bubble-collapse-induced shock wave are obtained. Besides, based on characteristic signals, both the maximum and the minimum bubble radii at each oscillation cycle are determined, as are the corresponding oscillating periods.

Experimental investigation of starting characteristics and wave propagation from a shallow open cavity and its acoustic emission at supersonic speed

NASA Astrophysics Data System (ADS)

Pandian, S.; Desikan, S. L. N.; Niranjan, Sahoo

2018-01-01

Experiments were carried out on a shallow open cavity (L/D = 5) at a supersonic Mach number (M = 1.8) to understand its transient starting characteristics, wave propagation (inside and outside the cavity) during one vortex shedding cycle, and acoustic emission. Starting characteristics and wave propagation were visualized through time resolved schlieren images, while acoustic emissions were captured through unsteady pressure measurements. Results showed a complex shock system during the starting process which includes characteristics of the bifurcated shock system, shock train, flow separation, and shock wave boundary layer interaction. In one vortex shedding cycle, vortex convection from cavity leading edge to cavity trailing edge was observed. Flow features outside the cavity demonstrated the formation and downstream movement of a λ-shock due to the interaction of shock from the cavity leading edge and shock due to vortex and generation of waves on account of shear layer impingement at the cavity trailing edge. On the other hand, interesting wave structures and its propagation were monitored inside the cavity. In one vortex shedding cycle, two waves such as a reflected compression wave from a cavity leading edge in the previous vortex shedding cycle and a compression wave due to the reflection of Mach wave at the cavity trailing edge corner in the current vortex shedding cycle were visualized. The acoustic emission from the cavity indicated that the 2nd to 4th modes/tones are dominant, whereas the 1st mode contains broadband spectrum. In the present studies, the cavity feedback mechanism was demonstrated through a derived parameter coherence coefficient.

Kink-bands: Shock deformation of biotite resulting from a nuclear explosion

USGS Publications Warehouse

Cummings, D.

1965-01-01

Microscopic examination of granodiorite samples from the shock region around a nuclear explosion reveals sharply folded lens-shaped zones (kink-bands) in the mineral biotite. Fifty percent of these zones are oriented approximately 90?? to the direction of shock-wave propagation, but other zones are symmetrically concentrated at shear angles of 50?? and 70?? to the direction of shock-wave propagation.

Asymmetric Shock Wave Generation in a Microwave Rocket Using a Magnetic Field

NASA Astrophysics Data System (ADS)

Takahashi, Masayuki

2017-10-01

A plasma pattern is reproduced by coupling simulations between a particle-in- cell with Monte Carlo collisions model and a finite-difference time-domain simulation for an electromagnetic wave propagation when an external magnetic field is applied to the breakdown volume inside a microwave-rocket nozzle. The propagation speed and energy-absorption rate of the plasma are estimated based on the breakdown simulation, and these are utilized to reproduce shock wave propagation, which provides impulsive thrust for the microwave rocket. The shock wave propagation is numerically reproduced by solving the compressible Euler equation with an energy source of the microwave heating. The shock wave is asymmetrically generated inside the nozzle when the electron cyclotron resonance region has a lateral offset, which generates lateral and angular impulses for postural control of the vehicle. It is possible to develop an integrated device to maintain beaming ight of the microwave rocket, achieving both axial thrust improvement and postural control, by controlling the spatial distribution of the external magnetic field.

Shock Waves in a Bose-Einstein Condensate

NASA Technical Reports Server (NTRS)

Kulikov, Igor; Zak, Michail

2005-01-01

A paper presents a theoretical study of shock waves in a trapped Bose-Einstein condensate (BEC). The mathematical model of the BEC in this study is a nonlinear Schroedinger equation (NLSE) in which (1) the role of the wave function of a single particle in the traditional Schroedinger equation is played by a space- and time-dependent complex order parameter (x,t) proportional to the square root of the density of atoms and (2) the atoms engage in a repulsive interaction characterized by a potential proportional to | (x,t)|2. Equations that describe macroscopic perturbations of the BEC at zero temperature are derived from the NLSE and simplifying assumptions are made, leading to equations for the propagation of sound waves and the transformation of sound waves into shock waves. Equations for the speeds of shock waves and the relationships between jumps of velocity and density across shock fronts are derived. Similarities and differences between this theory and the classical theory of sound waves and shocks in ordinary gases are noted. The present theory is illustrated by solving the equations for the example of a shock wave propagating in a cigar-shaped BEC.

Optodynamic characterization of shock waves after laser-induced breakdown in water.

PubMed

Petkovsek, Rok; Mozina, Janez; Mocnik, Grisa

2005-05-30

Plasma and a cavitation bubble develop at the site of laser-induced breakdown in water. Their formation and the propagation of the shock wave were monitored by a beam-deflection probe and an arm-compensated interferometer. The interferometer part of the setup was used to determine the relative position of the laser-induced breakdown. The time-of-flight data from the breakdown site to the probe beam yielded the velocity, and from the velocity the shock-wave pressure amplitudes were calculated. Two regions were found where the pressure decays with different exponents, pointing to a strong attenuation mechanism in the initial phase of the shock-wave propagation.

Laser Beam Propagation Through Inhomogeneous Media with Shock-Like Profiles: Modeling and Computing

NASA Technical Reports Server (NTRS)

Adamovsky, Grigory; Ida, Nathan

1997-01-01

Wave propagation in inhomogeneous media has been studied for such diverse applications as propagation of radiowaves in atmosphere, light propagation through thin films and in inhomogeneous waveguides, flow visualization, and others. In recent years an increased interest has been developed in wave propagation through shocks in supersonic flows. Results of experiments conducted in the past few years has shown such interesting phenomena as a laser beam splitting and spreading. The paper describes a model constructed to propagate a laser beam through shock-like inhomogeneous media. Numerical techniques are presented to compute the beam through such media. The results of computation are presented, discussed, and compared with experimental data.

Comparison of acoustic shock waves generated by micro and nanosecond lasers for a smart laser surgery system

NASA Astrophysics Data System (ADS)

Nguendon Kenhagho, Hervé K.; Rauter, Georg; Guzman, Raphael; C. Cattin, Philippe; Zam, Azhar

2018-02-01

Characterization of acoustic shock wave will guarantee efficient tissue differentiation as feedback to reduce the probability of undesirable damaging (i.e. cutting) of tissues in laser surgery applications. We ablated hard (bone) and soft (muscle) tissues using a nanosecond pulsed Nd:YAG laser at 532 nm and a microsecond pulsed Er:YAG laser at 2.94 μm. When the intense short ns-pulsed laser is applied to material, the energy gain causes locally a plasma at the ablated spot that expands and propagates as an acoustic shock wave with a rarefaction wave behind the shock front. However, when using a μs-pulsed Er:YAG laser for material ablation, the acoustic shock wave is generated during the explosion of the ablated material. We measured and compared the emitted acoustic shock wave generated by a ns-pulsed Nd:YAG laser and a μs-pulsed Er:YAG laser measured by a calibrated microphone. As the acoustic shock wave attenuates as it propagates through air, the distance between ablation spots and a calibrated microphone was at 5 cm. We present the measurements on the propagation characteristics of the laser generated acoustic shock wave by measuring the arrival time-of-flight with a calibrated microphone and the energy-dependent evolution of acoustic parameters such as peak-topeak pressure, the ratio of the peak-to-peak pressures for the laser induced breakdown in air, the ablated muscle and the bone, and the spectral energy.

Shock wave propagation in a magnetic flux tube

NASA Astrophysics Data System (ADS)

Ferriz-Mas, A.; Moreno-Insertis, F.

1992-12-01

The propagation of a shock wave in a magnetic flux tube is studied within the framework of the Brinkley-Kirkwood theory adapted to a radiating gas. Simplified thermodynamic paths along which the compressed plasma returns to its initial state are considered. It is assumed that the undisturbed medium is uniform and that the flux tube is optically thin. The shock waves investigated, which are described with the aid of the thin flux-tube approximation, are essentially slow magnetohydrodynamic shocks modified by the constraint of lateral pressure balance between the flux tube and the surrounding field-free fluid; the confining external pressure must be balanced by the internal gas plus magnetic pressures. Exact analytical solutions giving the evolution of the shock wave are obtained for the case of weak shocks.

Numerical Simulation of Focused Shock Shear Waves in Soft Solids and a Two-Dimensional Nonlinear Homogeneous Model of the Brain

PubMed Central

Giammarinaro, B.; Coulouvrat, F.; Pinton, G.

2016-01-01

Shear waves that propagate in soft solids, such as the brain, are strongly nonlinear and can develop into shock waves in less than one wavelength. We hypothesize that these shear shock waves could be responsible for certain types of traumatic brain injuries (TBI) and that the spherical geometry of the skull bone could focus shear waves deep in the brain, generating diffuse axonal injuries. Theoretical models and numerical methods that describe nonlinear polarized shear waves in soft solids such as the brain are presented. They include the cubic nonlinearities that are characteristic of soft solids and the specific types of nonclassical attenuation and dispersion observed in soft tissues and the brain. The numerical methods are validated with analytical solutions, where possible, and with self-similar scaling laws where no known solutions exist. Initial conditions based on a human head X-ray microtomography (CT) were used to simulate focused shear shock waves in the brain. Three regimes are investigated with shock wave formation distances of 2.54 m, 0.018 m, and 0.0064 m. We demonstrate that under realistic loading scenarios, with nonlinear properties consistent with measurements in the brain, and when the shock wave propagation distance and focal distance coincide, nonlinear propagation can easily overcome attenuation to generate shear shocks deep inside the brain. Due to these effects, the accelerations in the focal are larger by a factor of 15 compared to acceleration at the skull surface. These results suggest that shock wave focusing could be responsible for diffuse axonal injuries. PMID:26833489

Shock Waves Propagation in Scope of the Nonlocal Theory of Dynamical Plasticity

NASA Astrophysics Data System (ADS)

Khantuleva, Tatyana A.

2004-07-01

From the point of view of the modern statistical mechanics the problems on shock compression of solids require a reformulation in terms of highly nonequilibrium effects arising inside the wave front. The self-organization during the multiscale and multistage momentum and energy exchange are originated by the correlation function. The theory of dynamic plasticity has been developed by the author on the base of the self-consistent nonlocal hydrodynamic approach had been applied to the shock wave propagation in solids. Nonlocal balance equations describe both the reversible wave type transport at the initial stage and the diffusive (dissipative) one in the end. The involved inverse influence of the mesoeffects on the wave propagation makes the formulation of problems self-consistent and involves a concept of the cybernetic control close-loop.

Current topics in shock waves; Proceedings of the International Symposium on Shock Waves and Shock Tubes, 17th, Lehigh University, Bethlehem, PA, July 17-21, 1989

NASA Astrophysics Data System (ADS)

Kim, Yong W.

Various papers on shock waves are presented. The general topics addressed include: shock formation, focusing, and implosion; shock reflection and diffraction; turbulence; laser-produced plasmas and waves; ionization and shock-plasma interaction; chemical kinetics, pyrolysis, and soot formation; experimental facilities, techniques, and applications; ignition of detonation and combustion; particle entrainment and shock propagation through particle suspension; boundary layers and blast simulation; computational methods and numerical simulation.

Detonative propagation and accelerative expansion of the Crab Nebula shock front.

PubMed

Gao, Yang; Law, Chung K

2011-10-21

The accelerative expansion of the Crab Nebula's outer envelope is a mystery in dynamics, as a conventional expanding blast wave decelerates when bumping into the surrounding interstellar medium. Here we show that the strong relativistic pulsar wind bumping into its surrounding nebula induces energy-generating processes and initiates a detonation wave that propagates outward to form the current outer edge, namely, the shock front, of the nebula. The resulting detonation wave, with a reactive downstream, then provides the needed power to maintain propagation of the shock front. Furthermore, relaxation of the curvature-induced reduction of the propagation velocity from the initial state of formation to the asymptotic, planar state of Chapman-Jouguet propagation explains the observed accelerative expansion. Potential richness in incorporating reactive fronts in the description of various astronomical phenomena is expected. © 2011 American Physical Society

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

NASA Technical Reports Server (NTRS)

Cambier, Jean-Luc

1991-01-01

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Two-zone elastic-plastic single shock waves in solids.

PubMed

Zhakhovsky, Vasily V; Budzevich, Mikalai M; Inogamov, Nail A; Oleynik, Ivan I; White, Carter T

2011-09-23

By decoupling time and length scales in moving window molecular dynamics shock-wave simulations, a new regime of shock-wave propagation is uncovered characterized by a two-zone elastic-plastic shock-wave structure consisting of a leading elastic front followed by a plastic front, both moving with the same average speed and having a fixed net thickness that can extend to microns. The material in the elastic zone is in a metastable state that supports a pressure that can substantially exceed the critical pressure characteristic of the onset of the well-known split-elastic-plastic, two-wave propagation. The two-zone elastic-plastic wave is a general phenomenon observed in simulations of a broad class of crystalline materials and is within the reach of current experimental techniques.

On A Problem Of Propagation Of Shock Waves Generated By Explosive Volcanic Eruptions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gusev, V. A.; Sobissevitch, A. L.

2008-06-24

Interdisciplinary study of flows of matter and energy in geospheres has become one of the most significant advances in Earth sciences. It is carried out by means of direct quantitative estimations based on detailed analysis of geological and geophysical observations and experimental data. The actual contribution is the interdisciplinary study of nonlinear acoustics and physical volcanology dedicated to shock wave propagation in a viscous and inhomogeneous medium. The equations governing evolution of shock waves with an arbitrary initial profile and an arbitrary cross-section of a beam are obtained. For the case of low viscous medium, the asymptotic solution meant tomore » calculate a profile of a shock wave in an arbitrary point has been derived. The analytical solution of the problem on propagation of shock pulses from atmosphere into a two-phase fluid-saturated geophysical medium is analysed. Quantitative estimations were carried out with respect to experimental results obtained in the course of real explosive volcanic eruptions.« less

Propagation of nonlinear shock waves for the generalised Oskolkov equation and its dynamic motions in the presence of an external periodic perturbation

NASA Astrophysics Data System (ADS)

Ak, Turgut; Aydemir, Tugba; Saha, Asit; Kara, Abdul Hamid

2018-06-01

Propagation of nonlinear shock waves for the generalised Oskolkov equation and dynamic motions of the perturbed Oskolkov equation are investigated. Employing the unified method, a collection of exact shock wave solutions for the generalised Oskolkov equations is presented. Collocation finite element method is applied to the generalised Oskolkov equation for checking the accuracy of the proposed method by two test problems including the motion of shock wave and evolution of waves with Gaussian and undular bore initial conditions. Considering an external periodic perturbation, the dynamic motions of the perturbed generalised Oskolkov equation are studied depending on the system parameters with the help of phase portrait and time series plot. The perturbed generalised Oskolkov equation exhibits period-3, quasiperiodic and chaotic motions for some special values of the system parameters, whereas the generalised Oskolkov equation presents shock waves in the absence of external periodic perturbation.

On the response of rubbers at high strain rates.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Niemczura, Johnathan Greenberg

In this report, we examine the propagation of tensile waves of finite deformation in rubbers through experiments and analysis. Attention is focused on the propagation of one-dimensional dispersive and shock waves in strips of latex and nitrile rubber. Tensile wave propagation experiments were conducted at high strain-rates by holding one end fixed and displacing the other end at a constant velocity. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in the rubber strips. Analysis of the response through the theory of finite waves and quantitative matching between themore » experimental observations and analytical predictions was used to determine an appropriate instantaneous elastic response for the rubbers. This analysis also yields the tensile shock adiabat for rubber. Dispersive waves as well as shock waves are also observed in free-retraction experiments; these are used to quantify hysteretic effects in rubber.« less

Propagation of acoustic shock waves between parallel rigid boundaries and into shadow zones

DOE Office of Scientific and Technical Information (OSTI.GOV)

Desjouy, C., E-mail: cyril.desjouy@gmail.com; Ollivier, S.; Dragna, D.

2015-10-28

The study of acoustic shock propagation in complex environments is of great interest for urban acoustics, but also for source localization, an underlying problematic in military applications. To give a better understanding of the phenomenon taking place during the propagation of acoustic shocks, laboratory-scale experiments and numerical simulations were performed to study the propagation of weak shock waves between parallel rigid boundaries, and into shadow zones created by corners. In particular, this work focuses on the study of the local interactions taking place between incident, reflected, and diffracted waves according to the geometry in both regular or irregular – alsomore » called Von Neumann – regimes of reflection. In this latter case, an irregular reflection can lead to the formation of a Mach stem that can modify the spatial distribution of the acoustic pressure. Short duration acoustic shock waves were produced by a 20 kilovolts electric spark source and a schlieren optical method was used to visualize the incident shockfront and the reflection/diffraction patterns. Experimental results are compared to numerical simulations based on the high-order finite difference solution of the two dimensional Navier-Stokes equations.« less

The propagation of the shock wave from a strong explosion in a plane-parallel stratified medium: the Kompaneets approximation

NASA Astrophysics Data System (ADS)

Olano, C. A.

2009-11-01

Context: Using certain simplifications, Kompaneets derived a partial differential equation that states the local geometrical and kinematical conditions that each surface element of a shock wave, created by a point blast in a stratified gaseous medium, must satisfy. Kompaneets could solve his equation analytically for the case of a wave propagating in an exponentially stratified medium, obtaining the form of the shock front at progressive evolutionary stages. Complete analytical solutions of the Kompaneets equation for shock wave motion in further plane-parallel stratified media were not found, except for radially stratified media. Aims: We aim to analytically solve the Kompaneets equation for the motion of a shock wave in different plane-parallel stratified media that can reflect a wide variety of astrophysical contexts. We were particularly interested in solving the Kompaneets equation for a strong explosion in the interstellar medium of the Galactic disk, in which, due to intense winds and explosions of stars, gigantic gaseous structures known as superbubbles and supershells are formed. Methods: Using the Kompaneets approximation, we derived a pair of equations that we call adapted Kompaneets equations, that govern the propagation of a shock wave in a stratified medium and that permit us to obtain solutions in parametric form. The solutions provided by the system of adapted Kompaneets equations are equivalent to those of the Kompaneets equation. We solved the adapted Kompaneets equations for shock wave propagation in a generic stratified medium by means of a power-series method. Results: Using the series solution for a shock wave in a generic medium, we obtained the series solutions for four specific media whose respective density distributions in the direction perpendicular to the stratification plane are of an exponential, power-law type (one with exponent k=-1 and the other with k =-2) and a quadratic hyperbolic-secant. From these series solutions, we deduced exact solutions for the four media in terms of elemental functions. The exact solution for shock wave propagation in a medium of quadratic hyperbolic-secant density distribution is very appropriate to describe the growth of superbubbles in the Galactic disk. Member of the Carrera del Investigador Científico del CONICET, Argentina.

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading

NASA Technical Reports Server (NTRS)

Gentz, Steven J.; Ordway, David O; Parsons, David S.; Garrison, Craig M.; Rodgers, C. Steven; Collins, Brian W.

2015-01-01

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approx. 9 inches from the source) dominated by direct wave propagation, mid-field environment (approx. 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This report documents the outcome of the assessment.

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading. [Appendices

NASA Technical Reports Server (NTRS)

Gentz, Steven J.; Ordway, David O.; Parsons, David S.; Garrison, Craig M.; Rodgers, C. Steven; Collins, Brian W.

2015-01-01

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (9 inches from the source) dominated by direct wave propagation, mid-field environment (approximately 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading: Appendices

NASA Technical Reports Server (NTRS)

Gentz, Steven J.; Ordway, David O.; Parsons, David S.; Garrison, Craig M.; Rodgers, C. Steven; Collins, Brian W.

2015-01-01

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approx. 9 inches from the source) dominated by direct wave propagation, mid-field environment (approx. 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

CYLINDRICAL WAVES OF FINITE AMPLITUDE IN DISSIPATIVE MEDIUM (in Russian)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Naugol'nykh, K.A.; Soluyan, S.I.; Khokhlov, R.V.

1962-07-01

Propagation of diverging and converging cylindrical waves in a nonlinear, viscous, heat conducting medium is analyzed using approximation methods. The KrylovBogolyubov method was used for small Raynold's numbers, and the method of S. I. Soluyan et al. (Vest. Mosk. Univ. ser. phys. and astronomy 3, 52-81, 1981), was used for large Raynold's numbers. The formation and dissipation of shock fronts and spatial dimensions of shock phenomena were analyzed. It is shown that the problem of finiteamplitude cylindrical wave propagation is identical to the problem of plane wave propagations in a medium with variable viscosity. (tr-auth)

Molecular Dynamics Simulations of Shock Wave Propagation across the Nitromethane Crystal-Melt Interface

NASA Astrophysics Data System (ADS)

Jiang, Shan; Sewell, Thomas D.; Thompson, Donald L.

2015-06-01

We are interested in understanding the fundamental processes that occur during propagation of shock waves across the crystal-melt interface in molecular substances. We have carried out molecular dynamics simulations of shock passage from the nitromethane (100)-oriented crystal into the melt and vice versa using the fully flexible, non-reactive Sorescu, Rice, and Thompson force field. A stable interface was established for a temperature near the melting point by using a combination of isobaric-isothermal (NPT) and isochoric-isothermal (NVT) simulations. The equilibrium bulk and interfacial regions were characterized using spatial-temporal distributions of molecular number density, kinetic and potential energy, and C-N bond orientations. Those same properties were calculated as functions of time during shock propagation. As expected, the local temperatures (intermolecular, intramolecular, and total) and stress states differed significantly between the liquid and crystal regions and depending on the direction of shock propagation. Substantial differences in the spatial distribution of shock-induced defect structures in the crystalline region were observed depending on the direction of shock propagation. Research supported by the U.S. Army Research Office.

Observation of laser-driven shock propagation by nanosecond time-resolved Raman spectroscopy

NASA Astrophysics Data System (ADS)

Yu, Guoyang; Zheng, Xianxu; Song, Yunfei; Zeng, Yangyang; Guo, Wencan; Zhao, Jun; Yang, Yanqiang

2015-01-01

An improved nanosecond time-resolved Raman spectroscopy is performed to observe laser-driven shock propagation in the anthracene/epoxy glue layer. The digital delay instead of optical delay line is introduced for sake of unlimited time range of detection, which enables the ability to observe both shock loading and shock unloading that always lasts several hundred nanoseconds. In this experiment, the peak pressure of shock wave, the pressure distribution, and the position of shock front in gauge layer were determined by fitting Raman spectra of anthracene using the Raman peak shift simulation. And, the velocity of shock wave was calculated by the time-dependent position of shock front.

The Universal Role of Tubulence in the Propagation of Strong Shocks and Detonation Waves

NASA Astrophysics Data System (ADS)

Lee, John H.

2001-06-01

The passage of a strong shock wave usually results in irreversible physical and chemical changes in the medium. If the chemical reactions are sufficiently exothermic, the shock wave can be self-propagating, i.e., sustained by the chemical energy release via the expansion work of the reaction products. Although shocks and detonations can be globally stable and propagate at constant velocities (in the direction of motion), their structure may be highly unstable and exhibit large hydrodynamic fluctuations, i.e., turbulence. Recent investigations on plastic deformation of polycrystalline material behind shock waves have revealed particle velocity dispersion at the mesoscopic level, a result of vortical rotational motion similar to that of turbulent fluid flows at high Reynolds number.1 Strong ionizing shocks in noble gases2, as well as dissociating shock waves in carbon dioxide,3 also demonstrate a turbulent density fluctuation in the non-equilibrium shock transition zone. Perhaps the most thoroughly investigated unstable structure is that of detonation waves in gaseous explosives.4 Detonation waves in liquid explosives such as nitromethane also take on similar unstable structure as gaseous detonations.5 There are also indications that detonations in solid explosives have a similar unsteady structure under certain conditions. Thus, it appears that it is more of a rule than an exception that the structure of strong shocks and detonations are unstable and exhibit turbulent-like fluctuations as improved diagnostics now permit us to look more closely at the meso- and micro-levels. Increasing attention is now devoted to the understanding of the shock waves at the micro-scale level in recent years. This is motivated by the need to formulate physical and chemical models that contain the correct physics capable of describing quantitatively the shock transition process. It should be noted that, in spite of its unstable 3-D structure, the steady 1-D conservation laws (in the direction of propagation) apply across the shock transition zone if the downstream equilibrium plane is taken far enough away to ensure the decay of the turbulent fluctuations. Thus, the Hugoniot properties of one-dimensional propagation of shock and detonation waves remain valid. However, the conservation laws do not describe the important propagation mechanisms (i.e., the physical and chemical processes that effect the transition from initial to the final state) in the wave structure. Since gaseous detonations enjoy the advantage of being able to be observed experimentally in great detail, its complex turbulent structure is now quite well established. Furthermore, the equation of state for perfect gases is well known and the chemistry of most gas phase reactions is also sufficiently understood quantitatively to permit detailed numerical simulation of the complex detonation structure. Thus, a good database of information exists for gaseous detonation, and in this paper we shall explore the turbulent structure of gaseous detonation with the aim of answering the question as to "why nature prefers to evoke such a complicated manner to effect its propagation." We will then attempt to generalize the discussion to the "terra incognita" of condensed phase materials where the structure is much less understood. 1. Meshcheryakov, Yu.I., and Atroshenko, S.A., Izv. Vyssh. Uchebn. Zaved. Fiz., 4, 105-123 (1992). 2. Glass, I.I, and Liu, W.S., J. Fluid Mech., 84(1), 55-77 (1978). 3. Griffiths, R.W., Sanderman, R.J., and Hornung, H.G., J. Phys. D., 8, 1681-1691 (1975). 4. Lee, J.H.S., Ann. Rev. Fluid Mech., 16, 311-336 (1984). 5. Mallory, H.D., J. Appl. Physics, 38, 5302-5306 (1967).

Shock waves in weakly compressed granular media.

PubMed

van den Wildenberg, Siet; van Loo, Rogier; van Hecke, Martin

2013-11-22

We experimentally probe nonlinear wave propagation in weakly compressed granular media and observe a crossover from quasilinear sound waves at low impact to shock waves at high impact. We show that this crossover impact grows with the confining pressure P0, whereas the shock wave speed is independent of P0-two hallmarks of granular shocks predicted recently. The shocks exhibit surprising power law attenuation, which we model with a logarithmic law implying that shock dissipation is weak and qualitatively different from other granular dissipation mechanisms. We show that elastic and potential energy balance in the leading part of the shocks.

Proceedings of the 15th International Symposium on Shock Waves and Shock Tubes

NASA Astrophysics Data System (ADS)

Bershader, Daniel; Hanson, Ronald

1986-09-01

One hundred ten papers were presented in 32 sessions. Topics included: The application of Hook-method spectroscopy to the diagnosis of shock-heated gases. The nonintrusive destruction of kidney stones by underwater focused shock waves. Several of the papers reflect the recent and continuing interest in shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive configurations. The major subject areas were: shock propagation and interactions; shock-general chemical kinetics; shock computation, modeling, and stability problems; shock wave aerodynamics; experimental methods; shocks in multiphase and heterogeneous media; high energy gas excitation and wave phenomena; and technical applications and shocks in condensed matter.

Propagation characteristics of pulverized coal and gas two-phase flow during an outburst.

PubMed

Zhou, Aitao; Wang, Kai; Fan, Lingpeng; Tao, Bo

2017-01-01

Coal and gas outbursts are dynamic failures that can involve the ejection of thousands tons of pulverized coal, as well as considerable volumes of gas, into a limited working space within a short period. The two-phase flow of gas and pulverized coal that occurs during an outburst can lead to fatalities and destroy underground equipment. This article examines the interaction mechanism between pulverized coal and gas flow. Based on the role of gas expansion energy in the development stage of outbursts, a numerical simulation method is proposed for investigating the propagation characteristics of the two-phase flow. This simulation method was verified by a shock tube experiment involving pulverized coal and gas flow. The experimental and simulated results both demonstrate that the instantaneous ejection of pulverized coal and gas flow can form outburst shock waves. These are attenuated along the propagation direction, and the volume fraction of pulverized coal in the two-phase flow has significant influence on attenuation of the outburst shock wave. As a whole, pulverized coal flow has a negative impact on gas flow, which makes a great loss of large amounts of initial energy, blocking the propagation of gas flow. According to comparison of numerical results for different roadway types, the attenuation effect of T-type roadways is best. In the propagation of shock wave, reflection and diffraction of shock wave interact through the complex roadway types.

Propagation characteristics of pulverized coal and gas two-phase flow during an outburst

PubMed Central

Zhou, Aitao; Wang, Kai; Fan, Lingpeng; Tao, Bo

2017-01-01

Coal and gas outbursts are dynamic failures that can involve the ejection of thousands tons of pulverized coal, as well as considerable volumes of gas, into a limited working space within a short period. The two-phase flow of gas and pulverized coal that occurs during an outburst can lead to fatalities and destroy underground equipment. This article examines the interaction mechanism between pulverized coal and gas flow. Based on the role of gas expansion energy in the development stage of outbursts, a numerical simulation method is proposed for investigating the propagation characteristics of the two-phase flow. This simulation method was verified by a shock tube experiment involving pulverized coal and gas flow. The experimental and simulated results both demonstrate that the instantaneous ejection of pulverized coal and gas flow can form outburst shock waves. These are attenuated along the propagation direction, and the volume fraction of pulverized coal in the two-phase flow has significant influence on attenuation of the outburst shock wave. As a whole, pulverized coal flow has a negative impact on gas flow, which makes a great loss of large amounts of initial energy, blocking the propagation of gas flow. According to comparison of numerical results for different roadway types, the attenuation effect of T-type roadways is best. In the propagation of shock wave, reflection and diffraction of shock wave interact through the complex roadway types. PMID:28727738

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading. Volume 2, Part 1; Appendices

NASA Technical Reports Server (NTRS)

Gentz, Steven J.; Ordway, David O.; Parsons, David S.; Garrison, Craig M.; Rodgers, C. Steven; Collins, Brian W.

2015-01-01

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approximately 9 inches from the source) dominated by direct wave propagation, mid-field environment (approximately 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

Collisionless slow shocks in magnetotail reconnection

NASA Astrophysics Data System (ADS)

Cremer, Michael; Scholer, Manfred

The kinetic structure of collisionless slow shocks in the magnetotail is studied by solving the Riemann problem of the collapse of a current sheet with a normal magnetic field component using 2-D hybrid simulations. The collapse results in a current layer with a hot isotropic distribution and backstreaming ions in a boundary layer. The lobe plasma outside and within the boundary layer exhibits a large perpendicular to parallel temperature anisotropy. Waves in both regions propagate parallel to the magnetic field. In a second experiment a spatially limited high density beam is injected into a low beta background plasma and the subsequent wave excitation is studied. A model for slow shocks bounding the reconnection layer in the magnetotail is proposed where backstreaming ions first excite obliquely propagating waves by the electromagnetic ion/ion cyclotron instability, which lead to perpendicular heating. The T⊥/T∥ temperature anisotropy subsequently excites parallel propagating Alfvén ion cyclotron waves, which are convected into the slow shock and are refracted in the downstream region.

Shock tubes and waves; Proceedings of the Sixteenth International Symposium, Rheinisch-Westfaelische Technische Hochschule, Aachen, Federal Republic of Germany, July 26-31, 1987

NASA Astrophysics Data System (ADS)

Groenig, Hans

Topics discussed in this volume include shock wave structure, propagation, and interaction; shocks in condensed matter, dusty gases, and multiphase media; chemical processes and related combustion and detonation phenomena; shock wave reflection, diffraction, and focusing; computational fluid dynamic code development and shock wave application; blast and detonation waves; advanced shock tube technology and measuring technique; and shock wave applications. Papers are presented on dust explosions, the dynamics of shock waves in certain dense gases, studies of condensation kinetics behind incident shock waves, the autoignition mechanism of n-butane behind a reflected shock wave, and a numerical simulation of the focusing process of reflected shock waves. Attention is also given to the equilibrium shock tube flow of real gases, blast waves generated by planar detonations, modern diagnostic methods for high-speed flows, and interaction between induced waves and electric discharge in a very high repetition rate excimer laser.

Evidence for self-refraction in a convergence zone: NPE (Nonlinear progressive wave equation) model results

NASA Technical Reports Server (NTRS)

Mcdonald, B. Edward; Plante, Daniel R.

1989-01-01

The nonlinear progressive wave equation (NPE) model was developed by the Naval Ocean Research and Development Activity during 1982 to 1987 to study nonlinear effects in long range oceanic propagation of finite amplitude acoustic waves, including weak shocks. The NPE model was applied to propagation of a generic shock wave (initial condition provided by Sandia Division 1533) in a few illustrative environments. The following consequences of nonlinearity are seen by comparing linear and nonlinear NPE results: (1) a decrease in shock strength versus range (a well-known result of entropy increases at the shock front); (2) an increase in the convergence zone range; and (3) a vertical meandering of the energy path about the corresponding linear ray path. Items (2) and (3) are manifestations of self-refraction.

CONTRIBUTION OF VELOCITY VORTICES AND FAST SHOCK REFLECTION AND REFRACTION TO THE FORMATION OF EUV WAVES IN SOLAR ERUPTIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Hongjuan; Liu, Siqing; Gong, Jiancun

2015-06-01

We numerically study the detailed evolutionary features of the wave-like disturbance and its propagation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode shock at the front of the rope, which then produces a type II radio burst. The expansion of the fast shock, whichmore » is associated with outward motion, takes place in various directions, and the downward expansion shows the reflection and the refraction as a result of the non-uniform background plasma. The reflected component of the fast shock propagates upward and the refracted component propagates downward. As the refracted component reaches the boundary surface, a weak echo is excited. The Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may arise from various processes that are revealed in the new simulations.« less

Shock wave propagation within a confined multi-chamber system

NASA Astrophysics Data System (ADS)

Julien, B.; Sochet, I.; Tadini, P.; Vaillant, T.

2018-07-01

The influence of a variation of the opening ratios of rooms and side walls on the propagation of a shock wave within a confined multi-chamber system is analyzed through the evolution of some of the shock parameters (maximum overpressure and positive impulse). The shock wave is generated by the detonation of a hemispherical gaseous charge in one of the rooms. Several small-scale experiments have been carried out using an adjustable model representative of a pyrotechnic workshop. Using the same approach as for a previous article dealing with the impact of the volume of the rooms, we were able to link the evolution of the arrival time of the shock wave within the building with the reference obtained in the free field. Moreover, using a new parameter taking into account the opening ratios of the rooms and side walls, a predictive law was developed to model the maximal overpressure in the rooms.

Observation of shock transverse waves in elastic media.

PubMed

Catheline, S; Gennisson, J-L; Tanter, M; Fink, M

2003-10-17

We report the first experimental observation of a shock transverse wave propagating in an elastic medium. This observation was possible because the propagation medium, a soft solid, allows one to reach a very high Mach number. In this extreme configuration, the shock formation is observed over a distance of less than a few wavelengths, thanks to a prototype of an ultrafast scanner (that acquires 5000 frames per second). A comparison of these new experimental data with theoretical predictions, based on a modified Burger's equation, shows good agreement.

Prediction of Sound Waves Propagating Through a Nozzle Without/With a Shock Wave Using the Space-Time CE/SE Method

NASA Technical Reports Server (NTRS)

Wang, Xiao-Yen; Chang, Sin-Chung; Jorgenson, Philip C. E.

2000-01-01

The benchmark problems in Category 1 (Internal Propagation) of the third Computational Aeroacoustics (CAA) Work-shop sponsored by NASA Glenn Research Center are solved using the space-time conservation element and solution element (CE/SE) method. The first problem addresses the propagation of sound waves through a nearly choked transonic nozzle. The second one concerns shock-sound interaction in a supersonic nozzle. A quasi one-dimension CE/SE Euler solver for a nonuniform mesh is developed and employed to solve both problems. Numerical solutions are compared with the analytical solution for both problems. It is demonstrated that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple way. Furthermore, the simple nonreflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well.

Longitudinal waves in a perpendicular collisionless plasma shock. IV - Gradient B.

NASA Technical Reports Server (NTRS)

Gary, S. P.

1972-01-01

The consideration of elastic waves in a Vlasov plasma of unmagnetized ions and magnetized electrons undergoing E x B electron drift and gradient B drift, pursued in the earlier three parts, is brought to conclusion in this last part of the longitudinal wave study in a collisionless plasma shock. Detailed calculations of the effects of the beta sub e dimensionless parameter on the E x B electron drift instability are presented. It is shown that the range of propagation of the elastic waves about the perpendicular remains quite narrow, and that, for oblique propagation, the already narrow angular range of unstable waves is decreased by increases in the value of the beta sub e dimensionless parameter. Also, increases in wave number generally reduce the growth rate and the angular range of propagation.

Behavior of dusty real gas on adiabatic propagation of cylindrical imploding strong shock waves

NASA Astrophysics Data System (ADS)

Gangwar, P. K.

2018-05-01

In this paper, CCW method has been used to study the behavior of dusty real gas on adiabatic propagation of cylindrical imploding strong shock waves. The strength of overtaking waves is estimated under the assumption that both C+ and C- disturbances propagate in non-uniform region of same density distribution. It is assumed that the dusty gas is the mixture of a real gas and a large number of small spherical solid particles of uniform size. The solid particles are uniformly distributed in the medium. Maintaining equilibrium flow conditions, the expressions for shock strength has been derived both for freely propagation as well as under the effect of overtaking disturbances. The variation of all flow variables with propagation distance, mass concentration of solid particles in the mixture and the ratio of solid particles to the initial density of gas have been computed and discussed through graphs. It is found that the presence of dust particles in the gases medium has significant effects on the variation of flow variables and the shock is strengthened under the influence of overtaking disturbances. The results accomplished here been compared with those for ideal gas.

Modeling Propagation of Shock Waves in Metals

DOE Office of Scientific and Technical Information (OSTI.GOV)

Howard, W M; Molitoris, J D

2005-08-19

We present modeling results for the propagation of strong shock waves in metals. In particular, we use an arbitrary Lagrange Eulerian (ALE3D) code to model the propagation of strong pressure waves (P {approx} 300 to 400 kbars) generated with high explosives in contact with aluminum cylinders. The aluminum cylinders are assumed to be both flat-topped and have large-amplitude curved surfaces. We use 3D Lagrange mechanics. For the aluminum we use a rate-independent Steinberg-Guinan model, where the yield strength and shear modulus depend on pressure, density and temperature. The calculation of the melt temperature is based on the Lindermann law. Atmore » melt the yield strength and shear modulus is set to zero. The pressure is represented as a seven-term polynomial as a function of density. For the HMX-based high explosive, we use a JWL, with a program burn model that give the correct detonation velocity and C-J pressure (P {approx} 390 kbars). For the case of the large-amplitude curved surface, we discuss the evolving shock structure in terms of the early shock propagation experiments by Sakharov.« less

Modeling Propagation of Shock Waves in Metals

NASA Astrophysics Data System (ADS)

Howard, W. M.; Molitoris, J. D.

2006-07-01

We present modeling results for the propagation of strong shock waves in metals. In particular, we use an arbitrary Lagrange Eulerian (ALE3D) code to model the propagation of strong pressure waves (P ˜ 300 to 400 kbars) generated with high explosives in contact with aluminum cylinders. The aluminum cylinders are assumed to be both flat-topped and have large-amplitude curved surfaces. We use 3D Lagrange mechanics. For the aluminum we use a rate-independent Steinberg-Guinan model, where the yield strength and shear modulus depend on pressure, density and temperature. The calculation of the melt temperature is based on the Lindermann law. At melt the yield strength and shear modulus is set to zero. The pressure is represented as a seven-term polynomial as a function of density. For the HMX-based high explosive, we use a JWL, with a program burn model that give the correct detonation velocity and C-J pressure (P ˜ 390 kbars). For the case of the large-amplitude curved surface, we discuss the evolving shock structure in terms of the early shock propagation experiments by Sakharov.

Propagation of shock waves in elastic solids caused by cavitation microjet impact. II: Application in extracorporeal shock wave lithotripsy.

PubMed

Zhong, P; Chuong, C J; Preminger, G M

1993-07-01

To better understand the mechanism of stone fragmentation during extracorporeal shock wave lithotripsy (ESWL), the model developed in Part I [P. Zhong and C.J. Chuong, J. Acoust. Soc. Am. 94, 19-28 (1993)] is applied to study cavitation microjet impingement and its resultant shock wave propagation in renal calculi. Impact pressure at the stone boundary and stress, strain at the propagating shock fronts in the stone were calculated for typical ESWL loading conditions. At the anterior surface of the stone, the jet induced compressive stress can vary from 0.82 approximately 4 times that of the water hammer pressure depending on the contact angles; whereas the jet-induced shear stress can achieve its maximum, with a magnitude of 30% approximately 54% of the water hammer pressure, near the detachment of the longitudinal (or P) wave in the solid. Comparison of model predictions with material failure strengths of renal calculi suggests that jet impact can lead to stone surface erosion by combined compressive and shear loadings at the jet impacting surface, and spalling failure by tensile forces at the distal surface of the stone. Comparing responses from four different stone types suggests that cystine is the most difficult stone to fragment in ESWL, as observed from clinical experience.

Development of relativistic shock waves in viscous gluon matter

NASA Astrophysics Data System (ADS)

Bouras, I.; Molnár, E.; Niemi, H.; Xu, Z.; El, A.; Fochler, O.; Greiner, C.; Rischke, D. H.

2009-11-01

To investigate the formation and the propagation of relativistic shock waves in viscous gluon matter we solve the relativistic Riemann problem using a microscopic parton cascade. We demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to entropy density ratio η/s. We show that an η/s ratio larger than 0.2 prevents the development of well-defined shock waves on time scales typical for ultrarelativistic heavy-ion collisions. These findings are confirmed by viscous hydrodynamic calculations.

Molecular dynamics simulations of shock waves in oriented nitromethane single crystals.

PubMed

He, Lan; Sewell, Thomas D; Thompson, Donald L

2011-03-28

The structural relaxation of crystalline nitromethane initially at T = 200 K subjected to moderate (~15 GPa) supported shocks on the (100), (010), and (001) crystal planes has been studied using microcanonical molecular dynamics with the nonreactive Sorescu-Rice-Thompson force field [D. C. Sorescu, B. M. Rice, and D. L. Thompson, J. Phys. Chem. B 104, 8406 (2000)]. The responses to the shocks were determined by monitoring the mass density, the intermolecular, intramolecular, and total temperatures (average kinetic energies), the partitioning of total kinetic energy among Cartesian directions, the radial distribution functions for directions perpendicular to those of shock propagation, the mean-square displacements in directions perpendicular to those of shock propagation, and the time dependence of molecular rotational relaxation as a function of time. The results show that the mechanical response of crystalline nitromethane strongly depends on the orientation of the shock wave. Shocks propagating along [100] and [001] result in translational disordering in some crystal planes but not in others, a phenomenon that we refer to as plane-specific disordering; whereas for [010] the shock-induced stresses are relieved by a complicated structural rearrangement that leads to a paracrystalline structure. The plane-specific translational disordering is more complete by the end of the simulations (~6 ps) for shock propagation along [001] than along [100]. Transient excitation of the intermolecular degrees of freedom occurs in the immediate vicinity of the shock front for all three orientations; the effect is most pronounced for the [010] shock. In all three cases excitation of molecular vibrations occurs more slowly than the intermolecular excitation. The intermolecular and intramolecular temperatures are nearly equal by the end of the simulations, with 400-500 K of net shock heating. Results for two-dimensional mean-square molecular center-of-mass displacements, calculated as a function of time since shock wave passage in planes perpendicular to the direction of shock propagation, show that the molecular translational mobility in the picoseconds following shock wave passage is greatest for [001] and least for the [010] case. In all cases the root-mean-square center-of-mass displacement is small compared to the molecular diameter of nitromethane on the time scale of the simulations. The calculated time scales for the approach to thermal equilibrium are generally consistent with the predictions of a recent theoretical analysis due to Hooper [J. Chem. Phys. 132, 014507 (2010)].

Modeling shock waves in an ideal gas: combining the Burnett approximation and Holian's conjecture.

PubMed

He, Yi-Guang; Tang, Xiu-Zhang; Pu, Yi-Kang

2008-07-01

We model a shock wave in an ideal gas by combining the Burnett approximation and Holian's conjecture. We use the temperature in the direction of shock propagation rather than the average temperature in the Burnett transport coefficients. The shock wave profiles and shock thickness are compared with other theories. The results are found to agree better with the nonequilibrium molecular dynamics (NEMD) and direct simulation Monte Carlo (DSMC) data than the Burnett equations and the modified Navier-Stokes theory.

Effect of nonthermal electrons on the propagation characteristics and stability of two-dimensional nonlinear electrostatic coherent structures in relativistic electron positron ion plasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masood, W.; National Centre for Physics; Rizvi, H.

2011-06-15

Two-dimensional propagation of nonlinear ion acoustic shock and solitary waves in an unmagnetized plasma consisting of nonthermal electrons, Boltzmannian positrons, and singly charged hot ions streaming with relativistic velocities are investigated. The system of fluid equations is reduced to Kadomtsev-Petviashvili-Burgers and Kadomtsev-Petviashvili (KP) equations in the limit of small amplitude perturbation. The dependence of the ion acoustic shock and solitary waves on various plasma parameters are explored in detail. Interestingly, it is observed that increasing the nonthermal electron population increases the wave dispersion which enervates the strength of the ion acoustic shock wave; however, the same effect leads to anmore » enhancement of the soliton amplitude due to the absence of dissipation in the KP equation. The present investigation may be useful to understand the two-dimensional propagation characteristics of small but finite amplitude localized shock and solitary structures in planetary magnetospheres and auroral plasmas where nonthermal populations of electrons have been observed by several satellite missions.« less

Experimental Investigations on Microshock Waves and Contact Surfaces

NASA Astrophysics Data System (ADS)

Kai, Yun; Garen, Walter; Teubner, Ulrich

2018-02-01

The present work reports on progress in the research of a microshock wave. Because of the lack of a good understanding of the propagation mechanism of the microshock flow system (shock wave, contact surface, and boundary layer), the current work concentrates on measuring microshock flows with special attention paid to the contact surface. A novel setup involving a glass capillary (with a 200 or 300 μ m hydraulic diameter D ) and a high-speed magnetic valve is applied to generate a shock wave with a maximum initial Mach number of 1.3. The current work applies a laser differential interferometer to perform noncontact measurements of the microshock flow's trajectory, velocity, and density. The current work presents microscale measurements of the shock-contact distance L that solves the problem of calculating the scaling factor Sc =Re ×D /(4 L ) (introduced by Brouillette), which is a parameter characterizing the scaling effects of shock waves. The results show that in contrast to macroscopic shock waves, shock waves at the microscale have a different propagation or attenuation mechanism (key issue of this Letter) which cannot be described by the conventional "leaky piston" model. The main attenuation mechanism of microshock flow may be the ever slower moving contact surface, which drives the shock wave. Different from other measurements using pressure transducers, the current setup for density measurements resolves the whole microshock flow system.

Obliquely propagating low frequency electromagnetic shock waves in two dimensional quantum magnetoplasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masood, W.

2009-04-15

Linear and nonlinear propagation characteristics of low frequency magnetoacoustic waves in quantum magnetoplasmas are studied employing the quantum magnetohydrodynamic model. In this regard, a quantum Kadomtsev-Petviashvili-Burgers (KPB) equation is derived using the small amplitude expansion method. The dissipation is introduced by taking into account the kinematic viscosity among the plasma constituents. Furthermore, the solution of KPB equation is presented using the tangent hyperbolic (tanh) method. The variation in the fast and slow magnetoacoustic shock profiles with the quantum Bohm potential via increasing number density, obliqueness angle {theta}, magnetic field, and the resistivity are also investigated. It is observed that themore » aforementioned plasma parameters significantly modify the propagation characteristics of nonlinear magnetoacoustic shock waves in quantum magnetoplasmas. The relevance of the present investigation with regard to dense astrophysical environments is also pointed out.« less

Detonation suppression in hydrogen-air mixtures using porous coatings on the walls

NASA Astrophysics Data System (ADS)

Bivol, G. Yu.; Golovastov, S. V.; Golub, V. V.

2018-05-01

We considered the problem of detonation suppression and weakening of blast wave effects occurring during the combustion of hydrogen-air mixtures in confined spaces. The gasdynamic processes during combustion of hydrogen, an alternative environmentally friendly fuel, were also considered. Detonation decay and flame propagation in hydrogen-air mixtures were experimentally investigated in rectangular cross-section channels with solid walls and two types of porous coatings: steel wool and polyurethane foam. Shock wave pressure dynamics inside the section with porous coating were studied using pressure sensors; flame front propagation was studied using photodiodes and high-speed camera visualization. For all mixtures, the detonation wave formed before entering the section with porous coating. For both porous materials, the steady detonation wave decoupled in the porous section of the channel into a shock wave and flame front propagating with a velocity around the Chapman-Jouguet acoustic velocity. By the end of the porous section, shock wave pressure reductions of 70 and 85% were achieved for the polyurethane foam and steel wool, respectively. The dependence of the flame velocity on the mixture composition (equivalence ratio) is presented.

On the Generation of Hydrodynamic Shocks by Mixed Beams and Occurrence of Sunquakes in Flares

NASA Astrophysics Data System (ADS)

Zharkova, Valentina; Zharkov, Sergei

2015-11-01

Observations of solar flares with sunquakes by space- and ground-based instruments reveal essentially different dynamics of seismic events in different flares. Some sunquakes are found to be closely associated with the locations of hard X-ray (HXR) and white-light (WL) emission, while others are located outside either of them. In this article we investigate possible sources causing a seismic response in a form of hydrodynamic shocks produced by the injection of mixed (electron plus proton) beams, discuss the velocities of these shocks, and the depths where they deposit the bulk of their energy and momentum. The simulation of hydrodynamic shocks in flaring atmospheres induced by electron-rich and proton-rich beams reveals that the linear depth of the shock termination is shifted beneath the level of the quiet solar photosphere on a distance from 200 to 5000 km. The parameters of these atmospheric hydrodynamic shocks are used as initial condition for another hydrodynamic model developed for acoustic-wave propagation in the solar interior (Zharkov, Mon. Not. Roy. Astron. Soc. 431, 3414, 2013). The model reveals that the depth of energy and momentum deposition by the atmospheric shocks strongly affects the propagation velocity of the acoustic-wave packet in the interior. The locations of the first bounces from the photosphere of acoustic waves generated in the vicinity of a flare are seen as ripples on the solar surface, or sunquakes. Mixed proton-dominated beams are found to produce a strong supersonic shock at depths 200 - 300 km under the level of the quiet-Sun photosphere and in this way produce well-observable acoustic waves, while electron-dominated beams create a slightly supersonic shock propagating down to 5000 km under the photosphere. This shock can only generate acoustic waves at the top layers beneath the photosphere since the shock velocity very quickly drops below the local sound speed. The distance Δ of the first bounce of the generated acoustic waves is discussed in relation to the minimal phase velocities of wave packets defined by the acoustic cutoff frequency and the parameters of atmospheric shock termination beneath the photosphere.

The existence of electron-acoustic shock waves and their interactions in a non-Maxwellian plasma with q-nonextensive distributed electrons

DOE Office of Scientific and Technical Information (OSTI.GOV)

Han, Jiu-Ning; He, Yong-Lin; Han, Zhen-Hai

2013-07-15

We present a theoretical investigation for the nonlinear interaction between electron-acoustic shock waves in a nonextensive two-electron plasma. The interaction is governed by a pair of Korteweg-de Vries-Burgers equations. We focus on studying the colliding effects on the propagation of shock waves, more specifically, we have studied the effects of plasma parameters, i.e., the nonextensive parameter q, the “hot” to “cold” electron number density ratio α, and the normalized electron kinematic viscosity η{sub 0} on the trajectory changes (phase shifts) of shock waves. It is found that there are trajectory changes (phase shifts) for both colliding shock waves in themore » present plasma system. We also noted that the nonlinearity has no decisive effect on the trajectory changes, the occurrence of trajectory changes may be due to the combined role played by the dispersion and dissipation of the nonlinear structure. Our theoretical study may be beneficial to understand the propagation and interaction of nonlinear electrostatic waves and may brings a possibility to develop the nonlinear theory of electron-acoustic waves in astrophysical plasma systems.« less

Shock wave polarizations and optical metrics in the Born and the Born–Infeld electrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Minz, Christoph, E-mail: christoph.minz@alumni.tu-berlin.de; Borzeszkowski, Horst-Heino von, E-mail: borzeszk@mailbox.tu-berlin.de; Chrobok, Thoralf, E-mail: tchrobok@mailbox.tu-berlin.de

We analyze the behavior of shock waves in nonlinear theories of electrodynamics. For this, by use of generalized Hadamard step functions of increasing order, the electromagnetic potential is developed in a series expansion near the shock wave front. This brings about a corresponding expansion of the respective electromagnetic field equations which allows for deriving relations that determine the jump coefficients in the expansion series of the potential. We compute the components of a suitable gauge-normalized version of the jump coefficients given for a prescribed tetrad compatible with the shock front foliation. The solution of the first-order jump relations shows that,more » in contrast to linear Maxwell’s electrodynamics, in general the propagation of shock waves in nonlinear theories is governed by optical metrics and polarization conditions describing the propagation of two differently polarized waves (leading to a possible appearance of birefringence). In detail, shock waves are analyzed in the Born and Born–Infeld theories verifying that the Born–Infeld model exhibits no birefringence and the Born model does. The obtained results are compared to those ones found in literature. New results for the polarization of the two different waves are derived for Born-type electrodynamics.« less

Properties of a Laser Shock Wave in Al-Cu Alloy under Elevated Temperatures: A Molecular Dynamics Simulation Study

PubMed Central

Meng, Xiankai; Zhou, Jianzhong; Huang, Shu; Su, Chun; Sheng, Jie

2017-01-01

The laser shock wave (LSW) generated by the interaction between a laser and a material has been widely used in laser manufacturing, such as laser shock peening and laser shock forming. However, due to the high strain rate, the propagation of LSW in materials, especially LSW at elevated temperatures, is difficult to study through experimental methods. A molecular dynamics simulation was used in this study to investigate the propagation of LSW in an Al-Cu alloy. The Hugoniot relations of LSW were obtained at different temperatures and the effects of elevated temperatures on shock velocity and shock pressure were analyzed. Then the elastic and plastic wave of the LSW was researched. Finally, the evolution of dislocations induced by LSW and its mechanism under elevated temperatures was explored. The results indicate that the shock velocity and shock pressure induced by LSW both decrease with the increasing temperatures. Moreover, the velocity of elastic wave and plastic wave both decrease with the increasing treatment temperature, while their difference decreases as the temperature increases. Moreover, the dislocation atoms increases with the increasing temperatures before 2 ps, while it decreases with the increasing temperatures after 2 ps. The reason for the results is related to the formation and evolution of extended dislocations. PMID:28772433

Properties of a Laser Shock Wave in Al-Cu Alloy under Elevated Temperatures: A Molecular Dynamics Simulation Study.

PubMed

Meng, Xiankai; Zhou, Jianzhong; Huang, Shu; Su, Chun; Sheng, Jie

2017-01-18

The laser shock wave (LSW) generated by the interaction between a laser and a material has been widely used in laser manufacturing, such as laser shock peening and laser shock forming. However, due to the high strain rate, the propagation of LSW in materials, especially LSW at elevated temperatures, is difficult to study through experimental methods. A molecular dynamics simulation was used in this study to investigate the propagation of LSW in an Al-Cu alloy. The Hugoniot relations of LSW were obtained at different temperatures and the effects of elevated temperatures on shock velocity and shock pressure were analyzed. Then the elastic and plastic wave of the LSW was researched. Finally, the evolution of dislocations induced by LSW and its mechanism under elevated temperatures was explored. The results indicate that the shock velocity and shock pressure induced by LSW both decrease with the increasing temperatures. Moreover, the velocity of elastic wave and plastic wave both decrease with the increasing treatment temperature, while their difference decreases as the temperature increases. Moreover, the dislocation atoms increases with the increasing temperatures before 2 ps, while it decreases with the increasing temperatures after 2 ps. The reason for the results is related to the formation and evolution of extended dislocations.

Role of helmet in the mechanics of shock wave propagation under blast loading conditions.

PubMed

Ganpule, S; Gu, L; Alai, A; Chandra, N

2012-01-01

The effectiveness of helmets in extenuating the primary shock waves generated by the explosions of improvised explosive devices is not clearly understood. In this work, the role of helmet on the overpressurisation and impulse experienced by the head were examined. The shock wave-head interactions were studied under three different cases: (i) unprotected head, (ii) head with helmet but with varying head-helmet gaps and (iii) head covered with helmet and tightly fitting foam pads. The intensification effect was discussed by examining the shock wave flow pattern and verified with experiments. A helmet with a better protection against shock wave is suggested.

Biodamage via shock waves initiated by irradiation with ions.

PubMed

Surdutovich, Eugene; Yakubovich, Alexander V; Solov'yov, Andrey V

2013-01-01

Radiation damage following the ionising radiation of tissue has different scenarios and mechanisms depending on the projectiles or radiation modality. We investigate the radiation damage effects due to shock waves produced by ions. We analyse the strength of the shock wave capable of directly producing DNA strand breaks and, depending on the ion's linear energy transfer, estimate the radius from the ion's path, within which DNA damage by the shock wave mechanism is dominant. At much smaller values of linear energy transfer, the shock waves turn out to be instrumental in propagating reactive species formed close to the ion's path to large distances, successfully competing with diffusion.

Photoacoustic shock wave emission and cavitation from structured optical fiber tips

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mohammadzadeh, M.; Gonzalez-Avila, S. R.; Ohl, C. D., E-mail: cdohl@ntu.edu.sg

Photoacoustic waves generated at the tip of an optical fiber consist of a compressive shock wave followed by tensile diffraction waves. These tensile waves overlap along the fiber axis and form a cloud of cavitation bubbles. We demonstrate that shaping the fiber tip through micromachining alters the number and direction of the emitted waves and cavitation clouds. Shock wave emission and cavitation patterns from five distinctively shaped fiber tips have been studied experimentally and compared to a linear wave propagation model. In particular, multiple shock wave emission and generation of strong tension away from the fiber axis have been realizedmore » using modified fiber tips. These altered waveforms may be applied for novel microsurgery protocols, such as fiber-based histotripsy, by utilizing bubble-shock wave interaction.« less

Propagation of Electron Acoustic Soliton, Periodic and Shock Waves in Dissipative Plasma with a q-Nonextensive Electron Velocity Distribution

NASA Astrophysics Data System (ADS)

El-Hanbaly, A. M.; El-Shewy, E. K.; Elgarayhi, A.; Kassem, A. I.

2015-11-01

The nonlinear properties of small amplitude electron-acoustic (EA) solitary and shock waves in a homogeneous system of unmagnetized collisionless plasma with nonextensive distribution for hot electrons have been investigated. A reductive perturbation method used to obtain the Kadomstev-Petviashvili-Burgers equation. Bifurcation analysis has been discussed for non-dissipative system in the absence of Burgers term and reveals different classes of the traveling wave solutions. The obtained solutions are related to periodic and soliton waves and their behavior are shown graphically. In the presence of the Burgers term, the EXP-function method is used to solve the Kadomstev-Petviashvili-Burgers equation and the obtained solution is related to shock wave. The obtained results may be helpful in better conception of waves propagation in various space plasma environments as well as in inertial confinement fusion laboratory plasmas.

Fourth-power law structure of the shock wave fronts in metals and ceramics

NASA Astrophysics Data System (ADS)

Bayandin, Yuriy; Naimark, Oleg; Saveleva, Natalia

2017-06-01

The plate impact experiments were performed for solids during last fifty years. It was established that the dependence between the strain rate and the shock wave amplitude for metals and ceramics expressed by a fourth-power law. Present study is focused on the theoretical investigation and numerical simulation of plane shock wave propagation in metals and ceramics. Statistically based constitutive model of solid with defects (microcracks and microshears) was developed to provide the relation between damage induced mechanisms of structural relaxation, thermally activated plastic flow and material reactions for extreme loading conditions. Original approach based on the wide range constitutive equations was proposed for the numerical simulation of multiscale damage-failure transition mechanisms and plane shock wave propagation in solids with defects in the range of strain rate 103 -108s-1 . It was shown that mechanisms of plastic relaxation and damage-failure transitions are linked to the multiscale kinetics of defects leading to the self-similar nature of shock wave fronts in metals and ceramics. The work was supported by the Russian Science Foundation (Project No. 14-19-01173).

Shock waves in aviation security and safety

NASA Astrophysics Data System (ADS)

Settles, G. S.; Keane, B. T.; Anderson, B. W.; Gatto, J. A.

Accident investigations such as of Pan Am 103 and TWA 800 reveal the key role of shock-wave propagation in destroying the aircraft when an on-board explosion occurs. This paper surveys shock wave propagation inside an aircraft fuselage, caused either by a terrorist device or by accident, and provides some new experimental results. While aircraft-hardening research has been under way for more than a decade, no such experiments to date have used the crucial tool of high-speed optical imaging to visualize shock motion. Here, Penn State's Full-Scale Schlieren flow visualization facility yields the first shock-motion images in aviation security scenarios: 1) Explosions beneath full-size aircraft seats occupied by mannequins, 2) Explosions inside partially-filled luggage containers, and 3) Luggage-container explosions resulting in hull-holing. Both single-frame images and drum-camera movies are obtained. The implications of these results are discussed, though the overall topic must still be considered in its infancy.

Nonplanar dust-ion acoustic shock waves with transverse perturbation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xue Jukui

2005-01-01

The nonlinear dust-ion acoustic shock waves in dusty plasmas with the combined effects of bounded cylindrical/spherical geometry, the transverse perturbation, the dust charge fluctuation, and the nonthermal electrons are studied. Using the perturbation method, a cylindrical/spherical Kadomtsev-Petviashvili Burgers equation that describes the dust-ion acoustic shock waves is deduced. A particular solution of the cylindrical/spherical Kadomtsev-Petviashvili Burgers equation is also obtained. It is shown that the dust-ion acoustic shock wave propagating in cylindrical/spherical geometry with transverse perturbation will be slightly deformed as time goes on.

Experimental Study of the Shock Waves Produced by Condenser Discharge in a Gas Tube (thesis); ETUDE EXPERIMENTALE DES ONDES DE CHOC PRODUITES PAR DECHARGES D'UN CONDENSATEUR DANS UN TUBE A GAZ (thesis)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Der Agobian, R.

1964-10-31

The shock waves Droduced by condenser discharge in a gas tube were investigated. The study was limited to wave velocities less than five times the speed of sound, propagated in gas at low pressure (several mm Hg). A method was designed and perfected for the detection of the shock waves that are insufficiently rapid to produce gas ionization. This method consisted of the creation of an autonomous plasma, before the arrival of the wave, which was then modified by the wave passage. two methods were used for the detection of phenomena accompanying the passage of the shock waves, an opticalmore » method and a radioelectric method. The qualitative study of the modifications produced on the wave passage showed the remarkable correlation existing between the results obtained by the two methods. The experimental results on the propagation laws for shock waves in a low-diameter tube agreed with theory. The variations of the coefficient oi recombination were determined as a iunction of the electron temperature, and the results were in good agreement with the Bates theory. It was shown that the electron gas of the plasma had the same increase of density as a neutral gas during the passage of a shock wave. The variations of the frequency of electron collisions on passage of the shock wave could be explained by considering the electron--ion collisions with respect to electron-- atom collisions. (J.S.R.)« less

Drift ion acoustic shock waves in an inhomogeneous two-dimensional quantum magnetoplasma

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masood, W.; Siddiq, M.; Karim, S.

2009-04-15

Linear and nonlinear propagation characteristics of drift ion acoustic waves are investigated in an inhomogeneous quantum plasma with neutrals in the background employing the quantum hydrodynamics (QHD) model. In this regard, a quantum Kadomtsev-Petviashvili-Burgers (KPB) equation is derived for the first time. It is shown that the ion acoustic wave couples with the drift wave if the parallel motion of ions is taken into account. Discrepancies in the earlier works on drift solitons and shocks in inhomogeneous plasmas are also pointed out and a correct theoretical framework is presented to study the one-dimensional as well as the two-dimensional propagation ofmore » shock waves in an inhomogeneous quantum plasma. Furthermore, the solution of KPB equation is presented using the tangent hyperbolic (tanh) method. The variation of the shock profile with the quantum Bohm potential, collision frequency, and ratio of drift to shock velocity in the comoving frame, v{sub *}/u, are also investigated. It is found that increasing the number density and collision frequency enhances the strength of the shock. It is also shown that the fast drift shock (i.e., v{sub *}/u>0) increases, whereas the slow drift shock (i.e., v{sub *}/u<0) decreases the strength of the shock. The relevance of the present investigation with regard to dense astrophysical environments is also pointed out.« less

Modeling the Propagation of Shock Waves in Metals

NASA Astrophysics Data System (ADS)

Howard, W. Michael

2005-07-01

We present modeling results for the propagation of strong shock waves in metals. In particular, we use an arbitrary Lagrange Eulerian (ALE3D) code to model the propagation of strong pressure waves (P ˜300 to 400 kbars) generated with high explosives in contact with aluminum cylinders. The aluminum cylinders are assumed to be both flat-topped and have large-amplitude curved surfaces. We use 3D Lagrange mechanics. For the aluminum we use a rate-independent Steinberg-Guinan model, where the yield strength and bulk modulus depends on pressure, density and temperature. The calculation of the melt temperature is based on the Lindermann law. At melt the yield strength and bulk modulus is set to zero. The pressure is represented as a seven-term polynomial as a function of density. For the HMX-based high explosive, we use a JWL, with a program burn model that gives the correct detonation velocity and C-J pressure (P ˜ 390 kbars). For the case of the large-amplitude curved surface, we discuss the evolving shock structure in terms of the early shock propagation experiments by Sakharov. We also discuss the dependence of our results upon our material model for aluminum.

The structure of steady shock waves in porous metals

NASA Astrophysics Data System (ADS)

Czarnota, Christophe; Molinari, Alain; Mercier, Sébastien

2017-10-01

The paper aims at developing an understanding of steady shock wave propagation in a ductile metallic material containing voids. Porosity is assumed to be less than 0.3 and voids are not connected (foams are not considered). As the shock wave is traveling in the porous medium, the voids are facing a rapid collapse. During this dynamic compaction process, material particles are subjected to very high acceleration in the vicinity of voids, thus generating acceleration forces at the microscale that influence the overall response of the porous material. Analyzing how stationary shocks are influenced by these micro-inertia effects is the main goal of this work. The focus is essentially on the shock structure, ignoring oscillatory motion of pores prevailing at the tail of the shock wave. Following the constitutive framework developed by Molinari and Ravichandran (2004) for the analysis of steady shock waves in dense metals, an analytical approach of steady state propagation of plastic shocks in porous metals is proposed. The initial void size appears as a characteristic internal length that scales the overall dynamic response, thereby contributing to the structuring of the shock front. This key feature is not captured by standard damage models where the porosity stands for the single damage parameter with no contribution of the void size. The results obtained in this work provide a new insight in the fundamental understanding of shock waves in porous media. In particular, a new scaling law relating the shock width to the initial void radius is obtained when micro-inertia effects are significant.

Ionospheric response to the shock and acoustic waves excited by the launch of the Shenzhou 10 spacecraft

NASA Astrophysics Data System (ADS)

Ding, Feng; Wan, Weixing; Mao, Tian; Wang, Min; Ning, Baiqi; Zhao, Biqiang; Xiong, Bo

2014-05-01

We used a dense GPS network in China to track the ionospheric response to waves excited by the launch of the rocket that carried Shenzhou 10 spacecraft on 11 June 2013. The long-distance propagation of shock and acoustic waves were observed on both sides of the rocket's trajectory. On the southern side, the wave structures (characterized by a horizontal extension of ~1400 km and initial amplitudes of 0.3 total electron content unit (TECU) and 0.1 TECU for the shock and acoustic waves, respectively), traveled southwestward a distance of ~1500 km at mean velocities of 1011 m s-1 and 709 m s-1, respectively. On the northern side, northward propagating waves were seen to travel a distance of ~600 km with much smaller amplitudes of less than 0.05 TECU. Subsequent waves with amplitudes of less than 0.02 TECU could also be seen. Clear wave structures were found at a distance of ~600-2000 km from launch site.

Shock-activated electrochemical power supplies

DOEpatents

Benedick, William B.; Graham, Robert A.; Morosin, Bruno

1988-01-01

A shock-activated electrochemical power supply is provided which is initiated extremely rapidly and which has a long shelf life. Electrochemical power supplies of this invention are initiated much faster than conventional thermal batteries. Power supplies of this invention comprise an inactive electrolyte and means for generating a high-pressure shock wave such that the shock wave is propagated through the electrolytes rendering the electrolyte electrochemically active.

Shock-activated electrochemical power supplies

DOEpatents

Benedick, W.B.; Graham, R.A.; Morosin, B.

1988-11-08

A shock-activated electrochemical power supply is provided which is initiated extremely rapidly and which has a long shelf life. Electrochemical power supplies of this invention are initiated much faster than conventional thermal batteries. Power supplies of this invention comprise an inactive electrolyte and means for generating a high-pressure shock wave such that the shock wave is propagated through the electrolytes rendering the electrolyte electrochemically active. 2 figs.

High-Speed Photography of Detonation Propagation in Dynamically Precompressed Liquid Explosives

NASA Astrophysics Data System (ADS)

Petel, Oren; Higgins, Andrew; Yoshinaka, Akio; Zhang, Fan

2007-06-01

The propagation of detonation in shock compressed nitromethane was observed with a high speed framing camera. The test explosive, nitromethane, was compressed by a reverberating shock wave to pressures on the order of 10 GPa prior to being detonated by a secondary detonation event. The pressure and density in the test explosive prior to detonation was determined using two methods: manganin strain gauge measurements and LS-DYNA simulations. The velocity of the detonation front was determined from consecutive frames and correlated to the density of the explosive post-reverberating shock wave and prior to being detonated. Observing detonation propagation under these non-ambient conditions provides data which can be useful in the validation of equation of state models.

An in vivo study of low back pain and shock absorption in the human locomotor system.

PubMed

Voloshin, A; Wosk, J

1982-01-01

In this second of three papers, the principles of a non-invasive in vivo method to quantitatively evaluate the shock absorbing capacity of the human musculoskeletal system and the correlation of this shock absorbing capacity with low back pain (LPB) symptoms are presented. The experiments involved patients suffering from low back pain (as well as other degenerative joint diseases) and healthy patients. The obtained results reveal that low back pain correlates with the reduced capacity of the human musculoskeletal system between the femoral condyle and the forehead to attenuate incoming shock waves. Examination of the absolute values of the amplitude of the propagated waves leads to the conclusion that the human locomotor system, which possesses reduced attenuation capacity, tries to prevent overloading of the head from insufficiently attenuated shock waves. Results of the present investigation support the idea that the repetitive loading resulting from gait generates intermittent waves that propagate through the entire human musculoskeletal system from the heel up to the head. These waves are gradually attenuated along this course by the natural shock absorbers (bone and soft tissues). Contemporary methods for examination of the human musculoskeletal system may by improved by using the proposed non-invasive in vivo technique for quantitative characterization of the locomotor system's shock absorbing capacity.

Double shock front formation in cylindrical radiative blast waves produced by laser irradiation of krypton gas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kim, I.; Quevedo, H. J.; Feldman, S.

2013-12-15

Radiative blast waves were created by irradiating a krypton cluster source from a supersonic jet with a high intensity femtosecond laser pulse. It was found that the radiation from the shock surface is absorbed in the optically thick upstream medium creating a radiative heat wave that travels supersonically ahead of the main shock. As the blast wave propagates into the heated medium, it slows and loses energy, and the radiative heat wave also slows down. When the radiative heat wave slows down to the transonic regime, a secondary shock in the ionization precursor is produced. This paper presents experimental datamore » characterizing both the initial and secondary shocks and numerical simulations to analyze the double-shock dynamics.« less

Spherical shock-wave propagation in three-dimensional granular packings.

PubMed

Xue, Kun; Bai, Chun-Hua

2011-02-01

We investigate numerically the spherical shock-wave propagation in an open dense granular packing perturbed by the sudden expansion of a spherical intruder in the interior of the pack, focusing on the correlation between geometrical fabrics and propagating properties. The measurements of the temporal and spatial variations in a variety of propagating properties define a consistent serrated wave substructure with characteristic length on the orders of particle diameters. Further inspection of particle packing reveals a well-defined particle layering that persists several particle diameters away from the intruder, although its dominant effects are only within one to two diameters. This interface-induced layering not only exactly coincides with the serrated wave profile, but also highlights the competition between two energy transmission mechanisms involving distinct transport speeds. The alternating dominances between these two mechanisms contribute to the nonlinear wave propagation on the particle scale. Moreover, the proliferation of intricate three-dimensional contact force networks suggests the anisotropic stress transmission, which is found to also arise from the localized packing structure in the vicinity of the intruder.

Ultrafast dynamic response of single-crystal β-HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)

NASA Astrophysics Data System (ADS)

Zaug, Joseph M.; Austin, Ryan A.; Armstrong, Michael R.; Crowhurst, Jonathan C.; Goldman, Nir; Ferranti, Louis; Saw, Cheng K.; Swan, Raymond A.; Gross, Richard; Fried, Laurence E.

2018-05-01

We report experimental and computational studies of shock wave dynamics in single-crystal β-HMX on an ultrafast time scale. Here, a laser-based compression drive (˜1 ns in duration; stresses of up to ˜40 GPa) is used to propagate shock waves normal to the (110) and (010) lattice planes. Ultrafast time-domain interferometry measurements reveal distinct, time-dependent relationships between the shock wave velocity and particle velocity for each crystal orientation, which suggest evolving physical processes on a sub-nanosecond time scale. To help interpret the experimental data, elastic shock wave response was simulated using a finite-strain model of crystal thermoelasticity. At early propagation times (<500 ps), the model is in agreement with the data, which indicates that the mechanical response is dominated by thermoelastic deformation. The model agreement depends on the inclusion of nonlinear elastic effects in both the spherical and deviatoric stress-strain responses. This is achieved by employing an equation-of-state and a pressure-dependent stiffness tensor, which was computed via atomistic simulation. At later times (>500 ps), the crystal samples exhibit signatures of inelastic deformation, structural phase transformation, or chemical reaction, depending on the direction of wave propagation.

Stationary Shock Waves with Oscillating Front in Dislocation Systems of Semiconductors

NASA Astrophysics Data System (ADS)

Gestrin, S. G.; Shchukina, E. V.

2018-05-01

The paper presents a study of weakly nonlinear wave processes in the cylindrical region of a hole gas surrounding a negatively charged dislocation in an n-type semiconductor crystal. It is shown that shock waves propagating along the dislocation are the solutions of the Korteweg-de Vries-Burgers equation when the dispersion and dissipation of medium are taken into account. Estimates are obtained for the basic physical parameters characterizing the shock wave and the region inside the Reed cylinder.

Demonstration of Dispersive Rarefaction Shocks in Hollow Elliptical Cylinder Chains

NASA Astrophysics Data System (ADS)

Kim, H.; Kim, E.; Chong, C.; Kevrekidis, P. G.; Yang, J.

2018-05-01

We report an experimental and numerical demonstration of dispersive rarefaction shocks (DRS) in a 3D-printed soft chain of hollow elliptical cylinders. We find that, in contrast to conventional nonlinear waves, these DRS have their lower amplitude components travel faster, while the higher amplitude ones propagate slower. This results in the backward-tilted shape of the front of the wave (the rarefaction segment) and the breakage of wave tails into a modulated waveform (the dispersive shock segment). Examining the DRS under various impact conditions, we find the counterintuitive feature that the higher striker velocity causes the slower propagation of the DRS. These unique features can be useful for mitigating impact controllably and efficiently without relying on material damping or plasticity effects.

Demonstration of Dispersive Rarefaction Shocks in Hollow Elliptical Cylinder Chains.

PubMed

Kim, H; Kim, E; Chong, C; Kevrekidis, P G; Yang, J

2018-05-11

We report an experimental and numerical demonstration of dispersive rarefaction shocks (DRS) in a 3D-printed soft chain of hollow elliptical cylinders. We find that, in contrast to conventional nonlinear waves, these DRS have their lower amplitude components travel faster, while the higher amplitude ones propagate slower. This results in the backward-tilted shape of the front of the wave (the rarefaction segment) and the breakage of wave tails into a modulated waveform (the dispersive shock segment). Examining the DRS under various impact conditions, we find the counterintuitive feature that the higher striker velocity causes the slower propagation of the DRS. These unique features can be useful for mitigating impact controllably and efficiently without relying on material damping or plasticity effects.

Parabolic equation for nonlinear acoustic wave propagation in inhomogeneous moving media

NASA Astrophysics Data System (ADS)

Aver'yanov, M. V.; Khokhlova, V. A.; Sapozhnikov, O. A.; Blanc-Benon, Ph.; Cleveland, R. O.

2006-12-01

A new parabolic equation is derived to describe the propagation of nonlinear sound waves in inhomogeneous moving media. The equation accounts for diffraction, nonlinearity, absorption, scalar inhomogeneities (density and sound speed), and vectorial inhomogeneities (flow). A numerical algorithm employed earlier to solve the KZK equation is adapted to this more general case. A two-dimensional version of the algorithm is used to investigate the propagation of nonlinear periodic waves in media with random inhomogeneities. For the case of scalar inhomogeneities, including the case of a flow parallel to the wave propagation direction, a complex acoustic field structure with multiple caustics is obtained. Inclusion of the transverse component of vectorial random inhomogeneities has little effect on the acoustic field. However, when a uniform transverse flow is present, the field structure is shifted without changing its morphology. The impact of nonlinearity is twofold: it produces strong shock waves in focal regions, while, outside the caustics, it produces higher harmonics without any shocks. When the intensity is averaged across the beam propagating through a random medium, it evolves similarly to the intensity of a plane nonlinear wave, indicating that the transverse redistribution of acoustic energy gives no considerable contribution to nonlinear absorption.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mehta, Y.; Neal, C.; Salari, K.

Propagation of a strong shock through a bed of particles results in complex wave dynamics such as a reflected shock, a transmitted shock, and highly unsteady flow inside the particle bed. In this paper we present three-dimensional numerical simulations of shock propagation in air over a random bed of particles. We assume the flow is inviscid and governed by the Euler equations of gas dynamics. Simulations are carried out by varying the volume fraction of the particle bed at a fixed shock Mach number. We compute the unsteady inviscid streamwise and transverse drag coefficients as a function of time formore » each particle in the random bed as a function of volume fraction. We show that (i) there are significant variations in the peak drag for the particles in the bed, (ii) the mean peak drag as a function of streamwise distance through the bed decreases with a slope that increases as the volume fraction increases, and (iii) the deviation from the mean peak drag does not correlate with local volume fraction. We also present the local Mach number and pressure contours for the different volume fractions to explain the various observed complex physical mechanisms occurring during the shock-particle interactions. Since the shock interaction with the random bed of particles leads to transmitted and reflected waves, we compute the average flow properties to characterize the strength of the transmitted and reflected shock waves and quantify the energy dissipation inside the particle bed. Finally, to better understand the complex wave dynamics in a random bed, we consider a simpler approximation of a planar shock propagating in a duct with a sudden area change. We obtain Riemann solutions to this problem, which are used to compare with fully resolved numerical simulations.« less

Modeling of High-Velocity Flows in ITAM Impulse Facilities

DTIC Science & Technology

2010-04-01

up to 150 ms; Adiabatic compression wind tunnels up to 100 ms; Shock tubes... shock tubes. Basic and applied aerodynamic research has been performed in these wind tunnels in the range of Mach numbers М = 6 20 for many years...passage of a shock wave propagating over a cold rarefied gas filling the wind tunnel . When the gas heated in the shock wave (plug) passes around the

Potential applications of low-energy shock waves in functional urology.

PubMed

Wang, Hung-Jen; Cheng, Jai-Hong; Chuang, Yao-Chi

2017-08-01

A shock wave, which carries energy and can propagate through a medium, is a type of continuous transmitted sonic wave with a frequency of 16 Hz-20 MHz. It is accompanied by processes involving rapid energy transformations. The energy associated with shock waves has been harnessed and used for various applications in medical science. High-energy extracorporeal shock wave therapy is the most successful application of shock waves, and has been used to disintegrate urolithiasis for 30 years. At lower energy levels, however, shock waves have enhanced expression of vascular endothelial growth factor, endothelial nitric oxide synthase, proliferating cell nuclear antigen, chemoattractant factors and recruitment of progenitor cells; shock waves have also improved tissue regeneration. Low-energy shock wave therapy has been used clinically with musculoskeletal disorders, ischemic cardiovascular disorders and erectile dysfunction, through the mechanisms of neovascularization, anti-inflammation and tissue regeneration. Furthermore, low-energy shock waves have been proposed to temporarily increase tissue permeability and facilitate intravesical drug delivery. The present review article provides information on the basics of shock wave physics, mechanisms of action on the biological system and potential applications in functional urology. © 2017 The Japanese Urological Association.

Detecting Fragmentation of Kidney Stones in Lithotripsy by Means of Shock Wave Scattering

NASA Astrophysics Data System (ADS)

Sapozhnikov, Oleg A.; Trusov, Leonid A.; Owen, Neil R.; Bailey, Michael R.; Cleveland, Robin O.

2006-05-01

Although extracorporeal shock wave lithotripsy (a procedure of kidney stone comminution using focused shock waves) has been used clinically for many years, a proper monitoring of the stone fragmentation is still undeveloped. A method considered here is based on recording shock wave scattering signals with a focused receiver placed far from the stone, outside the patient body. When a fracture occurs in the stone or the stone becomes smaller, the elastic waves in the stone will propagate differently (e.g. shear waves will not cross a fracture) which, in turn, will change the scattered acoustic wave in the surrounding medium. Theoretical studies of the scattering phenomenon are based on a linear elastic model to predict shock wave scattering by a stone, with and without crack present in it. The elastic waves in the stone and the nearby liquid were modeled using a finite difference time domain approach. The subsequent acoustic propagation of the scattered waves into the far-field was calculated using the Helmholtz-Kirchhoff integral. Experimental studies were conducted using a research electrohydraulic lithotripter that produced the same acoustic output as an unmodified Dornier HM3 clinical lithotripter. Artificial stones, made from Ultracal-30 gypsum and acrylic, were used as targets. The stones had cylindrical shape and were positioned co-axially with the lithotripter axis. The scattered wave was measured by focused broadband PVDF hydrophone. It was shown that the size of the stone noticeably changed the signature of the reflected wave.

Acceleration Wave Propagation in Hyperelastic Rods of Variable Cross-Section.

DTIC Science & Technology

1981-07-01

direction of propagation. Many authors have considered both static and dynamic problems for such materials, of whom we mention only Antman [2] and... Antman and Jordan [3] who studied the Kirchhoff problem for nonlinearly elastic rods and qualitative properties in general, Jeffrey and Teymur [4] and...Jeffrey and Suhubi [5] who considered shock wave formation and acceleration wave propagation through periodically layered media, and Antman and Liu [6

Shock wave attenuation by grids and orifice plates

NASA Astrophysics Data System (ADS)

Britan, A.; Igra, O.; Ben-Dor, G.; Shapiro, H.

2006-11-01

The interaction of weak shock waves with porous barriers of different geometries and porosities is examined. Installing a barrier inside the shock tube test section will cause the development of the following wave pattern upon a head-on collision between the incident shock wave and the barrier: a reflected shock from the barrier and a transmitted shock propagating towards the shock tube end wall. Once the transmitted shock wave reaches the end wall it is reflected back towards the barrier. This is the beginning of multiple reflections between the barrier and the end wall. This full cycle of shock reflections/interactions resulting from the incident shock wave collision with the barrier can be studied in a single shock tube test. A one-dimensional (1D), inviscid flow model was proposed for simulating the flow resulting from the initial collision of the incident shock wave with the barrier. Fairly good agreement is found between experimental findings and simulations based on a 1D flow model. Based on obtained numerical and experimental findings an optimal design procedure for shock wave attenuator is suggested. The suggested attenuator may ensure the safety of the shelter’s ventilation systems.

Modeling of an electrohydraulic lithotripter with the KZK equation.

PubMed

Averkiou, M A; Cleveland, R O

1999-07-01

The acoustic pressure field of an electrohydraulic extracorporeal shock wave lithotripter is modeled with a nonlinear parabolic wave equation (the KZK equation). The model accounts for diffraction, nonlinearity, and thermoviscous absorption. A numerical algorithm for solving the KZK equation in the time domain is used to model sound propagation from the mouth of the ellipsoidal reflector of the lithotripter. Propagation within the reflector is modeled with geometrical acoustics. It is shown that nonlinear distortion within the ellipsoidal reflector can play an important role for certain parameters. Calculated waveforms are compared with waveforms measured in a clinical lithotripter and good agreement is found. It is shown that the spatial location of the maximum negative pressure occurs pre-focally which suggests that the strongest cavitation activity will also be in front of the focus. Propagation of shock waves from a lithotripter with a pressure release reflector is considered and because of nonlinear propagation the focal waveform is not the inverse of the rigid reflector. Results from propagation through tissue are presented; waveforms are similar to those predicted in water except that the higher absorption in the tissue decreases the peak amplitude and lengthens the rise time of the shock.

Experimental and numerical investigations of shock wave propagation through a bifurcation

NASA Astrophysics Data System (ADS)

Marty, A.; Daniel, E.; Massoni, J.; Biamino, L.; Houas, L.; Leriche, D.; Jourdan, G.

2018-02-01

The propagation of a planar shock wave through a split channel is both experimentally and numerically studied. Experiments were conducted in a square cross-sectional shock tube having a main channel which splits into two symmetric secondary channels, for three different shock wave Mach numbers ranging from about 1.1 to 1.7. High-speed schlieren visualizations were used along with pressure measurements to analyze the main physical mechanisms that govern shock wave diffraction. It is shown that the flow behind the transmitted shock wave through the bifurcation resulted in a highly two-dimensional unsteady and non-uniform flow accompanied with significant pressure loss. In parallel, numerical simulations based on the solution of the Euler equations with a second-order Godunov scheme confirmed the experimental results with good agreement. Finally, a parametric study was carried out using numerical analysis where the angular displacement of the two channels that define the bifurcation was changed from 90° , 45° , 20° , and 0° . We found that the angular displacement does not significantly affect the overpressure experience in either of the two channels and that the area of the expansion region is the important variable affecting overpressure, the effect being, in the present case, a decrease of almost one half.

Normal shock wave reflection on porous compressible material

NASA Astrophysics Data System (ADS)

Gvozdeva, L. G.; Faresov, Iu. M.; Brossard, J.; Charpentier, N.

The present experimental investigation of the interaction of plane shock waves in air and a rigid wall coated with flat layers of expanded polymers was conducted in a standard shock tube and a diaphragm with an initial test section pressure of 100,000 Pa. The Mach number of the incident shock wave was varied from 1.1 to 2.7; the peak pressures measured on the wall behind polyurethane at various incident wave Mach numbers are compared with calculated values, with the ideal model of propagation, and with the reflection of shock waves in a porous material that is understood as a homogeneous mixture. The effect of elasticity and permeability of the porous material structure on the rigid wall's pressure pulse parameters is qualitatively studied.

Impact of the volume of rooms on shock wave propagation within a multi-chamber system

NASA Astrophysics Data System (ADS)

Julien, B.; Sochet, I.; Vaillant, T.

2016-03-01

The behavior of a shock wave generated by a hemispherical gaseous charge and propagating within a confined multi-chamber system is analyzed through the evolution of some of the shock parameters (maximum overpressure and positive impulse). The influence of a variation in the volume of the rooms on the pressure history inside the building is also studied. Several small-scale experiments have been carried out using an adjustable model representative of a pyrotechnic workshop. The experimental results show that the pressure histories are very complex. Yet, using a global approach, we were able to link the evolution of the arrival time of the shock wave within the building with the reference obtained in the free field. New parameters were developed to best fit the experimental maximal overpressure in the cells and in the corridor leading to two predictive laws used to estimate the maximal overpressure in the model.

Chromospheric heating by acoustic shock waves

NASA Technical Reports Server (NTRS)

Jordan, Stuart D.

1993-01-01

Work by Anderson & Athay (1989) suggests that the mechanical energy required to heat the quiet solar chromosphere might be due to the dissipation of weak acoustic shocks. The calculations reported here demonstrate that a simple picture of chromospheric shock heating by acoustic waves propagating upward through a model solar atmosphere, free of both magnetic fields and local inhomogeneities, cannot reproduce their chromospheric model. The primary reason is the tendency for vertically propagating acoustic waves in the range of allowed periods to dissipate too low in the atmosphere, providing insufficient residual energy for the middle chromosphere. The effect of diverging magnetic fields and the corresponding expanding acoustic wavefronts on the mechanical dissipation length is then discussed as a means of preserving a quasi-acoustic heating hypothesis. It is argued that this effect, in a canopy that overlies the low chromosphere, might preserve the acoustic shock hypothesis consistent with the chromospheric radiation losses computed by Anderson & Athay.

Explosively Driven Shock Induced Damage in OFHC Copper

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koller, D. D.; Hixson, R. S.; Gray, G. T. III

OFHC Cu samples were subjected to shock loading using plane wave HE lenses to produce a uniaxial Taylor wave profile (shock followed by immediate release). Upon arrival of the shock wave at the free surface of the sample, the wave is reflected and propagates back into the sample as a release wave. It is the interaction of initial and reflected release waves that place the material in a localized state of tension which can ultimately result in damage and possibly complete failure of the material. The peak tensile stress and its location in the material are determined by the wavemore » shape. Damage evolution processes and localized behavior are discussed based on results from time-resolved free surface velocity (VISAR) interferometry and post shock metallurgical analysis of the soft recovered samples.« less

Head-on collision of normal shock waves with rigid porous materials

NASA Astrophysics Data System (ADS)

Levy, A.; Ben-Dor, G.; Skews, B. W.; Sorek, S.

1993-08-01

The head-on collision of a planar shock wave with a rigid porous material has been investigated experimentally in a 75 mm × 75 mm shock tube. The experimental study indicated that unlike the reflection from a flexible porous material (e.g., polyurethane foam) where the transmitted compression waves do not converge to a sharp shock wave, in the case of a rigid porous material (e.g., alumina) the transmitted compression waves do converge to a sharp shock wave, which decays as it propagates along the porous material. In addition to this major difference, many other differences were observed. They are outlined in the following sections. Based on these observations a suggestion modifying the phenomenology of the reflection/interaction process in the case a porous material with large permeability is proposed.

The Simultaneous Combination of Phase Contrast Imaging with In Situ X-ray diffraction from Shock Compressed Matter

NASA Astrophysics Data System (ADS)

McBride, Emma Elizabeth; Seiboth, Frank; Cooper, Leora; Frost, Mungo; Goede, Sebastian; Harmand, Marion; Levitan, Abe; McGonegle, David; Miyanishi, Kohei; Ozaki, Norimasa; Roedel, Melanie; Sun, Peihao; Wark, Justin; Hastings, Jerry; Glenzer, Siegfried; Fletcher, Luke

2017-10-01

Here, we present the simultaneous combination of phase contrast imaging (PCI) techniques with in situ X-ray diffraction to investigate multiple-wave features in laser-driven shock-compressed germanium. Experiments were conducted at the Matter at Extreme Conditions end station at the LCLS, and measurements were made perpendicular to the shock propagation direction. PCI allows one to take femtosecond snapshots of magnified real-space images of shock waves as they progress though matter. X-ray diffraction perpendicular to the shock propagation direction provides the opportunity to isolate and identify different waves and determine the crystal structure unambiguously. Here, we combine these two powerful techniques simultaneously, by using the same Be lens setup to focus the fundamental beam at 8.2 keV to a size of 1.5 mm on target for PCI and the 3rd harmonic at 24.6 keV to a spot size of 2 um on target for diffraction.

Effect of Pressure Gradients on the Initiation of PBX-9502 via Irregular (Mach) Reflection of Low Pressure Curved Shock Waves

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hull, Lawrence Mark; Miller, Phillip Isaac; Moro, Erik Allan

In the instance of multiple fragment impact on cased explosive, isolated curved shocks are generated in the explosive. These curved shocks propagate and may interact and form irregular or Mach reflections along the interaction loci, thereby producing a single shock that may be sufficient to initiate PBX-9501. However, the incident shocks are divergent and their intensity generally decreases as they expand, and the regions behind the Mach stem interaction loci are generally unsupported and allow release waves to rapidly affect the flow. The effects of release waves and divergent shocks may be considered theoretically through a “Shock Change Equation”.

Nonlinear Korteweg-de Vries-Burger equation for ion acoustic shock waves in a weakly relativistic electron-positron-ion plasma with thermal ions

NASA Astrophysics Data System (ADS)

Saeed, R.; Shah, Asif

2010-03-01

The nonlinear propagation of ion acoustic waves in electron-positron-ion plasma comprising of Boltzmannian electrons, positrons, and relativistic thermal ions has been examined. The Korteweg-de Vries-Burger equation has been derived by reductive perturbation technique, and its shock like solution is determined analytically through tangent hyperbolic method. The effect of various plasma parameters on strength and structure of shock wave is investigated. The pert graphical view of the results has been presented for illustration. It is observed that strength and steepness of the shock wave enervate with an increase in the ion temperature, relativistic streaming factor, positron concentrations, electron temperature and they accrue with an increase in coefficient of kinematic viscosity. The convective, dispersive, and dissipative properties of the plasma are also discussed. It is determined that the electron temperature has remarkable influence on the propagation and structure of nonlinear wave in such relativistic plasmas. The numerical analysis has been done based on the typical numerical data from a pulsar magnetosphere.

Detonation Propagation through Nitromethane Embedded Metal Foam

NASA Astrophysics Data System (ADS)

Lieberthal, Brandon; Maines, Warren R.; Stewart, D. Scott

2015-11-01

There is considerable interest in developing a better understanding of dynamic behaviors of multicomponent systems. We report results of Eulerian hydrodynamic simulations of shock waves propagating through metal foam at approximately 20% relative density and various porosities using a reactive flow model in the ALE3D software package. We investigate the applied pressure and energy of the shock wave and its effects on the fluid and the inert material interface. By varying pore sizes, as well as metal impedance, we predict the overall effects of heterogeneous material systems at the mesoscale. In addition, we observe a radially expanding blast front in these heterogeneous models and apply the theory of Detonation Shock Dynamics to the convergence behavior of the lead shock.

Enhanced kidney stone fragmentation by short delay tandem conventional and modified lithotriptor shock waves: a numerical analysis.

PubMed

Tham, Leung-Mun; Lee, Heow Pueh; Lu, Chun

2007-07-01

We evaluated the effectiveness of modified lithotriptor shock waves using computer models. Finite element models were used to simulate the propagation of lithotriptor shock waves in human renal calculi in vivo. Kidney stones were assumed to be spherical, homogeneous, isotropic and linearly elastic, and immersed in a continuum fluid. Single and tandem shock wave pulses modified to intensify the collapse of cavitation bubbles near the stone surface to increase fragmentation efficiency and suppress the expansion of intraluminal bubbles for decreased vascular injury were analyzed. The effectiveness of the modified shock waves was assessed by comparing the states of loading in the renal calculi induced by these shock waves to those produced by conventional shock waves. Our numerical simulations revealed that modified shock waves produced marginally lower stresses in spherical renal calculi than those produced by conventional shock waves. Tandem pulses of conventional or modified shock waves produced peak stresses in the front and back halves of the renal calculi. However, the single shock wave pulses generated significant peak stresses in only the back halves of the renal calculi. Our numerical simulations suggest that for direct stress wave induced fragmentation modified shock waves should be as effective as conventional shock waves for fragmenting kidney stones. Also, with a small interval of 20 microseconds between the pulses tandem pulse lithotripsy using modified or conventional shock waves could be considerably more effective than single pulse lithotripsy for fragmenting kidney stones.

FLARE-GENERATED SHOCK WAVE PROPAGATION THROUGH SOLAR CORONAL ARCADE LOOPS AND AN ASSOCIATED TYPE II RADIO BURST

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kumar, Pankaj; Cho, Kyung-Suk; Innes, D. E., E-mail: pankaj@kasi.re.kr

2016-09-01

This paper presents multiwavelength observations of a flare-generated type II radio burst. The kinematics of the shock derived from the type II burst closely match a fast extreme ultraviolet (EUV) wave seen propagating through coronal arcade loops. The EUV wave was closely associated with an impulsive M1.0 flare without a related coronal mass ejection, and was triggered at one of the footpoints of the arcade loops in active region NOAA 12035. It was initially observed in the 335 Å images from the Atmospheric Image Assembly with a speed of ∼800 km s{sup −1} and it accelerated to ∼1490 km s{supmore » −1} after passing through the arcade loops. A fan–spine magnetic topology was revealed at the flare site. A small, confined filament eruption (∼340 km s{sup −1}) was also observed moving in the opposite direction to the EUV wave. We suggest that breakout reconnection in the fan–spine topology triggered the flare and associated EUV wave that propagated as a fast shock through the arcade loops.« less

Evolutions of elastic-plastic shock compression waves in different materials

NASA Astrophysics Data System (ADS)

Kanel, G. I.; Zaretsky, E. B.; Razorenov, S. V.; Savinykh, A. S.; Garkushin, G. V.

2017-01-01

In the paper, we discuss such unexpected features in the wave evolution in solids as a departure from self-similar development of the wave process which is accompanied with apparent sub-sonic wave propagation, changes of shape of elastic precursor wave as a result of variations in the material structure and the temperature, unexpected peculiarities of reflection of elastic-plastic waves from free surface, effects of internal friction at shock compression of glasses and some other effects.

A Search for Early High-Energy Afterglows in BATSE Gamma-Ray Bursts

NASA Technical Reports Server (NTRS)

Giblin, Timothy W.

2003-01-01

The scope of this project was to perform a detailed search for the early high-energy afterglow component of gamma-ray bursts (GRBs) in the BATSE GRB data archive. GRBs are believed to be the product of shock waves generated in a relativistic outflow from the demise of extremely massive stars and/or binary neutron star mergers. The outflow undeniably encounters the ambient medium of the progenitor object and another shock wave is set up. A forward shock propagates into the medium and a reverse shock propagates through the ejecta. This "external" shock dissipates the kinetic energy of the ejecta in the form of radiation via synchrotron losses and slows the outflow eventually to a non-relativistic state. Radiation from the forward external shock is therefore expected to be long-lived, lasting days, weeks, and even months. This radiation is referred to as the 'afterglow'.

Finite element modeling of acoustic wave propagation and energy deposition in bone during extracorporeal shock wave treatment

NASA Astrophysics Data System (ADS)

Wang, Xiaofeng; Matula, Thomas J.; Ma, Yong; Liu, Zheng; Tu, Juan; Guo, Xiasheng; Zhang, Dong

2013-06-01

It is well known that extracorporeal shock wave treatment is capable of providing a non-surgical and relatively pain free alternative treatment modality for patients suffering from musculoskeletal disorders but do not respond well to conservative treatments. The major objective of current work is to investigate how the shock wave (SW) field would change if a bony structure exists in the path of the acoustic wave. Here, a model of finite element method (FEM) was developed based on linear elasticity and acoustic propagation equations to examine SW propagation and deflection near a mimic musculoskeletal bone. High-speed photography experiments were performed to record cavitation bubbles generated in SW field with the presence of mimic bone. By comparing experimental and simulated results, the effectiveness of FEM model could be verified and strain energy distributions in the bone were also predicted according to numerical simulations. The results show that (1) the SW field will be deflected with the presence of bony structure and varying deflection angles can be observed as the bone shifted up in the z-direction relative to SW geometric focus (F2 focus); (2) SW deflection angels predicted by the FEM model agree well with experimental results obtained from high-speed photographs; and (3) temporal evolutions of strain energy distribution in the bone can also be evaluated based on FEM model, with varied vertical distance between F2 focus and intended target point on the bone surface. The present studies indicate that, by combining MRI/CT scans and FEM modeling work, it is possible to better understand SW propagation characteristics and energy deposition in musculoskeletal structure during extracorporeal shock wave treatment, which is important for standardizing the treatment dosage, optimizing treatment protocols, and even providing patient-specific treatment guidance in clinic.

Potential for the Vishniac instability in ionizing shock waves propagating into cold gases

NASA Astrophysics Data System (ADS)

Robinson, A. P. L.; Pasley, J.

2018-05-01

The Vishniac instability was posited as an instability that could affect supernova remnants in their late stage of evolution when subject to strong radiative cooling, which can drive the effective ratio of specific heats below 1.3. The potential importance of this instability to these astrophysical objects has motivated a number of laser-driven laboratory studies. However, the Vishniac instability is essentially a dynamical instability that should operate independently of whatever physical processes happen to reduce the ratio of specific heats. In this paper, we examine the possibility that ionization and molecular dissociation processes can achieve this, and we show that this is possible for a certain range of shock wave Mach numbers for ionizing/dissociating shock waves propagating into cold atomic and molecular gases.

Shock Wave Dynamics in Weakly Ionized Plasmas

NASA Technical Reports Server (NTRS)

Johnson, Joseph A., III

1999-01-01

An investigation of the dynamics of shock waves in weakly ionized argon plasmas has been performed using a pressure ruptured shock tube. The velocity of the shock is observed to increase when the shock traverses the plasma. The observed increases cannot be accounted for by thermal effects alone. Possible mechanisms that could explain the anomalous behavior include a vibrational/translational relaxation in the nonequilibrium plasma, electron diffusion across the shock front resulting from high electron mobility, and the propagation of ion-acoustic waves generated at the shock front. Using a turbulence model based on reduced kinetic theory, analysis of the observed results suggest a role for turbulence in anomalous shock dynamics in weakly ionized media and plasma-induced hypersonic drag reduction.

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

NASA Technical Reports Server (NTRS)

Bershader, D. (Editor); Hanson, R. (Editor)

1986-01-01

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

NASA Astrophysics Data System (ADS)

Bershader, D.; Hanson, R.

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Nonlinear Acoustic Propagation into the Seafloor

NASA Astrophysics Data System (ADS)

McDonald, B. Edward

2006-05-01

Explosions near the seafloor result in shock waves entering a much more complicated medium than water or air. Nonlinearities may be increased by two processes inherent to granular media: (1) a poroelastic nonlinearity comparable to the addition of bubbles to water, and (2) the Hertz force resulting from elastic deformation of grains, proportional to the Youngs modulus of the grains times the strain rate to the power 3/2. These two types of nonlinearity for shock propagation into the seafloor are investigated using a variant of the NPE model. The traditional Taylor series expansion of the equation of state (pressure as a function of density) is not appropriate to the Hertz force in the limit of small strain. We present a simple nonlinear wave equation model for compressional waves in marine sediments that retains the Hertz force explicitly with overdensity to the power 3/2. Numerical results for shock propagation are compared with similarity solutions for quadratic nonlinearity and for the fractional nonlinearity of the Hertz force.

Initial Results from the Variable Intensity Sonic Boom Propagation Database

NASA Technical Reports Server (NTRS)

Haering, Edward A., Jr.; Cliatt, Larry J., II; Bunce, Thomas J.; Gabrielson, Thomas B.; Sparrow, Victor W.; Locey, Lance L.

2008-01-01

An extensive sonic boom propagation database with low- to normal-intensity booms (overpressures of 0.08 lbf/sq ft to 2.20 lbf/sq ft) was collected for propagation code validation, and initial results and flight research techniques are presented. Several arrays of microphones were used, including a 10 m tall tower to measure shock wave directionality and the effect of height above ground on acoustic level. A sailplane was employed to measure sonic booms above and within the atmospheric turbulent boundary layer, and the sailplane was positioned to intercept the shock waves between the supersonic airplane and the ground sensors. Sailplane and ground-level sonic boom recordings were used to generate atmospheric turbulence filter functions showing excellent agreement with ground measurements. The sonic boom prediction software PCBoom4 was employed as a preflight planning tool using preflight weather data. The measured data of shock wave directionality, arrival time, and overpressure gave excellent agreement with the PCBoom4-calculated results using the measured aircraft and atmospheric data as inputs. C-weighted acoustic levels generally decreased with increasing height above the ground. A-weighted and perceived levels usually were at a minimum for a height where the elevated microphone pressure rise time history was the straightest, which is a result of incident and ground-reflected shock waves interacting.

Chaos and ion heating in a slow shock

NASA Technical Reports Server (NTRS)

Lin, Y.; Lee, L. C.

1991-01-01

An ion heating mechanism is proposed of slow shocks, which is associated with the chaotic motion of particles in the downstream wave field. For a coherent electromagnetic wave propagating along the downstream magnetic field, corresponding to switch-off shocks, the particle motions are not chaotic. For an oblique wave, the interaction between the particles and the wave field may lead to chaotic particle motions. Such particles may be greatly thermalized within one wavelength after they are incident into the downstream wave field. The results can be used to explain the existence of the critical intermediate Mach number observed in the hybrid simulations.

In vitro study of the mechanical effects of shock-wave lithotripsy.

PubMed

Howard, D; Sturtevant, B

1997-01-01

Impulsive stress in repeated shock waves administered during extracorporeal shock-wave lithotripsy (ESWL) causes injury to kidney tissue. In a study of the mechanical input of ESWL, the effects of focused shock waves on thin planar polymeric membranes immersed in a variety of tissue-mimicking fluids have been examined. A direct mechanism of failure by shock compression and an indirect mechanism by bubble collapse have been observed. Thin membranes are easily damaged by bubble collapse. After propagating through cavitation-free acoustically heterogeneous media (liquids mixed with hollow glass spheres, and tissue) shock waves cause membranes to fail in fatigue by a shearing mechanism. As is characteristic of dynamic fatigue, the failure stress increases with strain rate, determined by the amplitude and rise time of the attenuated shock wave. Shocks with large amplitude and short rise time (i.e., in uniform media) cause no damage. Thus the inhomogeneity of tissue is likely to contribute to injury in ESWL. A definition of dose is proposed which yields a criterion for damage based on measurable shock wave properties.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ro, Stephen; Matzner, Christopher D., E-mail: ro@astro.utoronto.ca

Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, luminous blue variable eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing on nonlinear wave propagation, shock formation, and an early phase of the development of a weak shock. We identify the shock formation radius, showing that a heuristic estimate based on crossing characteristics matches an exact expansion around the wave front and verifying both with numerical experiments. Our general analytical condition for shock formation applies to one-dimensionalmore » motions within any static environment, including both eruptions and implosions. We also consider the early phase of shock energy dissipation. We find that waves of super-Eddington acoustic luminosity always create shocks, rather than damping by radiative diffusion. Therefore, shock formation is integral to super-Eddington outbursts.« less

Laser Light Scattering by Shock Waves

NASA Technical Reports Server (NTRS)

Panda, J.; Adamovsky, G.

1995-01-01

Scattering of coherent light as it propagates parallel to a shock wave, formed in front of a bluff cylindrical body placed in a supersonic stream, is studied experimentally and numerically. Two incident optical fields are considered. First, a large diameter collimated beam is allowed to pass through the shock containing flow. The light intensity distribution in the resultant shadowgraph image, measured by a low light CCD camera, shows well-defined fringes upstream and downstream of the shadow cast by the shock. In the second situation, a narrow laser beam is brought to a grazing incidence on the shock and the scattered light, which appears as a diverging sheet from the point of interaction, is visualized and measured on a screen placed normal to the laser path. Experiments are conducted on shocks formed at various free-stream Mach numbers, M, and total pressures, P(sub 0). It is found that the widths of the shock shadows in a shadowgraph image become independent of M and P(sub 0) when plotted against the jump in the refractive index, (Delta)n, created across the shock. The total scattered light measured from the narrow laser beam and shock interaction also follows the same trend. In the numerical part of the study, the shock is assumed to be a 'phase object', which introduces phase difference between the upstream and downstream propagating parts of the light disturbances. For a given shape and (Delta)n of the bow shock the phase and amplitude modulations are first calculated by ray tracing. The wave front is then propagated to the screen using the Fresnet diffraction equation. The calculated intensity distribution, for both of the incident optical fields, shows good agreement with the experimental data.

Shock wave interaction with laser-generated single bubbles.

PubMed

Sankin, G N; Simmons, W N; Zhu, S L; Zhong, P

2005-07-15

The interaction of a lithotripter shock wave (LSW) with laser-generated single vapor bubbles in water is investigated using high-speed photography and pressure measurement via a fiber-optic probe hydrophone. The interaction leads to nonspherical collapse of the bubble with secondary shock wave emission and microjet formation along the LSW propagation direction. The maximum pressure amplification is produced during the collapse phase of the bubble oscillation when the compressive pulse duration of the LSW matches with the forced collapse time of the bubble.

Optical distortion in the field of a lithotripter shock wave

NASA Astrophysics Data System (ADS)

Carnell, M. T.; Emmony, D. C.

1995-10-01

The schlieren observation of cavitation phenomena produced in the tail of a lithotripter shock wave has indicated the presence of some interesting features. The images produced appear to indicate that cavitation transients in the field of a shock wave propagate nonsymmetrically; this is not the case. The apparent lack of symmetry exhibited by the primary cavitation transients is due to a complex optical lensing effect, which is brought about by the change in refractive index associated with the pressure profile of the shock wave. Objects seen through or immersed in the shock-wave field of an electromagnetic acoustic transducer, such as cavitation, appear highly distorted because of the strong positive and negative lensing effects of the compression and rarefaction cycles of the shock wave. A modification of the schlieren technique called the scale method has been used to model the distortion introduced by the shock wave and consequently explain the cavitation distortion. The technique has also been used to quantitatively analyze and partially reconstruct the lithotripter shock wave. The combination of schlieren and scale imaging gives more information about the refractive index field and therefore the shock-wave structure itself.

Portable fiber optic coupled Doppler interferometer system for detonation and shock wave diagnostics

NASA Technical Reports Server (NTRS)

Fleming, Kevin J.

1993-01-01

Testing and analysis of shock wave characteristics such as detonators and ground shock propagation frequently require a method of measuring velocity and displacement of the surface of interest. One method of measurement is Doppler interferometry. The VISAR (Velocity Interferometer System for Any Reflector) uses Doppler interferometry and has gained wide acceptance as the preferred tool for shock measurement. An important asset of VISAR is that it measures velocity and displacement nonintrusively.

Propagation of Pressure Waves, Caused by a Thermal Shock, in Liquid Metals Containing Gas Bubbles

NASA Astrophysics Data System (ADS)

Okita, Kohei; Takagi, Shu; Matsumoto, Yoichiro

The propagation of pressure waves caused by a thermal shock in liquid mercury containing micro gas bubbles has been simulated numerically. In the present study, we clarify the influences of the introduced bubble size and void fraction on the absorption of thermal expansion of liquid mercury and attenuation of pressure waves. The mass, momentum and energy conservation equations for both bubbly mixture and gas inside each bubble are solved, in which the bubble dynamics is represented by the Keller equation. The results show that when the initial void fraction is larger than the rate of the thermal expansion of liquid mercury, the pressure rise caused by the thermal expansion decreases with decreasing the bubble radius, because of the increase of the natural frequency of bubbly mixture. On the other hand, as the bubble radius increases, the peak of pressure waves which propagate at the sound speed of mixture decreases gradually due to the dispersion effect of mixture. When the natural frequency of the mixture with large bubbles is lower than that of the thremal shock, the peak pressure at the wall increases because the pressure waves propagate through the mixture at the sound speed of liquid mercury. The comparison of the results with and without heat transfer through the gas liquid interface shows that the pressure waves are attenuated greatly by the thermal damping effect with the decrease of the void fraction which enhances the nonlinearity of bubble oscillation.

Theoretical study of the effect of the size of a high-energy proton beam of the Large Hadron Collider on the formation and propagation of shock waves in copper irradiated by 450-GeV proton beams

NASA Astrophysics Data System (ADS)

Ryazanov, A. I.; Stepakov, A. V.; Vasilyev, Ya. S.; Ferrari, A.

2014-02-01

The interaction of 450-GeV protons with copper, which is the material of the collimators of the Large Hadron Collider, has been theoretically studied. A theoretical model for the formation and propagation of shock waves has been proposed on the basis of the analysis of the energy released by a proton beam in the electronic subsystem of the material owing to the deceleration of secondary particles appearing in nuclear reactions induced by this beam on the electronic subsystem of the material. The subsequent transfer of the energy from the excited electronic subsystem to the crystal lattice through the electron-phonon interaction has been described within the thermal spike model [I.M. Lifshitz, M.I. Kaganov, and L.V. Tanatarov, Sov. Phys. JETP 4, 173 (1957); I.M. Lifshitz, M.I. Kaganov, and L.V. Tanatarov, At. Energ. 6, 391 (1959); K. Yasui, Nucl. Instrum. Methods Phys. Res., Sect. B 90, 409 (1994)]. The model of the formation of shock waves involves energy exchange processes between excited electronic and ionic subsystems of the irradiated material and is based on the hydrodynamic approximation proposed by Zel'dovich [Ya.B. Zel'dovich and Yu.P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena (Nauka, Moscow, 1966; Dover, New York, 2002)]. This model makes it possible to obtain the space-time distributions of the main physical characteristics (temperatures of the ionic and electronic subsystems, density, pressure, etc.) in materials irradiated by high-energy proton beams and to analyze the formation and propagation of shock waves in them. The nonlinear differential equations describing the conservation laws of mass, energy, and momentum of electrons and ions in the Euler variables in the case of the propagation of shock waves has been solved with the Godunov scheme [S. K. Godunov, A.V. Zabrodin, M.Ya. Ivanov, A.N. Kraiko, and G.P. Prokopov, Numerical Solution of Multidimensional Problems in Gas Dynamics (Nauka, Moscow, 1976) [in Russian

The relationship between elastic constants and structure of shock waves in a zinc single crystal

NASA Astrophysics Data System (ADS)

Krivosheina, M. N.; Kobenko, S. V.; Tuch, E. V.

2017-12-01

The paper provides a 3D finite element simulation of shock-loaded anisotropic single crystals on the example of a Zn plate under impact using a mathematical model, which allows for anisotropy in hydrostatic stress and wave velocities in elastic and plastic ranges. The simulation results agree with experimental data, showing the absence of shock wave splitting into an elastic precursor and a plastic wave in Zn single crystals impacted in the [0001] direction. It is assumed that the absence of an elastic precursor under impact loading of a zinc single crystal along the [0001] direction is determined by the anomalously large ratio of the c/a-axes and close values of the propagation velocities of longitudinal and bulk elastic waves. It is shown that an increase in only one elastic constant along the [0001] direction results in shock wave splitting into an elastic precursor and a shock wave of "plastic" compression.

Application of the wavenumber jump condition to the normal and oblique interaction of a plane acoustic wave and a plane shock

NASA Technical Reports Server (NTRS)

Kleinstein, G. G.; Gunzburger, M. D.

1977-01-01

The kinematics of normal and oblique interactions between a plane acoustic wave and a plane shock wave are investigated separately using an approach whereby the shock is considered as a sharp discontinuity surface separating two half-spaces, so that the dispersion relation on either side of the shock and the wavenumber jump condition across a discontinuity surface completely specify the kinematics of the problem in the whole space independently of the acoustic-field dynamics. The normal interaction is analyzed for a stationary shock, and the spectral change of the incident wave is investigated. The normal interaction is then examined for the case of a shock wave traveling into an ambient region where an acoustic disturbance is propagating in the opposite direction. Detailed attention is given to the consequences of the existence of a critical shock speed above which the frequency of the transmitted wave becomes negative. Finally, the oblique interaction with a fixed shock is considered, and the existence and nature of the transmitted wave is investigated, particularly as a function of the angle of incidence.

Off-center blast in a shocked medium

DOE PAGES

Duncan-Miller, Gabrielle Christiane; Stone, William D.

2017-11-16

When multiple blasts occur at different times, the situation arises in which a blast wave is propagating into a medium that has already been shocked. Determining the evolution in shape of the second shock is not trivial, as it is propagating into air that is not only non-uniform, but also non-stationary. To accomplish this task, we employ the method of Kompaneets to determine the shape of a shock in a non-uniform media. We also draw from the work of Korycansky [1] on an off-center explosion in a medium with radially varying density. Extending this to treat non-stationary flow, and makingmore » use of approximations to the Sedov solution for the point blast problem, we are able to determine an analytic expression for the evolving shape of the second shock. Specifically, we consider the case of a shock in air at standard ambient temperature and pressure, with the second shock occurring shortly after the original blast wave reaches it, as in a sympathetic detonation.« less

Off-center blast in a shocked medium

DOE Office of Scientific and Technical Information (OSTI.GOV)

Duncan-Miller, Gabrielle Christiane; Stone, William D.

When multiple blasts occur at different times, the situation arises in which a blast wave is propagating into a medium that has already been shocked. Determining the evolution in shape of the second shock is not trivial, as it is propagating into air that is not only non-uniform, but also non-stationary. To accomplish this task, we employ the method of Kompaneets to determine the shape of a shock in a non-uniform media. We also draw from the work of Korycansky [1] on an off-center explosion in a medium with radially varying density. Extending this to treat non-stationary flow, and makingmore » use of approximations to the Sedov solution for the point blast problem, we are able to determine an analytic expression for the evolving shape of the second shock. Specifically, we consider the case of a shock in air at standard ambient temperature and pressure, with the second shock occurring shortly after the original blast wave reaches it, as in a sympathetic detonation.« less

Modeling multiscale evolution of numerous voids in shocked brittle material.

PubMed

Yu, Yin; Wang, Wenqiang; He, Hongliang; Lu, Tiecheng

2014-04-01

The influence of the evolution of numerous voids on macroscopic properties of materials is a multiscale problem that challenges computational research. A shock-wave compression model for brittle material, which can obtain both microscopic evolution and macroscopic shock properties, was developed using discrete element methods (lattice model). Using a model interaction-parameter-mapping procedure, qualitative features, as well as trends in the calculated shock-wave profiles, are shown to agree with experimental results. The shock wave splits into an elastic wave and a deformation wave in porous brittle materials, indicating significant shock plasticity. Void collapses in the deformation wave were the natural reason for volume shrinkage and deformation. However, media slippage and rotation deformations indicated by complex vortex patterns composed of relative velocity vectors were also confirmed as an important source of shock plasticity. With increasing pressure, the contribution from slippage deformation to the final plastic strain increased. Porosity was found to determine the amplitude of the elastic wave; porosity and shock stress together determine propagation speed of the deformation wave, as well as stress and strain on the final equilibrium state. Thus, shock behaviors of porous brittle material can be systematically designed for specific applications.

LASER APPLICATIONS AND OTHER TOPICS IN QUANTUM ELECTRONICS: Mechanism for shock wave merging in magnetised plasma: criteria and efficiency of formation of low-frequency magnetosonic waves

NASA Astrophysics Data System (ADS)

Tishchenko, V. N.; Shaikhislamov, I. F.

2010-08-01

The mechanism of merging of shock waves produced by a pulsating energy source is considered for magnetised plasma. The criteria for the emergence of this mechanism are found and its high efficiency for producing low-frequency magnetosonic waves, which have the form of a jet and propagate at large distances without attenuation, is shown.

Simulations of Shock-induced Bubble Collapse near Hard and Soft Objects

NASA Astrophysics Data System (ADS)

Rodriguez, Mauro; Johnsen, Eric

2016-11-01

Understanding the dynamics of cavitation bubbles and shock waves in and near hard and soft objects is important particularly in various naval and medical applications. Two examples are therapeutic ultrasound procedures, which utilize this phenomenon for breaking kidney stones (lithotripsy) and ablation of pathogenic tissue (histotripsy), and erosion to elastomeric coatings on propellers. Although not fully understood, the damage mechanism combines the effect of the incoming pulses and cavitation produced by the high tension of the pulses. To understand the damage mechanism, it is of key interest to quantifying the influence of the shock waves on the material and the response of the material to the shock waves. A novel Eulerian numerical approach for simulating shock and acoustic wave propagation in viscoelastic media is leveraged to understand this influence. High-fidelity simulations of the bubble collapse dynamics for various experimental configurations (i.e. the viscous or viscoelastic material surrounding the bubble and neighboring object's rigidity are varied) will be conducted. In particular, we will discuss the shock emission from collapse and its propagation in the neighboring object, including stresses thereby produced. This research was supported in part by ONR Grant N00014-12-1-0751 under Dr. Ki-Han Kim and by NSF Grant Number CBET 1253157.

Shock tube and shock wave research; Proceedings of the Eleventh International Symposium, University of Washington, Seattle, Wash., July 11-14, 1977

NASA Technical Reports Server (NTRS)

Ahlborn, B. (Editor); Hertzberg, A.; Russell, D.

1978-01-01

Papers are presented on the applications of shock-wave technology to the study of hydrodynamics, the use of the pressure-wave machine for charging diesel engines, and measurements of the heat-transfer rate in gas-turbine components. Consideration is given to shock propagation along 90-degree bends, the explosive dissemination of liquids, and rotational and vibrational relaxation behind weak shock waves in water vapor. Shock phenomena associated with expansion flows are described and stratospheric-related research using the shock tube is outlined. Attention is given to shock-wave ignition of magnesium powders, Mach reflection and boundary layers, and transition in the shock-induced unsteady boundary layer on a flat plate. Shock-tube measurements of induction and post-induction rates for low-Btu gas mixtures are presented and shock-initiated ignition in COS-N2O-Ar mixtures is described. Cluster growth rates in supersaturated lead vapor are presented and a study of laser-induced plasma motion in a solenoidal magnetic field is reviewed.

Magnetosonic shock wave in collisional pair-ion plasma

DOE Office of Scientific and Technical Information (OSTI.GOV)

Adak, Ashish, E-mail: ashish-adak@yahoo.com; Khan, Manoranjan, E-mail: mkhan.ju@gmail.com; Sikdar, Arnab, E-mail: arnabs.ju@gmail.com

2016-06-15

Nonlinear propagation of magnetosonic shock wave has been studied in collisional magnetized pair-ion plasma. The masses of both ions are same but the temperatures are slightly different. Two fluid model has been taken to describe the model. Two different modes of the magnetosonic wave have been obtained. The dynamics of the nonlinear magnetosonic wave is governed by the Korteweg-de Vries Burgers' equation. It has been shown that the ion-ion collision is the source of dissipation that causes the Burgers' term which is responsible for the shock structures in equal mass pair-ion plasma. The numerical investigations reveal that the magnetosonic wavemore » exhibits both oscillatory and monotonic shock structures depending on the strength of the dissipation. The nonlinear wave exhibited the oscillatory shock wave for strong magnetic field (weak dissipation) and monotonic shock wave for weak magnetic field (strong dissipation). The results have been discussed in the context of the fullerene pair-ion plasma experiments.« less

Experimental Insights into the Mechanisms of Particle Acceleration by Shock Waves

NASA Astrophysics Data System (ADS)

Scolamacchia, T.; Scheu, B.; Dingwell, D. B.

2011-12-01

The generation of shock waves is common during explosive volcanic eruptions. Particles acceleration following shock wave propagation has been experimentally observed suggesting the potential hazard related to this phenomenon. Experiments and numerical models focused on the dynamics of formation and propagation of different types of shock waves when overpressurized eruptive mixtures are suddenly released in the atmosphere, using a pseudo-gas approximation to model those mixtures. Nevertheless, the results of several studies indicated that the mechanism of coupling between a gas and solid particles is valid for a limited grain-size range, which at present is not well defined. We are investigating particle acceleration mechanisms using a vertical shock tube consisting of a high-pressure steel autoclave (450 mm long, 28 mm in diameter), pressurized with argon, and a low-pressure 140 mm long acrylic glass autoclave, with the same internal diameter of the HP reservoir. Shock waves are generated by Ar decompression at atmospheric pressures at Pres/Pamb 100:1 to 150:1, through the failure of a diaphragm. Experiments were performed either with empty autoclave or suspending solid analogue particles 150 μm in size inside the LP autoclave. Incident Mach number varied from 1.7 to 2.1. Absolute and relative pressure sensors monitored P histories during the entire process, and a high-speed camera recorded particles movement at 20,000 to 30,000 fps. Preliminary results indicate pressure multiplication at the contact between shock waves and the particles in a time lapse of 100s μs, suggesting a possible different mechanism with respect to gas-particle coupling for particle acceleration.

Shock initiated thermal and chemical responses of HMX crystal from ReaxFF molecular dynamics simulation.

PubMed

Zhou, Tingting; Song, Huajie; Liu, Yi; Huang, Fenglei

2014-07-21

To gain an atomistic-level understanding of the thermal and chemical responses of condensed energetic materials under thermal shock, we developed a thermal shock reactive dynamics (TS-RD) computational protocol using molecular dynamics simulation coupled with ReaxFF force field. β-Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) was selected as a a target explosive due to its wide usage in the military and industry. The results show that a thermal shock initiated by a large temperature gradient between the "hot" region and the "cold" region results in thermal expansion of the particles and induces a thermal-mechanical wave propagating back and forth in the system with an averaged velocity of 3.32 km s(-1). Heat propagating along the direction of thermal shock leads to a temperature increment of the system and thus chemical reaction initiation. Applying a continuum reactive heat conduction model combined with the temperature distribution obtained from the RD simulation, a heat conduction coefficient is derived as 0.80 W m(-1) K(-1). The chemical reaction mechanisms during thermal shock were analyzed, showing that the reaction is triggered by N-NO2 bond breaking followed by HONO elimination and ring fission. The propagation rates of the reaction front and reaction center are obtained to be 0.069 and 0.038 km s(-1), based on the time and spatial distribution of NO2. The pressure effect on the thermal shock was also investigated by employing uniaxial compression before the thermal shock. We find that compression significantly accelerates thermal-mechanical wave propagation and heat conduction, resulting in higher temperature and more excited molecules and thus earlier initiation and faster propagation of chemical reactions.

Fundamental structure of steady plastic shock waves in metals

NASA Astrophysics Data System (ADS)

Molinari, A.; Ravichandran, G.

2004-02-01

The propagation of steady plane shock waves in metallic materials is considered. Following the constitutive framework adopted by R. J. Clifton [Shock Waves and the Mechanical Properties of Solids, edited by J. J. Burke and V. Weiss (Syracuse University Press, Syracuse, N.Y., 1971), p. 73] for analyzing elastic-plastic transient waves, an analytical solution of the steady state propagation of plastic shocks is proposed. The problem is formulated in a Lagrangian setting appropriate for large deformations. The material response is characterized by a quasistatic tensile (compression) test (providing the isothermal strain hardening law). In addition the elastic response is determined up to second order elastic constants by ultrasonic measurements. Based on this simple information, it is shown that the shock kinetics can be quite well described for moderate shocks in aluminum with stress amplitude up to 10 GPa. Under the later assumption, the elastic response is assumed to be isentropic, and thermomechanical coupling is neglected. The model material considered here is aluminum, but the analysis is general and can be applied to any viscoplastic material subjected to moderate amplitude shocks. Comparisons with experimental data are made for the shock velocity, the particle velocity and the shock structure. The shock structure is obtained by quadrature of a first order differential equation, which provides analytical results under certain simplifying assumptions. The effects of material parameters and loading conditions on the shock kinetics and shock structure are discussed. The shock width is characterized by assuming an overstress formulation for the viscoplastic response. The effects on the shock structure of strain rate sensitivity are analyzed and the rationale for the J. W. Swegle and D. E. Grady [J. Appl. Phys. 58, 692 (1985)] universal scaling law for homogeneous materials is explored. Finally, the ability to deduce information on the viscoplastic response of materials subjected to very high strain rates from shock wave experiments is discussed.

Investigation on the cavitation effect of underwater shock near different boundaries

NASA Astrophysics Data System (ADS)

Xiao, Wei; Wei, Hai-peng; Feng, Liang

2017-08-01

When the shock wave of underwater explosion propagates to the surfaces of different boundaries, it gets reflected. Then, a negative pressure area is formed by the superposition of the incident wave and reflected wave. Cavitation occurs when the value of the negative pressure falls below the vapor pressure of water. An improved numerical model based on the spectral element method is applied to investigate the cavitation effect of underwater shock near different boundaries, mainly including the feature of cavitation effect near different boundaries and the influence of different parameters on cavitation effect. In the implementation of the improved numerical model, the bilinear equation of state is used to deal with the fluid field subjected to cavitation, and the field separation technique is employed to avoid the distortion of incident wave propagating through the mesh and the second-order doubly asymptotic approximation is applied to simulate the non-reflecting boundary. The main results are as follows. As the peak pressure and decay constant of shock wave increases, the range of cavitation domain increases, and the duration of cavitation increases. As the depth of water increases, the influence of cavitation on the dynamic response of spherical shell decreases.

Comparison of geometrical shock dynamics and kinematic models for shock-wave propagation

NASA Astrophysics Data System (ADS)

Ridoux, J.; Lardjane, N.; Monasse, L.; Coulouvrat, F.

2018-03-01

Geometrical shock dynamics (GSD) is a simplified model for nonlinear shock-wave propagation, based on the decomposition of the shock front into elementary ray tubes. Assuming small changes in the ray tube area, and neglecting the effect of the post-shock flow, a simple relation linking the local curvature and velocity of the front, known as the A{-}M rule, is obtained. More recently, a new simplified model, referred to as the kinematic model, was proposed. This model is obtained by combining the three-dimensional Euler equations and the Rankine-Hugoniot relations at the front, which leads to an equation for the normal variation of the shock Mach number at the wave front. In the same way as GSD, the kinematic model is closed by neglecting the post-shock flow effects. Although each model's approach is different, we prove their structural equivalence: the kinematic model can be rewritten under the form of GSD with a specific A{-}M relation. Both models are then compared through a wide variety of examples including experimental data or Eulerian simulation results when available. Attention is drawn to the simple cases of compression ramps and diffraction over convex corners. The analysis is completed by the more complex cases of the diffraction over a cylinder, a sphere, a mound, and a trough.

The effect of cosmic-ray acceleration on supernova blast wave dynamics

NASA Astrophysics Data System (ADS)

Pais, M.; Pfrommer, C.; Ehlert, K.; Pakmor, R.

2018-05-01

Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code AREPO for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.

A critical survey of wave propagation and impact in composite materials

NASA Technical Reports Server (NTRS)

Moon, F. C.

1973-01-01

A review of the field of stress waves in composite materials is presented covering the period up to December 1972. The major properties of waves in composites are discussed and a summary is made of the major experimental results in this field. Various theoretical models for analysis of wave propagation in laminated, fiber and particle reinforced composites are surveyed. The anisotropic, dispersive and dissipative properties of stress pulses and shock waves in such materials are reviewed. A review of the behavior of composites under impact loading is presented along with the application of wave propagation concepts to the determination of impact stresses in composite plates.

Elastic-plastic deformation of molybdenum single crystals shocked along [100

DOE PAGES

Mandal, A.; Gupta, Y. M.

2017-01-24

To understand the elastic-plastic deformation response of shock-compressed molybdenum (Mo) – a body-centered cubic (BCC) metal, single crystal samples were shocked along the [100] crystallographic orientation to an elastic impact stress of 12.5 GPa. Elastic-plastic wave profiles, measured at different propagation distances ranging between ~0.23 to 2.31 mm using laser interferometry, showed a time-dependent material response. Within experimental scatter, the measured elastic wave amplitudes were nearly constant over the propagation distances examined. These data point to a large and rapid elastic wave attenuation near the impact surface, before reaching a threshold value (elastic limit) of ~3.6 GPa. Numerical simulations ofmore » the measured wave profiles, performed using a dislocation-based continuum model, suggested that {110}<111> and/or {112}<111> slip systems are operative under shock loading. In contrast to shocked metal single crystals with close-packed structures, the measured wave profiles in Mo single crystals could not be explained in terms of dislocation multiplication alone. A dislocation generation mechanism, operative for shear stresses larger than that at the elastic limit, was required to model the rapid elastic wave attenuation and to provide a good overall match to the measured wave profiles. However, the physical basis for this mechanism was not established for the high-purity single crystal samples used in this study. As a result, the numerical simulations also suggested that Mo single crystals do not work harden significantly under shock loading in contrast to the behavior observed under quasi-static loading.« less

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

NASA Technical Reports Server (NTRS)

Przekwas, A. J.; Yang, H. Q.

1989-01-01

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

A Physical Model for Mass Ejection in Failed Supernovae

NASA Astrophysics Data System (ADS)

Coughlin, Eric Robert; Quataert, Eliot; Fernandez, Rodrigo; Kasen, Daniel

2018-01-01

During the core collapse of a massive star, the formation of the protoneutron star is accompanied by the emission of a significant amount of mass-energy (a few tenths of a Solar mass) in the form of neutrinos. This mass-energy loss generates an outward-propagating pressure wave that steepens into a shock near the stellar surface, potentially powering a weak transient associated with an otherwise-failed supernova -- where the shock associated with the original core collapse cannot unbind the envelope in a successful explosion. We provide both rough estimates of the energy contained in the shock that powers the transient and a general formalism for analyzing the propagation and steepening of the pressure wave, and we apply this formalism to polytropic stellar models. We compare our results to simulations, and we find excellent agreement in both the early evolution of the pressure wave and in the energy contained in the shock. Our estimates provide important constraints on the observational implications of failed supernovae.

An Extreme-ultraviolet Wave Generating Upward Secondary Waves in a Streamer-like Solar Structure

NASA Astrophysics Data System (ADS)

Zheng, Ruisheng; Chen, Yao; Feng, Shiwei; Wang, Bing; Song, Hongqiang

2018-05-01

Extreme-ultraviolet (EUV) waves, spectacular horizontally propagating disturbances in the low solar corona, always trigger horizontal secondary waves (SWs) when they encounter the ambient coronal structure. We present the first example of upward SWs in a streamer-like structure after the passing of an EUV wave. This event occurred on 2017 June 1. The EUV wave happened during a typical solar eruption including a filament eruption, a coronal mass ejection (CME), and a C6.6 flare. The EUV wave was associated with quasi-periodic fast propagating (QFP) wave trains and a type II radio burst that represented the existence of a coronal shock. The EUV wave had a fast initial velocity of ∼1000 km s‑1, comparable to high speeds of the shock and the QFP wave trains. Intriguingly, upward SWs rose slowly (∼80 km s‑1) in the streamer-like structure after the sweeping of the EUV wave. The upward SWs seemed to originate from limb brightenings that were caused by the EUV wave. All of the results show that the EUV wave is a fast-mode magnetohydrodynamic (MHD) shock wave, likely triggered by the flare impulses. We suggest that part of the EUV wave was probably trapped in the closed magnetic fields of the streamer-like structure, and upward SWs possibly resulted from the release of slow-mode trapped waves. It is believed that the interplay of the strong compression of the coronal shock and the configuration of the streamer-like structure is crucial for the formation of upward SWs.

Research on Formation Mechanism of Dynamic Response and Residual Stress of Sheet Metal Induced by Laser Shock Wave

NASA Astrophysics Data System (ADS)

Feng, Aixin; Cao, Yupeng; Wang, Heng; Zhang, Zhengang

2018-01-01

In order to reveal the quantitative control of the residual stress on the surface of metal materials, the relevant theoretical and experimental studies were carried out to investigate the dynamic response of metal thin plates and the formation mechanism of residual stress induced by laser shock wave. In this paper, the latest research trends on the surface residual stress of laser shock processing technology were elaborated. The main progress of laser shock wave propagation mechanism and dynamic response, laser shock, and surface residual stress were discussed. It is pointed out that the multi-scale characterization of laser and material, surface residual stress and microstructure change is a new hotspot in laser shock strengthening technology.

Prediction system of the 1-AU arrival times of CME-associated interplanetary shocks using three-dimensional simulations

NASA Astrophysics Data System (ADS)

den, Mitsue; Amo, Hiroyoshi; Sugihara, Kohta; Takei, Toshifumi; Ogawa, Tomoya; Tanaka, Takashi; Watari, Shinichi

We describe prediction system of the 1-AU arrival times of interplanetary shock waves associated with coromal mass ejections (CMEs). The system is based on modeling of the shock propagation using a three-dimensional adaptive mesh refinement (AMR) code. Once a CME is observed by LASCO/SOHO, firstly ambient solar wind is obtained by numerical simulation, which reproduces the solar wind parameters at that time observed by ACE spacecraft. Then we input the expansion speed and occurrence position data of that CME as initial condtions for an CME model, and 3D simulation of the CME and the shock propagation is perfomed until the shock wave passes the 1-AU. Input the parameters, execution of simulation and output of the result are available on Web, so a person who is not familiar with operation of computer or simulations or is not a researcher can use this system to predict the shock passage time. Simulated CME and shock evolution is visuallized at the same time with simulation and snap shots appear on the web automatically, so that user can follow the propagation. This system is expected to be useful for forecasters of space weather. We will describe the system and simulation model in detail.

Modeling elastic wave propagation in kidney stones with application to shock wave lithotripsy.

PubMed

Cleveland, Robin O; Sapozhnikov, Oleg A

2005-10-01

A time-domain finite-difference solution to the equations of linear elasticity was used to model the propagation of lithotripsy waves in kidney stones. The model was used to determine the loading on the stone (principal stresses and strains and maximum shear stresses and strains) due to the impact of lithotripsy shock waves. The simulations show that the peak loading induced in kidney stones is generated by constructive interference from shear waves launched from the outer edge of the stone with other waves in the stone. Notably the shear wave induced loads were significantly larger than the loads generated by the classic Hopkinson or spall effect. For simulations where the diameter of the focal spot of the lithotripter was smaller than that of the stone the loading decreased by more than 50%. The constructive interference was also sensitive to shock rise time and it was found that the peak tensile stress reduced by 30% as rise time increased from 25 to 150 ns. These results demonstrate that shear waves likely play a critical role in stone comminution and that lithotripters with large focal widths and short rise times should be effective at generating high stresses inside kidney stones.

DOE Office of Scientific and Technical Information (OSTI.GOV)

EL-Shamy, E. F., E-mail: emadel-shamy@hotmail.com; Department of Physics, College of Science, King Khalid University, P.O. 9004, Abha; Al-Asbali, A. M., E-mail: aliaa-ma@hotmail.com

A theoretical investigation is carried out to study the propagation and the head-on collision of dust-acoustic (DA) shock waves in a strongly coupled dusty plasma consisting of negative dust fluid, Maxwellian distributed electrons and ions. Applying the extended Poincaré–Lighthill–Kuo method, a couple of Korteweg–deVries–Burgers equations for describing DA shock waves are derived. This study is a first attempt to deduce the analytical phase shifts of DA shock waves after collision. The impacts of physical parameters such as the kinematic viscosity, the unperturbed electron-to-dust density ratio, parameter determining the effect of polarization force, the ion-to-electron temperature ratio, and the effective dustmore » temperature-to-ion temperature ratio on the structure and the collision of DA shock waves are examined. In addition, the results reveal the increase of the strength and the steepness of DA shock waves as the above mentioned parameters increase, which in turn leads to the increase of the phase shifts of DA shock waves after collision. The present model may be useful to describe the structure and the collision of DA shock waves in space and laboratory dusty plasmas.« less

On improvement to the Shock Propagation Model (SPM) applied to interplanetary shock transit time forecasting

NASA Astrophysics Data System (ADS)

Li, H. J.; Wei, F. S.; Feng, X. S.; Xie, Y. Q.

2008-09-01

This paper investigates methods to improve the predictions of Shock Arrival Time (SAT) of the original Shock Propagation Model (SPM). According to the classical blast wave theory adopted in the SPM, the shock propagating speed is determined by the total energy of the original explosion together with the background solar wind speed. Noting that there exists an intrinsic limit to the transit times computed by the SPM predictions for a specified ambient solar wind, we present a statistical analysis on the forecasting capability of the SPM using this intrinsic property. Two facts about SPM are found: (1) the error in shock energy estimation is not the only cause of the prediction errors and we should not expect that the accuracy of SPM to be improved drastically by an exact shock energy input; and (2) there are systematic differences in prediction results both for the strong shocks propagating into a slow ambient solar wind and for the weak shocks into a fast medium. Statistical analyses indicate the physical details of shock propagation and thus clearly point out directions of the future improvement of the SPM. A simple modification is presented here, which shows that there is room for improvement of SPM and thus that the original SPM is worthy of further development.

Dynamics of Solar Energetic Particles in the Presence of a Shock Wave

NASA Astrophysics Data System (ADS)

Timofeev, V. E.; Petukhov, Ivan; Petukhov, Stanislav; Starodubtsev, Sergei

2003-07-01

From the analysis of problem solutions on the solar energetic particle propagation in the presence of a plane shock wave described by the diffusion convective transport equation, the condition and manifestations for the influence of a shock wave on the SEP propagation in the solar wind have been determined. Solar energetic particles (SEP) in gradual events are generated by shock waves (see, for example, [1] and references there). The SEP generation region is limited, on the whole, by the solar corona. Proton fluxes of 470 MeV to 21 GeV energies, a maximum of which occur at a time when the shock in the atmosphere of the Sun reaches heights equal to 5 10 solar radii [2] indicate to it. It is also confirmed by the significant advancing of the occurrence time of maximum in the SEP intensity with kinetic energies more than 10 MeV relative to the shock front arrival moment to Earth's orbit. model calculations for the particles acceleration by the diffusive mechanism in conditions, typical for the solar corona, show that the time taken to pass the solar atmosphere by the shock is quite sufficient to form the particle spectrum corresponding to the SEP characteristics observed [3,4]. Lee and Ryan [5] investigated the problem of SEP gradual event generation, propagation and confirmed the close association between the diffusive acceleration mechanism and SEP events. The absence of depending of particle diffusion coefficients on the energy is a lack of this model. As an extension of preceding investigations, in this work the temporal dynamics of the particle spectrum in the presence of a plane shock for diffusion coefficients depending on the particle energy and also their change in time is studied. The SEP event from a moment of arising of a shock to a moment of it's arrival on the Earth's orbit can be divided on two stages: the first stage (duration is ˜ 1 hour) is a generation of SEP in the solar corona, the second stage (duration is ˜ 1 day) is a propagation in interplanetary space in the presence of a shock. Here we consider the second stage only which as believed to be began with the injection of the particle spectrum formed during the first stage.

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

DOE Office of Scientific and Technical Information (OSTI.GOV)

Luquet, David; Marchiano, Régis; Coulouvrat, François, E-mail: francois.coulouvrat@upmc.fr

2015-10-28

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength,more » the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.« less

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

NASA Astrophysics Data System (ADS)

Luquet, David; Marchiano, Régis; Coulouvrat, François

2015-10-01

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mandal, A.; Gupta, Y. M.

To understand the elastic-plastic deformation response of shock-compressed molybdenum (Mo) – a body-centered cubic (BCC) metal, single crystal samples were shocked along the [100] crystallographic orientation to an elastic impact stress of 12.5 GPa. Elastic-plastic wave profiles, measured at different propagation distances ranging between ~0.23 to 2.31 mm using laser interferometry, showed a time-dependent material response. Within experimental scatter, the measured elastic wave amplitudes were nearly constant over the propagation distances examined. These data point to a large and rapid elastic wave attenuation near the impact surface, before reaching a threshold value (elastic limit) of ~3.6 GPa. Numerical simulations ofmore » the measured wave profiles, performed using a dislocation-based continuum model, suggested that {110}<111> and/or {112}<111> slip systems are operative under shock loading. In contrast to shocked metal single crystals with close-packed structures, the measured wave profiles in Mo single crystals could not be explained in terms of dislocation multiplication alone. A dislocation generation mechanism, operative for shear stresses larger than that at the elastic limit, was required to model the rapid elastic wave attenuation and to provide a good overall match to the measured wave profiles. However, the physical basis for this mechanism was not established for the high-purity single crystal samples used in this study. As a result, the numerical simulations also suggested that Mo single crystals do not work harden significantly under shock loading in contrast to the behavior observed under quasi-static loading.« less

Cometary kilometric radio waves and plasma waves correlated with ion pick-up effect at Comet Halley

NASA Technical Reports Server (NTRS)

Oya, H.; Morioka, A.; Miyake, W.; Smith, E. J.; Tsurutani, B. T.

1985-01-01

Bow-shock movements at Comet Halley are inferred from the discrete spectra of the cometary kilometric radiation (30-195 kHz); the observed emissions can be interpreted as being generated and propagating from the moving shock. The shock motion is possibly associated with the time variation of the solar wind and cometary outgassing. It is concluded that these plasma wave phenomena are manifestations of ion pick-up processes, which occur even in a remote region 7 million to 10 million km from the cometary nucleus.

Structural Changes in Alloys of the Al-Cu-Mg System Under Ion Bombardment and Shock-Wave Loading

NASA Astrophysics Data System (ADS)

Ovchinnikov, V. V.; Gushchina, N. V.; Romanov, I. Yu.; Kaigorodova, L. I.; Grigor'ev, A. N.; Pavlenko, A. V.; Plokhoi, V. V.

2017-02-01

To confirm the hypothesis on the shock-wave nature of long-range effects upon corpuscular irradiation of condensed media presumably caused by emission and propagation of post-cascade shock waves, comparative experiments on ion beam modification and mechanical shock-wave loading of specimens of VD1 and D16 alloys of the Al-Cu-Mg system are performed. Direct analogy between the processes of microstructural change of cold-deformed VD1 and D16 alloys under mechanical shock loading and irradiation by beams of accelerated Ar+ ions (E = 20-40 keV) with low fluences (1015-1016 cm-2) is established. This demonstrates the important role of the dynamic long-range effects that have not yet been considered in classical radiation physics of solids.

Shock waves in strongly coupled plasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Khlebnikov, Sergei; Kruczenski, Martin; Michalogiorgakis, Georgios

2010-12-15

Shock waves are supersonic disturbances propagating in a fluid and giving rise to dissipation and drag. Weak shocks, i.e., those of small amplitude, can be well described within the hydrodynamic approximation. On the other hand, strong shocks are discontinuous within hydrodynamics and therefore probe the microscopics of the theory. In this paper, we consider the case of the strongly coupled N=4 plasma whose microscopic description, applicable for scales smaller than the inverse temperature, is given in terms of gravity in an asymptotically AdS{sub 5} space. In the gravity approximation, weak and strong shocks should be described by smooth metrics withmore » no discontinuities. For weak shocks, we find the dual metric in a derivative expansion, and for strong shocks we use linearized gravity to find the exponential tail that determines the width of the shock. In particular, we find that, when the velocity of the fluid relative to the shock approaches the speed of light v{yields}1 the penetration depth l scales as l{approx}(1-v{sup 2}){sup 1/4}. We compare the results with second-order hydrodynamics and the Israel-Stewart approximation. Although they all agree in the hydrodynamic regime of weak shocks, we show that there is not even qualitative agreement for strong shocks. For the gravity side, the existence of shock waves implies that there are disturbances of constant shape propagating on the horizon of the dual black holes.« less

Radiating dispersive shock waves in non-local optical media

PubMed Central

El, Gennady A.

2016-01-01

We consider the step Riemann problem for the system of equations describing the propagation of a coherent light beam in nematic liquid crystals, which is a general system describing nonlinear wave propagation in a number of different physical applications. While the equation governing the light beam is of defocusing nonlinear Schrödinger (NLS) equation type, the dispersive shock wave (DSW) generated from this initial condition has major differences from the standard DSW solution of the defocusing NLS equation. In particular, it is found that the DSW has positive polarity and generates resonant radiation which propagates ahead of it. Remarkably, the velocity of the lead soliton of the DSW is determined by the classical shock velocity. The solution for the radiative wavetrain is obtained using the Wentzel–Kramers–Brillouin approximation. It is shown that for sufficiently small initial jumps the nematic DSW is asymptotically governed by a Korteweg–de Vries equation with the fifth-order dispersion, which explicitly shows the resonance generating the radiation ahead of the DSW. The constructed asymptotic theory is shown to be in good agreement with the results of direct numerical simulations. PMID:27118911

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

NASA Technical Reports Server (NTRS)

Derkevorkian, Armen; Peterson, Lee; Kolaini, Ali R.; Hendricks, Terry J.; Nesmith, Bill J.

2016-01-01

An analytic approach is demonstrated to reveal potential pyroshock -driven dynamic effects causing power losses in the Thermo -Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi -Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high- fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large -amplitude suddenly -applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees -of-freedom (DOF). First, a quasi -static preload was applied on top of the TE module bar, then transient tri- axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA -AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high -fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Laser Light Scattering Diagnostic for Measurement of Flow Velocity in Vicinity of Propagating Shock Waves

NASA Technical Reports Server (NTRS)

Seasholtz, Richard G.; Buggele, Alvin E.

2002-01-01

A laser light scattering diagnostic for measurement of dynamic flow velocity at a point is described. The instrument is being developed for use in the study of propagating shock waves and detonation waves in pulse detonation engines under development at the NASA Glenn Research Center (GRC). The approach uses a Fabry-Perot interferometer to measure the Doppler shift of laser light scattered from small (submicron) particles in the flow. The high-speed detection system required to resolve the transient response as a shock wave crosses the probe volume uses fast response photodetectors, and a PC based data acquisition system. Preliminary results of measurements made in the GRC Mach 4, 10 by 25 cm supersonic wind tunnel are presented. Spontaneous condensation of water vapor in the flow is used as seed. The tunnel is supplied with continuous air flow at up to 45 psia and the flow is exhausted into the GRC laboratory-wide altitude exhaust system at pressures down to 0.3 psia.

Propagation of exponential shock wave in an axisymmetric rotating non-ideal dusty gas

NASA Astrophysics Data System (ADS)

Nath, G.

2016-09-01

One-dimensional unsteady isothermal and adiabatic flow behind a strong exponential shock wave propagating in a rotational axisymmetric mixture of non-ideal gas and small solid particles, which has variable azimuthal and axial fluid velocities, is analyzed. The shock wave is driven out by a piston moving with time according to exponential law. The azimuthal and axial components of the fluid velocity in the ambient medium are assumed to be varying and obeying exponential laws. In the present work, small solid particles are considered as pseudo-fluid with the assumption that the equilibrium flow-conditions are maintained in the flow-field, and the viscous-stress and heat conduction of the mixture are negligible. Solutions are obtained in both the cases, when the flow between the shock and the piston is isothermal or adiabatic by taking into account the components of vorticity vector and compressibility. It is found that the assumption of zero temperature gradient brings a profound change in the density, axial component of vorticity vector and compressibility distributions as compared to that of the adiabatic case. To investigate the behavior of the flow variables and the influence on the shock wave propagation by the parameter of non-idealness of the gas overline{b} in the mixture as well as by the mass concentration of solid particles in the mixture Kp and by the ratio of the density of solid particles to the initial density of the gas G1 are worked out in detail. It is interesting to note that the shock strength increases with an increase in G1 ; whereas it decreases with an increase in overline{b} . Also, a comparison between the solutions in the cases of isothermal and adiabatic flows is made.

Short wavelength ion waves upstream of the earth's bow shock

NASA Technical Reports Server (NTRS)

Fuselier, S. A.; Gurnett, D. A.

1984-01-01

The identification and explanation of short wavelength antenna interference effects observed in spacecraft plasma wave data have provided an important new method of determining limits on the wavelength, direction of propagation, and Doppler shift of short wavelength electrostatic waves. Using the ISEE-1 wideband electric field data, antenna interference effects have been identified in the ion waves upstream of the earth's bow shock. This identification implies that wavelengths of the upstream ion waves are shorter than the antenna length. The interference effects also provide new measurements of the direction of propagation of the ion waves. The new measurements show that the wave vectors of the ion waves are not parallel to the interplanetary magnetic field (IMF) as previously reported. The direction of propagation does not appear to be controlled by the IMF. In addition, analysis of the Doppler shift of the short wavelength ion waves has provided a measurement of the dispersion relation. The upper limit of the rest frame frequency was found to be on the order of the ion plasma frequency. At this frequency, the wavelength is on the order of a few times the Debye length. The results of this study now provide strong evidence that the ion waves in the upstream region are Doppler-shifted ion acoustic waves. Previously announced in STAR as N83-36328

On the fundamental unsteady fluid dynamics of shock-induced flows through ducts

NASA Astrophysics Data System (ADS)

Mendoza, Nicole Renee

Unsteady shock wave propagation through ducts has many applications, ranging from blast wave shelter design to advanced high-speed propulsion systems. The research objective of this study was improved fundamental understanding of the transient flow structures during unsteady shock wave propagation through rectangular ducts with varying cross-sectional area. This research focused on the fluid dynamics of the unsteady shock-induced flow fields, with an emphasis placed on understanding and characterizing the mechanisms behind flow compression (wave structures), flow induction (via shock waves), and enhanced mixing (via shock-induced viscous shear layers). A theoretical and numerical (CFD) parametric study was performed, in which the effects of these parameters on the unsteady flow fields were examined: incident shock strength, area ratio, and viscous mode (inviscid, laminar, and turbulent). Two geometries were considered: the backward-facing step (BFS) geometry, which provided a benchmark and conceptual framework, and the splitter plate (SP) geometry, which was a canonical representation of the engine flow path. The theoretical analysis was inviscid, quasi-1 D and quasi-steady; and the computational analysis was fully 2D, time-accurate, and VISCOUS. The theory provided the wave patterns and primary wave strengths for the BFS geometry, and the simulations verified the wave pattems and quantified the effects of geometry and viscosity. It was shown that the theoretical wave patterns on the BFS geometry can be used to systematically analyze the transient, 20, viscous flows on the SP geometry. This work also highlighted the importance and the role of oscillating shock and expansion waves in the development of these unsteady flows. The potential for both upstream and downstream flow induction was addressed. Positive upstream flow induction was not found in this study due to the persistent formation of an upstream-moving shock wave. Enhanced mixing was addressed by examining the evolution of the unsteady shear layer, its instability, and their effects on the flow field. The instability always appeared after the reflected shock interaction, and was exacerbated in the laminar cases and damped out in the turbulent cases. This research provided new understanding of the long-term evolution of these confined flows. Lastly, the turbulent work is one of the few turbulent studies on these flows.

Nonlinear Acoustics: Periodic Waveguide, Finite-Amplitude Propagation in a Medium Having a Distribution of Relaxation Processes, and Production of an Isolated Negative Pulse in Water

DTIC Science & Technology

1993-08-24

T. Blackstock, "Shock wave propagation and shape of the waveform," Conference on Lithotripsy (Extra-Corporeal Shock Wave Applications - Technical and...83, S5 (1988). 0574 0 b4 . D. T. Blackstock, "Physical aspects of lithotripsy ," Paper GG1, 115th Meeting, Acoustical Society of America, Seattle, 16...1991). kAlso supported in part by Grant NAG-1-1204 and University of Southampton , Eng- land. 23 1992 ONR Contract Code 1109 0 𔃻. James A. Ten Cate

Detonator Performance Characterization using Multi-Frame Laser Schlieren Imaging

NASA Astrophysics Data System (ADS)

Clarke, Steven; Landon, Colin; Murphy, Michael; Martinez, Michael; Mason, Thomas; Thomas, Keith

2009-06-01

Multi-frame Laser Schlieren Imaging of shock waves produced by detonators in transparent witness materials can be used to evaluate detonator performance. We use inverse calculations of the 2D propagation of shock waves in the EPIC finite element model computer code to calculate a temporal-spatial-pressure profile on the surface of the detonator that is consistent with the experimental shock waves from the schlieren imaging. Examples of calculated 2D temporal-spatial-pressure profiles from a range of detonator types (EFI --exploding foil initiators, DOI -- direct optical initiation, EBW -- exploding bridge wire, hotwire), detonator HE materials (PETN, HMX, etc), and HE densities. Also pressure interaction profiles from the interaction of multiple shock waves will be shown. LA-UR-09-00909.

Electromagnetic Whistler Precursors at Supercritical Interplanetary Shocks

NASA Technical Reports Server (NTRS)

Wilson, L. B., III

2012-01-01

We present observations of electromagnetic precursor waves, identified as whistler mode waves, at supercritical interplanetary shocks using the Wind search coil magnetometer. The precursors propagate obliquely with respect to the local magnetic field, shock normal vector, solar wind velocity, and they are not phase standing structures. All are right-hand polarized with respect to the magnetic field (spacecraft frame), and all but one are right-hand polarized with respect to the shock normal vector in the normal incidence frame. Particle distributions show signatures of specularly reflected gyrating ions, which may be a source of free energy for the observed modes. In one event, we simultaneously observe perpendicular ion heating and parallel electron acceleration, consistent with wave heating/acceleration due to these waves.

The measurement of shock waves following heel strike while running.

PubMed

Dickinson, J A; Cook, S D; Leinhardt, T M

1985-01-01

A non-invasive method for demonstrating the shock wave which propagates through the skeletal system following heel strike is described. This wave was not seen in force plate studies where adequate shock absorption was provided by running shoes. In the present study six subjects ran across a force plate without shoes before and after they were fatigued on a treadmill to demonstrate possible changes in the heel strike transient. Most of the parameters measured were not altered by fatigue, and a relationship between the shock wave and height, but not the weight of the runner was demonstrated. The different mechanisms leading to this phenomenon, and its implication in the areas of osteoarthritic degeneration and running mechanics are discussed.

Quantitative evaluation of the mechanical strength of titanium/composite bonding using laser-generated shock waves

NASA Astrophysics Data System (ADS)

Ducousso, M.; Bardy, S.; Rouchausse, Y.; Bergara, T.; Jenson, F.; Berthe, L.; Videau, L.; Cuvillier, N.

2018-03-01

Intense acoustic shock waves were applied to evaluate the mechanical strength of structural epoxy bonds between a TA6V4 titanium alloy and a 3D woven carbon/epoxy composite material. Two bond types with different mechanical strengths were obtained from two different adhesive reticulations, at 50% and 90% of conversion, resulting in longitudinal static strengths of 10 and 39 MPa and transverse strengths of 15 and 35 MPa, respectively. The GPa shock waves were generated using ns-scale intense laser pulses and reaction principles to a confined plasma expansion. Simulations taking into account the laser-matter interaction, plasma relaxation, and non-linear shock wave propagation were conducted to aid interpretation of the experiments. Good correlations were obtained between the experiments and the simulation and between different measurement methods of the mechanical strength (normalized tests vs laser-generated shock waves). Such results open the door toward certification of structural bonding.

Interaction of Impulsive Pressures of Cavitation Bubbles with Cell Membranes during Sonoporation

NASA Astrophysics Data System (ADS)

Kodama, Tetsuya; Koshiyama, Ken-ichiro; Tomita, Yukio; Suzuki, Maiko; Yano, Takeru; Fujikawa, Shigeo

2006-05-01

Ultrasound contrast agents (UCAs), are capable of enhancing non-invasive cytoplasmic molecular delivery in the presence of ultrasound. Collapse of UCAs may generate nano-scale cavitation bubbles, resulting in the transient permeabilization of the cell membrane. In the present study, we investigated the interaction of a cavitation bubble-induced shock wave with a cell membrane using acoustic theory and molecular dynamics (MD) simulation. From the theory, we obtained the shock wave propagation distance from the center of a cavitation bubble that would induce membrane damage. The MD simulation determined the relationship between the uptake of water molecules into the lipid bilayer and the shock wave. The interaction of the shock wave induced a structural change of the bilayer and subsequently increased the fluidity of each molecule. These changes in the bilayer due to shock waves may be an important factor in the use of UCAs to produce the transient membrane permeability during sonoporation.

Generation of Pc 1 waves by the ion temperature anisotropy associated with fast shocks caused by sudden impulses

NASA Technical Reports Server (NTRS)

Mandt, M. E.; Lee, L. C.

1991-01-01

The high correlation of Pc 1 events with magnetospheric compressions is known. A mechanism is proposed which leads to the generation of Pc 1 waves. The interaction of a dynamic pressure pulse with the earth's bow shock leads to the formation of a weak fast-mode shock propagating into the magnetoshealth. The shock wave can pass right through a tangential discontinuity (magnetopause) and into the magnetosphere, without disturbing either of the structures. In a quasiperpendicular geometry, the shock wave exhibits anisotropic heating. This anisotropy drives unstable ion-cyclotron waves which can contribute to the generation of the Pc 1 waves which are detected. The viability of the mechanism is demonstrated with simulations. This mechanism could explain the peak in the occurrence of observed Pc 1 waves in the postnoon sector where a field-aligned discontinuity in the solar wind would most often be parallel to the magnetopause surface due to the average Parker-spiral magnetic-field configuration.

Analysis of the Giacobini-Zinner bow wave

NASA Technical Reports Server (NTRS)

Smith, E. J.; Slavin, J. A.; Bame, S. J.; Thomsen, M. F.; Cowley, S. W. H.; Richardson, I. G.; Hovestadt, D.; Ipavich, F. M.; Ogilvie, K. W.; Coplan, M. A.

1986-01-01

The cometary bow wave of P/Giacobini-Zinner has been analyzed using the complete set of ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream have been analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity has then been compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The various measurements have also been compared with values derived from a Rankine-Hugoniot analysis. The results indicate that, inbound, the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained, however, arguments are presented that the bow wave is also likely to be a shock at this location.

Numerical investigation of the interaction between a planar shock wave with square and triangular bubbles containing different gases

NASA Astrophysics Data System (ADS)

Igra, Dan; Igra, Ozer

2018-05-01

The interaction between a planar shock wave and square and triangular bubbles containing either SF6, He, Ar, or CO2 is studied numerically. It is shown that, due to the existing large differences in the molecular weight, the specific heat ratio, and the acoustic impedance between these gases, different wave patterns and pressure distribution inside the bubbles are developed during the interaction process. In the case of heavy gases, the velocity of the shock wave propagating along the bubble inner surface is always less than that of the incident shock wave and higher than that of the transmitted shock wave. However, in the case of the light gas (He), the fastest one is the transmitted shock wave and the slowest one is the incident shock wave. The largest pressure jump is witnessed in the SF6 case, while the smallest pressure jump is seen in the helium case. There are also pronounced differences in the deformation of the investigated bubbles; while triangular bubbles filled with either Ar, CO2, or SF6 were deformed to a crescent shape, the helium bubble is deformed to a trapezoidal shape with three pairs of vortices emanating from its surface.

Coupled nonlinear drift and ion acoustic waves in dense dissipative electron-positron-ion magnetoplasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masood, W.; Siddiq, M.; Karim, S.

2009-11-15

Linear and nonlinear propagation characteristics of drift ion acoustic waves are investigated in an inhomogeneous electron-positron-ion (e-p-i) quantum magnetoplasma with neutrals in the background using the well known quantum hydrodynamic model. In this regard, Korteweg-de Vries-Burgers (KdVB) and Kadomtsev-Petviashvili-Burgers (KPB) equations are obtained. Furthermore, the solutions of KdVB and KPB equations are presented by using the tangent hyperbolic (tanh) method. The variation in the shock profile with the quantum Bohm potential, collision frequency, and the ratio of drift to shock velocity in the comoving frame, v{sub *}/u, is also investigated. It is found that increasing the positron concentration and collisionmore » frequency decreases the strength of the shock. It is also shown that when the localized structure propagates with velocity greater than the diamagnetic drift velocity (i.e., u>v{sub *}), the shock strength decreases. However, the shock strength is observed to increase when the localized structure propagates with velocity less than that of drift velocity (i.e., u

Comparison of CME/Shock Propagation Models with Heliospheric Imaging and In Situ Observations

NASA Astrophysics Data System (ADS)

Zhao, Xinhua; Liu, Ying D.; Inhester, Bernd; Feng, Xueshang; Wiegelmann, Thomas; Lu, Lei

2016-10-01

The prediction of the arrival time for fast coronal mass ejections (CMEs) and their associated shocks is highly desirable in space weather studies. In this paper, we use two shock propagation models, I.e., Data Guided Shock Time Of Arrival (DGSTOA) and Data Guided Shock Propagation Model (DGSPM), to predict the kinematical evolution of interplanetary shocks associated with fast CMEs. DGSTOA is based on the similarity theory of shock waves in the solar wind reference frame, and DGSPM is based on the non-similarity theory in the stationary reference frame. The inputs are the kinematics of the CME front at the maximum speed moment obtained from the geometric triangulation method applied to STEREO imaging observations together with the Harmonic Mean approximation. The outputs provide the subsequent propagation of the associated shock. We apply these models to the CMEs on 2012 January 19, January 23, and March 7. We find that the shock models predict reasonably well the shock’s propagation after the impulsive acceleration. The shock’s arrival time and local propagation speed at Earth predicted by these models are consistent with in situ measurements of WIND. We also employ the Drag-Based Model (DBM) as a comparison, and find that it predicts a steeper deceleration than the shock models after the rapid deceleration phase. The predictions of DBM at 1 au agree with the following ICME or sheath structure, not the preceding shock. These results demonstrate the applicability of the shock models used here for future arrival time prediction of interplanetary shocks associated with fast CMEs.

Shock Wave Propagation in Layered Planetary Interiors: Revisited

NASA Astrophysics Data System (ADS)

Arkani-Hamed, J.; Monteux, J.

2017-12-01

The end of the terrestrial planet accretion is characterized by numerous large impacts. About 90% of the mass of a large planet is accreted while the core mantle separation is occurring, because of the accretionary and the short-lived radio-isotope heating. The characteristics of the shockwave propagation, hence the existing scaling laws are poorly known within the layered planets. Here, we use iSALE-2D hydrocode simulations to calculate shock pressure in a differentiated Mars type body for impact velocities of 5-20 km/s, and impactor sizes of 100-400 km. We use two different rheologies for the target interior, an inviscid model ("no-stress model") and a pressure and damage-dependent strength model ("elaborated model"). To better characterize the shock pressure within the whole mantle as a function of distance from the impact site, we propose the following distribution: (1) a near field zone larger than the isobaric core that extends to 7-15 times the projectile radius into the target, where the peak shock pressure decays exponentially with increasing distance, (2) a far field zone where the pressure decays with distance following a power law. The shock pressure decreases more rapidly with distance in the near field for the elaborated model than for the no-stress model because of the influence of acoustic fluidization and damage. However to better illustrate the influence of the rheology on the shock propagation, we use the same expressions to fit the shock pressure with distance for both models. At the core-mantle boundary, CMB, the peak shock pressure jumps as the shock wave enters the core. We derived the boundary condition at CMB for the peak shock pressure. It is less sensitive to the impact velocity or the impactor size, but strongly depends on the rheology of the planet's mantle. Because of the lower shock wave velocity in the core compared to that in the mantle, the refracted shockwave propagates toward the symmetry axis of the planet, and the shock pressure in the core decreases following a second power law. In this study, we express the output obtained from iSALE hydrocodes by scaling laws to illustrate the influence of the ray angle relative to the axis of symmetry, the target rheology, the impactor size and the impact velocity. We use these shock-pressure scaling laws to determine the impact heating of terrestrial planets.

Burnett-Cattaneo continuum theory for shock waves.

PubMed

Holian, Brad Lee; Mareschal, Michel; Ravelo, Ramon

2011-02-01

We model strong shock-wave propagation, both in the ideal gas and in the dense Lennard-Jones fluid, using a refinement of earlier work, which accounts for the cold compression in the early stages of the shock rise by a nonlinear, Burnett-like, strain-rate dependence of the thermal conductivity, and relaxation of kinetic-temperature components on the hot, compressed side of the shock front. The relaxation of the disequilibrium among the three components of the kinetic temperature, namely, the difference between the component in the direction of a planar shock wave and those in the transverse directions, particularly in the region near the shock front, is accomplished at a much more quantitative level by a rigorous application of the Cattaneo-Maxwell relaxation equation to a reference solution, namely, the steady shock-wave solution of linear Navier-Stokes-Fourier theory, along with the nonlinear Burnett heat-flux term. Our new continuum theory is in nearly quantitative agreement with nonequilibrium molecular-dynamics simulations under strong shock-wave conditions, using relaxation parameters obtained from the reference solution. ©2011 American Physical Society

Ultrafast Fiber Bragg Grating Interrogation for Sensing in Detonation and Shock Wave Experiments.

PubMed

Rodriguez, George; Gilbertson, Steve M

2017-01-27

Chirped fiber Bragg grating (CFBG) sensors coupled to high speed interrogation systems are described as robust diagnostic approaches to monitoring shock wave and detonation front propagation tracking events for use in high energy density shock physics applications. Taking advantage of the linear distributed spatial encoding of the spectral band in single-mode CFBGs, embedded fiber systems and associated photonic interrogation methodologies are shown as an effective approach to sensing shock and detonation-driven loading processes along the CFBG length. Two approaches, one that detects spectral changes in the integrated spectrum of the CFBG and another coherent pulse interrogation approach that fully resolves its spectral response, shows that 100-MHz-1-GHz interrogation rates are possible with spatial resolution along the CFBG in the 50 µm to sub-millimeter range depending on the combination of CFBG parameters (i.e., length, chirp rate, spectrum) and interrogator design specifics. Results from several dynamic tests are used to demonstrate the performance of these high speed systems for shock and detonation propagation tracking under strong and weak shock pressure loading: (1) linear detonation front tracking in the plastic bonded explosive (PBX) PBX-9501; (2) tracking of radial decaying shock with crossover to non-destructive CFBG response; (3) shock wave tracking along an aluminum cylinder wall under weak loading accompanied by dynamic strain effects in the CFBG sensor.

Ultrafast Fiber Bragg Grating Interrogation for Sensing in Detonation and Shock Wave Experiments

PubMed Central

Rodriguez, George; Gilbertson, Steve M.

2017-01-01

Chirped fiber Bragg grating (CFBG) sensors coupled to high speed interrogation systems are described as robust diagnostic approaches to monitoring shock wave and detonation front propagation tracking events for use in high energy density shock physics applications. Taking advantage of the linear distributed spatial encoding of the spectral band in single-mode CFBGs, embedded fiber systems and associated photonic interrogation methodologies are shown as an effective approach to sensing shock and detonation-driven loading processes along the CFBG length. Two approaches, one that detects spectral changes in the integrated spectrum of the CFBG and another coherent pulse interrogation approach that fully resolves its spectral response, shows that 100-MHz–1-GHz interrogation rates are possible with spatial resolution along the CFBG in the 50 μm to sub-millimeter range depending on the combination of CFBG parameters (i.e., length, chirp rate, spectrum) and interrogator design specifics. Results from several dynamic tests are used to demonstrate the performance of these high speed systems for shock and detonation propagation tracking under strong and weak shock pressure loading: (1) linear detonation front tracking in the plastic bonded explosive (PBX) PBX-9501; (2) tracking of radial decaying shock with crossover to non-destructive CFBG response; (3) shock wave tracking along an aluminum cylinder wall under weak loading accompanied by dynamic strain effects in the CFBG sensor. PMID:28134819

Ultrafast Fiber Bragg Grating Interrogation for Sensing in Detonation and Shock Wave Experiments

DOE PAGES

Rodriguez, George; Gilbertson, Steve Michael

2017-01-27

Chirped fiber Bragg grating (CFBG) sensors coupled to high speed interrogation systems are described as robust diagnostic approaches to monitoring shock wave and detonation front propagation tracking events for use in high energy density shock physics applications. Taking advantage of the linear distributed spatial encoding of the spectral band in single-mode CFBGs, embedded fiber systems and associated photonic interrogation methodologies are shown as an effective approach to sensing shock and detonation-driven loading processes along the CFBG length. Two approaches, one that detects spectral changes in the integrated spectrum of the CFBG and another coherent pulse interrogation approach that fully resolvesmore » its spectral response, shows that 100-MHz–1-GHz interrogation rates are possible with spatial resolution along the CFBG in the 50 µm to sub-millimeter range depending on the combination of CFBG parameters (i.e., length, chirp rate, spectrum) and interrogator design specifics. In conclusion, results from several dynamic tests are used to demonstrate the performance of these high speed systems for shock and detonation propagation tracking under strong and weak shock pressure loading: (1) linear detonation front tracking in the plastic bonded explosive (PBX) PBX-9501; (2) tracking of radial decaying shock with crossover to non-destructive CFBG response; (3) shock wave tracking along an aluminum cylinder wall under weak loading accompanied by dynamic strain effects in the CFBG sensor.« less

Investigation of Dalton and Amagat’s laws for gas mixtures with shock propagation

DOE PAGES

Wayne, Patrick; Cooper, Sean; Simons, Dylan; ...

2017-06-20

Dalton's and Amagat's laws (also known as the law of partial pressures and the law of partial volumes respectively) are two well-known thermodynamic models describing gas mixtures. We focus our current research on determining the suitability of these models in predicting effects of shock propagation through gas mixtures. Experiments are conducted at the Shock Tube Facility at the University of New Mexico (UNM). The gas mixture used in these experiments consists of approximately 50% sulfur hexafluoride (SF6) and 50% helium (He) by mass. Fast response pressure transducers are used to obtain pressure readings both before and after the shock wave;more » these data are then used to determine the velocity of the shock wave. Temperature readings are obtained using an ultra-fast mercury cadmium telluride (MCT) infrared (IR) detector, with a response time on the order of nanoseconds. Coupled with a stabilized broadband infrared light source (operating at 1500 K), the detector provides pre- and post-shock line-of-sight readings of average temperature within the shock tube, which are used to determine the speed of sound in the gas mixture. Paired with the velocity of the shock wave, this information allows us to determine the Mach number. Our experimental results are compared with theoretical predictions of Dalton's and Amagat's laws to determine which one is more suitable.« less

Some aspects of the aeroacoustics of high-speed jets

NASA Technical Reports Server (NTRS)

Lighthill, James

1993-01-01

Some of the background to contemporary jet aeroacoustics is addressed. Then scaling laws for noise generation by low-Mach-number airflows and by turbulence convected at 'not so low' Mach number is reviewed. These laws take into account the influence of Doppler effects associated with the convection of aeroacoustic sources. Next, a uniformly valid Doppler-effect approximation exhibits the transition, with increasing Mach number of convection, from compact-source radiation at low Mach numbers to a statistical assemblage of conical shock waves radiated by eddies convected at supersonic speed. In jets, for example, supersonic eddy convection is typically found for jet exit speeds exceeding twice the atmospheric speed of sound. The Lecture continues by describing a new dynamical theory of the nonlinear propagation of such statistically random assemblages of conical shock waves. It is shown, both by a general theoretical analysis and by an illustrative computational study, how their propagation is dominated by a characteristic 'bunching' process. That process associated with a tendency for shock waves that have already formed unions with other shock waves to acquire an increased proneness to form further unions - acts so as to enhance the high-frequency part of the spectrum of noise emission from jets at these high exit speeds.

Effect of target-fixture geometry on shock-wave compacted copper powders

NASA Astrophysics Data System (ADS)

Kim, Wooyeol; Ahn, Dong-Hyun; Yoon, Jae Ik; Park, Lee Ju; Kim, Hyoung Seop

2018-01-01

In shock compaction with a single gas gun system, a target fixture is used to safely recover a powder compact processed by shock-wave dynamic impact. However, no standard fixture geometry exists, and its effect on the processed compact is not well studied. In this study, two types of fixture are used for the dynamic compaction of hydrogen-reduced copper powders, and the mechanical properties and microstructures are investigated using the Vickers microhardness test and electron backscatter diffraction, respectively. With the assistance of finite element method simulations, we analyze several shock parameters that are experimentally hard to control. The results of the simulations indicate that the target geometry clearly affects the characteristics of incident and reflected shock waves. The hardness distribution and the microstructure of the compacts also show their dependence on the geometry. With the results of the simulations and the experiment, it is concluded that the target geometry affects the shock wave propagation and wave interaction in the specimen.

Dynamics of Laser-Driven Shock Waves in Solid Targets

NASA Astrophysics Data System (ADS)

Aglitskiy, Y.; Karasik, M.; Velikovich, A. L.; Serlin, V.; Weaver, J.; Schmitt, A. J.; Obenschain, S. P.; Grun, J.; Metzler, N.; Zalesak, S. T.; Gardner, J. H.; Oh, J.; Harding, E. C.

2009-11-01

Accurate shock timing is a key issue of both indirect- and direct-drive laser fusions. The experiments on the Nike laser at NRL presented here were made possible by improvements in the imaging capability of our monochromatic x-ray diagnostics based on Bragg reflection from spherically curved crystals. Side-on imaging implemented on Nike makes it possible to observe dynamics of the shock wave and ablation front in laser-driven solid targets. We can choose to observe a sequence of 2D images or a continuous time evolution of an image resolved in one spatial dimension. A sequence of 300 ps snapshots taken using vanadium backlighter at 5.2 keV reveals propagation of a shock wave in a solid plastic target. The shape of the shock wave reflects the intensity distribution in the Nike beam. The streak records with continuous time resolution show the x-t trajectory of a laser-driven shock wave in a 10% solid density DVB foam.

Bandwidth Limitations in Characterization of High Intensity Focused Ultrasound Fields in the Presence of Shocks

NASA Astrophysics Data System (ADS)

Khokhlova, V. A.; Bessonova, O. V.; Soneson, J. E.; Canney, M. S.; Bailey, M. R.; Crum, L. A.

2010-03-01

Nonlinear propagation effects result in the formation of weak shocks in high intensity focused ultrasound (HIFU) fields. When shocks are present, the wave spectrum consists of hundreds of harmonics. In practice, shock waves are modeled using a finite number of harmonics and measured with hydrophones that have limited bandwidths. The goal of this work was to determine how many harmonics are necessary to model or measure peak pressures, intensity, and heat deposition rates of the HIFU fields. Numerical solutions of the Khokhlov-Zabolotskaya-Kuznetzov-type (KZK) nonlinear parabolic equation were obtained using two independent algorithms, compared, and analyzed for nonlinear propagation in water, in gel phantom, and in tissue. Measurements were performed in the focus of the HIFU field in the same media using fiber optic probe hydrophones of various bandwidths. Experimental data were compared to the simulation results.

An FDTD-based computer simulation platform for shock wave propagation in electrohydraulic lithotripsy.

PubMed

Yılmaz, Bülent; Çiftçi, Emre

2013-06-01

Extracorporeal Shock Wave Lithotripsy (ESWL) is based on disintegration of the kidney stone by delivering high-energy shock waves that are created outside the body and transmitted through the skin and body tissues. Nowadays high-energy shock waves are also used in orthopedic operations and investigated to be used in the treatment of myocardial infarction and cancer. Because of these new application areas novel lithotriptor designs are needed for different kinds of treatment strategies. In this study our aim was to develop a versatile computer simulation environment which would give the device designers working on various medical applications that use shock wave principle a substantial amount of flexibility while testing the effects of new parameters such as reflector size, material properties of the medium, water temperature, and different clinical scenarios. For this purpose, we created a finite-difference time-domain (FDTD)-based computational model in which most of the physical system parameters were defined as an input and/or as a variable in the simulations. We constructed a realistic computational model of a commercial electrohydraulic lithotriptor and optimized our simulation program using the results that were obtained by the manufacturer in an experimental setup. We, then, compared the simulation results with the results from an experimental setup in which oxygen level in water was varied. Finally, we studied the effects of changing the input parameters like ellipsoid size and material, temperature change in the wave propagation media, and shock wave source point misalignment. The simulation results were consistent with the experimental results and expected effects of variation in physical parameters of the system. The results of this study encourage further investigation and provide adequate evidence that the numerical modeling of a shock wave therapy system is feasible and can provide a practical means to test novel ideas in new device design procedures. Copyright © 2012 Elsevier Ireland Ltd. All rights reserved.

On the propagation of decaying planar shock and blast waves through non-uniform channels

NASA Astrophysics Data System (ADS)

Peace, J. T.; Lu, F. K.

2018-05-01

The propagation of planar decaying shock and blast waves in non-uniform channels is investigated with the use of a two-equation approximation of the generalized CCW theory. The effects of flow non-uniformity for the cases of an arbitrary strength decaying shock and blast wave in the strong shock limit are considered. Unlike the original CCW theory, the two-equation approximation takes into account the effects of initial temporal flow gradients in the flow properties behind the shock as the shock encounters an area change. A generalized order-of-magnitude analysis is carried out to analyze under which conditions the classical area-Mach (A-M) relation and two-equation approximation are valid given a time constant of decay for the flow properties behind the shock. It is shown that the two-equation approximation extends the applicability of the CCW theory to problems where flow non-uniformity behind the shock is orders of magnitude above that for appropriate use of the A-M relation. The behavior of the two-equation solution is presented for converging and diverging channels and compared against the A-M relation. It is shown that the second-order approximation and A-M relation have good agreement for converging geometries, such that the influence of flow non-uniformity behind the shock is negligible compared to the effects of changing area. Alternatively, the two-equation approximation is shown to be strongly dependent on the initial magnitude of flow non-uniformity in diverging geometries. Further, in diverging geometries, the inclusion of flow non-uniformity yields shock solutions that tend toward an acoustic wave faster than that predicted by the A-M relation.

Numerical techniques for solving nonlinear instability problems in smokeless tactical solid rocket motors. [finite difference technique

NASA Technical Reports Server (NTRS)

Baum, J. D.; Levine, J. N.

1980-01-01

The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.

Study of nonlinear electron-acoustic solitary and shock waves in a dissipative, nonplanar space plasma with superthermal hot electrons

DOE Office of Scientific and Technical Information (OSTI.GOV)

Han, Jiu-Ning, E-mail: hanjiuning@126.com; He, Yong-Lin; Luo, Jun-Hua

2014-01-15

With the consideration of the superthermal electron distribution, we present a theoretical investigation about the nonlinear propagation of electron-acoustic solitary and shock waves in a dissipative, nonplanar non-Maxwellian plasma comprised of cold electrons, superthermal hot electrons, and stationary ions. The reductive perturbation technique is used to obtain a modified Korteweg-de Vries Burgers equation for nonlinear waves in this plasma. We discuss the effects of various plasma parameters on the time evolution of nonplanar solitary waves, the profile of shock waves, and the nonlinear structure induced by the collision between planar solitary waves. It is found that these parameters have significantmore » effects on the properties of nonlinear waves and collision-induced nonlinear structure.« less

Numerical simulation of blast wave propagation in vicinity of standalone prism on flat plate

NASA Astrophysics Data System (ADS)

Valger, Svetlana; Fedorova, Natalya; Fedorov, Alexander

2018-03-01

In the paper, numerical simulation of shock wave propagation in the vicinity of a standalone prism and a prism with a cavity in front of it was carried out. The modeling was based on the solution of 3D Euler equations and Fluent software was used as a main computational tool. The algorithm for local dynamic mesh adaptation to high gradients of pressure was applied. The initial stage of the explosion of condensed explosive was described with the help of "Compressed balloon method". The research allowed describing the characteristic stages of the blast in a semi-closed space, the structure of secondary shock waves and their interaction with obstacles. The numerical approach in Fluent based on combining inviscid gas dynamics methods and "Compressed balloon method" was compared with the method which had been used by the authors earlier with the help of AUTODYN and which is based on the use of the hydrodynamic model of a material to describe state of detonation products. For the problem of shock wave propagation in the vicinity of standalone prism the comparison of the simulation results obtained using both the methods with the experimental data was performed on the dependence of static pressure and effective momentum on time for the characteristic points located on prism walls.

Ignition-and-Growth Modeling of NASA Standard Detonator and a Linear Shaped Charge

NASA Technical Reports Server (NTRS)

Oguz, Sirri

2010-01-01

The main objective of this study is to quantitatively investigate the ignition and shock sensitivity of NASA Standard Detonator (NSD) and the shock wave propagation of a linear shaped charge (LSC) after being shocked by NSD flyer plate. This combined explosive train was modeled as a coupled Arbitrary Lagrangian-Eulerian (ALE) model with LS-DYNA hydro code. An ignition-and-growth (I&G) reactive model based on unreacted and reacted Jones-Wilkins-Lee (JWL) equations of state was used to simulate the shock initiation. Various NSD-to-LSC stand-off distances were analyzed to calculate the shock initiation (or failure to initiate) and detonation wave propagation along the shaped charge. Simulation results were verified by experimental data which included VISAR tests for NSD flyer plate velocity measurement and an aluminum target severance test for LSC performance verification. Parameters used for the analysis were obtained from various published data or by using CHEETAH thermo-chemical code.

High-resolution imaging of a shock front in plastic by phase contrast imaging at LCLS

NASA Astrophysics Data System (ADS)

Beckwith, M.; Jiang, S.; Zhao, Y.; Schropp, A.; Fernandez-Panella, A.; Rinderknecht, H. G.; Wilks, S.; Fournier, K.; Galtier, E.; Xing, Z.; Granados, E.; Gamboa, E.; Glenzer, S. H.; Heimann, P.; Zastrau, U.; Cho, B. I.; Eggert, J. H.; Collins, G. W.; Ping, Y.

2017-10-01

Understanding the propagation of shock waves is important for many areas of high energy density physics, including inertial confinement fusion (ICF) and shock compression science. In order to probe the shock front structures in detail, a diagnostic capable of detecting both the small spatial and temporal changes in the material is required. Here we show the experiment using hard X-ray phase contrast imaging (PCI) to probe the shock wave propagation in polyimide with submicron spatial resolution. The experiment was performed at the Matter in Extreme Conditions (MEC) endstation of the Linac Coherent Lightsource (LCLS). PCI together with the femtosecond time scales of x-ray free electron lasers enables the imaging of optically opaque materials that undergo rapid temporal and spatial changes. The result reveals the evolution of the density profile with time. Work performed under DOE Contract No. DE-AC52-07NA27344 with support from OFES Early Career and LLNL LDRD program.

On the numerical calculation of hydrodynamic shock waves in atmospheres by an FCT method

NASA Astrophysics Data System (ADS)

Schmitz, F.; Fleck, B.

1993-11-01

The numerical calculation of vertically propagating hydrodynamic shock waves in a plane atmosphere by the ETBFCT-version of the Flux Corrected Transport (FCT) method by Boris and Book is discussed. The results are compared with results obtained by a characteristic method with shock fitting. We show that the use of the internal energy density as a dependent variable instead of the total energy density can give very inaccurate results. Consequent discretization rules for the gravitational source terms are derived. The improvement of the results by an additional iteration step is discussed. It appears that the FCT method is an excellent method for the accurate calculation of shock waves in an atmosphere.

Shear Shock Waves Observed in the Brain

NASA Astrophysics Data System (ADS)

Espíndola, David; Lee, Stephen; Pinton, Gianmarco

2017-10-01

The internal deformation of the brain is far more complex than the rigid motion of the skull. An ultrasound imaging technique that we have developed has a combination of penetration, frame-rate, and motion-detection accuracy required to directly observe the formation and evolution of shear shock waves in the brain. Experiments at low impacts on the traumatic-brain-injury scale demonstrate that they are spontaneously generated and propagate within the porcine brain. Compared to the initially smooth impact, the acceleration at the shock front is amplified up to a factor of 8.5. This highly localized increase in acceleration suggests that shear shock waves are a previously unappreciated mechanism that could play a significant role in traumatic brain injury.

Experimental characterization of plasma formation and shockwave propagation induced by high power pulsed underwater electrical discharge.

PubMed

Claverie, A; Deroy, J; Boustie, M; Avrillaud, G; Chuvatin, A; Mazanchenko, E; Demol, G; Dramane, B

2014-06-01

High power pulsed electrical discharges into liquids are investigated for new industrial applications based on the efficiency of controlled shock waves. We present here new experimental data obtained by combination of detailed high speed imaging equipments. It allows the visualization of the very first instants of plasma discharge formation, and then the pulsations of the gaseous bubble with an accurate timing of events. The time history of the expansion/compression of this bubble leads to an estimation of the energy effectively transferred to water during the discharge. Finally, the consecutive shock generation driven by this pulsating bubble is optically monitored by shadowgraphs and schlieren setup. These data provide essential information about the geometrical pattern and chronometry associated with the shock wave generation and propagation.

Experimental evidence of nonthermal acceleration of relativistic electrons by an intensive laser pulse

NASA Astrophysics Data System (ADS)

Kuramitsu, Y.; Nakanii, N.; Kondo, K.; Sakawa, Y.; Mori, Y.; Miura, E.; Tsuji, K.; Kimura, K.; Fukumochi, S.; Kashihara, M.; Tanimoto, T.; Nakamura, H.; Ishikura, T.; Takeda, K.; Tampo, M.; Kodama, R.; Kitagawa, Y.; Mima, K.; Tanaka, K. A.; Hoshino, M.; Takabe, H.

2011-02-01

Nonthermal acceleration of relativistic electrons is investigated with an intensive laser pulse. An energy distribution function of energetic particles in the universe or cosmic rays is well represented by a power-law spectrum, therefore, nonthermal acceleration is essential to understand the origin of cosmic rays. A possible candidate for the origin of cosmic rays is wakefield acceleration at relativistic astrophysical perpendicular shocks. The wakefield is considered to be excited by large-amplitude precursor light waves in the upstream of the shocks. Substituting an intensive laser pulse for the large amplitude light waves, we performed a model experiment of the shock environments in a laboratory plasma. An intensive laser pulse was propagated in a plasma tube created by imploding a hollow polystyrene cylinder, as the large amplitude light waves propagated in the upstream plasma at an astrophysical shock. Nonthermal electrons were generated, and the energy distribution functions of the electrons have a power-law component with an index of ~2. We described the detailed procedures to obtain the nonthermal components from data obtained by an electron spectrometer.

Nonlinear Propagation of Planet-Generated Tidal Waves

NASA Technical Reports Server (NTRS)

Rafikov, R. R.

2002-01-01

The propagation and evolution of planet-generated density waves in protoplanetary disks is considered. The evolution of waves, leading to shock formation and wake dissipation, is followed in the weakly nonlinear regime. The 2001 local approach of Goodman and Rafikov is extended to include the effects of surface density and temperature variations in the disk as well as the disk cylindrical geometry and nonuniform shear. Wave damping due to shocks is demonstrated to be a nonlocal process spanning a significant fraction of the disk. Torques induced by the planet could be significant drivers of disk evolution on timescales of approx. 10(exp 6)-10(exp 7) yr, even in the absence of strong background viscosity. A global prescription for angular momentum deposition is developed that could be incorporated into the study of gap formation in a gaseous disk around the planet.

Suspended liquid particle disturbance on laser-induced blast wave and low density distribution

NASA Astrophysics Data System (ADS)

Ukai, Takahiro; Zare-Behtash, Hossein; Kontis, Konstantinos

2017-12-01

The impurity effect of suspended liquid particles on the laser-induced gas breakdown was experimentally investigated in quiescent gas. The focus of this study is the investigation of the influence of the impurities on the shock wave structure as well as the low density distribution. A 532 nm Nd:YAG laser beam with an 188 mJ/pulse was focused on the chamber filled with suspended liquid particles 0.9 ± 0.63 μm in diameter. Several shock waves are generated by multiple gas breakdowns along the beam path in the breakdown with particles. Four types of shock wave structures can be observed: (1) the dual blast waves with a similar shock radius, (2) the dual blast waves with a large shock radius at the lower breakdown, (3) the dual blast waves with a large shock radius at the upper breakdown, and (4) the triple blast waves. The independent blast waves interact with each other and enhance the shock strength behind the shock front in the lateral direction. The triple blast waves lead to the strongest shock wave in all cases. The shock wave front that propagates toward the opposite laser focal spot impinges on one another, and thereafter a transmitted shock wave (TSW) appears. The TSW interacts with the low density core called a kernel; the kernel then longitudinally expands quickly due to a Richtmyer-Meshkov-like instability. The laser-particle interaction causes an increase in the kernel volume which is approximately five times as large as that in the gas breakdown without particles. In addition, the laser-particle interaction can improve the laser energy efficiency.

Re-forming supercritical quasi-parallel shocks. II - Mechanism for wave generation and front re-formation

NASA Technical Reports Server (NTRS)

Winske, D.; Thomas, V. A.; Omidi, N.; Quest, K. B.

1990-01-01

This paper continues the study of Thomas et al. (1990) in which hybrid simulations of quasi-parallel shocks were performed in one and two spatial dimensions. To identify the wave generation processes, the electromagnetic structure of the shock is examined by performing a number of one-dimensional hybrid simulations of quasi-parallel shocks for various upstream conditions. In addition, numerical experiments were carried out in which the backstreaming ions were removed from calculations to show their fundamental importance in reformation process. The calculations show that the waves are excited before ions can propagate far enough upstream to generate resonant modes. At some later times, the waves are regenerated at the leading edge of the interface, with properties like those of their initial interactions.

Shock wave attenuation in a micro-channel

NASA Astrophysics Data System (ADS)

Giordano, J.; Perrier, P.; Meister, L.; Brouillette, M.

2018-05-01

This work presents optical measurements of shock wave attenuation in a glass micro-channel. This transparent facility, with a cross section ranging from 1 mm× 150 μm to 1 mm× 500 μm, allowed for the use of high-speed schlieren videography to visualize the propagation of a shock wave within the entire micro-channel and to quantify velocity attenuation of the wave due to wall effects. In this paper, we present the experimental technique and the relevant data treatment we have used to increase the sensitivity of shock wave detection. Then, we compared our experimental results for different channel widths, lengths, and shock wave velocities with the analytical model for shock attenuation proposed by Russell (J Fluid Mech 27(2):305-314, 1967), which assumes laminar flow, and by Mirels (Attenuation in a shock tube due to unsteady-boundary-layer action, NACA Report 1333, 1957) for turbulent flow. We found that these models are inadequate to predict the observed data, owing to the presence of fully developed flow which violates the basic assumption of these models. The data are also compared with the empirical shock attenuation models proposed by Zeitoun (Phys Fluids 27(1):011701, 2015) and Deshpande and Puranik (Shock Waves 26(4):465-475, 2016), where better agreement is observed. Finally, we presented experimental data for the flow field behind the shock wave from measurements of the Mach wave angle which shows globally decreasing flow Mach numbers due to viscous wall effects.

Wave propagation characteristics of a magnetic granular chain

NASA Astrophysics Data System (ADS)

Leng, Dingxin; Liu, Guijie; Sun, Lingyu; Wang, Xiaojie

2017-10-01

We investigate the wave propagation characteristics of a horizontal alignment of magnetic grains under a non-uniform magnetic field. The magnetic force of each grain is obtained using Maxwell's principle. The contact interaction of grains is based on Hertz potential. The effects of magnetic field strength on the dynamic responses of a granular chain under strong, intermediate, and weak amplitudes of incident impulses in comparison with static precompression force are studied. Different wave propagation modes induced by the magnetic field are observed. The applied field strength demonstrably reinforces the granular-position-dependent behaviors of decreasing amplitude and increasing wave propagation velocity. The magnetic field-induced features of a magnetic granular chain have potential applications in adaptive structures for shock attenuation.

Kinetic structures of quasi-perpendicular shocks in global particle-in-cell simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Peng, Ivy Bo, E-mail: bopeng@kth.se; Markidis, Stefano; Laure, Erwin

2015-09-15

We carried out global Particle-in-Cell simulations of the interaction between the solar wind and a magnetosphere to study the kinetic collisionless physics in super-critical quasi-perpendicular shocks. After an initial simulation transient, a collisionless bow shock forms as a result of the interaction of the solar wind and a planet magnetic dipole. The shock ramp has a thickness of approximately one ion skin depth and is followed by a trailing wave train in the shock downstream. At the downstream edge of the bow shock, whistler waves propagate along the magnetic field lines and the presence of electron cyclotron waves has beenmore » identified. A small part of the solar wind ion population is specularly reflected by the shock while a larger part is deflected and heated by the shock. Solar wind ions and electrons are heated in the perpendicular directions. Ions are accelerated in the perpendicular direction in the trailing wave train region. This work is an initial effort to study the electron and ion kinetic effects developed near the bow shock in a realistic magnetic field configuration.« less

Cytoplasmic molecular delivery with shock waves: importance of impulse.

PubMed Central

Kodama, T; Hamblin, M R; Doukas, A G

2000-01-01

Cell permeabilization using shock waves may be a way of introducing macromolecules and small polar molecules into the cytoplasm, and may have applications in gene therapy and anticancer drug delivery. The pressure profile of a shock wave indicates its energy content, and shock-wave propagation in tissue is associated with cellular displacement, leading to the development of cell deformation. In the present study, three different shock-wave sources were investigated; argon fluoride excimer laser, ruby laser, and shock tube. The duration of the pressure pulse of the shock tube was 100 times longer than the lasers. The uptake of two fluorophores, calcein (molecular weight: 622) and fluorescein isothiocyanate-dextran (molecular weight: 71,600), into HL-60 human promyelocytic leukemia cells was investigated. The intracellular fluorescence was measured by a spectrofluorometer, and the cells were examined by confocal fluorescence microscopy. A single shock wave generated by the shock tube delivered both fluorophores into approximately 50% of the cells (p < 0.01), whereas shock waves from the lasers did not. The cell survival fraction was >0.95. Confocal microscopy showed that, in the case of calcein, there was a uniform fluorescence throughout the cell, whereas, in the case of FITC-dextran, the fluorescence was sometimes in the nucleus and at other times not. We conclude that the impulse of the shock wave (i.e., the pressure integrated over time), rather than the peak pressure, was a dominant factor for causing fluorophore uptake into living cells, and that shock waves might have changed the permeability of the nuclear membrane and transferred molecules directly into the nucleus. PMID:11023888

Charged Particle Distribution near the Shock Front in a Glow Discharge

NASA Astrophysics Data System (ADS)

Baryshnikov, A. S.; Basargin, I. V.; Bezverkhnii, N. O.; Bobashev, S. V.; Monakhov, N. A.; Popov, P. A.; Sakharov, V. A.; Chistyakova, M. V.

2018-02-01

The charged particle distribution near the front of a shock wave propagating in the glow discharge plasma has been investigated. It has been found that the ion concentration before the front varies nonmonotonically. Behind the shock front, the charged particle concentration varies smoothly in contrast to the neutral component density.

Tailoring the Blast Exposure Conditions in the Shock Tube for Generating Pure, Primary Shock Waves: The End Plate Facilitates Elimination of Secondary Loading of the Specimen

PubMed Central

Misistia, Anthony; Kahali, Sudeepto; Sundaramurthy, Aravind; Chandra, Namas

2016-01-01

The end plate mounted at the mouth of the shock tube is a versatile and effective implement to control and mitigate the end effects. We have performed a series of measurements of incident shock wave velocities and overpressures followed by quantification of impulse values (integral of pressure in time domain) for four different end plate configurations (0.625, 2, 4 inches, and an open end). Shock wave characteristics were monitored by high response rate pressure sensors allocated in six positions along the length of 6 meters long 229 mm square cross section shock tube. Tests were performed at three shock wave intensities, which was controlled by varying the Mylar membrane thickness (0.02, 0.04 and 0.06 inch). The end reflector plate installed at the exit of the shock tube allows precise control over the intensity of reflected waves penetrating into the shock tube. At the optimized distance of the tube to end plate gap the secondary waves were entirely eliminated from the test section, which was confirmed by pressure sensor at T4 location. This is pronounced finding for implementation of pure primary blast wave animal model. These data also suggest only deep in the shock tube experimental conditions allow exposure to a single shock wave free of artifacts. Our results provide detailed insight into spatiotemporal dynamics of shock waves with Friedlander waveform generated using helium as a driver gas and propagating in the air inside medium sized tube. Diffusion of driver gas (helium) inside the shock tube was responsible for velocity increase of reflected shock waves. Numerical simulations combined with experimental data suggest the shock wave attenuation mechanism is simply the expansion of the internal pressure. In the absence of any other postulated shock wave decay mechanisms, which were not implemented in the model the agreement between theory and experimental data is excellent. PMID:27603017

Tailoring the Blast Exposure Conditions in the Shock Tube for Generating Pure, Primary Shock Waves: The End Plate Facilitates Elimination of Secondary Loading of the Specimen.

PubMed

Kuriakose, Matthew; Skotak, Maciej; Misistia, Anthony; Kahali, Sudeepto; Sundaramurthy, Aravind; Chandra, Namas

2016-01-01

The end plate mounted at the mouth of the shock tube is a versatile and effective implement to control and mitigate the end effects. We have performed a series of measurements of incident shock wave velocities and overpressures followed by quantification of impulse values (integral of pressure in time domain) for four different end plate configurations (0.625, 2, 4 inches, and an open end). Shock wave characteristics were monitored by high response rate pressure sensors allocated in six positions along the length of 6 meters long 229 mm square cross section shock tube. Tests were performed at three shock wave intensities, which was controlled by varying the Mylar membrane thickness (0.02, 0.04 and 0.06 inch). The end reflector plate installed at the exit of the shock tube allows precise control over the intensity of reflected waves penetrating into the shock tube. At the optimized distance of the tube to end plate gap the secondary waves were entirely eliminated from the test section, which was confirmed by pressure sensor at T4 location. This is pronounced finding for implementation of pure primary blast wave animal model. These data also suggest only deep in the shock tube experimental conditions allow exposure to a single shock wave free of artifacts. Our results provide detailed insight into spatiotemporal dynamics of shock waves with Friedlander waveform generated using helium as a driver gas and propagating in the air inside medium sized tube. Diffusion of driver gas (helium) inside the shock tube was responsible for velocity increase of reflected shock waves. Numerical simulations combined with experimental data suggest the shock wave attenuation mechanism is simply the expansion of the internal pressure. In the absence of any other postulated shock wave decay mechanisms, which were not implemented in the model the agreement between theory and experimental data is excellent.

Shock Wave Propagation in Cementitious Materials at Micro/Meso Scales

DTIC Science & Technology

2015-08-31

ABSTRACT 16. SECURITY CLASSIFICATION OF: Shock wave response of heterogeneous materials like cement and concrete is greatly influenced by the...constituents and their statistical distributions. The microstructure of cement is complex due to the presence of unhydrated water, nano /micro pores, and other...heterogeneous materials like cement and concrete is greatly influenced by the constituents and their statistical distributions. The microstructure of cement

High-Speed Photography of Detonation Propagation in Dynamically Precompressed Liquid Explosives

NASA Astrophysics Data System (ADS)

Petel, O. E.; Higgins, A. J.; Yoshinaka, A. C.; Zhang, F.

2007-12-01

The propagation of detonation in shock-compressed nitromethane was observed with a high-speed framing camera. The test explosive, nitromethane, was compressed by a reverberating shock wave to pressures as high as 10 GPa prior to being detonated by a secondary detonation event. The pressure and density in the test explosive prior to detonation were determined using two methods: manganin stress gauge measurements and LS-DYNA simulations. The velocity of the detonation front was determined from consecutive frames and correlated to the density of the reverberating shock-compressed explosive prior to detonation. Observing detonation propagation under these non-ambient conditions provides data which can be useful in the validation of equation of state models.

Acceleration and propagation of energetic charged particles in the inner heliosphere.

NASA Astrophysics Data System (ADS)

Kallenrode, M. B.

1995-02-01

Both particle propagation and acceleration are intimately related to the strength of scattering. The author reviews some developments in our understanding of interplanetary propagation, in particular the dawn of a solution of the well-known discrepancy problem between mean free paths derived from quasi-linear theory and from fits to observational data. With this much improved understanding of particle scattering one can re-evaluate the understanding of particle acceleration at interplanetary shocks. Special attention is paid to the model of coupled hydrodynamic wave excitation and ion acceleration at shocks.

High frequency generation in the corona: Resonant cavities

NASA Astrophysics Data System (ADS)

Santamaria, I. C.; Van Doorsselaere, T.

2018-03-01

Aims: Null points are prominent magnetic field singularities in which the magnetic field strength strongly decreases in very small spatial scales. Around null points, predicted to be ubiquitous in the solar chromosphere and corona, the wave behavior changes considerably. Null points are also responsible for driving very energetic phenomena, and for contributing to chromospheric and coronal heating. In previous works we demonstrated that slow magneto-acoustic shock waves were generated in the chromosphere propagate through the null point, thereby producing a train of secondary shocks escaping along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 MHz. The present work aims to investigate this high frequency region around a coronal null point to give a plausible explanation to its generation at that particular frequency. Methods: We carried out a set of two-dimensional numerical simulations of wave propagation in the neighborhood of a null point located in the corona. We varied both the amplitude of the driver and the atmospheric properties to investigate the sensitivity of the high frequency waves to these parameters. Results: We demonstrate that the wave frequency is sensitive to the atmospheric parameters in the corona, but it is independent of the strength of the driver. Thus, the null point behaves as a resonant cavity generating waves at specific frequencies that depend on the background equilibrium model. Moreover, we conclude that the high frequency wave train generated at the null point is not necessarily a result of the interaction between the null point and a shock wave. This wave train can be also developed by the interaction between the null point and fast acoustic-like magneto-acoustic waves, that is, this interaction within the linear regime.

Liquid-surface entrainment induced by shocked air stream

NASA Astrophysics Data System (ADS)

Rodriguez, V.; Jourdan, G.; Marty, A.; Allou, A.; Parisse, J.-D.

2018-02-01

Recently, we experimentally studied, in a shock tube environment, shock waves propagating over horizontal free water layers having depths of 10, 20, and 30 mm for shock wave Mach numbers M_is equal to 1.1 and 1.4. The qualitative interaction process was observed by means of high-speed visualizations, and the pressures arising in the air and in the water layer were measured and interpreted in terms of the various incident and refracted shock waves in air and water; in particular, it was concluded that the compression wave in the water is driven by the planar shock wave in the air. Additional experiments have been conducted and the novel contributions of the present technical note are quantitative results regarding the liquid-surface entrainment. At low Mach number (M_is=1.1 ), we show that the velocity of the droplets ejected into the air is independent of the water depth, unlike the wavelength of initial ripples and the angle of ejection. When the shock wave strength increases (M_is=1.4 ), the dispersion of a very thin droplet mist and a single large wave take place. We show that the thickening of the water mist and the velocity of the subsequent large wave decreases with the water-layer depth.

Nonlinear reflection of shock shear waves in soft elastic media.

PubMed

Pinton, Gianmarco; Coulouvrat, François; Gennisson, Jean-Luc; Tanter, Mickaël

2010-02-01

For fluids, the theoretical investigation of shock wave reflection has a good agreement with experiments when the incident shock Mach number is large. But when it is small, theory predicts that Mach reflections are physically unrealistic, which contradicts experimental evidence. This von Neumann paradox is investigated for shear shock waves in soft elastic solids with theory and simulations. The nonlinear elastic wave equation is approximated by a paraxial wave equation with a cubic nonlinear term. This equation is solved numerically with finite differences and the Godunov scheme. Three reflection regimes are observed. Theory is developed for shock propagation by applying the Rankine-Hugoniot relations and entropic constraints. A characteristic parameter relating diffraction and non-linearity is introduced and its theoretical values are shown to match numerical observations. The numerical solution is then applied to von Neumann reflection, where curved reflected and Mach shocks are observed. Finally, the case of weak von Neumann reflection, where there is no reflected shock, is examined. The smooth but non-monotonic transition between these three reflection regimes, from linear Snell-Descartes to perfect grazing case, provides a solution to the acoustical von Neumann paradox for the shear wave equation. This transition is similar to the quadratic non-linearity in fluids.

Whistler mode waves observed by MGF search coil magnetometer -Polarization and wave normal features of upstream waves near the bow-shock

NASA Astrophysics Data System (ADS)

Hayashi, K.; Matsui, H.; Kawano, H.; Yamamoto, T.; Kokubun, S.

1994-12-01

Whistler mode waves observed in the upstream region very close to the bow-shock is focused from the initial survey for magnetic fed data in a frequency range between 1Hz and 50Hz observed by the search coil magnetometer on board the Geotail satellite. Based on the three component wave form data polarization and wave-normal characteristics of foreshock waves is first shown as dynamic spectra for the whole Fourier components of the 50 Hz band width. Intense whistler mode waves generated in the foot region of the bow-shock are found strongly controlled in the observed polarization dependent on the angle between directions of the wave propagation and the solar wind flow but not very dependent on frequency. Our simple scheme to derive the ware characteristics which is effective to survey large amount of data continuously growing is also introduced.

Nonlinear Alfvén wave propagating in ideal MHD plasmas

NASA Astrophysics Data System (ADS)

Zheng, Jugao; Chen, Yinhua; Yu, Mingyang

2016-01-01

The behavior of nonlinear Alfvén waves propagating in ideal MHD plasmas is investigated numerically. It is found that in a one-dimensional weakly nonlinear system an Alfvén wave train can excite two longitudinal disturbances, namely an acoustic wave and a ponderomotively driven disturbance, which behave differently for β \\gt 1 and β \\lt 1, where β is the ratio of plasma-to-magnetic pressures. In a strongly nonlinear system, the Alfvén wave train is modulated and can steepen to form shocks, leading to significant dissipation due to appearance of current sheets at magnetic-pressure minima. For periodic boundary condition, we find that the Alfvén wave transfers its energy to the plasma and heats it during the shock formation. In two-dimensional systems, fast magneto-acoustic wave generation due to Alfvén wave phase mixing is considered. It is found that the process depends on the amplitude and frequency of the Alfvén waves, as well as their speed gradients and the pressure of the background plasma.

Attenuation of the dynamic yield point of shocked aluminum using elastodynamic simulations of dislocation dynamics.

PubMed

Gurrutxaga-Lerma, Beñat; Balint, Daniel S; Dini, Daniele; Eakins, Daniel E; Sutton, Adrian P

2015-05-01

When a metal is subjected to extremely rapid compression, a shock wave is launched that generates dislocations as it propagates. The shock wave evolves into a characteristic two-wave structure, with an elastic wave preceding a plastic front. It has been known for more than six decades that the amplitude of the elastic wave decays the farther it travels into the metal: this is known as "the decay of the elastic precursor." The amplitude of the elastic precursor is a dynamic yield point because it marks the transition from elastic to plastic behavior. In this Letter we provide a full explanation of this attenuation using the first method of dislocation dynamics to treat the time dependence of the elastic fields of dislocations explicitly. We show that the decay of the elastic precursor is a result of the interference of the elastic shock wave with elastic waves emanating from dislocations nucleated in the shock front. Our simulations reproduce quantitatively recent experiments on the decay of the elastic precursor in aluminum and its dependence on strain rate.

Two-Dimensional Imaging Velocimetry of Heterogeneous Flow and Brittle Failure in Diamond

NASA Astrophysics Data System (ADS)

Ali, S. J.; Smith, R.; Erskine, D.; Eggert, J.; Celliers, P. M.; Collins, G. W.; Jeanloz, R.

2014-12-01

Understanding the nature and dynamics of heterogeneous flow in diamond subjected to shock compression is important for many fields of research, from inertial confinement fusion to the study of carbon rich planets. Waves propagating through a shocked material can be significantly altered by the various deformation mechanisms present in shocked materials, including anisotropic sound speeds, phase transformations, plastic/inelastic flow and brittle failure. Quantifying the spatial and temporal effects of these deformation mechanisms has been limited by a lack of diagnostics capable of obtaining simultaneous micron resolution spatial measurements and nanosecond resolution time measurements. We have utilized the 2D Janus High Resolution Velocimeter at LLNL to study the time and space dependence of fracture in shock-compressed diamond above the Hugoniot elastic limit. Previous work on the OMEGA laser facility (Rochester) has shown that the free-surface reflectivity of μm-grained diamond samples drops linearly with increasing sample pressure, whereas under the same conditions the reflectivity of nm-grained samples remains unaffected. These disparate observations can be understood by way of better documenting fracture in high-strain compression of diamond. To this end, we have imaged the development and evolution of elastic-wave propagation, plastic-wave propagation and fracture networks in the three primary orientations of single-crystal diamond, as well as in microcrystalline and nanocrystalline diamond, and find that the deformation behavior depends sensitively on the orientation and crystallinity of the diamonds.

Dissipative processes under the shock compression of glass

NASA Astrophysics Data System (ADS)

Savinykh, A. S.; Kanel, G. I.; Cherepanov, I. A.; Razorenov, S. V.

2016-03-01

New experimental data on the behavior of the K8 and TF1 glasses under shock-wave loading conditions are obtained. It is found that the propagation of shock waves is close to the self-similar one in the maximum compression stress range 4-12 GPa. Deviations from a general deformation diagram, which are related to viscous dissipation, take place when the final state of compression is approached. The parameter region in which failure waves form in glass is found not to be limited to the elastic compression stress range, as was thought earlier. The failure front velocity increases with the shock compression stress. Outside the region covered by a failure wave, the glasses demonstrate a high tensile dynamic strength (6-7 GPa) in the case of elastic compression, and this strength is still very high after transition through the elastic limit in a compression wave.

Shock wave propagation in layered planetary embryos

NASA Astrophysics Data System (ADS)

Arkani-Hamed, Jafar; Ivanov, Boris A.

2014-05-01

The propagation of impact-induced shock wave inside a planetary embryo is investigated using the Hugoniot equations and a new scaling law, governing the particle velocity variations along a shock ray inside a spherical body. The scaling law is adopted to determine the impact heating of a growing embryo in its early stage when it is an undifferentiated and uniform body. The new scaling law, similar to other existing scaling laws, is not suitable for a large differentiated embryo consisting of a silicate mantle overlying an iron core. An algorithm is developed in this study on the basis of the ray theory in a spherically symmetric body which relates the shock parameters at the top of the core to those at the base of the mantle, thus enabling the adoption of scaling laws to estimate the impact heating of both the mantle and the core. The algorithm is applied to two embryo models: a simple two-layered model with a uniform mantle overlying a uniform core, and a model where the pre-shock density and acoustic velocity of the embryo are radially dependent. The former illustrates details of the particle velocity, shock pressure, and temperature increase behind the shock front in a 2D axisymmetric geometry. The latter provides a means to compare the results with those obtained by a hydrocode simulation. The agreement between the results of the two techniques in revealing the effects of the core-mantle boundary on the shock wave transmission across the boundary is encouraging.

Relativistic shock waves and Mach cones in viscous gluon matter

NASA Astrophysics Data System (ADS)

Bouras, Ioannis; Molnár, Etele; Niemi, Harri; Xu, Zhe; El, Andrej; Fochler, Oliver; Lauciello, Francesco; Greiner, Carsten; Rischke, Dirk H.

2010-06-01

To investigate the formation and the propagation of relativistic shock waves in viscous gluon matter we solve the relativistic Riemann problem using a microscopic parton cascade. We demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to entropy density ratio η/s. Furthermore we compare our results with those obtained by solving the relativistic causal dissipative fluid equations of Israel and Stewart (IS), in order to show the validity of the IS hydrodynamics. Employing the parton cascade we also investigate the formation of Mach shocks induced by a high-energy gluon traversing viscous gluon matter. For η/s = 0.08 a Mach cone structure is observed, whereas the signal smears out for η/s >= 0.32.

Heat-flow equation motivated by the ideal-gas shock wave.

PubMed

Holian, Brad Lee; Mareschal, Michel

2010-08-01

We present an equation for the heat-flux vector that goes beyond Fourier's Law of heat conduction, in order to model shockwave propagation in gases. Our approach is motivated by the observation of a disequilibrium among the three components of temperature, namely, the difference between the temperature component in the direction of a planar shock wave, versus those in the transverse directions. This difference is most prominent near the shock front. We test our heat-flow equation for the case of strong shock waves in the ideal gas, which has been studied in the past and compared to Navier-Stokes solutions. The new heat-flow treatment improves the agreement with nonequilibrium molecular-dynamics simulations of hard spheres under strong shockwave conditions.

Shock wave attenuation by water droplets

NASA Astrophysics Data System (ADS)

Eliasson, Veronica; Wan, Qian; Deiterding, Ralf

2017-11-01

The ongoing research on shock wave attenuation is fueled by the desire to predict and avoid damage caused by shock and blast waves. For example, during an explosion in an underground mine or subway tunnel, the shock front is forced to propagate in the direction of the channel. In this work, numerical simulations using water droplets in a 2D channel are conducted to study shock wave attenuation. Four different droplet configurations (1x1, 2x2, 3x3, and 4x4) are considered, where the total volume of water is kept constant throughout all the cases. Meanwhile, the incident shock Mach number was varied from 1.1 to 1.4 with increments of 0.1. The physical motion of the water droplets, such as the center-of-mass drift and velocity, and the energy exchange between air and water are quantitatively studied. Results for center-of-mass velocity, maximum peak pressure and impulse will be presented for all different cases that were studied. NSF CBET-1437412.

Off-center blast in a shocked medium

NASA Astrophysics Data System (ADS)

Duncan-Miller, G. C.; Stone, W. D.

2018-07-01

When multiple blasts occur at different times, the situation arises in which a blast wave is propagating into a medium that has already been shocked. Determining the evolution in the shape of the second shock is not trivial, as it is propagating into air that is not only non-uniform, but also non-stationary. To accomplish this task, we employ the method of Kompaneets to determine the shape of a shock in a non-uniform media. We also draw from the work of Korycansky (Astrophys J 398:184-189. https://doi.org/10.1086/171847 , 1992) on an off-center explosion in a medium with radially varying density. Extending this to treat non-stationary flow, and making use of approximations to the Sedov solution for the point blast problem, we are able to determine an analytic expression for the evolving shape of the second shock. In particular, we consider the case of a shock in air at standard ambient temperature and pressure, with the second shock occurring shortly after the original blast wave reaches it, as in a sympathetic detonation.

Numerical Study of Interaction of a Vortical Density Inhomogeneity with Shock and Expansion Waves

NASA Technical Reports Server (NTRS)

Povitsky, A.; Ofengeim, D.

1998-01-01

We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.

Ionospheric research opportunity

NASA Astrophysics Data System (ADS)

Rickel, Dwight

1985-05-01

Ground-based explosions have been exploited successfully in the past as a relatively controlled source for producing ionospheric disturbances. On June 25, the Defense Nuclear Agency will conduct a high explosives test on the northern section of the White Sands Missile Range. Approximately 4,800 tons of ammonium nitrate and fuel oil (ANFO) will be detonated at ground level, producing an acoustic shock wave with a surface pressure change of approximately 20 mbar at a 6 km range. This shock front will have sufficient strength to propagate into the ionosphere with at least a 10% change in the ambient pressure across the disturbance front in the lower F region. Such an ionospheric perturbation will give ionospheric researchers an excellent opportunity to investigate acoustic propagation at ionospheric heights, shock dissipation effect, the ion-neutral coupling process, acoustic-gravity wave (traveling ionospheric disturbance) generation mechanisms, and associated RF phenomena.

Generation and Micro-scale Effects of Electrostatic Waves in an Oblique Shock

NASA Astrophysics Data System (ADS)

Goodrich, K.; Ergun, R.; Schwartz, S. J.; Newman, D.; Johlander, A.; Argall, M. R.; Wilder, F. D.; Torbert, R. B.; Khotyaintsev, Y. V.; Lindqvist, P. A.; Strangeway, R. J.; Russell, C. T.; Giles, B. L.; Gershman, D. J.; Burch, J. L.

2017-12-01

We present an analysis of large amplitude (>100 mV/m), high frequency (≤1 kHz), electrostatic waves observed by MMS during an oblique bow shock crossing event. The observed waves primarily consist of electrostatic solitary waves (ESWs) and oblique ion plasma waves (IPWs). ESWs typically include nonlinear structures such as double layers, ion phase-space holes, and electron phase-space holes. Oblique IPWs are observed to be similar to ion acoustic waves, but can propagate up to 70° from the ambient magnetic field direction. Both wave-modes, particularly IPWs, are observed to have very short wavelengths ( 100 m) and are highly localized. While such wave-modes have been previously observed in the terrestrial bow shock, instrumental constraints have limited detailed insight into their generation and their effect on their plasma shock environment. Analysis of this oblique shock event shows evidence that ESWs and oblique IPWs can be generated through field-aligned currents associated with magnetic turbulence and through a counterstreaming ion instability respectively. We also present evidence that this wave activity can facilitate momentum exchange between ion populations, resulting in deceleration of incoming solar wind, and localized electron heating.

Properties of ultra low frequency upstream waves at Venus and Saturn: A comparison

NASA Technical Reports Server (NTRS)

Orlowski, D. S.; Russell, C. T.; Krauss-Varban, D.; Omidi, N.

1995-01-01

The upstream regions of all planets, except Pluto, have been investigated, using in situ spacecraft measurements and a variety of analysis techniques. The detailed studies at Earth indicate that these waves are generated locally in the magnetically connected solar wind by the interaction with ions backstreaming from the shock. However, since the properties of the solar wind vary with heliocentric distance and since properties of planetary shocks depend on plasma beta, interplanetary magnetic field (IMF) spiral angle and Mach number, the amount of heating, acceleration efficiencies, etc. significantly change with heliocentric distance. In turn the waves seen at each planet propagate not in the same but different (physical) propagation modes. In this paper we compare the ULF wave observations at an outer and an inner planet. We use the results of the ratio, quantites easily derivable with sufficient accuracy at each planet. We use the full electromagnetic dispersion relation for comparison with theoretical predictions.

Shock-wave propagation and reflection in semicrystalline polyethylene: A molecular-level investigation

NASA Astrophysics Data System (ADS)

Elder, Robert M.; O'Connor, Thomas C.; Chantawansri, Tanya L.; Sliozberg, Yelena R.; Sirk, Timothy W.; Yeh, In-Chul; Robbins, Mark O.; Andzelm, Jan W.

2017-09-01

Semicrystalline polyethylene (PE) is attractive for a variety of mechanically demanding applications, where shock compression can occur. Although often highly crystalline, PE invariably contains nanoscale amorphous domains that influence shock propagation. Our objective in this work is to study the effects of such domains. To this end, we adopt a novel approach wherein we parametrize a simple continuum-level theory based on the shock impedance from molecular dynamics (MD) simulations. Using this theory, we predict how crystalline/amorphous interfaces attenuate shocks via energy reflection due to the impedance mismatch between the phases. The theory predicts that these interfaces attenuate weak shocks more effectively than strong shocks. We compare the theory to explicit nonequilibrium MD simulations of compressive shocks in semicrystalline PE containing nanometer-scale amorphous regions of varying size, where we analyze the pressure response and reflection of energy. The theory and simulations show good agreement for strong shocks (≥1.0 km /s ), but for weak shocks (<1.0 km /s ) the simulations show enhanced energy reflection relative to the continuum predictions. Furthermore, the simulations show an effect not captured by the continuum theory: the size of amorphous regions is important. The theory assumes a sharp (discontinuous) interface between two bulk phases and a sharp change in thermodynamic and hydrodynamic quantities at the shock front. However, the simulations show that when amorphous domains are narrow—with widths comparable to the shock front—reflection is reduced compared to the predictions. We identify several nanoscale mechanisms that reduce the impedance mismatch, and thus reduce reflection, at thin amorphous domains. First, the two-wave elastic-plastic structure of shocks in crystalline PE allows the faster-moving elastic precursor wave to compress small amorphous domains before the plastic wave arrives. Second, confinement between stiff, ordered crystalline domains increases the stiffness and chain ordering in small amorphous regions. Moreover, in terms of stiffness the interfaces are similar in width to the shock front, which may contribute to the underprediction of the theory for weak shocks, where the shock front is widest. We conclude by discussing the significance of these results, namely, how they can be applied to tune shock attenuation for particular applications.

Advanced numerical models and material characterisation techniques for composite materials subject to impact and shock wave loading

NASA Astrophysics Data System (ADS)

Clegg, R. A.; White, D. M.; Hayhurst, C.; Ridel, W.; Harwick, W.; Hiermaier, S.

2003-09-01

The development and validation of an advanced material model for orthotropic materials, such as fibre reinforced composites, is described. The model is specifically designed to facilitate the numerical simulation of impact and shock wave propagation through orthotropic materials and the prediction of subsequent material damage. Initial development of the model concentrated on correctly representing shock wave propagation in composite materials under high and hypervelocity impact conditions [1]. This work has now been extended to further concentrate on the development of improved numerical models and material characterisation techniques for the prediction of damage, including residual strength, in fibre reinforced composite materials. The work is focussed on Kevlar-epoxy however materials such as CFRP are also being considered. The paper describes our most recent activities in relation to the implementation of advanced material modelling options in this area. These enable refined non-liner directional characteristics of composite materials to be modelled, in addition to the correct thermodynamic response under shock wave loading. The numerical work is backed by an extensive experimental programme covering a wide range of static and dynamic tests to facilitate derivation of model input data and to validate the predicted material response. Finally, the capability of the developing composite material model is discussed in relation to a hypervelocity impact problem.

Solar Wind - Magnetosheath - Magnetopause Interactions in Global Hybrid-Vlasov Simulations

NASA Astrophysics Data System (ADS)

Hoilijoki, S.; Pfau-Kempf, Y.; Ganse, U.; Hietala, H.; Cassak, P.; Walsh, B.; Juusola, L.; Jarvinen, R.; von Alfthan, S.; Palmroth, M.

2017-12-01

We present results of interactions of solar wind and Earth's magnetosphere in global hybrid-Vlasov simulations carried out using the Vlasiator model. Vlasiator propagates ions as velocity distribution functions by solving the Vlasov equation and electrons are treated as charge-neutralizing massless fluid. Vlasiator simulations show a strong coupling between the ion scale and global scale physics. Global scale phenomena affect the local physics and the local phenomena impact the global system. Our results have shown that mirror mode waves growing in the quasi-perpendicular magnetosheath have an impact on the local reconnection rates at the dayside magnetopause. Furthermore, multiple X-line reconnection at the dayside magnetopause leads to the formation of magnetic islands (2D flux transfer events), which launch bow waves upstream propagating through the magnetosheath. These steep bow waves have the ability to accelerate ions in the magnetosheath. When the bow waves reach the bow shock they are able to bulge the shock locally. The bulge in the shock decreases the angle between the interplanetary magnetic field and the shock normal and allows ions to be reflected back to the solar wind along the magnetic field lines. Consequently, Vlasiator simulations show that magnetosheath fluctuations affect magnetopause reconnection and reconnection may influence particle acceleration and reflection in the magnetosheath and solar wind.

Global Effects of Transmitted Shock Wave Propagation Through the Earth's Inner Magnetosphere: First Results from 3-D Hybrid Kinetic Modeling

NASA Technical Reports Server (NTRS)

Lipatov, A. S.; Sibeck, D. G.

2016-01-01

We use a new hybrid kinetic model to simulate the response of ring current, outer radiation belt, and plasmaspheric particle populations to impulsive interplanetary shocks. Since particle distributions attending the interplanetary shock waves and in the ring current and radiation belts are non-Maxwellian, waveparticle interactions play a crucial role in energy transport within the inner magnetosphere. Finite gyroradius effects become important in mass loading the shock waves with the background plasma in the presence of higher energy ring current and radiation belt ions and electrons. Initial results show that shocks cause strong deformations in the global structure of the ring current, radiation belt, and plasmasphere. The ion velocity distribution functions at the shock front, in the ring current, and in the radiation belt help us determine energy transport through the Earth's inner magnetosphere.

Relativistic shock waves in an electron-positron plasma

NASA Astrophysics Data System (ADS)

Tsintsadze, Levan N.

1995-12-01

The equations describing the detailed structure of radiation electromagnetic hydrodynamics for a relativistically hot electron-positron plasma are derived. Various discontinuities are studied by these equations. It is shown that the dependence of the electron (positron) mass on the temperature changes the structure of discontinuities, including shock waves, both qualitatively and quantitatively. Steady radiative shocks are considered, which can arise in steady flows, and which also can be used to describe the propagation of shocks when the shock thickness is very small as compared to the characteristic length over which the ambient medium changes significantly. First, the magnetohydrodynamic shock wave is treated as a discontinuity and jump relations, which relate the equilibrium states of the upstream and downstream plasma far from the front, are derived. Then the structure of the front itself is considered and tangential, contact (or entropy) and rotational discontinuities are investigated.

Focusing of shock waves induced by optical breakdown in water

PubMed Central

Sankin, Georgy N.; Zhou, Yufeng; Zhong, Pei

2008-01-01

The focusing of laser-generated shock waves by a truncated ellipsoidal reflector was experimentally and numerically investigated. Pressure waveform and distribution around the first (F1) and second foci (F2) of the ellipsoidal reflector were measured. A neodymium doped yttrium aluminum garnet laser of 1046 nm wavelength and 5 ns pulse duration was used to create an optical breakdown at F1, which generates a spherically diverging shock wave with a peak pressure of 2.1–5.9 MPa at 1.1 mm stand-off distance and a pulse width at half maximum of 36–65 ns. Upon reflection, a converging shock wave is produced which, upon arriving at F2, has a leading compressive wave with a peak pressure of 26 MPa and a zero-crossing pulse duration of 0.1 μs, followed by a trailing tensile wave of −3.3 MPa peak pressure and 0.2 μs pulse duration. The −6 dB beam size of the focused shock wave field is 1.6×0.2 mm2 along and transverse to the shock wave propagation direction. Formation of elongated plasmas at high laser energy levels limits the increase in the peak pressure at F2. General features in the waveform profile of the converging shock wave are in qualitative agreement with numerical simulations based on the Hamilton model. PMID:18537359

Modeling of shock wave propagation in large amplitude ultrasound.

PubMed

Pinton, Gianmarco F; Trahey, Gregg E

2008-01-01

The Rankine-Hugoniot relation for shock wave propagation describes the shock speed of a nonlinear wave. This paper investigates time-domain numerical methods that solve the nonlinear parabolic wave equation, or the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, and the conditions they require to satisfy the Rankine-Hugoniot relation. Two numerical methods commonly used in hyperbolic conservation laws are adapted to solve the KZK equation: Godunov's method and the monotonic upwind scheme for conservation laws (MUSCL). It is shown that they satisfy the Rankine-Hugoniot relation regardless of attenuation. These two methods are compared with the current implicit solution based method. When the attenuation is small, such as in water, the current method requires a degree of grid refinement that is computationally impractical. All three numerical methods are compared in simulations for lithotripters and high intensity focused ultrasound (HIFU) where the attenuation is small compared to the nonlinearity because much of the propagation occurs in water. The simulations are performed on grid sizes that are consistent with present-day computational resources but are not sufficiently refined for the current method to satisfy the Rankine-Hugoniot condition. It is shown that satisfying the Rankine-Hugoniot conditions has a significant impact on metrics relevant to lithotripsy (such as peak pressures) and HIFU (intensity). Because the Godunov and MUSCL schemes satisfy the Rankine-Hugoniot conditions on coarse grids, they are particularly advantageous for three-dimensional simulations.

The Interaction of Turbulence with Parallel and Perpendicular Shocks: Theory and Observations at 1 au

NASA Astrophysics Data System (ADS)

Adhikari, L.; Zank, G. P.; Hunana, P.; Hu, Q.

2016-12-01

Shocks are thought to be responsible for the amplification of turbulence as well as for generating turbulence throughout the heliosphere. We study the interaction of turbulence with parallel and perpendicular shock waves using the six-coupled-equation turbulence transport model of Zank et al. We model a 1D stationary shock wave using a hyperbolic tangent function and the Rankine-Hugoniot conditions for both a reduced model with four coupled equations and the full model. Eight quasi-parallel and five quasi-perpendicular events in the WIND spacecraft data sets are identified, and we compute the fluctuating magnetic and kinetic energy, the energy in forward and backward propagating modes, the total turbulent energy, the normalized residual energy, and the normalized cross helicity upstream and downstream of the observed shocks. We compare the observed fitted values upstream and downstream of the shock with numerical solutions to our model equations. The comparison shows that our theoretical results are in reasonable agreement with observations for both quasi-parallel and perpendicular shocks. We find that (1) the total turbulent energy, the energy in forward and backward propagating modes, and the normalized residual energy increase across the shock, (2) the normalized cross helicity increases or decreases across the shock, and (3) the correlation length increases upstream and downstream of the shock, and slightly flattens or decreases across the shock.

THE INTERACTION OF TURBULENCE WITH PARALLEL AND PERPENDICULAR SHOCKS: THEORY AND OBSERVATIONS AT 1 au

DOE Office of Scientific and Technical Information (OSTI.GOV)

Adhikari, L.; Zank, G. P.; Hunana, P.

Shocks are thought to be responsible for the amplification of turbulence as well as for generating turbulence throughout the heliosphere. We study the interaction of turbulence with parallel and perpendicular shock waves using the six-coupled-equation turbulence transport model of Zank et al. We model a 1D stationary shock wave using a hyperbolic tangent function and the Rankine–Hugoniot conditions for both a reduced model with four coupled equations and the full model. Eight quasi-parallel and five quasi-perpendicular events in the WIND spacecraft data sets are identified, and we compute the fluctuating magnetic and kinetic energy, the energy in forward and backwardmore » propagating modes, the total turbulent energy, the normalized residual energy, and the normalized cross helicity upstream and downstream of the observed shocks. We compare the observed fitted values upstream and downstream of the shock with numerical solutions to our model equations. The comparison shows that our theoretical results are in reasonable agreement with observations for both quasi-parallel and perpendicular shocks. We find that (1) the total turbulent energy, the energy in forward and backward propagating modes, and the normalized residual energy increase across the shock, (2) the normalized cross helicity increases or decreases across the shock, and (3) the correlation length increases upstream and downstream of the shock, and slightly flattens or decreases across the shock.« less

First simultaneous measurements of waves generated at the bow shock in the solar wind, the magnetosphere and on the ground

NASA Astrophysics Data System (ADS)

Clausen, L. B. N.; Yeoman, T. K.; Fear, R. C.; Behlke, R.; Lucek, E. A.; Engebretson, M. J.

2009-01-01

On 5 September 2002 the Geotail satellite observed the cone angle of the Interplanetary Magnetic Field (IMF) change to values below 30° during a 56 min interval between 18:14 and 19:10 UT. This triggered the generation of upstream waves at the bow shock, 13 RE downstream of the position of Geotail. Upstream generated waves were subsequently observed by Geotail between 18:30 and 18:48 UT, during times the IMF cone angle dropped below values of 10°. At 18:24 UT all four Cluster satellites simultaneously observed a sudden increase in wave power in all three magnetic field components, independent of their position in the dayside magnetosphere. We show that the 10 min delay between the change in IMF direction as observed by Geotail and the increase in wave power observed by Cluster is consistent with the propagation of the IMF change from the Geotail position to the bow shock and the propagation of the generated waves through the bow shock, magnetosheath and magnetosphere towards the position of the Cluster satellites. We go on to show that the wave power recorded by the Cluster satellites in the component containing the poloidal and compressional pulsations was broadband and unstructured; the power in the component containing toroidal oscillations was structured and shows the existence of multi-harmonic Alfvénic continuum waves on field lines. Model predictions of these frequencies fit well with the observations. An increase in wave power associated with the change in IMF direction was also registered by ground based magnetometers which were magnetically conjunct with the Cluster satellites during the event. To the best of our knowledge we present the first simultaneous observations of waves created by backstreaming ions at the bow shock in the solar wind, the dayside magnetosphere and on the ground.

Driven waves in a two-fluid plasma

NASA Astrophysics Data System (ADS)

Roberge, W. G.; Ciolek, Glenn E.

2007-12-01

We study the physics of wave propagation in a weakly ionized plasma, as it applies to the formation of multifluid, magnetohydrodynamics (MHD) shock waves. We model the plasma as separate charged and neutral fluids which are coupled by ion-neutral friction. At times much less than the ion-neutral drag time, the fluids are decoupled and so evolve independently. At later times, the evolution is determined by the large inertial mismatch between the charged and neutral particles. The neutral flow continues to evolve independently; the charged flow is driven by and slaved to the neutral flow by friction. We calculate this driven flow analytically by considering the special but realistic case where the charged fluid obeys linearized equations of motion. We carry out an extensive analysis of linear, driven, MHD waves. The physics of driven MHD waves is embodied in certain Green functions which describe wave propagation on short time-scales, ambipolar diffusion on long time-scales and transitional behaviour at intermediate times. By way of illustration, we give an approximate solution for the formation of a multifluid shock during the collision of two identical interstellar clouds. The collision produces forward and reverse J shocks in the neutral fluid and a transient in the charged fluid. The latter rapidly evolves into a pair of magnetic precursors on the J shocks, wherein the ions undergo force-free motion and the magnetic field grows monotonically with time. The flow appears to be self-similar at the time when linear analysis ceases to be valid.

Shock wave acceleration of protons in inhomogeneous plasma interacting with ultrashort intense laser pulses

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lecz, Zs.; Andreev, A.; Max-Born Institute, Berlin

The acceleration of protons, triggered by solitary waves in expanded solid targets is investigated using particle-in-cell simulations. The near-critical density plasma is irradiated by ultrashort high power laser pulses, which generate the solitary wave. The transformation of this soliton into a shock wave during propagation in plasma with exponentially decreasing density profile is described analytically, which allows to obtain a scaling law for the proton energy. The high quality proton bunch with small energy spread is produced by reflection from the shock-front. According to the 2D simulations, the mechanism is stable only if the laser pulse duration is shorter thanmore » the characteristic development time of the parasitic Weibel instability.« less

Magnetogasdynamics shock waves in a rotational axisymmetric non-ideal gas with increasing energy and conductive and radiative heat-fluxes

NASA Astrophysics Data System (ADS)

Nath, Gorakh

2016-07-01

Self-similar solutions are obtained for one-dimensional adiabatic flow behind a magnetogasdynamics cylindrical shock wave propagating in a rotational axisymmetric non ideal gas with increasing energy and conductive and radiative heat fluxes in presence of an azimuthal magnetic field. The fluid velocities and the azimuthal magnetic field in the ambient medium are assume to be varying and obeying power laws. In order to find the similarity solutions the angular velocity of the ambient medium is taken to be decreasing as the distance from the axis increases. The heat conduction is expressed in terms of Fourier's law and the radiation is considered to be the diffusion type for an optically thick grey gas model. The thermal conductivity and the absorption coefficient are assumed to vary with temperature and density. The effects of the presence of radiation and conduction, the non-idealness of the gas and the magnetic field on the shock propagation and the flow behind the shock are investigated.

Refraction of dispersive shock waves

NASA Astrophysics Data System (ADS)

El, G. A.; Khodorovskii, V. V.; Leszczyszyn, A. M.

2012-09-01

We study a dispersive counterpart of the classical gas dynamics problem of the interaction of a shock wave with a counter-propagating simple rarefaction wave, often referred to as the shock wave refraction. The refraction of a one-dimensional dispersive shock wave (DSW) due to its head-on collision with the centred rarefaction wave (RW) is considered in the framework of the defocusing nonlinear Schrödinger (NLS) equation. For the integrable cubic nonlinearity case we present a full asymptotic description of the DSW refraction by constructing appropriate exact solutions of the Whitham modulation equations in Riemann invariants. For the NLS equation with saturable nonlinearity, whose modulation system does not possess Riemann invariants, we take advantage of the recently developed method for the DSW description in non-integrable dispersive systems to obtain main physical parameters of the DSW refraction. The key features of the DSW-RW interaction predicted by our modulation theory analysis are confirmed by direct numerical solutions of the full dispersive problem.

Nonlinear hyperbolic theory of thermal waves in metals

NASA Technical Reports Server (NTRS)

Wilhelm, H. E.; Choi, S. H.

1975-01-01

A closed-form solution for cylindrical thermal waves in metals is given based on the nonlinear hyperbolic system of energy-conservation and heat-flux relaxation equations. It is shown that heat released from a line source propagates radially outward with finite speed in the form of a thermal wave which exhibits a discontinuous wave front. Unique nonlinear thermal-wave solutions exist up to a critical amount of driving energy, i.e., for larger energy releases, the thermal flow becomes multivalued (occurrence of shock waves). By comparison, it is demonstrated that the parabolic thermal-wave theory gives, in general, a misleading picture of the profile and propagation of thermal waves and leads to physical (infinite speed of heat propagation) and mathematical (divergent energy integrals) difficulties. Attention is drawn to the importance of temporal heat-flux relaxation for the physical understanding of fast transient processes such as thermal waves and more general explosions and implosions.

Features in the Behavior of the Solar Wind behind the Bow Shock Front near the Boundary of the Earth's Magnetosphere

NASA Astrophysics Data System (ADS)

Grib, S. A.; Leora, S. N.

2017-12-01

Macroscopic discontinuous structures observed in the solar wind are considered in the framework of magnetic hydrodynamics. The interaction of strong discontinuities is studied based on the solution of the generalized Riemann-Kochin problem. The appearance of discontinuities inside the magnetosheath after the collision of the solar wind shock wave with the bow shock front is taken into account. The propagation of secondary waves appearing in the magnetosheath is considered in the approximation of one-dimensional ideal magnetohydrodynamics. The appearance of a compression wave reflected from the magnetopause is indicated. The wave can nonlinearly break with the formation of a backward shock wave and cause the motion of the bow shock towards the Sun. The interaction between shock waves is considered with the well-known trial calculation method. It is assumed that the velocity of discontinuities in the magnetosheath in the first approximation is constant on the average. All reasonings and calculations correspond to consideration of a flow region with a velocity less than the magnetosonic speed near the Earth-Sun line. It is indicated that the results agree with the data from observations carried out on the WIND and Cluster spacecrafts.

High-frequency waves in the corona due to null points

NASA Astrophysics Data System (ADS)

Santamaria, I. C.; Khomenko, E.; Collados, M.; de Vicente, A.

2017-06-01

This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have shown that high-frequency waves are generated in such a magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer. Movies attached to Figs 2 and 4 are available at http://www.aanda.org

Kinetic Properties of an Interplanetary Shock Propagating inside a Coronal Mass Ejection

NASA Astrophysics Data System (ADS)

Liu, Mingzhe; Liu, Ying D.; Yang, Zhongwei; Wilson, L. B., III; Hu, Huidong

2018-05-01

We investigate the kinetic properties of a typical fast-mode shock inside an interplanetary coronal mass ejection (ICME) observed on 1998 August 6 at 1 au, including particle distributions and wave analysis with the in situ measurements from Wind. Key results are obtained concerning the shock and the shock–ICME interaction at kinetic scales: (1) gyrating ions, which may provide energy dissipation at the shock in addition to wave-particle interactions, are observed around the shock ramp; (2) despite the enhanced proton temperature anisotropy of the shocked plasma, the low plasma β inside the ICME constrains the shocked plasma under the thresholds of the ion cyclotron and mirror-mode instabilities; (3) whistler heat flux instabilities, which can pitch-angle scatter halo electrons through a cyclotron resonance, are observed around the shock, and can explain the disappearance of bi-directional electrons (BDEs) inside the ICME together with normal betatron acceleration; (4) whistler waves near the shock are likely associated with the whistler heat flux instabilities excited at the shock ramp, which is consistent with the result that the waves may originate from the shock ramp; (5) the whistlers share a similar characteristic with the shocklet whistlers observed by Wilson et al., providing possible evidence that the shock is decaying because of the strong magnetic field inside the ICME.

Mesoscopic simulations of shock-to-detonation transition in reactive liquid high explosive

NASA Astrophysics Data System (ADS)

Maillet, J. B.; Bourasseau, E.; Desbiens, N.; Vallverdu, G.; Stoltz, G.

2011-12-01

An extension of the model described in a previous work (see Maillet J. B. et al., EPL, 78 (2007) 68001) based on Dissipative Particle Dynamics is presented and applied to a liquid high explosive (HE), with thermodynamic properties mimicking those of liquid nitromethane. Large scale nonequilibrium simulations of reacting liquid HE with model kinetic under sustained shock conditions allow a better understanding of the shock-to-detonation transition in homogeneous explosives. Moreover, the propagation of the reactive wave appears discontinuous since ignition points in the shocked material can be activated by the compressive waves emitted from the onset of chemical reactions.

Boundary identification and error analysis of shocked material images

NASA Astrophysics Data System (ADS)

Hock, Margaret; Howard, Marylesa; Cooper, Leora; Meehan, Bernard; Nelson, Keith

2017-06-01

To compute quantities such as pressure and velocity from laser-induced shock waves propagating through materials, high-speed images are captured and analyzed. Shock images typically display high noise and spatially-varying intensities, causing conventional analysis techniques to have difficulty identifying boundaries in the images without making significant assumptions about the data. We present a novel machine learning algorithm that efficiently segments, or partitions, images with high noise and spatially-varying intensities, and provides error maps that describe a level of uncertainty in the partitioning. The user trains the algorithm by providing locations of known materials within the image but no assumptions are made on the geometries in the image. The error maps are used to provide lower and upper bounds on quantities of interest, such as velocity and pressure, once boundaries have been identified and propagated through equations of state. This algorithm will be demonstrated on images of shock waves with noise and aberrations to quantify properties of the wave as it progresses. DOE/NV/25946-3126 This work was done by National Security Technologies, LLC, under Contract No. DE- AC52-06NA25946 with the U.S. Department of Energy and supported by the SDRD Program.

Comparative study of the expansion dynamics of laser-driven plasma and shock wave in in-air and underwater ablation regimes

NASA Astrophysics Data System (ADS)

Nguyen, Thao T. P.; Tanabe, Rie; Ito, Yoshiro

2018-03-01

We compared the expansion characteristics of the plasma plumes and shock waves generated in laser-induced shock process between the two ablation regimes: in air and under water. The observation was made from the initial moment when the laser pulse hit the target until 1.5 μs. The shock processes were driven by focusing a single laser pulse (1064 nm, FWHM = 13 ns) onto the surface of epoxy-resin blocks using a 40-mm focal length lens. The estimated laser intensity at the target plane is approximate to 9 ×109Wcm-2 . We used the fast-imaging technique to observe the expansion of the plasma plume and a custom-designed time-resolved photoelasticity imaging technique to observe the propagation of shock waves with the time resolution of nanoseconds. We found that at the same intensity of the laser beam, the plasma expansion during the laser pulse follows different mechanisms: the plasma plume that grows in air follows a radiation-wave model while a detonation-wave model can explain the expansion of the plasma plume induced in water. The ideal blast wave theory can be used to predict the decay of the shock wave in air but is not appropriate to describe the decay of the shock wave induced under water.

A convergent series expansion for hyperbolic systems of conservation laws

NASA Technical Reports Server (NTRS)

Harabetian, E.

1985-01-01

The discontinuities piecewise analytic initial value problem for a wide class of conservation laws is considered which includes the full three-dimensional Euler equations. The initial interaction at an arbitrary curved surface is resolved in time by a convergent series. Among other features the solution exhibits shock, contact, and expansion waves as well as sound waves propagating on characteristic surfaces. The expansion waves correspond to he one-dimensional rarefactions but have a more complicated structure. The sound waves are generated in place of zero strength shocks, and they are caused by mismatches in derivatives.

A numerical calculation of outward propagation of solar disturbances. [solar atmospheric model with shock wave propagation

NASA Technical Reports Server (NTRS)

Wu, S. T.

1974-01-01

The responses of the solar atmosphere due to an outward propagation shock are examined by employing the Lax-Wendroff method to solve the set of nonlinear partial differential equations in the model of the solar atmosphere. It is found that this theoretical model can be used to explain the solar phenomena of surge and spray. A criterion to discriminate the surge and spray is established and detailed information concerning the density, velocity, and temperature distribution with respect to the height and time is presented. The complete computer program is also included.

Nonthermal ions and associated magnetic field behavior at a quasi-parallel earth's bow shock

NASA Technical Reports Server (NTRS)

Wilkinson, W. P.; Pardaens, A. K.; Schwartz, S. J.; Burgess, D.; Luehr, H.; Kessel, R. L.; Dunlop, M.; Farrugia, C. J.

1993-01-01

Attention is given to ion and magnetic field measurements at the earth's bow shock from the AMPTE-UKS and -IRM spacecraft, which were examined in high time resolution during a 45-min interval when the field remained closely aligned with the model bow shock normal. Dense ion beams were detected almost exclusively in the midst of short-duration periods of turbulent magnetic field wave activity. Many examples of propagation at large elevation angles relative to the ecliptic plane, which is inconsistent with reflection in the standard model shock configuration, were discovered. The associated waves are elliptically polarized and are preferentially left-handed in the observer's frame of reference, but are less confined to the maximum variance plane than other previously studied foreshock waves. The association of the wave activity with the ion beams suggests that the former may be triggered by an ion-driven instability, and possible candidates are discussed.

Nonplanar electrostatic shock waves in dense plasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masood, W.; Rizvi, H.

2010-02-15

Two-dimensional quantum ion acoustic shock waves (QIASWs) are studied in an unmagnetized plasma consisting of electrons and ions. In this regard, a nonplanar quantum Kadomtsev-Petviashvili-Burgers (QKPB) equation is derived using the small amplitude perturbation expansion method. Using the tangent hyperbolic method, an analytical solution of the planar QKPB equation is obtained and subsequently used as the initial profile to numerically solve the nonplanar QKPB equation. It is observed that the increasing number density (and correspondingly the quantum Bohm potential) and kinematic viscosity affect the propagation characteristics of the QIASW. The temporal evolution of the nonplanar QIASW is investigated both inmore » Cartesian and polar planes and the results are discussed from the numerical stand point. The results of the present study may be applicable in the study of propagation of small amplitude localized electrostatic shock structures in dense astrophysical environments.« less

Test of a new heat-flow equation for dense-fluid shock waves.

PubMed

Holian, Brad Lee; Mareschal, Michel; Ravelo, Ramon

2010-09-21

Using a recently proposed equation for the heat-flux vector that goes beyond Fourier's Law of heat conduction, we model shockwave propagation in the dense Lennard-Jones fluid. Disequilibrium among the three components of temperature, namely, the difference between the kinetic temperature in the direction of a planar shock wave and those in the transverse directions, particularly in the region near the shock front, gives rise to a new transport (equilibration) mechanism not seen in usual one-dimensional heat-flow situations. The modification of the heat-flow equation was tested earlier for the case of strong shock waves in the ideal gas, which had been studied in the past and compared to Navier-Stokes-Fourier solutions. Now, the Lennard-Jones fluid, whose equation of state and transport properties have been determined from independent calculations, allows us to study the case where potential, as well as kinetic contributions are important. The new heat-flow treatment improves the agreement with nonequilibrium molecular-dynamics simulations under strong shock wave conditions, compared to Navier-Stokes.

The shock sensitivity of nitromethane/methanol mixtures

NASA Astrophysics Data System (ADS)

Bartram, Brian; Dattelbaum, Dana; Sheffield, Steve; Gibson, Lee

2013-06-01

The dilution of liquid explosives has multiple effects on detonation properties including an increase in critical diameter, spatiotemporal lengthening of the chemical reaction zone, and the development of propagating wave instabilities. Earlier detonation studies of NM/methanol mixtures have shown several effects of increasing dilution, including: 1) a continual increase in the critical diameter, 2) lowering of the Chapman-Jouguet detonation pressure, and 3) slowing of the steady detonation velocity (Koldunov et al., Comb. Expl. Shock Waves). Here, we present the results of a series of gas gun-driven plate-impact experiments to study the shock-to-detonation transition in NM/methanol mixtures. Embedded electromagnetic gauges were used to obtain in situ particle velocity wave profiles at multiple Lagrangian positions in the initiating explosive mixture. From the wave profiles obtained in each experiment, an unreacted Hugoniot locus, the initiation mechanism, and the overtake-time-to-detonation were obtained as a function of shock input condition for mixture concentrations from 100% NM to 50 wt%/50 wt% NM/methanol. Desensitization with dilution is less than expected. For example, little change in overtake time occurs in 80 wt%/20 wt% NM/methanol when compared with neat NM. Furthermore, the shock wave profiles from the gauges indicate that wave instabilities grow in as the overdriven detonation wave settles down following the shock-to-detonation transition.

Shock-induced bubble jetting into a viscous fluid with application to tissue injury in shock-wave lithotripsy.

PubMed

Freund, J B; Shukla, R K; Evan, A P

2009-11-01

Shock waves in liquids are known to cause spherical gas bubbles to rapidly collapse and form strong re-entrant jets in the direction of the propagating shock. The interaction of these jets with an adjacent viscous liquid is investigated using finite-volume simulation methods. This configuration serves as a model for tissue injury during shock-wave lithotripsy, a medical procedure to remove kidney stones. In this case, the viscous fluid provides a crude model for the tissue. It is found that for viscosities comparable to what might be expected in tissue, the jet that forms upon collapse of a small bubble fails to penetrate deeply into the viscous fluid "tissue." A simple model reproduces the penetration distance versus viscosity observed in the simulations and leads to a phenomenological model for the spreading of injury with multiple shocks. For a reasonable selection of a single efficiency parameter, this model is able to reproduce in vivo observations of an apparent 1000-shock threshold before wide-spread tissue injury occurs in targeted kidneys and the approximate extent of this injury after a typical clinical dose of 2000 shock waves.

Shock-induced bubble jetting into a viscous fluid with application to tissue injury in shock-wave lithotripsy

PubMed Central

Freund, J. B.; Shukla, R. K.; Evan, A. P.

2009-01-01

Shock waves in liquids are known to cause spherical gas bubbles to rapidly collapse and form strong re-entrant jets in the direction of the propagating shock. The interaction of these jets with an adjacent viscous liquid is investigated using finite-volume simulation methods. This configuration serves as a model for tissue injury during shock-wave lithotripsy, a medical procedure to remove kidney stones. In this case, the viscous fluid provides a crude model for the tissue. It is found that for viscosities comparable to what might be expected in tissue, the jet that forms upon collapse of a small bubble fails to penetrate deeply into the viscous fluid “tissue.” A simple model reproduces the penetration distance versus viscosity observed in the simulations and leads to a phenomenological model for the spreading of injury with multiple shocks. For a reasonable selection of a single efficiency parameter, this model is able to reproduce in vivo observations of an apparent 1000-shock threshold before wide-spread tissue injury occurs in targeted kidneys and the approximate extent of this injury after a typical clinical dose of 2000 shock waves. PMID:19894850

Hydrodynamics of a cold one-dimensional fluid: the problem of strong shock waves

NASA Astrophysics Data System (ADS)

Hurtado, Pablo I.

2005-03-01

We study a shock wave induced by an infinitely massive piston propagating into a one-dimensional cold gas. The cold gas is modelled as a collection of hard rods which are initially at rest, so the temperature is zero. Most of our results are based on simulations of a gas of rods with binary mass distribution, and we partcularly focus on the case of spatially alternating masses. We find that the properties of the resulting shock wave are in striking contrast with those predicted by hydrodynamic and kinetic approaches, e.g., the flow-field profiles relax algebraically toward their equilibrium values. In addition, most relevant observables characterizing local thermodynamic equilibrium and equipartition decay as a power law of the distance to the shock layer. The exponents of these power laws depend non-monotonously on the mass ratio. Similar interesting dependences on the mass ratio also characterize the shock width, density and temperature overshoots, etc.

CHROMOSPHERIC AND CORONAL WAVE GENERATION IN A MAGNETIC FLUX SHEATH

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kato, Yoshiaki; Hansteen, Viggo; Gudiksen, Boris

2016-08-10

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab “pump” the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field inmore » the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of ν ≈ 4 mHz. We conclude that this process of “magnetic pumping” is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.« less

Time resolved interferometric study of the plasma plume induced shock wave in confined geometry: Two-dimensional mapping of the ambient and plasma density

DOE Office of Scientific and Technical Information (OSTI.GOV)

Choudhury, Kaushik; Singh, R. K.; Kumar, Ajai, E-mail: ajai@ipr.res.in

2016-04-15

An experimental investigation of the laser produced plasma induced shock wave in the presence of confining walls placed along the axial as well as the lateral direction has been performed. A time resolved Mach Zehnder interferometer is set up to track the primary as well as the reflected shock waves and its effect on the evolving plasma plume has been studied. An attempt has been made to discriminate the electronic and medium density contributions towards the changes in the refractive index of the medium. Two dimensional spatial distributions for both ambient medium density and plasma density (electron density) have beenmore » obtained by employing customised inversion technique and algorithm on the recorded interferograms. The observed density pattern of the surrounding medium in the presence of confining walls is correlated with the reflected shock wave propagation in the medium. Further, the shock wave plasma interaction and the subsequent changes in the shape and density of the plasma plume in confined geometry are briefly described.« less

Measuring the flexoelectric coefficient of bulk barium titanate from a shock wave experiment

NASA Astrophysics Data System (ADS)

Hu, Taotao; Deng, Qian; Liang, Xu; Shen, Shengping

2017-08-01

In this paper, a phenomenon of polarization introduced by shock waves is experimentally studied. Although this phenomenon has been reported previously in the community of physics, this is the first time to link it to flexoelectricity, the coupling between electric polarization and strain gradients in dielectrics. As the shock waves propagate in a dielectric material, electric polarization is thought to be induced by the strain gradient at the shock front. First, we control the first-order hydrogen gas gun to impact and generate shock waves in unpolarized bulk barium titanate (BT) samples. Then, a high-precision oscilloscope is used to measure the voltage generated by the flexoelectric effect. Based on experimental results, strain elastic wave theory, and flexoelectric theory, a longitudinal flexoelectric coefficient of the bulk BT sample is calculated to be μ 11 = 17.33 × 10 - 6 C/m, which is in accord with the published transverse flexoelectric coefficient. This method effectively suppresses the majority of drawbacks in the quasi-static and low frequency dynamic techniques and provides more reliable results of flexoelectric behaviors.

Observations of Electromagnetic Whistler Precursors at Supercritical Interplanetary Shocks

NASA Technical Reports Server (NTRS)

Wilson, L. B., III; Koval, A.; Szabo, Adam; Breneman, A.; Cattell, C. A.; Goetz, K.; Kellogg, P. J.; Kersten, K.; Kasper, J. C.; Maruca, B. A.;

Density Shock Waves in Confined Microswimmers

NASA Astrophysics Data System (ADS)

Tsang, Alan Cheng Hou; Kanso, Eva

2016-01-01

Motile and driven particles confined in microfluidic channels exhibit interesting emergent behavior, from propagating density bands to density shock waves. A deeper understanding of the physical mechanisms responsible for these emergent structures is relevant to a number of physical and biomedical applications. Here, we study the formation of density shock waves in the context of an idealized model of microswimmers confined in a narrow channel and subject to a uniform external flow. Interestingly, these density shock waves exhibit a transition from "subsonic" with compression at the back to "supersonic" with compression at the front of the population as the intensity of the external flow increases. This behavior is the result of a nontrivial interplay between hydrodynamic interactions and geometric confinement, and it is confirmed by a novel quasilinear wave model that properly captures the dependence of the shock formation on the external flow. These findings can be used to guide the development of novel mechanisms for controlling the emergent density distribution and the average population speed, with potentially profound implications on various processes in industry and biotechnology, such as the transport and sorting of cells in flow channels.

Effects of Alfvénic Drift on Diffusive Shock Acceleration at Weak Cluster Shocks

NASA Astrophysics Data System (ADS)

Kang, Hyesung; Ryu, Dongsu

2018-03-01

Non-detection of γ-ray emission from galaxy clusters has challenged diffusive shock acceleration (DSA) of cosmic-ray (CR) protons at weak collisionless shocks that are expected to form in the intracluster medium. As an effort to address this problem, we here explore possible roles of Alfvén waves self-excited via resonant streaming instability during the CR acceleration at parallel shocks. The mean drift of Alfvén waves may either increase or decrease the scattering center compression ratio, depending on the postshock cross-helicity, leading to either flatter or steeper CR spectra. We first examine such effects at planar shocks, based on the transport of Alfvén waves in the small amplitude limit. For the shock parameters relevant to cluster shocks, Alfvénic drift flattens the CR spectrum slightly, resulting in a small increase of the CR acceleration efficiency, η. We then consider two additional, physically motivated cases: (1) postshock waves are isotropized via MHD and plasma processes across the shock transition, and (2) postshock waves contain only forward waves propagating along with the flow due to a possible gradient of CR pressure behind the shock. In these cases, Alfvénic drift could reduce η by as much as a factor of five for weak cluster shocks. For the canonical parameters adopted here, we suggest η ∼ 10‑4–10‑2 for shocks with sonic Mach number M s ≈ 2–3. The possible reduction of η may help ease the tension between non-detection of γ-rays from galaxy clusters and DSA predictions.

Modeling and simulation of pressure waves generated by nano-thermite reactions

NASA Astrophysics Data System (ADS)

Martirosyan, Karen S.; Zyskin, Maxim; Jenkins, Charles M.; (Yuki) Horie, Yasuyuki

2012-11-01

This paper reports the modeling of pressure waves from the explosive reaction of nano-thermites consisting of mixtures of nanosized aluminum and oxidizer granules. Such nanostructured thermites have higher energy density (up to 26 kJ/cm3) and can generate a transient pressure pulse four times larger than that from trinitrotoluene (TNT) based on volume equivalence. A plausible explanation for the high pressure generation is that the reaction times are much shorter than the time for a shock wave to propagate away from the reagents region so that all the reaction energy is dumped into the gaseous products almost instantaneously and thereby a strong shock wave is generated. The goal of the modeling is to characterize the gas dynamic behavior for thermite reactions in a cylindrical reaction chamber and to model the experimentally measured pressure histories. To simplify the details of the initial stage of the explosive reaction, it is assumed that the reaction generates a one dimensional shock wave into an air-filled cylinder and propagates down the tube in a self-similar mode. Experimental data for Al/Bi2O3 mixtures were used to validate the model with attention focused on the ratio of specific heats and the drag coefficient. Model predictions are in good agreement with the measured pressure histories.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rodriguez, George; Gilbertson, Steve Michael

Chirped fiber Bragg grating (CFBG) sensors coupled to high speed interrogation systems are described as robust diagnostic approaches to monitoring shock wave and detonation front propagation tracking events for use in high energy density shock physics applications. Taking advantage of the linear distributed spatial encoding of the spectral band in single-mode CFBGs, embedded fiber systems and associated photonic interrogation methodologies are shown as an effective approach to sensing shock and detonation-driven loading processes along the CFBG length. Two approaches, one that detects spectral changes in the integrated spectrum of the CFBG and another coherent pulse interrogation approach that fully resolvesmore » its spectral response, shows that 100-MHz–1-GHz interrogation rates are possible with spatial resolution along the CFBG in the 50 µm to sub-millimeter range depending on the combination of CFBG parameters (i.e., length, chirp rate, spectrum) and interrogator design specifics. In conclusion, results from several dynamic tests are used to demonstrate the performance of these high speed systems for shock and detonation propagation tracking under strong and weak shock pressure loading: (1) linear detonation front tracking in the plastic bonded explosive (PBX) PBX-9501; (2) tracking of radial decaying shock with crossover to non-destructive CFBG response; (3) shock wave tracking along an aluminum cylinder wall under weak loading accompanied by dynamic strain effects in the CFBG sensor.« less

Entropy jump across an inviscid shock wave

NASA Technical Reports Server (NTRS)

Salas, Manuel D.; Iollo, Angelo

1995-01-01

The shock jump conditions for the Euler equations in their primitive form are derived by using generalized functions. The shock profiles for specific volume, speed, and pressure are shown to be the same, however density has a different shock profile. Careful study of the equations that govern the entropy shows that the inviscid entropy profile has a local maximum within the shock layer. We demonstrate that because of this phenomenon, the entropy, propagation equation cannot be used as a conservation law.

High Strain Rate and Shock-Induced Deformation in Metals

NASA Astrophysics Data System (ADS)

Ravelo, Ramon

2012-02-01

Large-scale non-equilibrium molecular Dynamics (MD) simulations are now commonly used to study material deformation at high strain rates (10^9-10^12 s-1). They can provide detailed information-- such as defect morphology, dislocation densities, and temperature and stress profiles, unavailable or hard to measure experimentally. Computational studies of shock-induced plasticity and melting in fcc and bcc single, mono-crystal metals, exhibit generic characteristics: high elastic limits, large directional anisotropies in the yield stress and pre-melting much below the equilibrium melt temperature for shock wave propagation along specific crystallographic directions. These generic features in the response of single crystals subjected to high strain rates of deformation can be explained from the changes in the energy landscape of the uniaxially compressed crystal lattice. For time scales relevant to dynamic shock loading, the directional-dependence of the yield strength in single crystals is shown to be due to the onset of instabilities in elastic-wave propagation velocities. The elastic-plastic transition threshold can accurately be predicted by a wave-propagation stability analysis. These strain-induced instabilities create incipient defect structures, which can be quite different from the ones, which characterize the long-time, asymptotic state of the compressed solid. With increase compression and strain rate, plastic deformation via extended defects gives way to amorphization associated with the loss in shear rigidity along specific deformation paths. The hot amorphous or (super-cooled liquid) metal re-crystallizes at rates, which depend on the temperature difference between the amorphous solid and the equilibrium melt line. This plastic-amorphous transition threshold can be computed from shear-waves stability analyses. Examples from selected fcc and bcc metals will be presented employing semi-empirical potentials of the embedded atom method (EAM) type as well as results from density functional theory calculations.

Studies of Shock Wave Interaction with a Curtain of Massive Particles

NASA Astrophysics Data System (ADS)

Lingampally, Sumanth Reddy; Wayne, Patrick; Cooper, Sean; Izard, Ricardo Gonzalez; Jacobs, Gustaaf; Vorobieff, Peter

2017-11-01

Interaction of a shock wave with planar and perturbed curtain of massive particles is studied experimentally. To form the curtain, solid soda lime particles (30-50 micron diameter) are dropped from a hopper fitted with mesh sieves and vibrated with a motor. The curtain forms when the particles move through a rectangular slot in the top of the test section of the shock tube used in experiment. The curtain can be either planar or perturbed in the horizontal plane (parallel to the shock direction) based on the shape of the slot. This setup generates a particle curtain with a volume fraction varying between 2 and 8 percent along its vertical height. A laser illuminates the curtain in vertical and horizontal planes. When the diaphragm separating the driver and the driven section is ruptured, shock waves with Mach numbers ranging from 1 to 2, depending on the pressure, propagate down the driven section and into test section. The phenomena following the shock wave impingement on the particle curtain are captured using an Apogee Alta U42 camera. This work is supported by the National Science Foundation Grant 1603915/1603326.

Observations of Upward Propagating Waves in the Transition Region and Corona above Sunspots

NASA Astrophysics Data System (ADS)

Hou, Zhenyong; Huang, Zhenghua; Xia, Lidong; Li, Bo; Fu, Hui

2018-03-01

We present observations of persistent oscillations of some bright features in the upper-chromosphere/transition region above sunspots taken by IRIS SJ 1400 Å and upward propagating quasi-periodic disturbances along coronal loops rooted in the same region taken by the AIA 171 Å passband. The oscillations of the features are cyclic oscillatory motions without any obvious damping. The amplitudes of the spatial displacements of the oscillations are about 1″. The apparent velocities of the oscillations are comparable to the sound speed in the chromosphere, but the upward motions are slightly larger than that of the downward. The intensity variations can take 24%–53% of the background, suggesting nonlinearity of the oscillations. The FFT power spectra of the oscillations show a dominant peak at a period of about 3 minutes, which is consistent with the omnipresent 3 minute oscillations in sunspots. The amplitudes of the intensity variations of the upward propagating coronal disturbances are 10%–15% of the background. The coronal disturbances have a period of about 3 minutes, and propagate upward along the coronal loops with apparent velocities in a range of 30 ∼ 80 km s‑1. We propose a scenario in which the observed transition region oscillations are powered continuously by upward propagating shocks, and the upward propagating coronal disturbances can be the recurrent plasma flows driven by shocks or responses of degenerated shocks that become slow magnetic-acoustic waves after heating the plasma in the coronal loops at their transition-region bases.

Generalized Sagdeev potential theory for shock waves modeling

NASA Astrophysics Data System (ADS)

Akbari-Moghanjoughi, M.

2017-05-01

In this paper, we develop an innovative approach to study the shock wave propagation using the Sagdeev potential method. We also present an analytical solution for Korteweg de Vries Burgers (KdVB) and modified KdVB equation families with a generalized form of the nonlinearity term which agrees well with the numerical one. The novelty of the current approach is that it is based on a simple analogy of the particle in a classical potential with the variable particle energy providing one with a deeper physical insight into the problem and can easily be extended to more complex physical situations. We find that the current method well describes both monotonic and oscillatory natures of the dispersive-diffusive shock structures in different viscous fluid configurations. It is particularly important that all essential parameters of the shock structure can be deduced directly from the Sagdeev potential in small and large potential approximation regimes. Using the new method, we find that supercnoidal waves can decay into either compressive or rarefactive shock waves depending on the initial wave amplitude. Current investigation provides a general platform to study a wide range of phenomena related to nonlinear wave damping and interactions in diverse fluids including plasmas.

Numerical Investigation of Shock Wave Propagation in Bone-Like Tissue

NASA Astrophysics Data System (ADS)

Nelms, Matt; Rajendran, Arunachalam

In this investigation, the effects of shock wave propagation in bone-like biomineralized tissue was investigated. The Alligator gar (Atractosteus spatula) exoskeleton is comprised of many disparate scales that provide a biological analog for potential design of flexible protective material systems. The penetration resistant fish scale was modeled by simulating a plate impact test configuration using ABAQUS®finite element (FE) software. The gar scale is identified as a two-phase, (1) hydroxyapatite mineral and (2) collagen protein, biological composite with two distinct layers where a stiff, ceramic-like ganoine overlays a soft, highly ductile bone. The geometry and variation of elastic modulus were determined from high-resolution scanning electron microscopy and dynamic nanoindentation experimentation to develop an idealized computational model for RVE-based FE simulations. The numerical analysis shows the effects of different functional material property variations on the stress histories and energy dissipation generated by wave propagation. Given the constitutive behaviors of the two layers are distinctly different, a brittle tensile damage model was employed to describe the ganoine and Drucker-Prager plasticity was used for the nonlinear response of the bone.

Nonlinear waveform distortion and shock formation in the near field of a continuous wave piston source

NASA Astrophysics Data System (ADS)

Sapozhnikov, Oleg A.; Khokhlova, Vera A.; Cathignol, Dominique

2004-05-01

A classical effect of nonlinear acoustics is that a plane sinusoidal acoustic wave propagating in a nonlinear medium transforms to a sawtooth wave with one shock per cycle. However, the waveform evolution can be quite different in the near field of a plane source due to diffraction. Previous numerical simulations of nonlinear acoustic waves in the near field of a circular piston source predict the development of two shocks per wave cycle [Khokhlova et al., J. Acoust. Soc. Am. 110, 95-108 (2001)]. Moreover, at some locations the peak pressure may be up to 4 times the source amplitude. The motivation of this work was to experimentally verify and further explain the phenomena of the nonlinear waveform distortion. Measurements were conducted in water with a 47-mm-diameter unfocused transducer, working at 1-MHz frequency. For pressure amplitudes higher than 0.5 MPa, two shocks per cycle were observed in the waveform beyond the last minimum of the fundamental harmonic amplitude. With the increase of the observation distance, these two shocks collided and formed one shock (per cycle), i.e., the waveform developed into the classical sawtooth wave. The experimental results were in a very good agreement with the modeling based on the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation.

Kidney damage in extracorporeal shock wave lithotripsy: a numerical approach for different shock profiles.

PubMed

Weinberg, Kerstin; Ortiz, Michael

2009-08-01

In shock-wave lithotripsy--a medical procedure to fragment kidney stones--the patient is subjected to hypersonic waves focused at the kidney stone. Although this procedure is widely applied, the physics behind this medical treatment, in particular the question of how the injuries to the surrounding kidney tissue arise, is still under investigation. To contribute to the solution of this problem, two- and three-dimensional numerical simulations of a human kidney under shock-wave loading are presented. For this purpose a constitutive model of the bio-mechanical system kidney is introduced, which is able to map large visco-elastic deformations and, in particular, material damage. The specific phenomena of cavitation induced oscillating bubbles is modeled here as an evolution of spherical pores within the soft kidney tissue. By means of large scale finite element simulations, we study the shock-wave propagation into the kidney tissue, adapt unknown material parameters and analyze the resulting stress states. The simulations predict localized damage in the human kidney in the same regions as observed in animal experiments. Furthermore, the numerical results suggest that in first instance the pressure amplitude of the shock wave impulse (and not so much its exact time-pressure profile) is responsible for damaging the kidney tissue.

Pressure wave propagation in fluid-filled co-axial elastic tubes. Part 1: Basic theory.

PubMed

Berkouk, K; Carpenter, P W; Lucey, A D

2003-12-01

Our work is motivated by ideas about the pathogenesis of syringomyelia. This is a serious disease characterized by the appearance of longitudinal cavities within the spinal cord. Its causes are unknown, but pressure propagation is probably implicated. We have developed an inviscid theory for the propagation of pressure waves in co-axial, fluid-filled, elastic tubes. This is intended as a simple model of the intraspinal cerebrospinal-fluid system. Our approach is based on the classic theory for the propagation of longitudinal waves in single, fluid-filled, elastic tubes. We show that for small-amplitude waves the governing equations reduce to the classic wave equation. The wave speed is found to be a strong function of the ratio of the tubes' cross-sectional areas. It is found that the leading edge of a transmural pressure pulse tends to generate compressive waves with converging wave fronts. Consequently, the leading edge of the pressure pulse steepens to form a shock-like elastic jump. A weakly nonlinear theory is developed for such an elastic jump.

Investigate the shock focusing under a single vortex disturbance using 2D Saint-Venant equations with a shock-capturing scheme

NASA Astrophysics Data System (ADS)

Zhao, Jiaquan; Li, Renfu; Wu, Haiyan

2018-02-01

In order to characterize the flow structure and the effect of acoustic waves caused by the shock-vortex interaction on the performance of the shock focusing, the incident plane shock wave with a single disturbance vortex focusing in a parabolic cavity is simulated systematically through solving the two-dimensional, unsteady Saint-Venant equations with the two order HLL scheme of Riemann solvers. The simulations show that the dilatation effect to be dominant in the net vorticity generation, while the baroclinic effect is dominate in the absence of initial vortex disturbance. Moreover, the simulations show that the time evolution of maximum focusing pressure with initial vortex is more complicate than that without initial vortex, which has a lot of relevance with the presence of quadrupolar acoustic wave structure induced by shock-vortex interaction and its propagation in the cavity. Among shock and other disturbance parameters, the shock Mach number, vortex Mach number and the shape of parabolic reflector proved to play a critical role in the focusing of shock waves and the strength of viscous dissipation, which in turn govern the evolution of maximum focusing pressure due to the gas dynamic focus, the change in dissipation rate and the coincidence of motion disturbance vortex with aerodynamic focus point.

Intracellular calcium and vulnerability to fibrillation and defibrillation in Langendorff-perfused rabbit ventricles.

PubMed

Hwang, Gyo-Seung; Hayashi, Hideki; Tang, Liang; Ogawa, Masahiro; Hernandez, Heidy; Tan, Alex Y; Li, Hongmei; Karagueuzian, Hrayr S; Weiss, James N; Lin, Shien-Fong; Chen, Peng-Sheng

2006-12-12

The role of intracellular calcium (Ca(i)) in defibrillation and vulnerability is unclear. We simultaneously mapped epicardial membrane potential and Ca(i) during shock on T-wave episodes (n=104) and attempted defibrillation episodes (n=173) in 17 Langendorff-perfused rabbit ventricles. Unsuccessful and type B successful defibrillation shocks were followed by heterogeneous distribution of Ca(i), including regions of low Ca(i) surrounded by elevated Ca(i) ("Ca(i) sinkholes") 31+/-12 ms after shock. The first postshock activation then originated from the Ca(i) sinkhole 53+/-14 ms after the shock. No sinkholes were present in type A successful defibrillation. A Ca(i) sinkhole also was present 39+/-32 ms after a shock on T that induced ventricular fibrillation, followed 22+/-15 ms later by propagated wave fronts that arose from the same site. This wave propagated to form a spiral wave and initiated ventricular fibrillation. Thapsigargin and ryanodine significantly decreased the upper limit of vulnerability and defibrillation threshold. We studied an additional 7 rabbits after left ventricular endocardial cryoablation, resulting in a thin layer of surviving epicardium. Ca(i) sinkholes occurred 31+/-12 ms after the shock, followed in 19+/-7 ms by first postshock activation in 63 episodes of unsuccessful defibrillation. At the Ca(i) sinkhole, the rise of Ca(i) preceded the rise of epicardial membrane potential in 5 episodes. There is a heterogeneous postshock distribution of Ca(i). The first postshock activation always occurs from a Ca(i) sinkhole. The Ca(i) prefluorescence at the first postshock early site suggests that reverse excitation-contraction coupling might be responsible for the initiation of postshock activations that lead to ventricular fibrillation.

Predictions of lithium interactions with earth's bow shock in the presence of wave activity

NASA Technical Reports Server (NTRS)

Decker, R. B.; Lui, A. T. Y.; Vlahos, L.

1984-01-01

The results of a test-particle simulation studying the movement of a lithium tracer ion injected upstream of the bow shock are reported. Wave activity consists of parallel and antiparallel propagating Alfven waves characterized by a frequency power spectrum within a frequency or range of amplitudes defined separately in the upstream and downstream regions. The results show that even a moderate level of wave activity can substantially change the results obtained in the absence of waves. Among the effects observed are: (1) increased ion transmission; (2) both the average energy gain and spread about the average are increased for transmitted and reflected particles; (3) the average final pitch angle for transmitted particles tends to 90 deg, and the spread of reflected particles is reduced; and (4) the spatial dispersion of the ions on the bow shock after a single encounter is increased.

A generic efficient adaptive grid scheme for rocket propulsion modeling

NASA Technical Reports Server (NTRS)

Mo, J. D.; Chow, Alan S.

1993-01-01

The objective of this research is to develop an efficient, time-accurate numerical algorithm to discretize the Navier-Stokes equations for the predictions of internal one-, two-dimensional and axisymmetric flows. A generic, efficient, elliptic adaptive grid generator is implicitly coupled with the Lower-Upper factorization scheme in the development of ALUNS computer code. The calculations of one-dimensional shock tube wave propagation and two-dimensional shock wave capture, wave-wave interactions, shock wave-boundary interactions show that the developed scheme is stable, accurate and extremely robust. The adaptive grid generator produced a very favorable grid network by a grid speed technique. This generic adaptive grid generator is also applied in the PARC and FDNS codes and the computational results for solid rocket nozzle flowfield and crystal growth modeling by those codes will be presented in the conference, too. This research work is being supported by NASA/MSFC.

Evolution of wave patterns and temperature field in shock-tube flow

NASA Astrophysics Data System (ADS)

Kiverin, A. D.; Yakovenko, I. S.

2018-05-01

The paper is devoted to the numerical analysis of wave patterns behind a shock wave propagating in a tube filled with a gaseous mixture. It is shown that the flow inside the boundary layer behind the shock wave is unstable, and the way the instability develops fully corresponds to the solution obtained for the boundary layer over a flat plate. Vortical perturbations inside the boundary layer determine the nonuniformity of the temperature field. In turn, exactly these nonuniformities define the way the ignition kernels arise in the combustible mixture after the reflected shock interaction with the boundary layer. In particular, the temperature nonuniformity determines the spatial limitations of probable ignition kernel position relative to the end wall and side walls of the tube. In the case of low-intensity incident shocks the ignition could start not farther than the point of first interaction between the reflected shock wave and roller vortices formed in the process of boundary layer development. Proposed physical mechanisms are formulated in general terms and can be used for interpretation of the experimental data in any systems with a delayed exothermal reaction start. It is also shown that contact surface thickening occurs due to its interaction with Tollmien-Schlichting waves. This conclusion is of importance for understanding the features of ignition in shock tubes operating in the over-tailored regime.

Simultaneous Observations of a Large-scale Wave Event in the Solar Atmosphere: From Photosphere to Corona

NASA Astrophysics Data System (ADS)

Shen, Yuandeng; Liu, Yu

2012-06-01

For the first time, we report a large-scale wave that was observed simultaneously in the photosphere, chromosphere, transition region, and low corona layers of the solar atmosphere. Using the high temporal and high spatial resolution observations taken by the Solar Magnetic Activity Research Telescope at Hida Observatory and the Atmospheric Imaging Assembly (AIA) on board Solar Dynamic Observatory, we find that the wave evolved synchronously at different heights of the solar atmosphere, and it propagated at a speed of 605 km s-1 and showed a significant deceleration (-424 m s-2) in the extreme-ultraviolet (EUV) observations. During the initial stage, the wave speed in the EUV observations was 1000 km s-1, similar to those measured from the AIA 1700 Å (967 km s-1) and 1600 Å (893 km s-1) observations. The wave was reflected by a remote region with open fields, and a slower wave-like feature at a speed of 220 km s-1 was also identified following the primary fast wave. In addition, a type-II radio burst was observed to be associated with the wave. We conclude that this wave should be a fast magnetosonic shock wave, which was first driven by the associated coronal mass ejection and then propagated freely in the corona. As the shock wave propagated, its legs swept the solar surface and thereby resulted in the wave signatures observed in the lower layers of the solar atmosphere. The slower wave-like structure following the primary wave was probably caused by the reconfiguration of the low coronal magnetic fields, as predicted in the field-line stretching model.

A new shock-capturing numerical scheme for ideal hydrodynamics

NASA Astrophysics Data System (ADS)

Fecková, Z.; Tomášik, B.

2015-05-01

We present a new algorithm for solving ideal relativistic hydrodynamics based on Godunov method with an exact solution of Riemann problem for an arbitrary equation of state. Standard numerical tests are executed, such as the sound wave propagation and the shock tube problem. Low numerical viscosity and high precision are attained with proper discretization.

DETECTION OF SHOCK MERGING IN THE CHROMOSPHERE OF A SOLAR PORE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chae, Jongchul; Song, Donguk; Seo, Minju

2015-06-01

It was theoretically demonstrated that a shock propagating in the solar atmosphere can overtake another and merge with it. We provide clear observational evidence that shock merging does occur quite often in the chromosphere of sunspots. Using Hα imaging spectral data taken by the Fast Imaging Solar Spectrograph of the 1.6 m New Solar Telescope at the Big Bear Soar Observatory, we construct time–distance maps of line-of-sight velocities along two appropriately chosen cuts in a pore. The maps show a number of alternating redshift and blueshift ridges, and we identify each interface between a preceding redshift ridge and the followingmore » blueshift ridge as a shock ridge. The important finding of ours is that two successive shock ridges often merge with each other. This finding can be theoretically explained by the merging of magneto-acoustic shock waves propagating with lower speeds of about 10 km s{sup −1} and those propagating at higher speeds of about 16–22 km s{sup −1}. The shock merging is an important nonlinear dynamical process of the solar chromosphere that can bridge the gap between higher-frequency chromospheric oscillations and lower-frequency dynamic phenomena such as fibrils.« less

Dislocation pinning effects induced by nano-precipitates during warm laser shock peening: Dislocation dynamic simulation and experiments

NASA Astrophysics Data System (ADS)

Liao, Yiliang; Ye, Chang; Gao, Huang; Kim, Bong-Joong; Suslov, Sergey; Stach, Eric A.; Cheng, Gary J.

2011-07-01

Warm laser shock peening (WLSP) is a new high strain rate surface strengthening process that has been demonstrated to significantly improve the fatigue performance of metallic components. This improvement is mainly due to the interaction of dislocations with highly dense nanoscale precipitates, which are generated by dynamic precipitation during the WLSP process. In this paper, the dislocation pinning effects induced by the nanoscale precipitates during WLSP are systematically studied. Aluminum alloy 6061 and AISI 4140 steel are selected as the materials with which to conduct WLSP experiments. Multiscale discrete dislocation dynamics (MDDD) simulation is conducted in order to investigate the interaction of dislocations and precipitates during the shock wave propagation. The evolution of dislocation structures during the shock wave propagation is studied. The dislocation structures after WLSP are characterized via transmission electron microscopy and are compared with the results of the MDDD simulation. The results show that nano-precipitates facilitate the generation of highly dense and uniformly distributed dislocation structures. The dislocation pinning effect is strongly affected by the density, size, and space distribution of nano-precipitates.

Acoustic and mechanical properties of renal calculi: implications in shock wave lithotripsy.

PubMed

Chuong, C J; Zhong, P; Preminger, G M

1993-12-01

The acoustic and mechanical properties of renal calculi dictate how a stone interacts with the mechanical forces produced by shock wave lithotripsy; thus, these properties are directly related to the success of the treatment. Using an ultrasound pulse transmission technique, we measured both longitudinal and transverse (or shear) wave propagation speeds in nine groups of renal calculi with different chemical compositions. We also measured stone density using a pycnometer based on Archimedes' principle. From these measurements, we calculated wave impedance and dynamic mechanical properties of the renal stones. Calcium oxalate monohydrate and cystine stones had higher longitudinal and transverse wave speeds, wave impedances, and dynamic moduli (bulk modulus, Young's modulus, and shear modulus), suggesting that these stones are more difficult to fragment. Phosphate stones (carbonate apatite and magnesium ammonium phosphate hydrogen) were found to have lower values of these properties, suggesting they are more amenable to shock wave fragmentation. These data provide a physical explanation for the significant differences in stone fragility observed clinically.

On the dynamics of a shock-bubble interaction

NASA Technical Reports Server (NTRS)

Quirk, James J.; Karni, Smadar

1994-01-01

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Shock-wave studies of anomalous compressibility of glassy carbon

DOE Office of Scientific and Technical Information (OSTI.GOV)

Molodets, A. M., E-mail: molodets@icp.ac.ru; Golyshev, A. A.; Savinykh, A. S.

2016-02-15

The physico-mechanical properties of amorphous glassy carbon are investigated under shock compression up to 10 GPa. Experiments are carried out on the continuous recording of the mass velocity of compression pulses propagating in glassy carbon samples with initial densities of 1.502(5) g/cm{sup 3} and 1.55(2) g/cm{sup 3}. It is shown that, in both cases, a compression wave in glassy carbon contains a leading precursor with amplitude of 0.135(5) GPa. It is established that, in the range of pressures up to 2 GPa, a shock discontinuity in glassy carbon is transformed into a broadened compression wave, and shock waves are formedmore » in the release wave, which generally means the anomalous compressibility of the material in both the compression and release waves. It is shown that, at pressure higher than 3 GPa, anomalous behavior turns into normal behavior, accompanied by the formation of a shock compression wave. In the investigated area of pressure, possible structural changes in glassy carbon under shock compression have a reversible character. A physico-mechanical model of glassy carbon is proposed that involves the equation of state and a constitutive relation for Poisson’s ratio and allows the numerical simulation of physico-mechanical and thermophysical properties of glassy carbon of different densities in the region of its anomalous compressibility.« less

Photographic laboratory studies of explosions.

NASA Technical Reports Server (NTRS)

Kamel, M. M.; Oppenheim, A. K.

1973-01-01

Description of a series of cinematographic studies of explosions made with a high-speed rotating-mirror streak camera which uses a high-frequency stroboscopic ruby laser as the light source. The results obtained mainly concern explosions initiated by focused laser irradiation from a pulsed neodymium laser in a detonating gas consisting essentially of an equimolar mixture of acetylene and oxygen at an initial pressure of 100 torr at room temperature. Among the most significant observations were observations of a spherical blast wave preceded by a Chapman-Jouguet detonation which is stabilized immediately after initiation, the merging of a spherical flame with a shock front of the blast wave in which the flame is propagating, the division of a spherical detonation front into a shock wave and flame, and the generation of shock waves by a network of spherical flames.

Shock interactions with heterogeneous energetic materials

NASA Astrophysics Data System (ADS)

Yarrington, Cole D.; Wixom, Ryan R.; Damm, David L.

2018-03-01

The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet, the majority of computational studies aimed at predicting phenomena affected by these processes, such as the initiation and propagation of detonation waves in explosives or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed. The measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of the microstructure along with a fully dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on the dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide a clear insight into the nature of threshold behavior and are a way to understand complex physical phenomena.

Shock interactions with heterogeneous energetic materials

DOE PAGES

Yarrington, Cole D.; Wixom, Ryan R.; Damm, David L.

2018-03-14

The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet the majority of computational studies aimed at predicting phenomena affected by these processes, such as initiation and propagation of detonation waves in explosives, or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed.more » Measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics (DFT-MD) derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of microstructure along with a fully-dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide clear insight into the nature of threshold behavior, and are a way to understand complex physical phenomena.« less

On meteor-generated infrasound. [propagation characteristics during entry into earth atmosphere

NASA Technical Reports Server (NTRS)

Revelle, D. O.

1976-01-01

The characteristics of generation and propagation of infrasonic pressure waves excited during meteor entry into the earth's atmosphere are studied. Existing line source blast wave theory is applied to infrasonic airwave data from four bright fire-balls. It is shown that the strong shock behavior of the blast wave is confined to a cylinderical region with a radius proportional to the product of the meteor Mach number and its diameter. A description of the wave form far from the source is provided. Infrasonic data reported elsewhere are analyzed. All the results should be considered as preliminary, and additional work is under way to refine the estimates obtained.

THE COUPLED EVOLUTION OF ELECTRONS AND IONS IN CORONAL MASS EJECTION-DRIVEN SHOCKS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manchester IV, W. B.; Van der Holst, B.; Toth, G.

2012-09-01

We present simulations of coronal mass ejections (CMEs) performed with a new two-temperature coronal model developed at the University of Michigan, which is able to address the coupled thermodynamics of the electron and proton populations in the context of a single fluid. This model employs heat conduction for electrons, constant adiabatic index ({gamma} = 5/3), and includes Alfven wave pressure to accelerate the solar wind. The Wang-Sheeley-Arge empirical model is used to determine the Alfven wave pressure necessary to produce the observed bimodal solar wind speed. The Alfven waves are dissipated as they propagate from the Sun and heat protonsmore » on open magnetic field lines to temperatures above 2 MK. The model is driven by empirical boundary conditions that includes GONG magnetogram data to calculate the coronal field, and STEREO/EUVI observations to specify the density and temperature at the coronal boundary by the Differential Emission Measure Tomography method. With this model, we simulate the propagation of fast CMEs and study the thermodynamics of CME-driven shocks. Since the thermal speed of the electrons greatly exceeds the speed of the CME, only protons are directly heated by the shock. Coulomb collisions low in the corona couple the protons and electrons allowing heat exchange between the two species. However, the coupling is so brief that the electrons never achieve more than 10% of the maximum temperature of the protons. We find that heat is able to conduct on open magnetic field lines and rapidly propagates ahead of the CME to form a shock precursor of hot electrons.« less

Mechanism of Gaseous Detonation Propagation Through Reactant Layers Bounded by Inert Gas

NASA Astrophysics Data System (ADS)

Houim, Ryan

2017-11-01

Vapor cloud explosions and rotating detonation engines involve the propagation of gaseous detonations through a layer of reactants that is bounded by inert gas. Mechanistic understanding of how detonations propagate stably or fail in these scenarios is incomplete. Numerical simulations were used to investigate mechanisms of gaseous detonation propagation through reactant layers bounded by inert gas. The reactant layer was a stoichiometric mixture of C2H4/O2 at 1 atm and 300K and is 4 detonation cells in height. Cases where the inert gas temperature was 300, 1500, and 3500 K will be discussed. The detonation failed for the 300 K case and propagated marginally for the 1500 K case. Surprisingly, the detonation propagated stably for the 3500 K case. A shock structure forms that involves a detached shock in the inert gas and a series of oblique shocks in the reactants. A small local explosion is triggered when the Mach stem of a detonation cell interacts with the compressed reactants behind one of these oblique shocks. The resulting pressure wave produces a new Mach stem and a new triple point that leads to a stable detonation. Preliminary results on the influence of a deflagration at the inert/reactant interface on the stability of a layered detonation will be discussed.

Constraining the source location of the 30 May 2015 (Mw 7.9) Bonin deep-focus earthquake using seismogram envelopes of high-frequency P waveforms: Occurrence of deep-focus earthquake at the bottom of a subducting slab

NASA Astrophysics Data System (ADS)

Takemura, Shunsuke; Maeda, Takuto; Furumura, Takashi; Obara, Kazushige

2016-05-01

In this study, the source location of the 30 May 2015 (Mw 7.9) deep-focus Bonin earthquake was constrained using P wave seismograms recorded across Japan. We focus on propagation characteristics of high-frequency P wave. Deep-focus intraslab earthquakes typically show spindle-shaped seismogram envelopes with peak delays of several seconds and subsequent long-duration coda waves; however, both the main shock and aftershock of the 2015 Bonin event exhibited pulse-like P wave propagations with high apparent velocities (~12.2 km/s). Such P wave propagation features were reproduced by finite-difference method simulations of seismic wave propagation in the case of slab-bottom source. The pulse-like P wave seismogram envelopes observed from the 2015 Bonin earthquake show that its source was located at the bottom of the Pacific slab at a depth of ~680 km, rather than within its middle or upper regions.

Use of high-speed visualization for the study of shock-wave interactions with deformable porous materials

NASA Astrophysics Data System (ADS)

Skews, Beric W.; Glick, Gavin; Doyle, Graham K.; Lamond, Paul W.

1997-05-01

This paper describes the use of high-speed photography, and videography, in the study of material distortion and movement when a shock wave traverses a highly deformable porous structure, such as a blob of foam or a porous bed of particles. The effects of surface porosity can be significant in determining the nature of reflection of shock waves from surfaces. Not only are wave geometries substantially modified but the resulting wall pressures are also strongly affected. It, in addition, the surface is highly deformable by being made up of an elastic matrix or a collection of discrete particles, then the reflection geometry and loading can be even more complex. It is known, for example, that shock wave impact on open-cell polyurethane foam attached to a wall can cause a significant increase in pressure on the wall compared to reflection off a plane rigid wall without covering. The motion of the interface is an essential consideration in understanding the dynamics of these interactions. These studies could have application to the effects of blast wave propagation over complex surfaces such as forests, grasslands, and snow; as well as in establishing the efficacy of safety padding and attenuation materials under shock and impact loading conditions. Studies on an assortment of materials are presented, using a variety of visualization techniques. Recording methods used range from short duration flash photography (both shadow and schlieren), through multi-frame videography; to single frame, multi-exposure video capture with a camera capable of rates up to 1 million pictures per second. In the case of shock wave impact on specimens of polyurethane foam, the results clearly show the expulsion and reingestion of shock heated gas from within the foam body as the material collapses and then recovers, coupled with longitudinal and transverse oscillations of the body of the foam material. For blast wave propagation over porous beds, occurrence of particle lift off, bed fluidization, and the generation of surface dunes are evident. The recordings allow the calculation of the velocities and accelerations of the various interfaces and particles to be made, using suitable image processing techniques. Thus, estimates may be made of the unsteady drag forces acting on the individual particles.

Field experiments to determine wave propagation principles and mechanical properties of snow

NASA Astrophysics Data System (ADS)

Simioni, Stephan; Gebhard, Felix; Dual, Jürg; Schweizer, Jürg

2017-04-01

To understand the release of snow avalanches by explosions one needs to know how acoustic waves travel above and within the snowpack. Hitherto, wave propagation was investigated in the laboratory with small samples or in the field in the shock wave region. We developed a measurement system and layout to derive wave attenuation in snow, wave speeds and elastic moduli on small-scale (1-2 m) field experiments to close the gap between the lab scale (0.1 m) and the scale of artificial release (10-100 m). We used solid explosives and hammer blows to create the load and accelerometers to measure the resulting wave within the snowpack. The strong attenuation we observed indicates that we measured the second longitudinal wave which propagates through the pore space. The wave speeds, however, corresponded to the speeds of the first longitudinal wave within the ice skeleton. The elastic moduli were high on the order of several tens of MPa for lower densities (150 kg m-3) and agreed well with earlier lab studies, in particular for the higher densities 250-400 kg m-3). However, the scatter was rather large as expected for in-situ experiments in the layered snow cover. In addition, we measured accelerations during propagation saw test experiments. The propagation of cracks during this type of snow instability test has mainly been studied by analysing the bending of the slab (due to the saw cut) using particle tracking velocimetry. We used the accelerometers to measure crack propagation speeds. The wave speeds were slightly higher for most experiments than reported previously. Furthermore, in some experiments, we encountered to different wave types with one propagating at a higher speed. This finding may be interpreted as the actual crack propagation and the settling of the weak layer (collapse wave). Our results show that field measurements of propagation properties are feasible and that crack propagation as observed during propagation saw tests may involve different processes that need to be further investigated.

Experiment and analysis of shock waves radiated from pulse laser focusing in a gelatin gel

NASA Astrophysics Data System (ADS)

Nakamura, Nobuyuki; Ando, Keita

2017-11-01

A fundamental understanding of shock and bubble dynamics in human tissues is essential to laser application for medical purposes. Here, we experimentally study the dynamics of shock waves in viscoelastic media. A nanosecond laser pulse of wavelength at 532 nm and of energy up to 2.66 +/- 0.09 mJ was focused through a microscope objective lens (10 x, NA = 0.30) into a gel of gelatin concentration at 3 and 10 wt%; a shock wave and a bubble can be generated, respectively, by rapid expansion of the laser-induced plasma and local heat deposition after the plasma recombines. The shock propagation and the bubble growth were recorded by a ultra-high-speed camera at 100 Mfps. The shock evolution was determined by image analysis of the recording and the shock pressure in the near field was computed according to the Rankine-Hugoniot relation. The far-field pressure was measured by a hydrophone. In the poster, we will present the decay rate of the shock pressure in the near and far fields and examine viscous effects on the shock dynamics. The Research Grant of Keio Leading-edge Laboratory of Science & Technology.

Trajectory Control of Small Rotating Projectiles by Laser Sparks

NASA Astrophysics Data System (ADS)

Starikovskiy, Andrey; Limbach, Christopher; Miles, Richard

2015-09-01

The possibility of controlling the trajectory of the supersonic motion of a rotating axisymmetric projectile using a remotely generated laser spark was investigated. The dynamic images of the interaction of thermal inhomogeneity created by the laser spark with the bow shock in front of the projectile were obtained. The criterion for a strong shock wave interaction with the thermal inhomogeneity at different angles of a shock wave was derived. Significant changes in the configuration of the bow shock wave and changes in the pressure distribution over the surface of the rotating projectile can appear for laser spark temperature of T' = 2500-3000 K. The experiment showed that strong interaction takes place for both plane and oblique shock waves. The measurement of the velocity of the precession of the rotating projectile axis from the initial position in time showed that the angle of attack of the projectile deviates with a typical time of perturbation propagation along the projectile's surface. Thus the laser spark can change the trajectory of the rotating projectile, moving at supersonic speed, through the creation of thermal heterogeneity in front of it.

Modeling of plasticity and fracture of metals at shock loading

NASA Astrophysics Data System (ADS)

Mayer, A. E.; Khishchenko, K. V.; Levashov, P. R.; Mayer, P. N.

2013-05-01

In this paper, we present a model of dislocation plasticity and fracture of metals, which in combination with the wide-range equation of state and the continuum mechanics equations is a necessary component for simulation of the shock-wave loading. We take into account immobilization of dislocations and nucleation of micro-voids in weakened zones of substance; this is distinguished feature of the present version of the model. Accounting of the dislocations immobilization provides a better description of the unloading wave structure, while the detailed consideration of processes in the weakened zones expands the domain of applicability of fracture model to higher strain rates. We compare our results with the experimental data for the shock loading of aluminum, copper, and nickel samples; the comparison indicates satisfactory description of the elastic precursor, unloading wave, and spall pulse. Using the model, we investigate intently the early stage of the shock formation in solids; it is found out that the elastic precursor is formed even for a strong shock wave, and initially the precursor has very large amplitude and propagation velocity.

Development of a shock wave adhesion test for composite bonds by pulsed laser and mechanical impacts

NASA Astrophysics Data System (ADS)

Ecault, R.; Boustie, M.; Touchard, F.; Arrigoni, M.; Berthe, L.

2014-05-01

Evaluating the bonding quality of composite material is becoming one of the main challenges faced by aeronautic industries. This work aims to the development of a technique using shock wave, which would enable to quantify the bonding mechanical quality. Laser shock experiments were carried out. This technique enables high tensile stress generation in the thickness of composite bonds. The resulting damage has been quantified using different methods such as confocal microscopy, ultrasound and cross section observation. The discrimination between a correct bond and a weak bond was possible thanks to these experiments. Nevertheless, laser sources are not well adapted for optimization of such a test because of often fixed settings. That is why mechanical impacts on bonded composites were also performed in this work. By changing the thickness of aluminum projectiles, the generated tensile stresses by the shock wave propagation were moved toward the composite/bond interface. The made observations prove that the technique optimization is possible. The key parameters for the development of a bonding test using shock waves have been identified.

Development of a shock wave adhesion test for composite bonds by laser pulsed and mechanical impacts

NASA Astrophysics Data System (ADS)

Ecault, Romain; Boustie, Michel; Touchard, Fabienne; Arrigoni, Michel; Berthe, Laurent; CNRS Collaboration

2013-06-01

Evaluating the bonding quality of composite material is becoming one of the main challenges faced by aeronautic industries. This work aims the development of a technique using shock wave, which would enable to quantify the bonding mechanical quality. Laser shock experiments were carried out. This technique enables high tensile stress generation in the thickness of composite bond without any mechanical contact. The resulting damage has been quantified using different method such as confocal microscopy, ultrasound and cross section observation. The discrimination between a correct bond and a weak bond was possible thanks to these experiments. Nevertheless, laser sources are not well adapted for optimization of such a test since it has often fixed parameters. That is why mechanical impacts bonded composites were also performed in this work. By changing the thickness of aluminum projectiles, the tensile stresses generated by the shock wave propagation were moved toward the composite/bond interface. The observations made prove that the optimization of the technique is possible. The key parameters for the development of a bonding test using shock wave have been identified.

Simple and accurate theory for strong shock waves in a dense hard-sphere fluid.

PubMed

Montanero, J M; López de Haro, M; Santos, A; Garzó, V

1999-12-01

Following an earlier work by Holian et al. [Phys. Rev. E 47, R24 (1993)] for a dilute gas, we present a theory for strong shock waves in a hard-sphere fluid described by the Enskog equation. The idea is to use the Navier-Stokes hydrodynamic equations but taking the temperature in the direction of shock propagation rather than the actual temperature in the computation of the transport coefficients. In general, for finite densities, this theory agrees much better with Monte Carlo simulations than the Navier-Stokes and (linear) Burnett theories, in contrast to the well-known superiority of the Burnett theory for dilute gases.

ON THE CONNECTION BETWEEN PROPAGATING SOLAR CORONAL DISTURBANCES AND CHROMOSPHERIC FOOTPOINTS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bryans, P.; McIntosh, S. W.; Moortel, I. De

2016-09-20

The Interface Region Imaging Spectrograph ( IRIS ) provides an unparalleled opportunity to explore the (thermal) interface between the chromosphere, transition region, and the coronal plasma observed by the Atmospheric Imaging Assembly (AIA) of the Solar Dynamics Observatory ( SDO ). The SDO /AIA observations of coronal loop footpoints show strong recurring upward propagating signals—“propagating coronal disturbances” (PCDs) with apparent speeds of the order of 100–120 km s{sup −1}. That signal has a clear signature in the slit-jaw images of IRIS in addition to identifiable spectral signatures and diagnostics in the Mg iih (2803 Å) line. In analyzing the Mgmore » iih line, we are able to observe the presence of magnetoacoustic shock waves that are also present in the vicinity of the coronal loop footpoints. We see there is enough of a correspondence between the shock propagation in Mg iih, the evolution of the Si iv line profiles, and the PCD evolution to indicate that these waves are an important ingredient for PCDs. In addition, the strong flows in the jet-like features in the IRIS Si iv slit-jaw images are also associated with PCDs, such that waves and flows both appear to be contributing to the signals observed at the footpoints of PCDs.« less

Implicit approximate-factorization schemes for the low-frequency transonic equation

NASA Technical Reports Server (NTRS)

Ballhaus, W. F.; Steger, J. L.

1975-01-01

Two- and three-level implicit finite-difference algorithms for the low-frequency transonic small disturbance-equation are constructed using approximate factorization techniques. The schemes are unconditionally stable for the model linear problem. For nonlinear mixed flows, the schemes maintain stability by the use of conservatively switched difference operators for which stability is maintained only if shock propagation is restricted to be less than one spatial grid point per time step. The shock-capturing properties of the schemes were studied for various shock motions that might be encountered in problems of engineering interest. Computed results for a model airfoil problem that produces a flow field similar to that about a helicopter rotor in forward flight show the development of a shock wave and its subsequent propagation upstream off the front of the airfoil.

Simulations of viscous and compressible gas-gas flows using high-order finite difference schemes

NASA Astrophysics Data System (ADS)

Capuano, M.; Bogey, C.; Spelt, P. D. M.

2018-05-01

A computational method for the simulation of viscous and compressible gas-gas flows is presented. It consists in solving the Navier-Stokes equations associated with a convection equation governing the motion of the interface between two gases using high-order finite-difference schemes. A discontinuity-capturing methodology based on sensors and a spatial filter enables capturing shock waves and deformable interfaces. One-dimensional test cases are performed as validation and to justify choices in the numerical method. The results compare well with analytical solutions. Shock waves and interfaces are accurately propagated, and remain sharp. Subsequently, two-dimensional flows are considered including viscosity and thermal conductivity. In Richtmyer-Meshkov instability, generated on an air-SF6 interface, the influence of the mesh refinement on the instability shape is studied, and the temporal variations of the instability amplitude is compared with experimental data. Finally, for a plane shock wave propagating in air and impacting a cylindrical bubble filled with helium or R22, numerical Schlieren pictures obtained using different grid refinements are found to compare well with experimental shadow-photographs. The mass conservation is verified from the temporal variations of the mass of the bubble. The mean velocities of pressure waves and bubble interface are similar to those obtained experimentally.

The Shock and Vibration Digest. Volume 12, Number 3.

DTIC Science & Technology

1980-03-01

this problem by Mallik to design the fuselage so that it acts as a band pass and M-d [211. Two :.ipes of support were con- filter, filtering out the...370-373 (1975). 237-245 (1975). 11. Harari, A., "Wave Propagation in Cylindrical 21. Mallik , A.K. and Mead, D.J., "Free Vibration Shells with Finite...1183 (1973). 29 _ _ 23. Singh, K. and Mallik , A.K., "Wave Propagation sure Fields," J. Sound Vib., 28 (2), pp 247- and Vibration Response of a

A computational study of laser-supported detonation waves propagating up an oblique incident beam

NASA Astrophysics Data System (ADS)

Bohn, C. L.; Crawford, M. L.

1987-01-01

A series of numerical experiments was conducted to study the propagation of laser-supported detonation waves (LSDWs) in the case that a CO2 laser beam strikes an aluminum surface obliquely in air. A reflected shock formed at the aluminum surface was more prominent at higher angles of incidence theta of the beam, but otherwise the hydrodynamics of the plasma and the LSDW were insensitive to theta. Furthermore, the total impulse delivered to the aluminum varied approximately as 1/cos theta, a result that can be modeled with elementary blast-wave theory.

Tunable evolutions of shock absorption and energy partitioning in magnetic granular chains

NASA Astrophysics Data System (ADS)

Leng, Dingxin; Liu, Guijie; Sun, Lingyu

2018-01-01

In this paper, we investigate the tunable characteristics of shock waves propagating in one-dimensional magnetic granular chains at various chain lengths and magnetic flux densities. According to the Hertz contact theory and Maxwell principle, a discrete element model with coupling elastic and field-induced interaction potentials of adjacent magnetic grains is proposed. We also present hard-sphere approximation analysis to describe the energy partitioning features of magnetic granular chains. The results demonstrate that, for a fixed magnetic field strength, when the chain length is greater than two times of the wave width of the solitary wave, the chain length has little effect on the output energy of the system; for a fixed chain length, the shock absorption and energy partitioning features of magnetic granular chains are remarkably influenced by varying magnetic flux densities. This study implies that the magnetic granular chain is potential to construct adaptive shock absorption components for impulse mitigation.

Shock induced Richtmyer-Meshkov instability in the presence of a wall boundary layer

NASA Astrophysics Data System (ADS)

Jourdan, G.; Billiotte, M.; Houas, L.

1996-06-01

An experimental investigation on gaseous mixing zones originated from the Richtmyer-Meshkov instability has been undertaken in a square cross section shock tube. Mass concentration fields, of one of the two mixing constituents, have been determined within the mixing zone when the shock wave passes from the heavy gas to the light one, from one gas to an other of close density, and from the light gas to the heavy one. Results have been obtained before and after the coming back of the reflected shock wave. The diagnostic method is based on the infrared absorption of one of the two constituents of the mixing zone. It is shown that the mixing zone is strongly deformed by the wall boundary layer. The consequence is the presence of strong gradients of concentration in the direction perpendicular to the shock wave propagation. Finally, it is pointed out that the mixing goes more homogeneous when the Atwood number tends to zero.

Dynamic response of phenolic resin and its carbon-nanotube composites to shock wave loading

DOE PAGES

Arman, B.; An, Q.; Luo, S. N.; ...

2011-01-04

We investigate with nonreactive molecular dynamics simulations the dynamic response of phenolic resin and its carbon-nanotube (CNT) composites to shock wave compression. For phenolic resin, our simulations yield shock states in agreement with experiments on similar polymers except the “phase change” observed in experiments, indicating that such phase change is chemical in nature. The elastic–plastic transition is characterized by shear stress relaxation and atomic-level slip, and phenolic resin shows strong strain hardening. Shock loading of the CNT-resin composites is applied parallel or perpendicular to the CNT axis, and the composites demonstrate anisotropy in wave propagation, yield and CNT deformation. Themore » CNTs induce stress concentrations in the composites and may increase the yield strength. Our simulations indicate that the bulk shock response of the composites depends on the volume fraction, length ratio, impact cross-section, and geometry of the CNT components; the short CNTs in current simulations have insignificant effect on the bulk response of resin polymer.« less

Peculiarities of evolutions of elastic-plastic shock compression waves in different materials

NASA Astrophysics Data System (ADS)

Kanel, G. I.; Savinykh, A. S.; Garkushin, G. V.; Razorenov, S. V.; Ashitkov, S. I.; Zaretsky, E. B.

2016-11-01

In the paper, we discuss such unexpected features in the wave evolution in solids as strongly nonlinear uniaxial elastic compression in a picosecond time range, a departure from self-similar development of the wave process which is accompanied with apparent sub-sonic wave propagation, changes of shape of elastic precursor wave as a result of variations in the material structure and the temperature, unexpected peculiarities of reflection of elastic-plastic waves from free surface.

Shock Propagation In Crustal Rock

DTIC Science & Technology

1991-04-29

liquid produced above the melting point during unloading. Figure 5 displays calculations of release adiabats in the mixed phase regime based on the...muscovite [Bridgman, 1949]. The zero-pressure densities at points 1, 2 and 3 correspond to mixtures of orthoclase + A1203+H20, of wadeite + kyanite ...shocked tantalum and the high pressure melting point , in Shock Waves in Condensed Matter-1983, edited by J. R. Asay, R. A. Graham and G. K. Straub, pp. 91

Optical diagnostics of turbulent mixing in explosively-driven shock tube

NASA Astrophysics Data System (ADS)

Anderson, James; Hargather, Michael

2016-11-01

Explosively-driven shock tube experiments were performed to investigate the turbulent mixing of explosive product gases and ambient air. A small detonator initiated Al / I2O5 thermite, which produced a shock wave and expanding product gases. Schlieren and imaging spectroscopy were applied simultaneously along a common optical path to identify correlations between turbulent structures and spatially-resolved absorbance. The schlieren imaging identifies flow features including shock waves and turbulent structures while the imaging spectroscopy identifies regions of iodine gas presence in the product gases. Pressure transducers located before and after the optical diagnostic section measure time-resolved pressure. Shock speed is measured from tracking the leading edge of the shockwave in the schlieren images and from the pressure transducers. The turbulent mixing characteristics were determined using digital image processing. Results show changes in shock speed, product gas propagation, and species concentrations for varied explosive charge mass. Funded by DTRA Grant HDTRA1-14-1-0070.

Turbulent Water Coupling in Shock Wave Lithotripsy

PubMed Central

Lautz, Jaclyn; Sankin, Georgy; Zhong, Pei

2013-01-01

Previous studies have demonstrated that stone comminution decreases with increased pulse repetition frequency as a result of bubble proliferation in the cavitation field of a shock wave lithotripter (Pishchalnikov et al., 2011). If cavitation nuclei remain in the propagation path of successive lithotripter pulses, especially in the acoustic coupling cushion of the shock wave source, they will consume part of the incident wave energy, leading to reduced tensile pressure in the focal region and thus lower stone comminution efficiency. We introduce a method to remove cavitation nuclei from the coupling cushion between successive shock exposures using a jet of degassed water. As a result, pre-focal bubble nuclei lifetime quantified by B-mode ultrasound imaging was reduced from 7 s to 0.3 s by a jet with an exit velocity of 62 cm/s. Stone fragmentation (percent mass < 2 mm) after 250 shocks delivered at 1 Hz was enhanced from 22 ± 6% to 33 ± 5% (p = 0.007) in water without interposing tissue mimicking materials. Stone fragmentation after 500 shocks delivered at 2 Hz was increased from 18 ± 6% to 28 ± 8% (p = 0.04) with an interposing tissue phantom of 8 cm thick. These results demonstrate the critical influence of cavitation bubbles in the coupling cushion on stone comminution and suggest a potential strategy to improve the efficacy of contemporary shock wave lithotripters. PMID:23322027

Shock interaction behind a pair of cylindrical obstacles

NASA Astrophysics Data System (ADS)

Liu, Heng; Mazumdar, Raoul; Eliasson, Veronica

2014-11-01

The body of work focuses on two-dimensional numerical simulations of shock interaction with a pair of cylindrical obstacles, varying the obstacle separation and incident shock strength. With the shock waves propagating parallel to the center-line between the two cylindrical obstacles, the shock strengths simulated vary from a Mach of 1.4 to a Mach of 2.4, against a wide range of obstacle separation distance to their diameters. These cases are simulated via a software package called Overture, which is used to solve the inviscid Euler equations of gas dynamics on overlapping grids with adaptive mesh refinement. The goal of these cases is to find a so-called ``safe'' region for obstacle spacing and varying shock Mach numbers, such that the pressure in the ``safe'' region is reduced downstream of the obstacles. The benefits apply to both building and armor design for the purpose of shock wave mitigation to keep humans and equipment safe. The results obtained from the simulations confirm that the length of the ``safe'' region and the degree of shock wave attenuation depend on the ratio of obstacle separation distance to obstacle diameter. The influence of various Mach number is also discussed.

Preliminary airborne measurements for the SR-71 sonic boom propagation experiment

NASA Technical Reports Server (NTRS)

Haering, Edward A., Jr.; Ehernberger, L. J.; Whitmore, Stephen A.

1995-01-01

SR-71 sonic boom signatures were measured to validate sonic boom propagation prediction codes. An SR-71 aircraft generated sonic booms from Mach 1.25 to Mach 1.6, at altitudes of 31,000 to 48,000 ft, and at various gross weights. An F-16XL aircraft measured the SR-71 near-field shock waves from close to the aircraft to more than 8,000 ft below, gathering 105 signatures. A YO-3A aircraft measured the SR-71 sonic booms from 21,000 to 38,000 feet below, recording 17 passes. The sonic booms at ground level and atmospheric data were recorded for each flight. Data analysis is underway. Preliminary results show that shock wave patterns and coalescence vary with SR-71 gross weight, Mach number, and altitude. For example, noncoalesced shock wave signatures were measured by the YO-3A at 21,000 ft below the SR-71 aircraft while at a low gross weight, Mach 1.25, and 31,000-ft altitude. This paper describes the design and execution of the flight research experiment. Instrumentation and flight maneuvers of the SR-71, F-16XL, and YO-3A aircraft and sample sonic boom signatures are included.

Switch-on Shock and Nonlinear Kink Alfvén Waves in Solar Polar Jets

NASA Astrophysics Data System (ADS)

DeVore, C. Richard; Karpen, Judith T.; Antiochos, Spiro K.; Uritsky, Vadim

2016-05-01

It is widely accepted that solar polar jets are produced by fast magnetic reconnection in the low corona, whether driven directly by flux emergence from below or indirectly by instability onset above the photosphere. In either scenario, twisted flux on closed magnetic field lines reconnects with untwisted flux on nearby open field lines. Part of the twist is inherited by the newly reconnected open flux, which rapidly relaxes due to magnetic tension forces that transmit the twist impulsively into the outer corona and heliosphere. We propose that this transfer of twist launches switch-on MHD shock waves, which propagate parallel to the ambient coronal magnetic field ahead of the shock and convect a perpendicular component of magnetic field behind the shock. In the frame moving with the shock front, the post-shock flow is precisely Alfvénic in all three directions, whereas the pre-shock flow is super-Alfvénic along the ambient magnetic field, yielding a density enhancement at the shock front. Nonlinear kink Alfvén waves are exact solutions of the time-dependent MHD equations in the post-shock region when the ambient corona is uniform and the magnetic field is straight. We have performed and analyzed 3D Cartesian and spherical simulations of polar jets driven by instability onset in the corona. The results of both simulations are consistent with the generation of MHD switch-on shocks trailed predominantly by incompressible kink Alfvén waves. It is noteworthy that the kink waves are irrotational, in sharp contrast to the vorticity-bearing torsional waves reported from previous numerical studies. We will discuss the implications of the results for understanding solar polar jets and predicting their heliospheric signatures. Our research was supported by NASA’s LWS TR&T and H-SR programs.

The soft X-ray background as a supernova blast wave viewed from inside - Solar abundance models

NASA Technical Reports Server (NTRS)

Edgar, R. J.

1986-01-01

A model of the soft X-ray background is presented in which the sun is assumed to be inside an active supernova blast wave. The blast wave evolves in a preexisting cavity. The broad band surface brightnesses is explained by such a blast wave with an explosion energy of E sub approximately 5 x 10 to the 50th power ergs and radius 80 to 100 pc, using solar abundances. An approach to treating the problem of large anisotropies in the ambient medium is also explored, accommodating the observed anticorrelation between the soft X-ray surface brightness and the 21 cm column density. It is found that only for post shock temperatures below 10 to the 6 power K a shock propagating into a density enhancement will be dimmer than a similar shock in a lower density region.

An artificial nonlinear diffusivity method for supersonic reacting flows with shocks

NASA Astrophysics Data System (ADS)

Fiorina, B.; Lele, S. K.

2007-03-01

A computational approach for modeling interactions between shocks waves, contact discontinuities and reactions zones with a high-order compact scheme is investigated. To prevent the formation of spurious oscillations around shocks, artificial nonlinear viscosity [A.W. Cook, W.H. Cabot, A high-wavenumber viscosity for high resolution numerical method, J. Comput. Phys. 195 (2004) 594-601] based on high-order derivative of the strain rate tensor is used. To capture temperature and species discontinuities a nonlinear diffusivity based on the entropy gradient is added. It is shown that the damping of 'wiggles' is controlled by the model constants and is largely independent of the mesh size and the shock strength. The same holds for the numerical shock thickness and allows a determination of the L2 error. In the shock tube problem, with fluids of different initial entropy separated by the diaphragm, an artificial diffusivity is required to accurately capture the contact surface. Finally, the method is applied to a shock wave propagating into a medium with non-uniform density/entropy and to a CJ detonation wave. Multi-dimensional formulation of the model is presented and is illustrated by a 2D oblique wave reflection from an inviscid wall, by a 2D supersonic blunt body flow and by a Mach reflection problem.

Partial Reflection and Trapping of a Fast-mode Wave in Solar Coronal Arcade Loops

NASA Astrophysics Data System (ADS)

Kumar, Pankaj; Innes, D. E.

2015-04-01

We report on the first direct observation of a fast-mode wave propagating along and perpendicular to cool (171 Å) arcade loops observed by the Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA). The wave was associated with an impulsive/compact flare near the edge of a sunspot. The EUV wavefront expanded radially outward from the flare center and decelerated in the corona from 1060 to 760 km s-1 within ˜3-4 minutes. Part of the EUV wave propagated along a large-scale arcade of cool loops and was partially reflected back to the flare site. The phase speed of the wave was about 1450 km s-1, which is interpreted as a fast-mode wave. A second overlying loop arcade, orientated perpendicular to the cool arcade, is heated and becomes visible in the AIA hot channels. These hot loops sway in time with the EUV wave, as it propagated to and fro along the lower loop arcade. We suggest that an impulsive energy release at one of the footpoints of the arcade loops causes the onset of an EUV shock wave that propagates along and perpendicular to the magnetic field.

Formation of multiple energy dispersion of H+, He+, and O+ ions in the inner magnetosphere in response to interplanetary shock

NASA Astrophysics Data System (ADS)

Tsuji, H.; Ebihara, Y.; Tanaka, T.

2017-04-01

An interplanetary (IP) shock has a large impact on magnetospheric ions. Satellite observations have shown that soon after arrival of the IP shock, overall intensity of the ions rapidly increases and multiple energy dispersion appears in an energy-time spectrogram of the ions. In order to understand the response of the magnetospheric ions to IP shock, we have performed test particle simulation under the electric and magnetic fields provided by the global magnetohydrodynamic simulation. We reconstructed the differential flux of H+, He+, and O+ ions at (7, 0, 0) Re in GSM coordinates by means of the semi-Lagrangian (phase space mapping) method. Simulation results show that the ions respond to the IP shock in two different ways. First, overall intensity of the flux gradually increases at all pitch angles. As the compressional wave propagates tailward, the magnetic field increases, which accelerates the ions due to the gyrobetatron. Second, multiple energy-time dispersion appears in the reconstructed spectrograms of the ion flux. The energy-time dispersion is caused by the ion moving toward mirror point together with tailward propagating compressional wave at off-equator. The ions are primarily accelerated by the drift betatron under the strong electric field looking dawnward. The dispersion is absent in the spectrogram of equatorially mirroring ions. The dispersion appears at higher energy for heavier ions. These features are consistent with the satellite observations. Because the acceleration depends on bounce phase, the bounce-averaged approximation is probably invalid for the ions during the interval of geomagnetic sudden commencement.