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 1,603 records · Page 89

SHOCK BREAKOUT IN DENSE MASS LOSS: LUMINOUS SUPERNOVAE

We examine the case where a circumstellar medium around a supernova is sufficiently opaque that a radiation-dominated shock propagates in the circumstellar region. The initial propagation of the shock front into the circumstellar region can be approximated by a self-similar solution that determines the radiative energy in a shocked shell; the eventual escape of this energy gives the maximum luminosity of the supernova. If the circumstellar density is described by {rho} = Dr {sup -2} out to a radius R{sub w} , where D is a constant, the properties of the shock breakout radiation depend on R{sub w} and R{sub d} {identical_to} {kappa}Dv{sub sh}/c, where {kappa} is the opacity and v{sub sh} is the shock velocity. If R{sub w} >R{sub d} , the rise to maximum light begins at {approx}R{sub d} /v{sub sh}; the duration of the rise is also {approx}R{sub d} /v{sub sh}; the outer parts of the opaque medium are extended and at low velocity at the time of peak luminosity; and a dense shell forms whose continued interaction with the dense mass loss gives a characteristic flatter portion of the declining light curve. If R{sub w} < R{sub d} , the rise to maximum light begins at R{sub w} /v{sub sh}; the duration of the rise is R {sup 2}{sub w}/v{sub sh} R{sub d} ; the outer parts of the opaque medium are not extended and are accelerated to high velocity by radiation pressure at the time of maximum luminosity; and a dense shell forms but does not affect the light curve near maximum. We argue that SN 2006gy is an example of the first kind of event, while SN 2010gx and related supernovae are examples of the second.

79 ASTRONOMY AND ASTROPHYSICS↗

COSMIC-RAY-INDUCED FILAMENTATION INSTABILITY IN COLLISIONLESS SHOCKS

We used unprecedentedly large two-dimensional and three-dimensional hybrid (kinetic ions-fluid electrons) simulations of non-relativistic collisionless strong shocks in order to investigate the effects of self-consistently accelerated ions on the overall shock dynamics. The current driven by suprathermal particles streaming ahead of the shock excites modes transverse to the background magnetic field. The Lorentz force induced by these self-amplified fields tends to excavate tubular, underdense, magnetic-field-depleted cavities that are advected with the fluid and perturb the shock surface, triggering downstream turbulent motions. These motions further amplify the magnetic field, up to factors of 50-100 in knot-like structures. Once downstream, the cavities tend to be filled by hot plasma plumes that compress and stretch the magnetic fields in elongated filaments; this effect is particularly evident if the shock propagates parallel to the background field. Highly magnetized knots and filaments may provide explanations for the rapid X-ray variability observed in RX J1713.7-3946 and for the regular pattern of X-ray bright stripes detected in Tycho's supernova remnant.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗

Sound velocities in shock-compressed soda lime glass: Melting and liquid-state response

Longitudinal sound velocities (and elastic moduli) were determined in soda lime glass (SLG) at peak shock stresses ranging between 40 and 90 GPa. Laser interferometry was used to obtain particle velocity histories and sound velocities by impacting SLG samples on lithium fluoride (LiF) optical windows. In all experiments, the SLG response consists of a sharp jump to a constant state followed by a release wave. The measured longitudinal sound velocities and moduli showed a marked decrease between 52 and 58 GPa, providing experimental evidence for the transformation from an amorphous solid to a liquid in shock-compressed SLG. Furthermore, the stress threshold (~55 GPa) for melting in SLG is considerably lower than the threshold reported in shock-compressed fused silica (~72 GPa), showing the effect of network-modifying cations on the onset of melting. The relative values of sound velocities, shock velocities, and the Hugoniot slopes—between 58 and 90 GPa—are fully consistent with the thermodynamic response of a shock-compressed liquid. Using the experimental results, the Grüneisen parameter (Γ) values were determined for liquid SLG to 90 GPa and then used to provide the Mie-Grüneisen equation of state for liquid SLG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach

Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisionless Shocks Driven by Supersonic Plasma Flows with Self-Generated Magnetic Fields

Collisionless shocks are ubiquitous in the Universe as a consequence of supersonic plasma flows sweeping through interstellar and intergalactic media. These shocks are the cause of many observed astrophysical phenomena, but details of shock structure and behavior remain controversial because of the lack of ways to study them experimentally. Laboratory experiments reported here, with astrophysically relevant plasma parameters, demonstrate for the first time the formation of a quasiperpendicular magnetized collisionless shock. In the upstream it is fringed by a filamented turbulent region, a rudiment for a secondary Weibel-driven shock. Furthermore, this turbulent structure is found responsible for electron acceleration to energies exceeding the average energy by two orders of magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-Generated Magnetic and Electric Fields at a Mach-6 Shock Front in a Low Density Helium Gas by Dual-Angle Proton Radiography

Shocks are abundant both in astrophysical and laboratory systems. While the electric fields generated at shock fronts have recently attracted great attention, the associated self-generated magnetic field is rarely studied, despite its ability to significantly affect the shock profile in the nonideal geometry where density and temperature gradients are not parallel. We report here the observation of a magnetic field at the front of a Mach ~6 shock propagating in a low-density helium gas system. Proton radiography from different projection angles not only confirms the magnetic field’s existence, but also provides a quantitative measurement of the field strength in the range ~5 to 7 T. X-ray spectrometry allowed inference of the density and temperature at the shock front, constraining the plasma conditions under which the magnetic and electric fields are generated. Furthermore, simulations with the particle-in-cell code lsp attribute the self-generation of the magnetic field to the Biermann battery effect (∇n e × ∇T e ).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock-wave studies of anomalous compressibility of glassy carbon

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

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optoheterodyne Doppler measurements of the ballistic expansion of the products of the shock wave-induced surface destruction: Experiment and theory

The results of optoheterodyne Doppler measurements of the ballistic expansion of the products of surface destruction under shock-wave loading are presented. The possibility of determining the physical characteristics of a rapidly flying dust cloud, including the microparticle velocities, the microparticle sizes, and the areal density of the dust cloud, is shown. A compact stand for performing experiments on shock-wave loading of metallic samples is described. Shock-wave loading is performed by a 100-µm-thick tantalum flyer plate accelerated to a velocity of 2.8 km/s. As the samples, lead plates having various thicknesses and the same surface roughness are used. At a shock-wave pressure of 31.5 GPa, the destruction products are solid microparticles about 50 µm in size. At a pressure of 42 and 88 GPa, a liquid-drop dust cloud with a particle size of 10–15 µm is formed. To interpret the spectral data on the optoheterodyne Doppler measurements of the expansion of the surface destruction products (spalled fragments, dust microparticles), a transport equation for the function of mutual coherence of a multiply scattered field is used. The Doppler spectra of a backscattered signal are calculated with the model developed for the dust cloud that appears when a shock wave reaches the sample surface at the parameters that are typical of an experimental situation. Qualitative changes are found in the spectra, depending on the optical thickness of the dust cloud. The obtained theoretical results are in agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental evaluation of low-pass shock isolation performance of elastomers using frequency-based Kolsky bar analyses

Elastomeric materials are used as shock isolation materials in a variety of environments to dampen vibrations and/or absorb energy from external impact to minimize energy transfer between two objects or bodies. Some applications require the shock isolation materials to behave as a low-pass mechanical filter to mitigate the shock/impact at high frequencies but transmit the energy at low frequencies with minimal attenuation. To fulfill this requirement, a shock isolation material needs to be carefully evaluated and selected with proper experimental design, procedures, and analyses. In this study, a Kolsky bar was modified with precompression (up to 15.5 kN) and confinement capabilities to evaluate low-pass shock isolation performance in terms of acceleration attenuation through a variety of elastomers. Also investigated were the effects of preload and specimen geometry on the low-pass shock isolation response.

36 MATERIALS SCIENCE↗

Barrier experiment: Shock initiation under complex loading

The barrier experiments are a variant of the gap test; a detonation wave in a donor HE impacts a barrier and drives a shock wave into an acceptor HE. The question we ask is: What is the trade-off between the barrier material and threshold barrier thickness to prevent the acceptor from detonating. This can be viewed from the perspective of shock initiation of the acceptor subject to a complex pressure drive condition. Here we consider key factors which affect whether or not the acceptor undergoes a shock-to-detonation transition. These include the following: shock impedance matches for the donor detonation wave into the barrier and then the barrier shock into the acceptor, the pressure gradient behind the donor detonation wave, and the curvature of detonation front in the donor. Numerical simulations are used to illustrate how these factors affect the reaction in the acceptor.

36 MATERIALS SCIENCE↗

Magnetar-driven shock breakout and double-peaked supernova light curves

The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (∼10{sup 14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T{sub eff} ≈ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar-driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SN-LSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double-peaked light curve may therefore only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe how such a delay in heating may naturally result from inefficient dissipation and thermalization of the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.

79 ASTRONOMY AND ASTROPHYSICS↗

HIGHLY EXCITED H{sub 2} IN HERBIG–HARO 7: FORMATION PUMPING IN SHOCKED MOLECULAR GAS?

We have obtained K -band spectra at R ∼ 5000 and an angular resolution of 0.″3 of a section of the Herbig–Haro 7 (HH7) bow shock, using the Near-Infrared Integral Field Spectrograph at Gemini North. Present in the portion of the data cube corresponding to the brightest part of the bow shock are emission lines of H{sub 2} with upper state energies ranging from ∼6000 K to the dissociation energy of H{sub 2}, ∼50,000 K. Because of low signal-to-noise ratios, the highest excitation lines cannot be easily seen elsewhere in the observed region. However, excitation temperatures, measured throughout much of the observed region using lines from levels as high as 25,000 K, are a strong function of upper level energy, indicating that the very highest levels are populated throughout. The level populations in the brightest region are well fit by a two-temperature model, with 98.5% of the emitting gas at T = 1800 K and 1.5% at T = 5200 K. The bulk of the H{sub 2} line emission in HH7, from the 1800 K gas, has previously been well-modeled by a continuous shock, but the 5200 K cozmponent is inconsistent with standalone standard continuous shock models. We discuss various possible origins for the hot component and suggest that this component is H{sub 2} newly reformed on dust grains and then ejected from them, presumably following dissociation of some of the H{sub 2} by the shock.

79 ASTRONOMY AND ASTROPHYSICS↗

Global MHD simulations of accretion disks in cataclysmic variables. I. The importance of spiral shocks

We present results from the first global 3D MHD simulations of accretion disks in cataclysmic variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared with that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component of the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states, both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. (1) The linear wave dispersion relation fits the pitch angles of spiral density waves very well. (2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. (3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate α{sub eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. (4) Using the mass accretion rate measured in our simulations to directly measure α defined in standard thin-disk theory, we find 0.02≲α{sub eff}≲0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae. In this regime, the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime and, therefore, future models with more realistic thermodynamics and non-ideal MHD are warranted.

79 ASTRONOMY AND ASTROPHYSICS↗

COLLISIONLESS ELECTRON–ION SHOCKS IN RELATIVISTIC UNMAGNETIZED JET–AMBIENT INTERACTIONS: NON-THERMAL ELECTRON INJECTION BY DOUBLE LAYER

The course of non-thermal electron ejection in relativistic unmagnetized electron–ion shocks is investigated by performing self-consistent particle-in-cell simulations. The shocks are excited through the injection of a relativistic jet into ambient plasma, leading to two distinct shocks (referred to as the trailing shock and leading shock) and a contact discontinuity. The Weibel-like instabilities heat the electrons up to approximately half of the ion kinetic energy. The double layers formed in the trailing and leading edges then accelerate the electrons up to the ion kinetic energy. The electron distribution function in the leading edge shows a clear, non-thermal power-law tail which contains ∼1% of electrons and ∼8% of the electron energy. Its power-law index is −2.6. The acceleration efficiency is ∼23% by number and ∼50% by energy, and the power-law index is −1.8 for the electron distribution function in the trailing edge. The effect of the dimensionality is examined by comparing the results of three-dimensional simulations with those of two-dimensional simulations. The comparison demonstrates that electron acceleration is more efficient in two dimensions.

79 ASTRONOMY AND ASTROPHYSICS↗

FLARE VERSUS SHOCK ACCELERATION OF HIGH-ENERGY PROTONS IN SOLAR ENERGETIC PARTICLE EVENTS

Recent studies have presented evidence for a significant to dominant role for a flare-resident acceleration process for high-energy protons in large (“gradual”) solar energetic particle (SEP) events, contrary to the more generally held view that such protons are primarily accelerated at shock waves driven by coronal mass ejections (CMEs). The new support for this flare-centric view is provided by correlations between the sizes of X-ray and/or microwave bursts and associated SEP events. For one such study that considered >100 MeV proton events, we present evidence based on CME speeds and widths, shock associations, and electron-to-proton ratios that indicates that events omitted from that investigation’s analysis should have been included. Inclusion of these outlying events reverses the study’s qualitative result and supports shock acceleration of >100 MeV protons. Examination of the ratios of 0.5 MeV electron intensities to >100 MeV proton intensities for the Grechnev et al. event sample provides additional support for shock acceleration of high-energy protons. Simply scaling up a classic “impulsive” SEP event to produce a large >100 MeV proton event implies the existence of prompt 0.5 MeV electron events that are approximately two orders of magnitude larger than are observed. While classic “impulsive” SEP events attributed to flares have high electron-to-proton ratios (≳5 × 10{sup 5}) due to a near absence of >100 MeV protons, large poorly connected (≥W120) gradual SEP events, attributed to widespread shock acceleration, have electron-to-proton ratios of ∼2 × 10{sup 3}, similar to those of comparably sized well-connected (W20–W90) SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗

On Weibull's Spectrum of Nonrelativistic Energetic Particles at IP Shocks: Observations and Theoretical Interpretation

Some interplanetary shocks are associated with short-term and sharp particle flux enhancements near the shock front. Such intensity enhancements, known as shock-spike events (SSEs), represent a class of relatively energetic phenomena as they may extend to energies of some tens of MeV or even beyond. Here we present an SSE case study in order to shed light on the nature of the particle acceleration involved in this kind of event. Our observations refer to an SSE registered on 2011 October 3 at 22:23 UT, by STEREO B instrumentation when, at a heliocentric distance of 1.08 au, the spacecraft was swept by a perpendicular shock moving away from the Sun. The main finding from the data analysis is that a Weibull distribution represents a good fitting function to the measured particle spectrum over the energy range from 0.1 to 30 MeV. To interpret such an observational result, we provide a theoretical derivation of the Weibull spectrum in the framework of the acceleration by “killed” stochastic processes exhibiting power-law growth in time of the velocity expectation, such as the classical Fermi process. We find an overall coherence between the experimental values of the Weibull spectrum parameters and their physical meaning within the above scenario. Hence, our approach based on the Weibull distribution proves to be useful for understanding SSEs. With regard to the present event, we also provide an alternative explanation of the Weibull spectrum in terms of shock-surfing acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical Simulations of Particle Acceleration at Interplanetary Quasi-perpendicular Shocks

Using test particle simulations we study particle acceleration at highly perpendicular (θ{sub Bn}⩾75{sup ∘}) shocks under conditions of modeling magnetic turbulence. We adopt a backward-in-time method to solve the Newton–Lorentz equation using the observed shock parameters for quasi-perpendicular interplanetary shocks, and compare the simulation results with ACE/EPAM observations to obtain the injection energy and timescale of particle acceleration. With our modeling and observations, we find that a large upstream speed is responsible for efficient particle acceleration. Our results also show that the quasi-perpendicular shocks are capable of accelerating thermal particles to high energies of the order of MeV for both kappa and Maxwellian upstream distributions, which may originate from the fact that in our model, the local background magnetic field has a component parallel to the shock normal.

79 ASTRONOMY AND ASTROPHYSICS↗