Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock structure”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Laboratory Resolved Structure of Supercritical Perpendicular Shocks

Supermagnetosonic perpendicular flows are magnetically driven by a large radius theta-pinch experiment. Fine spatial resolution and macroscopic coverage allow the full structure of the plasma-piston coupling to be resolved in laboratory experiment for the first time. Here, a moving ambipolar potential is observed to reflect unmagnetized ions to twice the piston speed. Magnetized electrons balance the radial potential via Hall currents and generate signature quadrupolar magnetic fields. Electron heating in the reflected ion foot is adiabatic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of Topology on Transient Dynamic and Shock Response of Polymeric Lattice Structures

Architected cellular materials, such as lattice structures, offer potential for tunable mechanical properties for dynamic applications of energy absorption and impact mitigation. In this work, the static and dynamic behavior of polymeric lattice structures was investigated through experiments on octet-truss, Kelvin, and cubic topologies with relative densities around 8%. Here, dynamic testing was conducted via direct impact experiments (25–70 m/s) with high-speed imaging coupled with digital image correlation and a polycarbonate Hopkinson pressure bar. Mechanical properties such as elastic wave speed, deformation modes, failure properties, particle velocities, and stress histories were extracted from experimental results. At low impact velocities, a transient dynamic response was observed which was composed of a compaction front initiating at the impact surface and additional deformation bands whose characteristics matched low strain-rate behavior. For higher impact velocities, shock analysis was carried out using compaction wave velocity and Eulerian Rankine–Hugoniot jump conditions with parameters determined from full-field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Saturation of Vortex Rings Ejected from Shock-Accelerated Interfaces

Structures evoking vortex rings can be discerned in shock-accelerated flows ranging from astrophysics to inertial confinement fusion. By constructing an analogy between vortex rings produced in conventional propulsion systems and rings generated by a shock impinging upon a high-aspect-ratio protrusion along a material interface, we extend classical, constant-density vortex-ring theory to compressible multifluid flows. Here, we further demonstrate saturation of such vortex rings as the protrusion aspect ratio is increased, thus explaining morphological differences observed in practice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Influence of Pore Surface Structure and Contents on Shock-Induced Collapse and Energy Localization

The majority of computational efforts in the shockwave driven collapse of porosity focus on idealized pore geometries and structures. These give significant insight into the physical mechanisms in play, but omit some potentially important effects of real energetic materials. Here, to address this, we simulate pore collapse in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for four types of pore surfaces and fillings: an idealized cylindrical pore, a pore surface coated with polymer, a melted/amorphous pore surface, and a pore filled with gas molecules. We find that these perturbations to the pore structure result in only minor changes to the highly complex collapse mechanisms. However, significant differences in hotspot temperature are observed, with peak temperatures ranging over several hundred Kelvin. The polymer-coated system has the most effect, lowering the hotspot temperature by 300 K, and also limiting the manifestation of intramolecular strain energy that is known to induce mechanochemistry.

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

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗
DOE OSTI2020

Shock-induced kinetics and cellular structures of liquid nitromethane detonation

Using a combination of high-speed diagnostics, optical pyrometry, velocimetry and video photography, we examine the spatiotemporal reaction kinetics of a prototypical high-energy explosive, liquid nitromethane (NM), as we drive this high explosive into detonation. The detonations were initiated by powerful shock waves whose durations (4 ns) were shorter than the characteristic time associated with the reaction (7 ns). Simple optical spectroscopy alone cannot characterize the kinetics with high time resolution because the reactants and products are flowing at a high velocity of ~6 km/s (6 µm/ns). Additionally, every detection volume behind the shock front contains material initiated over a range of times. High spatial resolution can be obtained by probing interfaces where the shock enters or breaks out of the NM. In addition, we obtain cellular patterns imprinted on the luminous shock front by the two-stage explosion in NM. These spatiotemporal patterns arise naturally as a result of the asymmetry produced by the moving shock front. In this way the shock front serves as a thin moving intrinsic optical gauge that reports and characterizes the two-stage NM explosion behind the front.

42 ENGINEERING↗
DOE OSTI2021

High-Speed X-Ray Stereo Digital Image Correlation in a Shock Tube

X-ray stereo digital image correlation (DIC) measurements were performed at 10 kHz on the internal surface of a jointed structure in a shock tube at a shock Mach number of 1.42 and compared with optical stereo DIC measurements on the outer, visible surface of the structure. The shock tube environment introduces temperature and density gradients in the gas through which the structure was imaged, resulting in spatial and temporal index of refraction variations. These variations cause bias errors in optical DIC measurements due to beam-steering but have minimal influence on x-ray DIC measurements. These results demonstrate the utility of time-resolved x-ray DIC measurements in complicated environments where optical measurements suffer severe errors and/or are precluded by lack of optical access.

42 ENGINEERING↗
DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in Its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. Here, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of $R={744.0}_{-128.0}^{+183.0}\,{R}_{\odot }$ and an ejected mass of ${M}_{\mathrm{ej}}={9.3}_{-5.8}^{+0.4}\,{M}_{\odot }$. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Phase transitions in Zr at sub-nanosecond time scales

Here solid-solid phase transitions are investigated in Zr thin films using shock compression induced by a short laser pulse (<1 ns). Shock wave profiles are measured at free surfaces for films of different thicknesses (a few micrometers) using chirped-pulse line velocimetry with 10 ps time resolution. Experiments are performed at pressures up to ~50 GPa, which is sufficient to reach the ω and β phases under equilibrium conditions. The shock wave structures are analyzed using a general Lagrangian analysis method, which allows for the calculation of stress–strain paths and assessments of phase transition behavior. In agreement with recent short laser pulse experiments using ultra-fast x-ray diffraction, we do not find any clear evidence of the α–ω transition, though this would be expected from the phase diagram. Instead, we infer a direct transformation to a metastable β-phase at lower shock pressures (<20 GPa) and equilibrium β at higher pressures. Through the velocimetry analysis, we find α–β transformation onset times of less than ~100 ps and completion times of less than ~200 ps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Crossed beam energy transfer in the presence of laser speckle ponderomotive self-focusing and nonlinear sound waves

Crossed beam energy transfer, CBET, is investigated by taking into account the speckle structure of smoothed laser beams that overlap in a plasma with an inhomogeneous flow profile. Using the two-dimensional simulation code Harmony, it is shown how ponderomotive self-focusing of laser speckles in crossing beams can significantly affect the transfer of energy from one beam to the other. The role of plasma flow in speckle self-focusing is investigated and revisited, in particular its consequences in terms of redirection and increasing angular spread of the laser beams due to beam bending and plasma-induced smoothing, respectively. In close-to-sonic flow, the onset of self-focusing in the beam speckle structure occurs at considerably lower beam intensities than expected for the case without flow. CBET and speckle self-focusing can, hence, occur together when two crossed beams with equal frequency resonantly exchange energy via their ponderomotively driven density perturbations flowing with sound speed. From the simulations, it is found that consequences of ponderomotive self-focusing can be expected above an average intensity threshold scaling as $I_L~ 2 × 10^{14} W cm^{-2}(λ_0/1 μm)^{-2}(T_e/ keV)$, with an impact on the spatial and temporal coherence of the transmitted light. The density perturbations due to the ponderomotive force of the crossing beams can locally be enhanced in self-focusing speckles, partly leading to shock-like structures. These structures eventually increase the effect of plasma-induced smoothing and are at the origin of the stronger angular spread.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of a Fine-Scale Layered Structure of the Atmosphere on Infrasound Signals from Fragmenting Meteoroids

We investigate the influence of a fine-scale (FS) layered structure in the atmosphere on the propagation of infrasound signals generated by fragmenting meteoroids. Using a pseudo-differential parabolic equation (PPE) approach, we model broadband acoustic signals from point sources at altitudes of 35–100 km. The presence of FS fluctuations in the stratosphere (37–45 km) and the lower thermosphere (100–120 km) modifies ray trajectories, causing multiple arrivals and prolonged signal durations at ground stations. In particular, meteoroids fragmenting at 80–100 km can produce two distinct thermospheric arrivals beyond 150 km range, while meteoroids descending to 50 km or below yield weak, long-lived arrivals within the acoustic shadow zone via antiguiding propagation and diffraction. Comparison with observed infrasound data confirms that FS-layered inhomogeneities can account for multi-arrival “N-waves,” broadening potential interpretations of meteoroid signals. The results also apply to other atmospheric-entry objects, such as sample return capsules, emphasizing how FS structure impacts shock wave propagation. In conclusion, our findings advance understanding of wavefield evolution in a layered atmosphere and have broad relevance for global infrasound monitoring of diverse phenomena (e.g., re-entry capsules, rocket launches, and large-scale explosions).

Aeroacoustics↗
DOE OSTI2022

Magnetic Field Amplification by a Plasma Cavitation Instability in Relativistic Shock Precursors

Abstract Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗