Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Collisionless Weibel shocks: Full formation mechanism and timing

Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

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

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 been 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Ablative stabilization of Rayleigh-Taylor instabilities resulting from a laser-driven radiative shock

The Rayleigh-Taylor (RT) instability is a common occurrence in nature, notably in astrophysical systems like supernovae, where it serves to mix the dense layers of the interior of an exploding star with the low-density stellar wind surrounding it, and in inertial confinement fusion experiments, where it mixes cooler materials with the central hot spot in an imploding capsule and stifles the desired nuclear reactions. In both of these examples, the radiative flux generated by strong shocks in the system may play a role in partially stabilizing RT instabilities. We present experiments performed on the National Ignition Facility, designed to isolate and study the role of radiation and heat conduction from a shock front in the stabilization of hydrodynamic instabilities. By varying the laser power delivered to a shock-tube target with an embedded, unstable interface, the radiative fluxes generated at the shock front could be controlled. We observe decreased RT growth when the shock significantly heats the medium around it, in contrast to a system where the shock did not produce significant heating. Both systems are modeled with a modified set of buoyancy-drag equations accounting for ablative stabilization, and the experimental results are consistent with ablative stabilization when the shock is radiative. This result has important implications for our understanding of astrophysical radiative shocks and supernova radiative hydrodynamics [Kuranz et al., Nature Communications 9(1), 1564 (2018)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Influence of baroclinic vorticity production on unsteady drag coefficient in shock–particle interaction

The influence of baroclinic vorticity production on the unsteady drag coefficient in shock–particle interaction is numerically studied in this work. Numerical simulations are performed for shock–particle interaction utilizing a high–resolution axisymmetric solver for the Euler equations that allows for multi-material interface and shock propagation in both the particle and surrounding medium. We consider an aluminum particle in nitromethane and allow for particle deformation. We compute the vorticity production and unsteady drag coefficient as a function of time to explain the complex physical mechanisms that occur during shock–particle interaction. We observe baroclinic vorticity production as the shock propagates over the particle and find that the vorticity is primarily generated at the surface of the particle. After the passage of the shock over the particle, the generated vortex traverses downstream, thus creating a sharpened particle edge and low pressure on the downstream side of the particle, followed by the trapping of the vortex at the particle edge. These mechanisms lead to the generation of a quasi-steady drag force even after the passage of the shock, thus suggesting the importance of baroclinic vorticity production to the unsteady drag coefficient. Finally, we compute the unsteady drag coefficient for various shock Mach numbers and particle ellipticities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Simulations of ion acceleration at non-relativistic shocks. I. Acceleration efficiency

We use two-dimensional and three-dimensional hybrid (kinetic ions-fluid electrons) simulations to investigate particle acceleration and magnetic field amplification at non-relativistic astrophysical shocks. We show that diffusive shock acceleration operates for quasi-parallel configurations (i.e., when the background magnetic field is almost aligned with the shock normal) and, for large sonic and Alfvénic Mach numbers, produces universal power-law spectra ∝p {sup –4}, where p is the particle momentum. The maximum energy of accelerated ions increases with time, and it is only limited by finite box size and run time. Acceleration is mainly efficient for parallel and quasi-parallel strong shocks, where 10%-20% of the bulk kinetic energy can be converted to energetic particles and becomes ineffective for quasi-perpendicular shocks. Also, the generation of magnetic turbulence correlates with efficient ion acceleration and vanishes for quasi-perpendicular configurations. At very oblique shocks, ions can be accelerated via shock drift acceleration, but they only gain a factor of a few in momentum and their maximum energy does not increase with time. These findings are consistent with the degree of polarization and the morphology of the radio and X-ray synchrotron emission observed, for instance, in the remnant of SN 1006. We also discuss the transition from thermal to non-thermal particles in the ion spectrum (supra-thermal region) and we identify two dynamical signatures peculiar of efficient particle acceleration, namely, the formation of an upstream precursor and the alteration of standard shock jump conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Shocks and star formation in Stephan's Quintet. I. Gemini spectroscopy of Hα-bright knots

We present a Gemini-GMOS spectroscopic study of Hubble Space Telescope (HST)-selected Hα-emitting regions in Stephan's Quintet (HCG 92), a nearby compact galaxy group, with the aim of disentangling the processes of shock-induced heating and star formation in its intra-group medium. The ≈40 sources are distributed across the system, but most densely concentrated in the ∼kiloparsec-long shock region. Their spectra neatly divide them into narrow- and broad-line emitters, and we decompose the latter into three or more emission peaks corresponding to spatial elements discernible in HST imaging. The emission-line ratios of the two populations of Hα-emitters confirm their nature as H II regions (90% of the sample) or molecular gas heated by a shock front propagating at ≲300 km s{sup –1}. Their redshift distribution reveals interesting three-dimensional structure with respect to gas-phase baryons, with no H II regions associated with shocked gas, no shocked regions in the intruder galaxy NGC 7318B, and a sharp boundary between shocks and star formation. We conclude that star formation is inhibited substantially, if not entirely, in the shock region. Attributing those H II regions projected against the shock to the intruder, we find a lopsided distribution of star formation in this galaxy, reminiscent of pileup regions in models of interacting galaxies. The Hα luminosities imply mass outputs, star formation rates, and efficiencies similar to nearby star-forming regions. Two large knots are an exception to this, being comparable in stellar output to the prolific 30 Doradus region. We also examine Stephan's Quintet in the context of compact galaxy group evolution, as a paradigm for intermittent star formation histories in the presence of a rich, X-ray-emitting intra-group medium. All spectra are provided as supplemental materials.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Shock waves and cosmic ray acceleration in the outskirts of galaxy clusters

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within r {sub 200} in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside r {sub 200} has the power-law slope of 4.25-4.5, indicating that the average Mach number of the shocks of main CR production is in the range of {sub CR} ≈ 3-4. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with M{sub s} ≳ 3 induced in the cluster outskirts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Spectrum and light curve of a supernova shock breakout through a thick Wolf-Rayet wind

Wolf-Rayet stars are known to eject winds. Thus, when a Wolf-Rayet star explodes as a supernova, a fast (≳ 40, 000 km s{sup –1}) shock is expected to be driven through a wind. We study the signal expected from a fast supernova shock propagating through an optically thick wind and find that the electrons behind the shock driven into the wind are efficiently cooled by inverse Compton over soft photons that were deposited by the radiation-mediated shock that crossed the star. Therefore, the bolometric luminosity is comparable to the kinetic energy flux through the shock, and the spectrum is found to be a power law, whose slope and frequency range depend on the number flux of soft photons available for cooling. Wolf-Rayet supernovae that explode through a thick wind have a high flux of soft photons, producing a flat spectrum, νF {sub ν} = Const, in the X-ray range of 0.1 ≲ T ≲ 50 keV. As the shock expands into an optically thin wind, the soft photons are no longer able to cool the shock that plows through the wind, and the bulk of the emission takes the form of a standard core-collapse supernova (without a wind). However, a small fraction of the soft photons is upscattered by the shocked wind and produces a transient unique X-ray signature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Magnetic field amplification in nonlinear diffusive shock acceleration including resonant and non-resonant cosmic-ray driven instabilities

We present a nonlinear Monte Carlo model of efficient diffusive shock acceleration where the magnetic turbulence responsible for particle diffusion is calculated self-consistently from the resonant cosmic-ray (CR) streaming instability, together with non-resonant short- and long-wavelength CR-current-driven instabilities. We include the backpressure from CRs interacting with the strongly amplified magnetic turbulence which decelerates and heats the super-Alfvénic flow in the extended shock precursor. Uniquely, in our plane-parallel, steady-state, multi-scale model, the full range of particles, from thermal (∼eV) injected at the viscous subshock to the escape of the highest energy CRs (∼PeV) from the shock precursor, are calculated consistently with the shock structure, precursor heating, magnetic field amplification, and scattering center drift relative to the background plasma. In addition, we show how the cascade of turbulence to shorter wavelengths influences the total shock compression, the downstream proton temperature, the magnetic fluctuation spectra, and accelerated particle spectra. A parameter survey is included where we vary shock parameters, the mode of magnetic turbulence generation, and turbulence cascading. From our survey results, we obtain scaling relations for the maximum particle momentum and amplified magnetic field as functions of shock speed, ambient density, and shock size.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The effect of stochastic re-acceleration on the energy spectrum of shock-accelerated protons

The energy spectra of particles in gradual solar energetic particle (SEP) events do not always have a power-law form attributed to the diffusive shock acceleration mechanism. In particular, the observed spectra in major SEP events can take the form of a broken (double) power law. In this paper, we study the effect of a process that can modify the power-law spectral form produced by the diffusive shock acceleration: the stochastic re-acceleration of energetic protons by enhanced Alfvénic turbulence in the downstream region of a shock wave. There are arguments suggesting that this process can be important when the shock propagates in the corona. We consider a coronal magnetic loop traversed by a shock and perform Monte Carlo simulations of interactions of shock-accelerated protons with Alfvén waves in the loop. The wave-particle interactions are treated self-consistently, so the finiteness of the available turbulent energy is taken into account. The initial energy spectrum of particles is taken to be a power law. The simulations reveal that the stochastic re-acceleration leads either to the formation of a spectrum that is described in a wide energy range by a power law (although the resulting power-law index is different from the initial one) or to a broken power-law spectrum. The resulting spectral form is determined by the ratio of the energy density of shock-accelerated protons to the wave energy density in the shock's downstream region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Non-thermal electron acceleration in low Mach number collisionless shocks. I. Particle energy spectra and acceleration mechanism

Electron acceleration to non-thermal energies in low Mach number (M{sub s} ≲ 5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M{sub s} = 3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15% of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p ≅ 2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration (SDA). The accelerated electrons are then reflected back upstream where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock for further energization via SDA. In summary, the self-generated waves allow for repeated cycles of SDA, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A VERY DEEP CHANDRA OBSERVATION OF THE GALAXY GROUP NGC 5813: AGN SHOCKS, FEEDBACK, AND OUTBURST HISTORY

We present results from a very deep (650 ks) Chandra X-ray observation of the galaxy group NGC 5813, the deepest Chandra observation of a galaxy group to date. This system uniquely shows three pairs of collinear cavities, with each pair associated with an unambiguous active galactic nucleus (AGN) outburst shock front. The implied mean kinetic power is roughly the same for each outburst, demonstrating that the average AGN kinetic luminosity can remain stable over long timescales (∼50 Myr). The two older outbursts have larger, roughly equal total energies as compared with the youngest outburst, implying that the youngest outburst is ongoing. We find that the gas radiative cooling rate and mean shock heating rate are well balanced at each shock front, suggesting that shock heating alone is sufficient to offset cooling and establish AGN/intracluster medium (ICM) feedback within at least the central 30 kpc. This heating takes place roughly isotropically and most strongly at small radii, as is required for feedback to operate. We suggest that shock heating may play a significant role in AGN feedback at smaller radii in other systems, where weak shocks are more difficult to detect. We find non-zero shock front widths that are too large to be explained by particle diffusion. Instead, all measured widths are consistent with shock broadening due to propagation through a turbulent ICM with a mean turbulent speed of ∼70 km s{sup −1}. Finally, we place lower limits on the temperature of any volume-filling thermal gas within the cavities that would balance the internal cavity pressure with the external ICM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RECOLLIMATION SHOCKS IN MAGNETIZED RELATIVISTIC JETS

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal, and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, significantly modifying the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behavior depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PRECURSORS TO INTERSTELLAR SHOCKS OF SOLAR ORIGIN

On or about 2012 August 25, the Voyager 1 spacecraft crossed the heliopause into the nearby interstellar plasma. In the nearly three years that the spacecraft has been in interstellar space, three notable particle and field disturbances have been observed, each apparently associated with a shock wave propagating outward from the Sun. Here, we present a detailed analysis of the third and most impressive of these disturbances, with brief comparisons to the two previous events, both of which have been previously reported. The shock responsible for the third event was first detected on 2014 February 17 by the onset of narrowband radio emissions from the approaching shock, followed on 2014 May 13 by the abrupt appearance of intense electron plasma oscillations generated by electrons streaming outward ahead of the shock. Finally, the shock arrived on 2014 August 25, as indicated by a jump in the magnetic field strength and the plasma density. Various disturbances in the intensity and anisotropy of galactic cosmic rays were also observed ahead of the shock, some of which are believed to be caused by the reflection and acceleration of cosmic rays by the magnetic field jump at the shock, and/or by interactions with upstream plasma waves. Comparisons to the two previous weaker events show somewhat similar precursor effects, although differing in certain details. Many of these effects are very similar to those observed in the region called the “foreshock” that occurs upstream of planetary bow shocks, only on a vastly larger spatial scale.

79 ASTRONOMY AND ASTROPHYSICS↗