Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock hydrodynamics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Shocks in Dense Clouds in the Vela Supernova Remnant: FUSE

We have obtained 8 LWRS FUSE spectra to study a recently identified interaction of the Vela supernova remnant with a dense cloud region along its western edge. The goal is to quantify the temperature, ionization, density, and abundance characteristics associated with this shock/dense cloud interface by means of UV absorption line studies. Our detection of high-velocity absorption line C I at +90 to +130 km/s with IUE toward a narrow region interior to the Vela SNR strongly suggests the Vela supernova remnant is interacting with a dense ISM or molecular cloud. The shock/dense cloud interface is suggested by (1) the rarity of detection of high-velocity C I seen in IUE spectra, (2) its very limited spatial distribution in the remnant, and (3) a marked decrease in X-ray emission in the region immediately west of the position of these stars where one also finds a 100 micron emission ridge in IRAS images. We have investigated the shock physics and general properties of this interaction region through a focussed UV absorption line study using FUSE spectra. We have FUSE data on OVI absorption lines observed toward 8 stars behind the Vela supernova remnant (SNR). We compare the OVI observations with IUE observations of CIV absorption toward the same stars. Most of the stars, which are all B stars, have complex continua making the extraction of absorption lines difficult. Three of the stars, HD 72088, HD 72089 and HD 72350, however, are rapid rotators (v sin i less than 100 km/s) making the derivation of absorption column densities much easier. We have measured OVI and CIV column densities for the "main component" (i.e. the low velocity component) for these stars. In addition, by removing the H2 line at 1032.35A (121.6 km/s relative to OVI), we find high velocity components of OVI at approximately 150 km/s that we attribute to the shock in the Vela SNR. The column density ratios and magnitudes are compared to both steady shock models and results of hydrodynamical SNR modeling. We find that the models require the shock to be relatively slow (approximately 100 - 170 km/s) to match the FUSE data. We discuss the implications of our results for models of the evolution of the Vela SNR.

Nichols, Joy↗
NASA NTRS1994

Postexplosion hydrodynamics of supernovae in red supergiants

Shock propagation, mixing, and clumping are studied in the explosion of red supergiants as Type II supernovae using a two-dimensional smooth particle hydrodynamic (SPH) code. We show that extensive Rayleigh-Talor instabilities develop in the ejecta in the wake of the reverse shock wave. In all cases, the shell structure of the progenitor is obliterated to leave a clumpy, well-mixed supernova remnant. However, the occurrence of mass loss during the lifetime of the progenitor can significantly reduce the amount of mixing. These results are independent of the Type II supernova explosion mechanism.

Herant, Marc↗
DOE OSTI2023

Observation of subcritical shocks in a collisional laboratory plasma: scale dependence near the resistive length

We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.

Physics↗
NASA NTRS2004

Underexpanded Screeching Jets From Circular, Rectangular, and Elliptic Nozzles

The screech frequency and amplitude, the shock spacing, the hydrodynamic-acoustic standing wave spacing, and the convective velocity of large organized structures are measured in the nominal Mach number range of 1.1 less than or = Mj less that or = l0.9 for supersonic, underexpanded jets exhausting from a circular, a rectangular and an elliptic nozzle. This provides a carefully measured data set useful in comparing the importance of various physical parameters in the screech generation process. The hydrodynamic-acoustic standing wave is formed between the potential pressure field of large turbulent structures and the acoustic pressure field of the screech sound. It has been demonstrated earlier that in the currently available screech frequency prediction models replacement of the shock spacing by the standing wave spacing provides an exact expression. In view of this newly found evidence, a comparison is made between the average standing wavelength and the average shock spacing. It is found that there exists a small, yet important, difference, which is dependent on the azimuthal screech mode. For example, in the flapping modes of circular, rectangular, and elliptic jets, the standing wavelength is slightly longer than the shock spacing, while for the helical screech mode in a circular jet the opposite is true. This difference accounts for the departure of the existing models from predicting the exact screech frequency. Another important parameter, necessary in screech prediction, is the convective velocity of the large organized structures. It is demonstrated that the presence of the hydrodynamic-acoustic standing wave, even inside the jet shear layer, becomes a significant source of error in the convective velocity data obtained using the conventional methods. However, a new relationship, using the standing wavelength and screech frequency is shown to provide more accurate results.

Panda, J.↗
NASA NTRS1986

Cosmic-ray shock acceleration in oblique MHD shocks

A one-dimensional, steady-state hydrodynamical model of cosmic-ray acceleration at oblique MHD shocks is presented. Upstream of the shock the incoming thermal plasma is subject to the adverse pressure gradient of the accelerated particles, the J x B force, as well as the thermal gas pressure gradient. The efficiency of the acceleration of cosmic-rays at the shock as a function of the upstream magnetic field obliquity and upstream plasma beta is investigated. Astrophysical applications of the results are briefly discussed.

Webb, G. M.↗
DOE OSTI2025

Observation of ponderomotively driven bow shock using Thomson scattering

High-energy speckled lasers are known to exert ponderomotive forces on a plasma. This can reduce flow transverse to the propagation of the beams. When coupled with a supersonic flow, this has been shown to lead to the formation of a shock that travels against the flow. Experiments conducted on the OMEGA laser facility have used Thomson scattering to observe density and velocity changes consistent with this ponderomotively driven shock. Comparisons of the data with hydrodynamic simulations with the ponderomotive force, particle-in-cell simulations with a full Maxwell field solver, and hydrodynamic simulations without the ponderomotive force show that this shock feature is only reproduced when accounting for the ponderomotive force.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock formation in flowing plasmas by temporally and spatially smoothed laser beams

The cumulative impact of multiple laser speckles on a supersonic plasma flow across optically smoothed laser beams is investigated. The bending of laser beams caused by ponderomotive laser–plasma coupling, together with flow, leads to plasma a momentum-conserving response that results in a deceleration of the flow. Once the flow velocity decreases to a subsonic level, the action of the laser beams can generate a shock within the plasma. This scenario has been predicted theoretically and confirmed by hydrodynamic simulations. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size, and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging only in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser–plasma coupling in the inertial confinement fusion experiments are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Suppression of Richtmyer-Meshkov Instability via Special Pairs of Shocks and Phase Transitions

The classical Richtmyer-Meshkov instability (RMI) is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. Here, we demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances, a phase-transitioning material may self-suppress RMI.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗
DOE OSTI2020

Hotspot formation due to shock-induced pore collapse in 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX): Role of pore shape and shock strength in collapse mechanism and temperature

The shock to detonation transition in heterogeneous high energy density solids starts with the spatial localization of mechanical energy into so-called hotspots that form due to the interaction between the leading wave and microstructural features and defects. We used large-scale molecular dynamics to characterize the hotspots resulting from the shock-induced collapse of cylindrical voids and elongated cracks focusing on the effect of shock strength, defect shape, and size. The temperature fields resulting from the collapse of cracks elongated along the shock direction show significantly higher sensitivity to both shock strength and size than cylindrical voids. Cracks 80 nm in length result in temperatures almost three times higher than voids 80 nm in diameter, reaching values corresponding to the ideal case of isentropic recompression of a gas. The molecular dynamics trajectories reveal the atomic origin of this contrasting behavior. While circular voids undergo a transition from viscoelastic pore collapse to a hydrodynamic regime with increasing shock strength, shock focusing in elongated cracks results in jetting and vaporization which, upon recompression, leads to increased heating.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Model development of supersonic trough wind with shocks

The time dependent one dimensional hydrodynamic equations describe the evolution of the thermal plasma flow along closed magnetic field lines outside of the plasmasphere. The convection of the supersonic polar wind onto a closed fieldline results in the assumed formation of collisionless plasma shocks. These shocks move earthward as the field line with its frozen-in plasma remains fixed or contracts with time to smaller L coordinates. The high equatorial plasma temperature (of the order of electron volts) produced by the shock process decreases with time if the flow is isothermal but it will increase if the contraction is under adiabatic conditions. Assuming adiabaticity a peak in the temperature forms at the equator in conjunction with a depression in the ion density. After an initial contraction, if the flux tube drifts to higher L coordinates the direction of the shock motion can be reversed so that the supersonic region will expand along the field line towards the state characterizing the supersonic polar wind. A rapid expansion will lower the equatorial density while the temperature decreases with time under adiabatic but not isothermal conditions.

Grebowsky, J. M.↗
DOE OSTI2023

Structures of strong shocks in low-density helium and neon gases

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ~350 µm with the helium gas mixture but less than 30 µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Additionally, both effects explain the longer precursor layer in the helium shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Measurements of the inert Hugoniot and observation of mechanical ignition in Ni(V) + Al reactive multilayers via x-ray diffraction

This study investigates the inert Hugoniot response, mechanical ignition, and reaction dynamics of Ni(V)+Al multilayers during longitudinal, laser-driven shock compression experiments. Ni(V)+Al multilayers, known for their self-propagating exothermic reactions, were subjected to longitudinal stresses exceeding 50 GPa using the laser shock facility within the Dynamic Compression Sector (DCS) at the Advanced Photon Source (APS). In situ x-ray diffraction (XRD) revealed that Ni(V) and Al were not in equilibrium during compression, with stress discrepancies attributed to twinning, grain structure effects, and/or dislocation density. However, the measured inert Hugoniot closely matched prior experimental and computational studies, confirming the utility of XRD for measuring the equation of state of thin, complex materials. Additionally, reaction was observed at significantly higher stresses than reported previously using laser-launched flyers. This discrepancy suggests a strong influence of externally imposed shear stress on reaction thresholds, which likely arose from deviations in flyer planarity during past experiments. Full reaction of the multilayer occurred within 40 ns after shock-wave passage, evidenced by complete melting of the constituents. Eulerian hydrocode simulations replicated experimental conditions, providing insights into equilibrium dynamics and experimental artifacts. The results highlight how even small shear forces facilitate ignition in Ni(V)+Al multilayers at lower stresses.

Crystallography↗
DOE OSTI2022

Design optimization for Richtmyer–Meshkov instability suppression at shock-compressed material interfaces

We report the Richtmyer–Meshkov instability (RMI) is a phenomenon that occurs at the interface of two substances of different densities due to an impulsive acceleration, such as a shock wave passing through this interface. Under these conditions, the instability can be seen as interface perturbations begin to grow into narrow jets or spikes of one substance that propagate into the other. In some cases, this interface may involve an elastic–plastic material, which can play a significant role in the development and behavior of the RMI. The ability to effectively control RMI jetting and spike growth is one major limiting factor in technological challenges, such as inertial confinement fusion, that involve using high-pressure shock waves to implode a fuel target. The propagation of RMI growth can lead to increased asymmetry in this implosion process and significantly reduce the obtained energy yield. We use hydrodynamics simulations of impactor shock-compression experiments and methods based in design optimization to suppress RMI spike growth by altering the geometry and other properties of a shock-compressed elastic–plastic material target that shares an interface with atmospheric air. These hydrodynamics simulations use an arbitrary Lagrangian–Eulerian method with a high-order finite element approach. Our results demonstrate that RMI suppression can be achieved by intentionally creating a separate upstream interface instability to counteract the growth of long narrow RMI spikes at an interface with initial perturbations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Physics consistent machine learning framework for inverse modeling with applications to ICF capsule implosions

In high energy density physics (HEDP) and inertial confinement fusion (ICF), predictive modeling is complicated by uncertainty in parameters that characterize various aspects of the modeled system, such as those characterizing material properties, equation of state (EOS), opacities, and initial conditions. Typically, however, these parameters are not directly observable. What is observed instead is a time sequence of radiographic projections using X-rays. In this work, we define a set of sparse hydrodynamic features derived from the outgoing shock profile and outer material edge, which can be obtained from radiographic measurements, to directly infer such parameters. Our machine learning (ML)-based methodology involves a pipeline of two architectures, a radiograph-to-features network (R2FNet) and a features-to-parameters network (F2PNet), that are trained independently and later combined to approximate a posterior distribution for the parameters from radiographs. We show that the machine learning architectures are able to accurately infer initial conditions and EOS parameters, and that the estimated parameters can be used in a hydrodynamics code to obtain density fields, shocks, and material interfaces that satisfy thermodynamic and hydrodynamic consistency. Finally, we demonstrate that features resulting from an unknown EOS model can be successfully mapped onto parameters of a chosen analytical EOS model, implying that network predictions are learning physics, with a degree of invariance to the underlying choice of EOS model. To the best of our knowledge, our framework is the first demonstration of recovering both thermodynamic and hydrodynamic consistent density fields from noisy radiographs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗