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,567 records · Page 87

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing shocks interacting with radiation waves with the Radishock experiment

Both radiation flows and shocks have been extensively studied in the laboratory in the past few decades due to their critical roles in many astrophysical and high-energy density physics processes. In the Radishock experiment, a halfraum-powered radiation wave is driven into a low-density foam and interacts with an ablatively driven, counter-propagating shock. The interacting waves produce a spike in energy density with a temperature greater than the local temperature of the individual waves. As in the successful predecessor experiment, COAX, the primary diagnostic uses absorption spectroscopy at many locations down the cylindrical target, enabling a spatial temperature inference of the radiation wave and its interactions with the shock. Combined with a radiography diagnostic that is capable of imaging the shock and interaction features, we are able to study and inform model predictions of the interaction spike phenomenon. We describe the underlying physics behind the shock interactions with the radiation front and the implications of this experimental study for a broad range of astrophysical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗

Numerical analysis of experiments on the generation of shock waves in aluminium under indirect (X-ray) action on the Iskra-5 facility

The dynamics of laser and X-ray radiation fields in experiments with cylindrical converter boxes (illuminators), which had earlier been carried out on the Iskra-5 laser facility (the second harmonic of iodine laser radiation, {lambda} = 0.66 {mu}m) was investigated in a sector approximation using the SND-LIRA numerical technique. In these experiments, the X-ray radiation temperature in the box was determined by measuring the velocity of the shock wave generated in the sample under investigation, which was located at the end of the cylindrical illuminator. Through simulations were made using the SND-LIRA code, which took into account the absorption of laser driver radiation at the box walls, the production of quasithermal radiation, as well as the formation and propagation of the shock wave in the sample under investigation. An analysis of the experiments permits determining the electron thermal flux limiter f: for f = 0.03 it is possible to match the experimental scaling data for X-ray in-box radiation temperature to the data of our simulations. The shock velocities obtained from the simulations are also consistent with experimental data. In particular, in the experiment with six laser beams (and a laser energy E{sub L} = 1380 J introduced into the box) the velocity of the shock front (determined from the position of a laser mark) after passage through a 50-{mu}m thick base aluminium layer was equal to 35{+-}1.6 km s{sup -1}, and in simulations to 36 km s{sup -1}. In the experiment with four laser beams (for E{sub L} = 850 J) the shock velocity (measured from the difference of transit times through the base aluminium layer and an additional thin aluminium platelet) was equal to 30{+-}3.6 km s{sup -1}, and in simulations to 30 km s{sup -1}. (interaction of laser radiation with matter)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock synthesis of quasicrystals with implications for their origin in asteroid collisions

Here, we designed a plate impact shock recovery experiment to simulate the starting materials and shock conditions associated with the only known natural quasicrystals, in the Khatyrka meteorite. At the boundaries among CuAl5, (Mg 0.75 Fe 0.25 2+ ) 2 SiO 4 olivine, and the stainless steel chamber walls, the recovered specimen contains numerous micron-scale grains of a quasicrystalline phase displaying face-centered icosahedral symmetry and low phason strain. The compositional range of the icosahedral phase is Al 68-73 Fe 11-16 Cu 10-12 Cr 1-4 Ni 1-2 and extends toward higher Al/(Cu+Fe) and Fe/Cu ratios than those reported for natural icosahedrite or for any previously known synthetic quasicrystal in the Al-Cu-Fe system. The shock-induced synthesis demonstrated in this experiment reinforces the evidence that natural quasicrystals formed during a shock event but leaves open the question of whether this synthesis pathway is attributable to the expanded thermodynamic stability range of the quasicrystalline phase at high pressure, to a favorable kinetic pathway that exists under shock conditions, or to both thermodynamic and kinetic factors.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

First detection of multi-shocks in RR Lyrae stars from Antarctica : A possible explanation of the Blazhko effect

We present the first detection of multi-shocks propagating through the atmosphere of the Blazhko star S Arae using uninterrupted, accurate optical photometric data collected during one polar night, 150 days from Antarctica at Dome C, with the Photometer AntarctIca eXtinction (PAIX). We acquired 89,736 CCD frames during 323 pulsation cycles and 3 Blazhko cycles. We detected two new light curve properties in the PAIX light curve, jump and rump, which we associated with two new post-maximum shock waves Sh{sub PM1} and Sh{sub PM2}. jump, lump, rump, bump, and hump are induced by five shock waves, with different amplitudes and origins, Sh{sub PM1}, Sh{sub PM}, Sh{sub PM2}, Sh{sub PM3}, and the main shock Sh{sub H+He}. Correlations between the length of rise time and light amplitude and Sh{sub PM3} are monotonous during three Blazhko cycles, but the pulsation curve is double peaked. We discuss the physical mechanisms driving the modulation of these quantities. Finally, we hypothesize that the origin of the Blazhko effect is a dynamical interaction between a multi-shock structure and an outflowing wind in a coronal structure.

79 ASTRONOMY AND ASTROPHYSICS↗

GALAXY-WIDE SHOCKS IN LATE-MERGER STAGE LUMINOUS INFRARED GALAXIES

We present an integral field spectroscopic study of two nearby luminous infrared galaxies (LIRGs), IC 1623 and NGC 3256, which exhibit evidence of widespread shock excitation induced by ongoing merger activity. We show the importance of carefully separating excitation due to shocks versus excitation by H II regions and the usefulness of integral field unit data in interpreting the complex processes in LIRGs. Our analysis focuses primarily on the emission line gas, which is extensive in both systems and is a result of the abundant ongoing star formation as well as widespread LINER-like excitation from shocks. We use emission line ratio maps, line kinematics, line-ratio diagnostics, and new models as methods for distinguishing and analyzing shocked gas in these systems. We discuss how our results inform the merger sequence associated with local U/LIRGs and the impact that widespread shock excitation has on the interpretation of emission line spectra and derived quantities of both local and high-redshift galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

OBLIQUE SHOCKS AS THE ORIGIN OF RADIO TO GAMMA-RAY VARIABILITY IN ACTIVE GALACTIC NUCLEI

The 'shock in jet' model for centimeter-waveband blazar variability is revisited, allowing for arbitrary shock orientation with respect to the jet flow direction, and both random and ordered magnetic field. It is shown that oblique shocks can explain events with swings in polarization position angle much less than the 90{sup 0} associated with transverse structures, while retaining the general characteristics of outbursts, including spectral behavior and level of peak percentage polarization. Models dominated by a force-free, minimum energy magnetic field configuration (essentially helical) display a shallow rise in percentage polarization and frequency-dependent swing in polarization position angle not in agreement with the results of single-dish monitoring observations, implying that the field is predominantly random in the quiescent state. Outbursts well explained by the 'shock in jet' model are present during {gamma}-ray flaring in several sources, supporting the idea that shock events are responsible for activity from the radio to {gamma}-ray bands.

79 ASTRONOMY AND ASTROPHYSICS↗

TEMPERATURE EQUILIBRATION BEHIND THE SHOCK FRONT: AN OPTICAL AND X-RAY STUDY OF RCW 86

We study the electron-proton temperature equilibration behind several shocks of the RCW 86 supernova remnant. To measure the proton temperature, we use published and new optical spectra, all from different locations on the remnant. For each location, we determine the electron temperature from X-ray spectra, and correct for temperature equilibration between the shock front and the location of the X-ray spectrum. We confirm the result of previous studies that the electron and proton temperatures behind shock fronts are consistent with equilibration for slow shocks and deviate for faster shocks. However, we cannot confirm the previously reported trend of T{sub e}/T{sub p}{proportional_to}1/v{sup 2}{sub s}.

79 ASTRONOMY AND ASTROPHYSICS↗

RELATIVISTIC SHOCK BREAKOUTS-A VARIETY OF GAMMA-RAY FLARES: FROM LOW-LUMINOSITY GAMMA-RAY BURSTS TO TYPE Ia SUPERNOVAE

The light from a shock breakout of stellar explosions, which carries a wealth of information, strongly depends on the shock velocity at the time of the breakout. The emission from Newtonian breakouts, typical in regular core-collapse supernovae (SNe), has been explored extensively. However, a large variety of explosions result in mildly or ultrarelativistic breakouts, where the observed signature is unknown. Here we calculate the luminosity and spectrum produced by relativistic breakouts. In order to do so, we improve the analytic description of relativistic radiation-mediated shocks and follow the system from the breakout itself, through the planar phase and into the spherical phase. We limit our calculation to cases where the post-breakout acceleration of the gas ends during the planar phase (i.e., the final gas Lorentz factor {approx}< 30). We find that spherical relativistic breakouts produce a flash of gamma rays with energy, E{sub bo}, temperature, T{sub bo}, and duration, t{sup obs} b{sub o}, that provide the breakout radius ( Almost-Equal-To 5 R{sub Sun }(t{sup obs}{sub bo}/10 s)(T{sub bo}/50 keV){sup 2}) and the Lorentz factor ( Almost-Equal-To T{sub bo}/50 keV). They also always satisfy a relativistic breakout relation (t{sup obs}{sub bo}/20 s) {approx} (E{sub bo}/10{sup 46} erg){sup 1/2}(T{sub bo}/50 keV){sup -2.68}. The breakout flare is typically followed, on longer timescales, by X-rays that carry a comparable energy. We apply our model to a variety of explosions, including Type Ia and .Ia SNe, accretion-induced collapse, energetic SNe, and gamma-ray bursts (GRBs). We find that all these events produce detectable gamma-ray signals, some of which may have already been seen. Some particular examples are: (1) relativistic shock breakouts provide a natural explanation to the energy, temperature, and timescales of low-luminosity GRBs. Indeed, all observed low-luminosity GRBs satisfy the relativistic breakout relation. (2) Nearby broad-line Type Ib/c (like SN 2002ap) may produce a detectable {gamma}-ray signal. (3) Galactic Type Ia SNe may produce detectable {gamma}-ray flares. We conclude that relativistic shock breakouts provide a generic process for the production of gamma-ray flares.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE AMPLIFICATION OF MAGNETIC FIELD BY A SUPERNOVA BLAST SHOCK WAVE IN A TURBULENT MEDIUM

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the amplification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar plasma. The blast wave is driven by injecting high pressure in the simulation domain. The interstellar magnetic field can be amplified by two different processes, occurring in different regions. One is facilitated by the fluid vorticity generated by the 'rippled' shock front interacting with the background turbulence. The resulting turbulent flow keeps amplifying the magnetic field, consistent with earlier work. The other process is facilitated by the growth of the Rayleigh-Taylor instability at the contact discontinuity between the ejecta and the shocked medium. This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We investigate the dependence of the amplification on numerical parameters such as grid-cell size and on various physical parameters. We show that the magnetic field has a characteristic radial profile such that the downstream magnetic field gets progressively stronger away from the shock. This is because the downstream magnetic field needs a finite time to reach the efficient amplification, and will get further amplified in the Rayleigh-Taylor region. In our simulation, we do not observe a systematic strong magnetic field within a small distance to the shock. This indicates that if the magnetic-field amplification in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional instabilities, plasma kinetics, and/or cosmic ray effect may need to be considered to explain the strong magnetic field in supernova remnants.

79 ASTRONOMY AND ASTROPHYSICS↗

X-RAY EMISSION LINE PROFILES FROM WIND CLUMP BOW SHOCKS IN MASSIVE STARS

The consequences of structured flows continue to be a pressing topic in relating spectral data to physical processes occurring in massive star winds. In a preceding paper, our group reported on hydrodynamic simulations of hypersonic flow past a rigid spherical clump to explore the structure of bow shocks that can form around wind clumps. Here we report on profiles of emission lines that arise from such bow shock morphologies. To compute emission line profiles, we adopt a two-component flow structure of wind and clumps using two 'beta' velocity laws. While individual bow shocks tend to generate double-horned emission line profiles, a group of bow shocks can lead to line profiles with a range of shapes with blueshifted peak emission that depends on the degree of X-ray photoabsorption by the interclump wind medium, the number of clump structures in the flow, and the radial distribution of the clumps. Using the two beta law prescription, the theoretical emission measure and temperature distribution throughout the wind can be derived. The emission measure tends to be power law, and the temperature distribution is broad in terms of wind velocity. Although restricted to the case of adiabatic cooling, our models highlight the influence of bow shock effects for hot plasma temperature and emission measure distributions in stellar winds and their impact on X-ray line profile shapes. Previous models have focused on geometrical considerations of the clumps and their distribution in the wind. Our results represent the first time that the temperature distribution of wind clump structures are explicitly and self-consistently accounted for in modeling X-ray line profile shapes for massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗

PARTICLE ACCELERATION AT A FLARE TERMINATION SHOCK: EFFECT OF LARGE-SCALE MAGNETIC TURBULENCE

We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations-known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks-are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ({delta}B{sup 2}/B 2{sub 0} {approx} 0.3), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

EFFECTS OF LEAKAGE NEUTRAL PARTICLES ON SHOCKS

In this paper, we investigate effects of neutral particles on shocks propagating into the partially ionized medium. We find that for 120 km s{sup -1} < u {sub sh} < 3000 km s{sup -1} (u {sub sh} is the shock velocity), about 10% of upstream neutral particles leak into the upstream region from the downstream region. Moreover, we investigate how the leakage neutral particles affect the upstream structure of the shock and particle accelerations. Using four-fluid approximations (upstream ions, upstream neutral particles, leakage neutral particles, and pickup ions), we provide analytical solutions of the precursor structure due to leakage neutral particles. It is shown that the upstream flow is decelerated in the precursor region and the shock compression ratio becomes smaller than without leakage neutral particles, but the total compression ratio does not change. Even if leakage of neutral particles is small (a few percent of total upstream particles), this smaller compression ratio of the shock can explain steep gamma-ray spectra from young supernova remnants. Furthermore, leakage neutral particles could amplify the magnetic field and heat the upstream region.

79 ASTRONOMY AND ASTROPHYSICS↗