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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
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 OSTI2013

EVIDENCE FOR CO SHOCK EXCITATION IN NGC 6240 FROM HERSCHEL SPIRE SPECTROSCOPY

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J = 4 - 3 through J = 13 - 12, 6 H{sub 2}O rotational lines, and [C I] and [N II] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H{sub 2} v = 1-0 S(1) and v = 2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v{sub s} = 10 km s{sup -1} and a pre-shock density n{sub H} = 5 Multiplication-Sign 10{sup 4} cm{sup -3}. We find that the solution best fitting the H{sub 2} lines is degenerate. The shock velocities and number densities range between v{sub s} = 17-47 km s{sup -1} and n{sub H} = 10{sup 7}-5 Multiplication-Sign 10{sup 4} cm{sup -3}, respectively. The H{sub 2} lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s{sup -1}) shocks while only a small fraction ({approx}< 1%) of the interstellar medium is exposed to the high-velocity shocks. This implies that the gas is rapidly losing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2025

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

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

DNS of strong shock and turbulence interactions with thermochemical non-equilibrium effects

The underlying physics in shock and turbulence interaction is essential for a better understanding of many natural processes as well as scientific and engineering applications, such as volcanic eruption, supernova explosion, detonation, medical application of shock wave lithotripsy to break up kidney stones, and energy application of the implosion of a cryogenic fuel capsule for inertial confinement fusion. One of the fundamental building blocks in these complex processes and applications is the canonical problem of isotropic turbulence and a normal shock. Unfortunately, even this fundamental problem is not well understood for strong shocks. Recent direct numerical simulation (DNS) results of perfect gas flow showed some new trends in turbulent statistics as mean Mach number is increased. In this paper, we first conduct extensive DNS studies on canonical strong shock and turbulence interaction problem of perfect gas flow with mean Mach numbers ranging from 2 to 30. The results show that strong shock waves have a significant effect on turbulent statistics. However, gas temperature increases dramatically after strong shocks so that numerical simulations based on perfect gas flow may not be enough. The effects of thermochemical nonequilibrium flow including internal energy excitations, translation-vibration energy relaxation, and chemical reactions among different species need to be considered. Therefore, we have developed and tested a new high-order shock-fitting solver for non-equilibrium flow simulations based on the 5-species air chemistry and recently thermal non-equilibrium models.

Wang, Xiaowen↗
DOE OSTI2017

The Presence of Turbulent and Ordered Local Structure within the ICME Shock-sheath and Its Contribution to Forbush Decrease

The transient interplanetary disturbances evoke short-time cosmic-ray flux decrease, which is known as Forbush decrease. The traditional model and understanding of Forbush decrease suggest that the sub-structure of an interplanetary counterpart of coronal mass ejection (ICME) independently contributes to cosmic-ray flux decrease. These sub-structures, shock-sheath, and magnetic cloud (MC) manifest as classical two-step Forbush decrease. The recent work by Raghav et al. has shown multi-step decreases and recoveries within the shock-sheath. However, this cannot be explained by the ideal shock-sheath barrier model. Furthermore, they suggested that local structures within the ICME’s sub-structure (MC and shock-sheath) could explain this deviation of the FD profile from the classical FD. Therefore, the present study attempts to investigate the cause of multi-step cosmic-ray flux decrease and respective recovery within the shock-sheath in detail. A 3D-hodogram method is utilized to obtain more details regarding the local structures within the shock-sheath. This method unambiguously suggests the formation of small-scale local structures within the ICME (shock-sheath and even in MC). Moreover, the method could differentiate the turbulent and ordered interplanetary magnetic field (IMF) regions within the sub-structures of ICME. The study explicitly suggests that the turbulent and ordered IMF regions within the shock-sheath do influence cosmic-ray variations differently.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗
NASA NTRS1971

Shock-metamorphic effects in the Luna-16 soil sample from Mare Fecunditatis

The results of intensive studies indicate that shock-metamorphic effects, characteristic of meteorite impact and virtually identical to those observed in Apollo samples, are common in fragments of the Luna-16 soil sample from Mare Fecunditatis. Two types of shock effects are present: (1) deformation and partial melting features in rock and mineral fragments (1-2 percent of fragments); and (2) heterogeneous glasses and glassy breccias produced by shock melting (70-80 percent of fragments). Shock effects were observed in pyroxene (deformation twin lamellae, multiple planar shock lamellae, extreme mosaicism, partial isotropization); in plagioclase (planar shock lamellae, complete isotropization to form maskelynite); and in basalt fragments (plagioclase isotropization, selective partial melting). The glasses exhibit several characteristics of shock melting, especially: (1) diversity in chemical composition; (2) association with shock mineral fragments and Ni-Fe spherules; and (3) heterogeneous schlieren and incipient fusion of mineral inclusions. Two types of source rocks are present in the Luna-16 sample; basaltic (85-90 percent) and feldspathic (10-15 percent). The basaltic rocks are predominant and generally occur as unshocked fragments, indicating that they form the bedrock underlying Mare Fecunditatis.

French, B. M.↗
NASA NTRS1972

Observation of nonuniform structure of the earth's bow shock correlated with interplanetary field orientation.

Explorer 33 and 35 magnetometers, on the western and eastern flanks of the earth's bow shock, respectively, observed the boundary concurrently between 0130 and 0430 UT, Oct. 30, 1968. Contrasting shock structures were recorded. Explorer 35 saw a quiet abrupt shock, whereas Explorer 33 saw an irregular noisy boundary with much upstream wave activity. The interplanetary field was roughly in the average archimedean spiral angle and was therefore approximately tangent to the shock at Explorer 35 and normal to the shock at Explorer 33. Gross motions and variable tilting of the aberrated shock probably contributed to the peculiar sequence of shock crossings at the two spacecraft. The observations support a model of the shock in which perpendicular and oblique collisionless structures coexist and form a nonuniform magnetosheath outer boundary.

Greenstadt, E. W.↗