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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Time-resolved diffraction of shock-released SiO 2 and diaplectic glass formation

Understanding how rock-forming minerals transform under shock loading is critical for modeling collisions between planetary bodies, interpreting the significance of shock features in minerals and for using them as diagnostic indicators of impact conditions, such as shock pressure. To date, our understanding of the formation processes experienced by shocked materials is based exclusively on ex situ analyses of recovered samples. Formation mechanisms and origins of commonly observed mesoscale material features, such as diaplectic (i.e., shocked) glass, remain therefore controversial and unresolvable. Here in this paper we show in situ pump-probe X-ray diffraction measurements on fused silica crystallizing to stishovite on shock compression and then converting to an amorphous phase on shock release in only 2.4 ns from 33.6 GPa. Recovered glass fragments suggest permanent densification. These observations of real-time diaplectic glass formation attest that it is a back-transformation product of stishovite with implications for revising traditional shock metamorphism stages.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Dynamics and microinstabilities at perpendicular collisionless shock: A comparison of large-scale two-dimensional full particle simulations with different ion to electron mass ratio

Large-scale two-dimensional (2D) full particle-in-cell (PIC) simulations are carried out for studying the relationship between the dynamics of a perpendicular shock and microinstabilities generated at the shock foot. The structure and dynamics of collisionless shocks are generally determined by Alfven Mach number and plasma beta, while microinstabilities at the shock foot are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the ratio of the plasma-to-cyclotron frequency. With a fixed Alfven Mach number and plasma beta, the ratio of the upstream bulk velocity to the electron thermal velocity is given as a function of the ion-to-electron mass ratio. The present 2D full PIC simulations with a relatively low Alfven Mach number (M{sub A} ∼ 6) show that the modified two-stream instability is dominant with higher ion-to-electron mass ratios. It is also confirmed that waves propagating downstream are more enhanced at the shock foot near the shock ramp as the mass ratio becomes higher. The result suggests that these waves play a role in the modification of the dynamics of collisionless shocks through the interaction with shock front ripples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Shocks in unmagnetized plasma with a shear flow: Stability and magnetic field generation

A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Shock compression response of highly reactive Ni + Al multilayered thin foils

The shock-compression response of Ni + Al multilayered thin foils is investigated using laser-accelerated thin-foil plate-impact experiments over the pressure range of 2 to 11 GPa. The foils contain alternating Ni and Al layers (parallel but not flat) of nominally 50 nm bilayer spacing. The goal is to determine the equation of state and shock-induced reactivity of these highly reactive fully dense thin-foil materials. The laser-accelerated thin-foil impact set-up involved combined use of photon-doppler-velocimetry to monitor the acceleration and impact velocity of an aluminum flyer, and VISAR interferometry was used to monitor the back free-surface velocity of the impacted Ni + Al multilayered target. The shock-compression response of the Ni + Al target foils was determined using experimentally measured parameters and impedance matching approach, with error bars identified considering systematic and experimental errors. Meso-scale CTH shock simulations were performed using real imported microstructures of the cross-sections of the multilayered Ni + Al foils to compute the Hugoniot response (assuming no reaction) for correlation with their experimentally determined equation of state. It was observed that at particle velocities below ∼150 m/s, the experimentally determined equation of state trend matches the CTH-predicted inert response and is consistent with the observed unreacted state of the recovered Ni + Al target foils from this velocity regime. At higher particle velocities, the experimentally determined equation of state deviates from the CTH-predicted inert response. A complete and self-sustained reaction is also seen in targets recovered from experiments performed at these higher particle velocities. The deviation in the measured equation of state, to higher shock speeds and expanded volumes, combined with the observation of complete reaction in the recovered multilayered foils, confirmed via microstructure characterization, is indicative of the occurrence of shock-induced chemical reaction occurring in the time-scale of the high-pressure state. TEM characterization of recovered shock-compressed (unreacted) Ni + Al multilayered foils exhibits distinct features of constituent mixing revealing jetted layers and inter-mixed regions. These features were primarily observed in the proximity of the undulations present in the alternating layers of the Ni + Al starting foils, suggesting the important role of such instabilities in promoting shock-induced intermetallic-forming reactions in the fully dense highly exothermic multilayered thin foils.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Ion-acoustic shocks with self-regulated ion reflection and acceleration

Here, an analytic solution describing an ion-acoustic collisionless shock, self-consistently with the evolution of shock-reflected ions, is obtained. The solution extends the classic soliton solution beyond a critical Mach number, where the soliton ceases to exist because of the upstream ion reflection. The reflection transforms the soliton into a shock with a trailing wave and a foot populated by the reflected ions. The solution relates parameters of the entire shock structure, such as the maximum and minimum of the potential in the trailing wave, the height of the foot, as well as the shock Mach number, to the number of reflected ions. This relation is resolvable for any given distribution of the upstream ions. In this paper, we have resolved it for a simple “box” distribution. Two separate models of electron interaction with the shock are considered. The first model corresponds to the standard Boltzmannian electron distribution in which case the critical shock Mach number only insignificantly increases from M ≈ 1:6 (no ion reflection) to M ≈ 1:8 (substantial reflection). The second model corresponds to adiabatically trapped electrons. They produce a stronger increase, from M ≈ 3:1 to M ≈ 4:5 . The shock foot that is supported by the reflected ions also accelerates them somewhat further. A self-similar foot expansion into the upstream medium is described analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Plasma ion stratification by weak planar shocks

We derive fluid equations for describing steady-state planar shocks of a moderate strength (0 < M - 1 ≲ 1 with M the shock Mach number) propagating through an unmagnetized quasineutral collisional plasma comprising two separate ion species. In addition to the standard fluid shock quantities, such as the total mass density, mass-flow velocity, and electron and average ion temperatures, the equations describe shock stratification in terms of variations in the relative concentrations and temperatures of the two ion species along the shock propagation direction. We have solved these equations analytically for weak shocks (0 < M - 1 << 1), with the results depending on M, ratios of the ion masses and charges, and the upstream mass fraction of one of the ion species. These analytical results are instrumental for gaining understanding in the behavior of weak shocks, and they have been used to verify kinetic simulations of shocks in multi-ion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Ion species stratification within strong shocks in two-ion plasmas

We report strong collisional shocks in multi-ion plasmas are featured in many environments, with Inertial Confinement Fusion (ICF) experiments being one prominent example. Recent work [Keenan et al., Phys. Rev. E 96, 053203 (2017)] answered in detail a number of outstanding questions concerning the kinetic structure of steady-state, planar plasma shocks, e.g., the shock width scaling by the Mach number, M. However, it did not discuss shock-driven ion-species stratification (e.g., relative concentration modification and temperature separation). These are important effects since many recent ICF experiments have evaded explanation by standard, single-fluid, radiation-hydrodynamic (rad-hydro) numerical simulations, and shock-driven fuel stratification likely contributes to this discrepancy. Employing the state-of-the-art Vlasov-Fokker-Planck code, iFP, along with multi-ion hydro simulations and semi-analytics, we quantify the ion stratification by planar shocks with the arbitrary Mach number and the relative species concentration for two-ion plasmas in terms of ion mass and charge ratios. In particular, for strong shocks, we find that the structure of the ion temperature separation has a nearly universal character across ion mass and charge ratios. Lastly, we find that the shock fronts are enriched with the lighter ion species and the enrichment scales as M 4 for M»1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Phenomenology of reverse-shock emission in the optical afterglows of gamma-ray bursts

We use a parent sample of 118 gamma-ray burst (GRB) afterglows, with known redshift and host galaxy extinction, to separate afterglows with and without signatures of dominant reverse-shock (RS) emission and to determine which physical conditions lead to a prominent reverse-shock emission. We identify 10 GRBs with reverse-shock signatures: 990123, 021004, 021211, 060908, 061126, 080319B, 081007, 090102, 090424, and 130427A. By modeling their optical afterglows with reverse- and forward-shock analytic light curves and using Monte Carlo simulations, we estimate the parameter space of the physical quantities describing the ejecta and circumburst medium. We find that physical properties cover a wide parameter space and do not seem to cluster around any preferential values. Comparing the rest-frame optical, X-ray, and high-energy properties of the larger sample of non-RS-dominated GRBs, we show that the early-time (<1 ks) optical spectral luminosity, X-ray afterglow luminosity, and γ-ray energy output of our reverse-shock dominated sample do not differ significantly from the general population at early times. However, the GRBs with dominant reverse-shock emission have fainter than average optical forward-shock emission at late times (>10 ks). We find that GRBs with an identifiable reverse-shock component show a high magnetization parameter R {sub B} = ε{sub B,r}/ε{sub B,f} ∼ 2-10{sup 4}. Our results are in agreement with the mildly magnetized baryonic jet model of GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Shock generated vorticity in spark discharges

Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 μs), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the shock wave for a range of electrical energies to measure the shock properties and estimate shock velocity and curvature. These measurements are combined with a model to calculate the vorticity generated. The measurements show that the highest vorticity is generated near the peak shock curvature location, and the shock curvature and strength increase with electrical energy deposited. A comparison of the vorticity estimated from the model to vorticity measurements from stereoscopic particle image velocimetry shows the results to be statistically equivalent. Furthermore, this suggests that the shock curvature and velocity contribute to the vortex rings induced by spark discharges.

42 ENGINEERING↗
DOE OSTI2018

Collisional and collisionless shocks

Shock waves are one of the most common plasma phenomena. They play a significant role in astrophysical and magnetospheric environments, as well as in laboratory plasmas. The paper focuses on two issues that underline important similarities and dissimilarities between collisional and collisionless shocks. The two types of shocks are similar in that they can be formed by the overtaking mechanism. We demonstrate that geometrical structures appearing early in the process of shock formation from the generic smooth initial state are quite universal and similar in both cases. Characterization of these structures (dubbed here the “mussel-shell”) is presented. Alongside with these similarities, there exist also significant differences between the two types of shocks. The classical collisional shocks usually connect two well-defined equilibrium states (those before and after the shock transition). Each of these equilibrium states is characterized by the thermodynamic parameters of density, temperature, and pressure whose upstream and downstream values are related by the continuity of mass, momentum, and energy flux. In the collisionless plasmas, however, this description does not work: the final state can be any of the much broader class of states constrained only by the requirement of being stable with respect to collisionless plasma instabilities. To find a final state (which does not have to be Maxwellian) one now has to follow the evolution of the system through the whole transition, and a lot of universality is lost. In some situations even the separation of scales between the global flow and shock transition may be lost. Taken together, these phenomena reveal a fascinating interplay of hydrodynamics, statistical mechanics, and plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Shock-induced chemistry of phenylacetylene

Gas gun-driven shock compression experiments of phenylacetylene using embedded electromagnetic gauging were used to obtain in situ particle velocity wave profiles at multiple Lagrangian positions at several shock input conditions. At shock conditions above 6 GPa, the input shock wave evolved over time and distance into a complex multiple wave structure due to shock-driven chemical reactions. The 3-wave structure was marked by a fast risetime 2nd wave, slower risetime 3rd wave, and unusual wave dynamics in the 1st wave. From the measured shock and particle velocities, the 1st wave, and intermediate and final product states associated with the chemical reactions were determined. A thermodynamically complete unreacted equation of state was calibrated to estimate the temperature rise along the shock locus. Use of this EOS with the measured 2nd and 3rd wave risetimes yielded highly statesensitive global reaction rates as a function of the shock locus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

On the importance of progenitor asymmetry to shock revival in core-collapse supernovae

The progenitor stars of core-collapse supernovae (CCSNe) are asymmetrically fluctuating due to turbulent convections in the late stages of their lives. The progenitor asymmetry at the pre-supernova stage has recently caught the attention as a new ingredient to facilitate shock revival in the delayed neutrino-heating mechanism. In this paper, we investigate the importance of the progenitor asymmetries to shock revival with a semi-analytical approach. Free parameters were chosen such that the time evolution of shock radii and mass accretion rates are compatible with the results of detailed numerical simulations of CCSNe in spherical symmetry. We first estimate the amplitude of asymmetries required for the shock revival by the impulsive change of pre-shock flows in the context of neutrino-heating mechanism, and then convert the amplitude to the corresponding amplitude in the pre-supernova phase by taking into account the growth of asymmetries during infall. Further, we apply our model to various types of progenitors and find that the requisite amplitude of pre-supernova asymmetry is roughly three times larger than the prediction by current stellar evolution models unless other additional physical ingredients such as multidimensional fluid instabilities and turbulent convections in post-shock flows aid shock revival. We thus conclude that progenitor asymmetries cannot trigger the shock revival by the impulsive way but rather play a supplementary role in reality.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Probabilistic Assessment of Damage from Same Shock Response Spectra due to Variations in Damping

Interpretation of field data from shock tests and subsequent assessment of product safety margins via laboratory testing are based on the shock response spectra (SRS). The SRS capture how a single degree of freedom (SDOF) structure responds to the shock at differing frequencies and, therefore, no longer contain the duration or other temporal parameters pertaining to the shock. A single duration can often be included in the technical specification or in the recreation of acceleration vs. time history from the specified SRS; however, there is little basis for that beyond technical judgment. The loss of such temporal information can result in the recreated SRS being the same while its effect on a system or component can be different. This paper attempts to quantify this deficiency as well as propose a simple method of capturing damping from shock waves that can allow the original waveform to be more accurately reconstructed from the SRS. In this study the decay rate associated with various frequencies that comprise the overall shock was varied. This variation in the decay rate leads to a variation in the acceleration vs. time history, which can be correlated to a “Damage Index” that captures the fatigue damage imparted to the object under shock. Several waveforms that have the same SRS but varying rates of decay for either high- or low-frequency components of the shock were investigated. The resulting variation in stress cycles and Damage Index is discussed in the context of the lognormal distribution of fatigue failure data. It is proposed that, along with the SRS, the decay rate is also captured to minimize the discrepancy between field data and representative laboratory tests.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Shock compression modeling of metallic single crystals: comparison of finite difference, steady wave, and analytical solutions

Background: The shock response of metallic single crystals can be captured using a micro-mechanical description of the thermoelastic-viscoplastic material response; however, using a such a description within the context of traditional numerical methods may introduce a physical artifacts. Advantages and disadvantages of complex material descriptions, in particular the viscoplastic response, must be framed within approximations introduced by numerical methods. Methods: Three methods of modeling the shock response of metallic single crystals are summarized: finite difference simulations, steady wave simulations, and algebraic solutions of the Rankine-Hugoniot jump conditions. For the former two numerical techniques, a dislocation density based framework describes the rate- and temperature-dependent shear strength on each slip system. For the latter analytical technique, a simple (two-parameter) rate- and temperature-independent linear hardening description is necessarily invoked to enable simultaneous solution of the governing equations. For all models, the same nonlinear thermoelastic energy potential incorporating elastic constants of up to order 3 is applied. Results: Solutions are compared for plate impact of highly symmetric orientations (all three methods) and low symmetry orientations (numerical methods only) of aluminum single crystals shocked to 5 GPa (weak shock regime) and 25 GPa (overdriven regime). Conclusions: For weak shocks, results of the two numerical methods are very similar, regardless of crystallographic orientation. For strong shocks, artificial viscosity affects the finite difference solution, and effects of transverse waves for the lower symmetry orientations not captured by the steady wave method become important. The analytical solution, which can only be applied to highly symmetric orientations, provides reasonable accuracy with regards to prediction of most variables in the final shocked state but, by construction, does not provide insight into the shock structure afforded by the numerical methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Data Based Shock Test Specification

It is vital that avionic packages used for testing and certifying the reliability and safety of U.S. nuclear weapons with platform aircraft survive exposure to shock environments during transportation and delivery. The objective of this research was to characterize the response to these transportation shock environments delivering accurate shock test specifications in order to set laboratory programming material and device certification rigor. Responses to shock events were analyzed in the frequency domain via the Shock Response Spectrum (SRS). Shocks were then grouped based on respective behavior of maximum response accelerations which were pseudorandomly resampled and compared to test data to form test specifications based on the MinerPalmgren hypothesis. In addition to discovering significant over testing in current shock specifications, a new systematic, data-driven approach to designing shock specifications was formulated.

42 ENGINEERING↗