Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock wave experiments”

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

Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501

We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. Overall, we conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.

42 ENGINEERING↗
DOE OSTI2022

Effect of laser preheat in magnetized liner inertial fusion at OMEGA

Integrated magnetized liner inertial fusion (MagLIF) experiments on the OMEGA laser system are modeled to investigate the effects of laser preheat on implosion performance. In both simulations and experiments, preheating the fuel with a laser enhances the neutron yield, with the maximum yield occurring at an optimal preheat laser energy. Increasing the preheat energy past the optimal value reduces the neutron yield. In simulations, the rate of the drop in neutron yield depends on the inclusion of the Nernst effect. The Nernst effect advects the magnetic field away from the center of the fuel region during the preheat stage of MagLIF on OMEGA and weakens magnetic-flux compression. If the Nernst effect is not included, the simulated drop in yield past the optimal preheat laser energy is more gradual as opposed to the sharper drop seen in experiments. Simulations are able to model trend in measured ion temperatures seen in experiments. Here, a mix model is used to show that including mix from the wall in the simulated fuel region further degrades the yield and lowers the optimal preheat laser energy. Simulations with mix predict that increasing the initial axial magnetic field may still enhance yield performance in integrated implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition

The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Calibration and characterization of the line-VISAR diagnostic at the HED-HIBEF instrument at the European XFEL

In dynamic-compression experiments, the line-imaging Velocity Interferometer System for Any Reflector (VISAR) is a well-established diagnostic used to probe the velocity history, including wave profiles derived from dynamically compressed interfaces and wavefronts, depending on material optical properties. Knowledge of the velocity history allows for the determination of the pressure achieved during compression. Such a VISAR analysis is often based on Fourier transform techniques and assumes that the recorded interferograms are free from image distortions. In this paper, we describe the VISAR diagnostic installed at the HED-HIBEF instrument located at the European XFEL along with its calibration and characterization. It comprises a two-color (532, 1064 nm), three-arm (with three velocity sensitivities) line imaging system. We provide a procedure to correct VISAR images for geometric distortions and evaluate the performance of the system using Fourier analysis. We finally discuss the spatial and temporal calibrations of the diagnostic. As an example, we compare the pressure extracted from the VISAR analysis of shock-compressed polyimide and silicon.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Experimental characterization of hot-electron emission and shock dynamics in the context of the shock ignition approach to inertial confinement fusion

We report on planar target experiments conducted on the OMEGA-EP laser facility performed in the context of the shock ignition (SI) approach to inertial confinement fusion. The experiment aimed at characterizing the propagation of strong shock in matter and the generation of hot electrons (HEs), with laser parameters relevant to SI (1-ns UV laser beams with I ~10 16 W/cm 2 ). Time-resolved radiographs of the propagating shock front were performed in order to study the hydrodynamic evolution. The hot-electron source was characterized in terms of Maxwellian temperature, T h , and laser to hot-electron energy conversion efficiency η using data from different x-ray spectrometers. The post-processing of these data gives a range of the possible values for T h and η [i.e., T h [keV] ϵ (20, 50) and η ϵ (2%, 13%)]. These values are used as input in hydrodynamic simulations to reproduce the results obtained in radiographs, thus constraining the range for the HE measurements. According to this procedure, we found that the laser converts ~10% ± 4% of energy into hot electrons with T h = 27 ± 8 keV. The paper shows how the coupling of different diagnostics and numerical tools is required to sufficiently constrain the problem, solving the large ambiguity coming from the post-processing of spectrometers data. The effect of the hot electrons on the shock dynamics is then discussed, showing an increase in the pressure around the shock front. Furthermore, the low temperature found in this experiment without pre-compression laser pulses could be advantageous for the SI scheme, but the high conversion efficiency may lead to an increase in the shell adiabat, with detrimental effects on the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Argon equation of state data to 1 TPa: Shock compression experiments and simulations

Argon is the most abundant noble gas on Earth and its noble, atomic fluid nature makes it an excellent candidate for comparison of experiment and theory at extreme conditions. For this study, we performed a combined computational and experimental study on shock compressed cryogenic liquid argon. Using Sandia's Z machine, we shock compressed liquid argon to 600 GPa and reshock states up to 950 GPa. Laser shock experiments at the Omega Laser facility extend the principal Hugoniot to 1000 GPa and provided temperature data along the principal Hugoniot. The plate impact experiments and laser shock experiments used well-characterized impedance matching standards and demonstrate consistent results between the two platforms over a common range. Density functional theory based molecular dynamics simulations provided additional data on the Hugoniot to 600 GPa. The combined experimental data and simulation results provide constraints on the development of new equation of state models at extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct observations of shock-induced melting in a porous solid using time-resolved x-ray diffraction

Here in this work we provide direct evidence of shock-induced melting and associated kinetics in a porous solid (aluminum powder) using time-resolved x-ray diffraction. Unambiguous evidence of melting in 50% porous aluminum (Al) powder samples, shocked to peak pressures between ~13 - 19 GPa, was provided by the broadening of the Debye-Scherrer ring corresponding to the (111) peak. Shocked Al powder did not melt completely in any of our experiments within the durations of measurement. Incomplete (partial) melting of the powder, even after several hundreds of nanoseconds of shock loading, provides insights into thermal transport with Al powder particles under high-pressure dynamic loading. Such insights are quite valuable for developing well-constrained melting models and thermodynamic equations of state for porous Al and other porous solids relevant to planetary and materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Observations of twinning microstructure in iron ramp-compressed through the α – ε phase transition

Single crystal iron was ramp-compressed along the [100] axis at 4 GPa/ns through the α–ϵ phase transformation at 13 GPa. Unlike shock experiments, in these ramp compression experiments the bcc lattice is observed to be nearly isotropically relaxed before the phase transition. There was a mixed α/ϵ phase region starting at ~13 GPa which was observed for 0.75 ns until the peak stress of 18 GPa was reached. In situ x-ray diffraction measurements show the formation of new hcp orientations not reported in shock or quasistatic compression experiments. The new hcp orientations appear to be caused by {10¯11} <10¯1 ¯2> → {¯1011} <¯101¯2> sequential twins that occur during the phase transition. This twinning mechanism relieves the shear strain caused by the bcc-hcp phase transition from the isotropically relaxed bcc phase while subject to uniaxial compression. Furthermore, these results show that in iron the induced microstructure through a phase transition and the phase transition mechanism depend on the loading history.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The Seismic Signature of a High-Energy Density Physics Laboratory and Its Potential for Measuring Time-Dependent Velocity Structure

The Z Machine at Sandia National Laboratories is a pulsed power facility for high-energy density physics experiments that can shock materials to extreme temperatures and pressures through a focused energy release of up to ~ 25 MJ in < 100 nanoseconds. It has been in operation for more than two decades and conducts up to ~ 100 experiments, or “shots,” per year. Based on a set of 74 known shot times from 2018, we determined that Z Machine shots produce detectable ~ 3–17 Hz ground motion 12 km away at the Albuquerque Seismological Laboratory, New Mexico (ANMO), borehole seismograph, with peak signal at ~ 7 Hz. The known shot waveforms were used to create a three-component template, leading to the detection of 2339 Z Machine shots since 1998 through single-station cross-correlation. Local seismic magnitude estimates range from local magnitude (M L ) —2 to —1.3 and indicate that only a small fraction of the shot energy is transmitted by seismic phases observable at 12 km distance. The most recent major facility renovation, which was intended to decrease mechanical dissipation, is associated with an abrupt decrease in observed seismic amplitudes at ANMO despite stable maximum shot energy. The highly repetitive impulsive sources are well suited to coda-wave interferometry to investigate time-dependent velocity structures. Relative velocity variations (dv/v) show an annual cycle with amplitude of ~ 0.2%. Local minima are observed in the late spring, and dv/v increases through the summer monsoon rainfall, possibly reflecting patchy saturation as rainfall infiltrates near the eastern edge of the Albuquerque basin. Importantly, the cumulative results demonstrate that forensic seismology can provide insight into long-term operation of facilities such as pulsed-power laboratories, and that their recurring signals may be valuable for studies of time-dependent structure

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. As a result, we measure little heating across the shock (<10 % of the ion kinetic energy) which is consistent with an absence of viscous dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Detonation-induced transformation of graphite to hexagonal diamond

In this work, we explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in situ synchrotron x-ray diffraction in the ns timescale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under uniaxial compression, the energy barrier for the transition toward hexagonal diamond is lower than that for cubic diamond. Finally, no indication of cubic diamond formation was observed up to >70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Observation of shock-front separation in multi-ion-species collisional plasma shocks

We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the plasma liner experiment. Visible plasma emission near the He-I 587.6-nm and Ar-II 476.5–514.5-nm lines is simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈5×10 14 cm –3 , ion and electron temperatures of ≈2.2 eV, and relative plasma-merging speed of 22 km/s), the observationally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are <1 cm, which correspond to ~50 post-shock thermal ion–ion mean free paths. These experimental length scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the CHICAGO code. However, we note there are differences between the experimentally inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal (rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multiphase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(\sigma _{\rho /\langle \rho \rangle }^{2\Gamma }-1)^{1/2}(1-\sigma _{\rho /\langle \rho \rangle }^{-2})^{-1/2}\mathcal {M}^{-1}\Gamma ^{-1/2}$, with the density dispersion σρ/<ρ>, the turbulent Mach number $\mathcal {M}$, and the polytropic exponent Γ. Purely solenoidal and purely compressive driving correspond to b ~ 1/3 and b ~ 1, respectively. As a result, using simulations in which a shock is driven into a multiphase medium with structures of different sizes and Γ < 1, we find b ~ 1 for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

79 ASTRONOMY AND ASTROPHYSICS↗