Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗
DOE OSTI2026

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Measurement of the sound velocity and Grüneisen parameter of polystyrene at inertial confinement fusion conditions

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of ICF experiments. The Grüneisen parameter increases with density below ~3.1 g/cm 3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. Finally, the discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Elastic moduli of hexagonal diamond and cubic diamond formed under shock compression

The relative stiffnesses and strengths of hexagonal diamond and cubic diamond constitute long standing unresolved issues, because relevant experimental results are lacking. Laser interferometry was used to experimentally determine the longitudinal sound speeds and moduli in hexagonal diamond and cubic diamond formed during the shock compression of graphite. The hexagonal diamond longitudinal moduli are significantly larger than the cubic diamond longitudinal moduli, and even exceed averaged cubic diamond single crystal values. Here, the measured hexagonal diamond longitudinal moduli, combined with high-pressure bulk moduli for cubic diamond single crystals, show that shock-formed hexagonal diamond shear moduli are larger than the shear moduli for cubic diamond single crystals.

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 OSTI2021

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗
DOE OSTI2021

Electrical conductivity of Sn at high pressure and temperature

To date, temperature and conductivity have many outstanding implications in extreme environments but are yet to be fully understood under high pressure and temperature dynamic conditions. Here, we introduce an approach to provide high quality electrical conductivity results under dynamic loading conditions. Emphasis is given to address the skin depth effect's influence in a dynamic loading experiment by using thin films. The thin film samples in this study were at least 100 times thinner than previous samples in dynamic electrical conductivity experiments, increasing the current density to its full potential across the sample's entire cross section. Consideration of the skin depth accounts for at minimum a 4× scaling factor to the final electrical conductivity result that has been neglected in previous dynamic electrical conductivity studies. We also obtained improved signal-to-noise ratio with custom diagnostics optimized for better electrical impedance matching. Furthermore, these considerations were applied to Sn to assess electrical conductivity at elevated pressure and temperature. The high signal-to-noise ratio with reduced skin depth influence results in Sn allowing observation of the conductivity changes related to solid-to-solid and solid-to-liquid phase transitions. Additionally, we calculate the Sn thermal conductivity using the Wiedemann-Franz law for our experiments and compare against thermal transport dependent temperature measurements from previous work.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sound speed measurements in silver shock compressed to 300 GPa: Solid-state transition, melting, and liquid-state response

To gain insight into the thermodynamic response of shock compressed Ag states corresponding to the face-centered-cubic (fcc) to body-centered-cubic (bcc) transformation and melting observed using in-situ x-ray diffraction (XRD) measurements, we present longitudinal sound speed results and their analysis in Ag at peak stresses ranging from 60 to 300 GPa. The measured sound speeds increased linearly with density compression to 171 GPa, showing that the sound speeds (and longitudinal moduli) in the fcc and bcc phases are very similar. Between 171 and 218 GPa, the sound speed dropped significantly, consistent with the melting reported using XRD measurements in shock compressed Ag. From 218 to 300 GPa, the increasing sound speeds and Hugoniot states provide a determination of the liquid phase Ag response. In particular, determination of the Grüneisen parameter (Γ) showed that the density – Grüneisen parameter product (ρΓ) for liquid Ag is constant, but differs significantly from that for solid Ag at ambient conditions. Furthermore, the Mie- Grüneisen equation of state can be used to describe the Hugoniot and off-Hugoniot response of liquid Ag.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

X-ray diffraction study of phase transformation dynamics of Fe and Fe-Si alloys along the shock Hugoniot using an x-ray free electron laser

The x-ray free electron laser (XFEL) enables probing a highly compressed material response at the subnanosecond timescale. Here we exploit the ultrafast XFEL pulse to combine reflection x-ray diffraction and laser-driven shock compression to perform a study of the phase transformation and stability in Fe and Fe-Si alloys. Our approach enables us to observe that solid-solid phase transformations occur in Fe and Fe- Si 8.5 wt % in ≤ 130 ps at ~ 130 GPa; no transformation is observed in Fe- Si 16 wt % up to 110 GPa. Density functional theory calculations predict similar phase relations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

X-ray diffraction measurements and pressure determination in nanosecond compression of solids up to 600 GPa

X-ray diffraction measurements under laser-driven dynamic compression now allow us to investigate the atomic structure of matter at TPa pressures and thousands of degree temperatures with broad implications for condensed matter physics, planetary science and astronomy. Furthermore, pressure determination in these experiments often relies on velocimetry measurements coupled with modeling that requires accurate knowledge of the optical and thermo-mechanical properties of a window material, resulting in significant systematic uncertainty. Here we demonstrate different approach applicable to X-ray diffraction experiments under quasi-isentropic ramp-compression based on the use of in-situ pressure calibrants, similar to the methods often adopted in static-compression experiments with diamond anvil cells. Focusing on experiments using a diamond window, we discuss challenges and mitigation strategies for the novel approach. Our study, in addition to providing new structural information of 5 metals up to hundreds of GPa, provides validation to the currently used methods based on time-resolved measurement of the diamond free-surface velocity, and reveals that the use of in-situ calibrants enables a factor of four reduction in the pressure uncertainty in these experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 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

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗