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 523 records · Page 29

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗

INTERSTELLAR H{sub 2}O MASERS FROM J SHOCKS

We present a model in which the 22 GHz H{sub 2}O masers observed in star-forming regions occur behind shocks propagating in dense regions (preshock density n{sub 0} {approx} 10{sup 6}-10{sup 8} cm{sup -3}). We focus on high-velocity (v{sub s} {approx}> 30 km s{sup -1}) dissociative J shocks in which the heat of H{sub 2} re-formation maintains a large column of {approx}300-400 K gas; at these temperatures the chemistry drives a considerable fraction of the oxygen not in CO to form H{sub 2}O. The H{sub 2}O column densities, the hydrogen densities, and the warm temperatures produced by these shocks are sufficiently high to enable powerful maser action. The observed brightness temperatures (generally {approx} 10{sup 11}-10{sup 14} K) are the result of coherent velocity regions that have dimensions in the shock plane that are 10-100 times the shock thickness of {approx}10{sup 13} cm. The masers are therefore beamed toward the observer, who typically views the shock ''edge-on'', or perpendicular to the shock velocity; the brightest masers are then observed with the lowest line-of-sight velocities with respect to the ambient gas. We present numerical and analytic studies of the dependence of the maser inversion, the resultant brightness temperature, the maser spot size and shape, the isotropic luminosity, and the maser region magnetic field on the shock parameters and the coherence path length; the overall result is that in galactic H{sub 2}O 22 GHz masers, these observed parameters can be produced in J shocks with n{sub 0} {approx} 10{sup 6}-10{sup 8} cm{sup -3} and v{sub s} {approx} 30-200 km s{sup -1}. A number of key observables such as maser shape, brightness temperature, and global isotropic luminosity depend only on the particle flux into the shock, j = n{sub 0} v{sub s} , rather than on n{sub 0} and v{sub s} separately.

79 ASTRONOMY AND ASTROPHYSICS↗

HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

79 ASTRONOMY AND ASTROPHYSICS↗

PARTICLE ACCELERATION AND WAVE EXCITATION IN QUASI-PARALLEL HIGH-MACH-NUMBER COLLISIONLESS SHOCKS: PARTICLE-IN-CELL SIMULATION

We herein investigate shock formation and particle acceleration processes for both protons and electrons in a quasi-parallel high-Mach-number collisionless shock through a long-term, large-scale, particle-in-cell simulation. We show that both protons and electrons are accelerated in the shock and that these accelerated particles generate large-amplitude Alfvénic waves in the upstream region of the shock. After the upstream waves have grown sufficiently, the local structure of the collisionless shock becomes substantially similar to that of a quasi-perpendicular shock due to the large transverse magnetic field of the waves. A fraction of protons are accelerated in the shock with a power-law-like energy distribution. The rate of proton injection to the acceleration process is approximately constant, and in the injection process, the phase-trapping mechanism for the protons by the upstream waves can play an important role. The dominant acceleration process is a Fermi-like process through repeated shock crossings of the protons. This process is a “fast” process in the sense that the time required for most of the accelerated protons to complete one cycle of the acceleration process is much shorter than the diffusion time. A fraction of the electrons are also accelerated by the same mechanism, and have a power-law-like energy distribution. However, the injection does not enter a steady state during the simulation, which may be related to the intermittent activity of the upstream waves. Upstream of the shock, a fraction of the electrons are pre-accelerated before reaching the shock, which may contribute to steady electron injection at a later time.

79 ASTRONOMY AND ASTROPHYSICS↗

The acceleration of electrons at collisionless shocks moving through a turbulent magnetic field

We perform a numerical-simulation study of the acceleration of electrons at shocks that propagate through a prespecified, kinematically defined turbulent magnetic field. The turbulence consists of broadband magnetic fluctuations that are embedded in the plasma and cover a range of wavelengths, the smallest of which is larger than the gyroradii of electrons that are initially injected into the system. We find that when the variance of the turbulent component of the upstream magnetic field is sufficiently large—σ{sup 2} ∼ 10 B{sub 0}{sup 2}, where B{sub 0} is the strength of the background magnetic field—electrons can be efficiently accelerated at a collisionless shock regardless of the orientation of the mean upstream magnetic field relative to the shock-normal direction. Since the local angle between the incident magnetic-field vector and the shock-normal vector can be quite large, electrons can be accelerated through shock-drift acceleration at the shock front. In the upstream region, electrons are mirrored back to the shock front leading to multiple shock encounters. Eventually the accelerated electrons are energetic enough that their gyroradii are of the same order as the wavelength of waves that are included in our description of the turbulent magnetic field. Our results are consistent with recent in situ observations at Saturn’s bow shock. This study may help us to understand the acceleration of electrons at shocks in space and astrophysical systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisionless Weibel shocks and electron acceleration in gamma-ray bursts

A study of collisionless external shocks in gamma-ray bursts is presented. The shock structure, electromagnetic field, and process of electron acceleration are assessed by performing a self-consistent 3D particle-in-cell simulation. In accordance with hydrodynamic shock systems, the shock consists of a reverse shock (RS) and forward shock separated by a contact discontinuity. The development and structure are controlled by the ion Weibel instability. The ion filaments are sources of strong transverse electromagnetic fields at both sides of the double shock structure over a length of 30–100 ion skin depths. Electrons are heated up to a maximum energy ϵ{sub ele}≈√(ϵ{sub b}), where ϵ is the energy normalized to the total incoming energy. Jet electrons are trapped in the RS transition region due to the presence of an ambipolar electric field and reflection by the strong transverse magnetic fields in the shocked region. In a process similar to shock surfing acceleration for ions, electrons experience drift motion and acceleration by ion filament transverse electric fields in the plane perpendicular to the shock propagation direction. Ultimately, accelerated jet electrons are convected back into the upstream.

79 ASTRONOMY AND ASTROPHYSICS↗

PHYSICAL CONDITIONS OF CORONAL PLASMA AT THE TRANSIT OF A SHOCK DRIVEN BY A CORONAL MASS EJECTION

We report here on the determination of plasma physical parameters across a shock driven by a coronal mass ejection using white light (WL) coronagraphic images and radio dynamic spectra (RDS). The event analyzed here is the spectacular eruption that occurred on 2011 June 7, a fast CME followed by the ejection of columns of chromospheric plasma, part of them falling back to the solar surface, associated with a M2.5 flare and a type-II radio burst. Images acquired by the Solar and Heliospheric Observatory/LASCO coronagraphs (C2 and C3) were employed to track the CME-driven shock in the corona between 2–12 R{sub ⊙} in an angular interval of about 110°. In this interval we derived two-dimensional (2D) maps of electron density, shock velocity, and shock compression ratio, and we measured the shock inclination angle with respect to the radial direction. Under plausible assumptions, these quantities were used to infer 2D maps of shock Mach number M{sub A} and strength of coronal magnetic fields at the shock's heights. We found that in the early phases (2–4 R{sub ⊙}) the whole shock surface is super-Alfvénic, while later on (i.e., higher up) it becomes super-Alfvénic only at the nose. This is in agreement with the location for the source of the observed type-II burst, as inferred from RDS combined with the shock kinematic and coronal densities derived from WL. For the first time, a coronal shock is used to derive a 2D map of the coronal magnetic field strength over intervals of 10 R{sub ⊙} altitude and ∼110° latitude.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMULATIONS AND THEORY OF ION INJECTION AT NON-RELATIVISTIC COLLISIONLESS SHOCKS

We use kinetic hybrid simulations (kinetic ions-fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also, by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗

Shock Propagation and Associated Particle Acceleration in the Presence of Ambient Solar-Wind Turbulence

The topic of this review paper is on the influence of solar wind turbulence on shock propagation and its consequence on the acceleration and transport of energetic particles at shocks. As the interplanetary shocks sweep through the turbulent solar wind, the shock surfaces fluctuate and ripple in a range of different scales. We discuss particle acceleration at rippled shocks in the presence of ambient solar-wind turbulence. This strongly affects particle acceleration and transport of energetic particles (both ions and electrons) at shock fronts. In particular, we point out that the effects of upstream turbulence is critical for understanding the variability of energetic particles at shocks. Moreover, the presence of pre-existing upstream turbulence significantly enhances the trapping near the shock of low-energy charged particles, including those near the thermal energy of the incident plasma, even when the shock propagates normal to the average magnetic field. Pre-existing turbulence, always present in space plasmas, provides a means for the efficient acceleration of low-energy particles and overcoming the well known injection problem at shocks.

, solar wind↗

Diffusive cosmic ray acceleration at relativistic shock waves with magnetostatic turbulence. II. Influence of a finite downstream medium

The diffusive acceleration of relativistic cosmic rays at parallel shock waves with magnetostatic turbulence and a finite size of the downstream medium is investigated. For ultrarelativistic shock speeds with Lorentz factor Γ{sub 1}≫1, both the differential momentum spectrum at the shock and the volume-integrated momentum spectrum are power-law distribution functions with different spectral indices as compared to the case of an infinitely extended downstream medium. However, the spectral differences are only modest as compared to the case of nonrelativistic shocks. The behavior of the momentum spectrum of shock-accelerated particles depends sensitively on the relativistic shock wave Peclet number G(p)=τ{sub D}(p)/τ{sub C}, i.e., the ratio between the diffusion and convection timescales of cosmic rays to propagate from the shock position to the downstream boundary z {sub 0}. For large values of G(p)≫1 the free-escape boundary has no influence on the effectiveness of particle acceleration, still providing a flat momentum power-law spectrum of the accelerated particles. In the opposite case of small Peclet numbers G(p)≪1 at all momenta, the momentum spectrum at the shock steepens to the greater spectral index ξ{sub 0}=3−s+(3.18/Γ{sub 1}), whereas the volume-integrated momentum spectrum flattens by the same factor 2−s for its power-law spectral index, where s denotes the spectral index of the downstream power spectrum of magnetostatic turbulence. This effectiveness of relativistic shocks in generating flat power-law momentum spectra irrespective of the Peclet number G(p) differs completely from the behavior of nonrelativistic shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

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

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.

79 ASTRONOMY AND ASTROPHYSICS↗

STRUCTURE, PROPAGATION, AND EXPANSION OF A CME-DRIVEN SHOCK IN THE HELIOSPHERE: A REVISIT OF THE 2012 JULY 23 EXTREME STORM

We examine the structure, propagation, and expansion of the shock associated with the 2012 July 23 extreme coronal mass ejection. Characteristics of the shock determined from multi-point imaging observations are compared to in situ measurements at different locations and a complex radio type II burst, which according to our definition has multiple branches that may not all be fundamental-harmonic related. The white-light shock signature can be modeled reasonably well by a spherical structure and was expanding backward even on the opposite side of the Sun. The expansion of the shock, which was roughly self-similar after the first ∼1.5 hr from launch, largely dominated over the translation of the shock center for the time period of interest. Our study also suggests a bow-shock morphology around the nose at later times due to the outward motion in combination with the expansion of the ejecta. The shock decayed and failed to reach Mercury in the backward direction and the Solar Terrestrial Relations Observatory B ( STEREO B ) and Venus in the lateral directions, as indicated by the imaging and in situ observations. The shock in the nose direction, however, may have persisted to the far outer heliosphere, with predicted impact on Dawn around 06:00 UT on July 25 and on Jupiter around 23:30 UT on July 27 by a magnetohydrodynamic model. The type II burst shows properties generally consistent with the spatial/temporal variations of the shock deduced from imaging and in situ observations. In particular, the low-frequency bands agree well with the in situ measurements of a very low density ahead of the shock at STEREO A .

79 ASTRONOMY AND ASTROPHYSICS↗

Particle Acceleration in Two Converging Shocks

Observations by spacecraft such as ACE , STEREO , and others show that there are proton spectral “breaks” with energy E {sub br} at 1–10 MeV in some large CME-driven shocks. Generally, a single shock with the diffusive acceleration mechanism would not predict the “broken” energy spectrum. The present paper focuses on two converging shocks to identify this energy spectral feature. In this case, the converging shocks comprise one forward CME-driven shock on 2006 December 13 and another backward Earth bow shock. We simulate the detailed particle acceleration processes in the region of the converging shocks using the Monte Carlo method. As a result, we not only obtain an extended energy spectrum with an energy “tail” up to a few 10 MeV higher than that in previous single shock model, but also we find an energy spectral “break” occurring on ∼5.5 MeV. The predicted energy spectral shape is consistent with observations from multiple spacecraft. The spectral “break,” then, in this case is caused by the interaction between the CME shock and Earth’s bow shock, and otherwise would not be present if Earth were not in the path of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗

The Dynamics of Very High Alfvén Mach Number Shocks in Space Plasmas

Astrophysical shocks, such as planetary bow shocks or supernova remnant shocks, are often in the high or very-high Mach number regime, and the structure of such shocks is crucial for understanding particle acceleration and plasma heating, as well inherently interesting. Recent magnetic field observations at Saturn’s bow shock, for Alfvén Mach numbers greater than about 25, have provided evidence for periodic non-stationarity, although the details of the ion- and electron-scale processes remain unclear due to limited plasma data. High-resolution, multi-spacecraft data are available for the terrestrial bow shock, but here the very high Mach number regime is only attained on extremely rare occasions. Here we present magnetic field and particle data from three such quasi-perpendicular shock crossings observed by the four-spacecraft Cluster mission. Although both ion reflection and the shock profile are modulated at the upstream ion gyroperiod timescale, the dominant wave growth in the foot takes place at sub-proton length scales and is consistent with being driven by the ion Weibel instability. The observed large-scale behavior depends strongly on cross-scale coupling between ion and electron processes, with ion reflection never fully suppressed, and this suggests a model of the shock dynamics that is in conflict with previous models of non-stationarity. Thus, the observations offer insight into the conditions prevalent in many inaccessible astrophysical environments, and provide important constraints for acceleration processes at such shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-Spacecraft Observations of Shocklets at an Interplanetary Shock

Interplanetary (IP) shocks are fundamental building blocks of the heliosphere, and the possibility to observe them in situ is crucial to address important aspects of energy conversion for a variety of astrophysical systems. Steepened waves known as shocklets are known to be important structures of planetary bow shocks, but they are very rarely observed related to IP shocks. We present here the first multi-spacecraft observations of shocklets observed by upstream of an unusually strong IP shock observed on 3 No v ember 2021 by several spacecraft at L1 and near-Earth solar wind. The same shock was detected also by radially aligned Solar Orbiter at 0.8 AU from the Sun, but no shocklets were identified from its data, introducing the possibility to study the environment in which shocklets developed. The Wind spacecraft has been used to characterize the shocklets, associated with pre-conditioning of the shock upstream by decelerating incoming plasma in the shock normal direction. Finally, using the Wind observations together with ACE and DSCOVR spacecraft at L1, as well as THEMIS B and THEMIS C in the near-Earth solar wind, the portion of interplanetary space filled with shocklets is addressed, and a lower limit for its extent is estimated to be of about 110 R E in the shock normal direction and 25 R E in the directions transverse to the shock normal. Using multiple spacecraft also reveals that for this strong IP shock, shocklets are observed for a large range of local obliquity estimates (9° –64°).

plasmas↗

Thermophysical properties of multi-shock compressed dense argon

In contrast to the single shock compression state that can be obtained directly via experimental measurements, the multi-shock compression states, however, have to be calculated with the aid of theoretical models. In order to determine experimentally the multiple shock states, a diagnostic approach with the Doppler pins system (DPS) and the pyrometer was used to probe multiple shocks in dense argon plasmas. Plasma was generated by a shock reverberation technique. The shock was produced using the flyer plate impact accelerated up to ∼6.1 km/s by a two-stage light gas gun and introduced into the plenum argon gas sample, which was pre-compressed from the environmental pressure to about 20 MPa. The time-resolved optical radiation histories were determined using a multi-wavelength channel optical transience radiance pyrometer. Simultaneously, the particle velocity profiles of the LiF window was measured with multi-DPS. The states of multi-shock compression argon plasma were determined from the measured shock velocities combining the particle velocity profiles. We performed the experiments on dense argon plasmas to determine the principal Hugonoit up to 21 GPa, the re-shock pressure up to 73 GPa, and the maximum measure pressure of the fourth shock up to 158 GPa. The results are used to validate the existing self-consistent variational theory model in the partial ionization region and create new theoretical models.

74 ATOMIC AND MOLECULAR PHYSICS↗