Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK WAVES”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ion Acceleration in Fermi-LAT Behind-the-Limb Solar Flares: The Role of Coronal Shock Waves

Here, we investigate the relationship between the gamma-ray emission measured with Fermi-LAT and radio signatures of coronal shock waves in four behind-the-limb (BTL) solar flares. All events were associated with metric type II radio burst. Both start and end times of the radio bursts were synchronized with the gamma-ray emission. The type II bursts associated with the BTL gamma-ray flares had higher speeds and lower formation heights than those of an average sample. These findings support the notion that the highly relativistic ions that produce the gamma-rays in BTL flares are accelerated at CME-driven propagating coronal shock waves rather than in large-scale coronal loops.

Corona↗
DOE OSTI2022

Shock wave formation from head-on collision of two subsonic vortex rings

Vortex ring collisions have attracted intense interest in both water and air studies (Baird in Proc R Soc Lond Ser Math Phys Sci 409:59–65, 1987, Poudel et al. in Phys Fluids 33:096105, 2021, Lim and Nickels in Nature 357:225, 1992, New et al. in Exp Fluids 57:109, 2016, Suzuki et al. in Geophys Res Lett 34, 2007, Yan et al. in J Fluids Eng 140:054502, 2018, New et al. in J Fluid Mech 899, 2020, Cheng et al. in Phys Fluids 31:067107, 2019, Hernández and Reyes in 29:103604, 2017, Mishra et al. in Phys Rev Fluids, 2021, Zednikova et al. in Chem Eng Technol 42:843–850, 2019, Kwon et al. in Nature 600:64–69, 2021). These toroidal structures spin around a central axis and travel in the original direction of impulse while spinning around the core until inertial forces become predominant causing the vortex flow to spontaneously decay to turbulence (Vortex Rings, https://projects.iq.harvard.edu/smrlab/vortex-rings). Previous studies have shown the collision of subsonic vortex rings resulting in reconnected vortex rings, but the production of a shock wave from the collision has not been demonstrated visibly (Lim and Nickels in Nature 357:225, 1992, Cheng et al. in Phys Fluids 31:067107, 2019). Here we present the formation of a shock wave due to the collision of explosively formed subsonic vortex rings. As the vortex rings travel at Mach 0.66 toward the collision point, they begin to trap high pressure air between them. Upon collision, high pressure air was imploded and released radially away from the axis of the collision, generating a visible shock wave traveling through and away from the colliding vortices at Mach 1.22. Our results demonstrate a pressure gradient with high pressure release creating a shock wave. We anticipate our study to be a starting point for more explosively formed vortex collisions. For example, explosives with different velocities of detonation could be tested to produce vortex rings of varying velocities.

42 ENGINEERING↗
DOE OSTI2021

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An Efficient B-Spline Lagrangian/Eulerian Method for Compressible Flow, Shock Waves, and Fracturing Solids

This study presents a new method for modeling the interaction between compressible flow, shock waves, and deformable structures, emphasizing destructive dynamics. Extending advances in time-splitting compressible flow and the Material Point Methods (MPM), we develop a hybrid Eulerian and Lagrangian/Eulerian scheme for monolithic flow-structure interactions. We adopt the second-order WENO scheme to advance the continuity equation. To stably resolve deforming boundaries with sub-cell particles, we propose a blending treatment of reflective and passable boundary conditions inspired by the theory of porous media. The strongly coupled velocity-pressure system is discretized with a new mixed-order finite element formulation employing B-spline shape functions. Shock wave propagation, temperature/density-induced buoyancy effects, and topology changes in solids are unitedly captured.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics of Three-Dimensional Shock-Wave/Boundary-Layer Interactions

Advances in measuring and understanding separated, nominally two-dimensional (2D) shock-wave/turbulent-boundary-layer interactions (STBLI) have triggered recent campaigns focused on three-dimensional (3D) STBLI, which display far greater configuration diversity. Nonetheless, unifying properties emerge for semi-infinite interactions, taking the form of conical asymptotic behavior where shock-generator specifics become insignificant. The contrast between 2D and 3D separation is substantial; the skewed vortical structure of 3D STBLI reflects the essentially 2D influence of the boundary layer on the 3D character of the swept shock. As with 2D STBLI, conical interactions engender prominent spectral content below that of the turbulent boundary layer. However, the uniform separation length scale, which is crucial to normalizing the lowest-frequency dynamics in 2D STBLI, is absent. Comparatively, the spectra of 3D STBLI are more representative of the mid-frequency, convective, shear-layer dynamics in 2D, while phenomena associated with 2D separation-shock breathing are muted. Asymptotic behavior breaks down in many regions important to 3D-STBLI dynamics, occurring in a configuration-dependent manner. Aspects of inceptive regions near shock generators and symmetry planes are reviewed. Focused efforts toward 3D modal and nonmodalanalyses, moving-shock/boundary-layer interactions, fluid/structure interactions, and flow control are suggested as directions for future work.

42 ENGINEERING↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2020

Airfoil leading edge blowing to control bow shock waves

Abstract This manuscript presents a detailed characterization of active control of bow shock waves via leading edge injection, including subsonic coolant ejection and the appearance of Coanda effects. The flow phenomena occurring at steady and pulsating flow injection regimes were analyzed using steady and unsteady two-dimensional Reynolds-Averaged Navier Stokes, leading to a precise evaluation of the thermal load and drag reductions. Steady supersonic injection yields the largest abatement in thermal load and aerodynamic drag, while subsonic or fluctuating ones can also provide significant improvements at reduced cooling mass flow rates. Furthermore, a Coanda effect, causing a non-symmetric flow topology, was observed and analyzed for reduced injection port size. This Coanda effect is due to the sudden expansion happening from the injection port to the main flow and it causes the flow topology at the leading edge to become non-symmetric despite the complete symmetry of the problem. This is the first time in the literature such a phenomenon is documented for a supersonic airfoil leading edge injection. Furthermore, it enables the design of novel flow control strategies for the leading edge shock topology and flow structures in supersonic flows.

42 ENGINEERING↗
DOE OSTI2026

Time-resolved unsteady shock motion and spectral analysis of a hypersonic shock-wave/boundary-layer interaction over a double-cone

Direct experimental measurements of double-cone shock-wave/boundary-layer interaction dynamics at Mach 7.2 were collected using a schlieren-based shock-detection algorithm at the University of Texas at San Antonio. The points of flow separation and reattachment were spatially and spectrally evaluated at 300 kHz for time-resolved characterization of the unsteadiness of this flow field. The separation bubble was found to rapidly expand and contract about a large mean length. Spectral analysis of the velocity fluctuations of this motion revealed this behavior to be a turbulent phenomenon, with spectral energy trends consistent with the Kolmogorov theory for the turbulent energy cascade.

Glasby, Ryan [ORNL]↗
DOE OSTI2022

Hybrid interferometric diagnostic for the refractive index measurement across nearly discrete shock waves

The measurement of high temperature gas properties is key for characterizing high speed flows. Nearly discrete changes in density across shock waves, in particular, are difficult to resolve through traditional fringe-counting interferometric methods. Existing techniques for estimating large fringe jumps are either resolution limited or require specialized window configurations. In this Letter, we describe a unique hybrid interferometric technique that combines narrowband fringes for high resolution and broadband fringes as an absolute reference to measure changes in refractive index with a resolution of up to 7 × 10 −8 across nearly discrete index changes of up to 1.5 × 10 −4 . By capturing fringes with an ultrahigh-speed camera, the refractive index changes across discrete shock fronts can be estimated inside a shock tube with high accuracy and time resolution. First, a novel hybrid calibration method for tracking finite fringes is discussed. Next, this technique is used to measure the post-initial-shock refractive indices for Mach 2.7 to 4.2 flows (pressures from 90.4 to 228.4 kPa). Results are then compared with theoretical values showing agreement within 2%.

Wang, Gwendolyn T. (ORCID:000000016271893X)↗
DOE OSTI2021

Observations on the fourth-power scaling of high-pressure shock waves in solids

This paper is concerned with the bounds and approximations that are found in the original derivation of the fourth-power scaling of shock waves in solids [D. Grady, Appl. Phys. Lett. 38, 825 (1981)]. The analysis is focused on the framework of the derivation and is independent of constitutive assumption, such as visco-plastic behavior. Results include an upper bound for the shock pressure and a restriction on the range of power coefficient for materials having a power-type function for the shock velocity–particle velocity relation. Relaxation of this restriction is proposed based on the idea that the rise time along the Rayleigh line is strain dependent. The idea led to the application of the Onsager relation to strain rate calculation that in turn resulted in a quadratic function for the strain rate–shock pressure relation. The fourth-power relation is obtained by generalizing the Onsager relation through the introduction of a quadratic dissipative potential in analogy to Rayleigh's dissipative potential.

Horie, Yasuyuki (ORCID:0000000245207246)↗
DOE OSTI2022

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shock-Wave Proprieties of High-Strength Ceramics

The shock and release response of a suite of high-strength ceramics has been investigated with controlled planar impact methods and velocity interferometry diagnostics. The study provides Hugoniot, yield and post-yield compressive strength, phase transformation and tensile spall material properties data for these ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spontaneous shock waves in pulse-stimulated flocks of Quincke rollers

Active matter demonstrates complex spatiotemporal self-organization not accessible at equilibrium and the emergence of collective behavior. Fluids comprised of microscopic Quincke rollers represent a popular realization of synthetic active matter. Temporal activity modulations, realized by modulated external electric fields, represent an effective tool to expand the variety of accessible dynamic states in active ensembles. Here, we report on the emergence of shockwave patterns composed of coherently moving particles energized by a pulsed electric field. The shockwaves emerge spontaneously and move faster than the average particle speed. Combining experiments, theory, and simulations, we demonstrate that the shockwaves originate from intermittent spontaneous vortex cores due to a vortex meandering instability. They occur when the rollers’ translational and rotational decoherence times, regulated by the electric pulse durations, become comparable. The phenomenon does not rely on the presence of confinement, and multiple shock waves continuously arise and vanish in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗