Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 OSTI2026

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗
DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sound speed determination in copper shock compressed to 190 GPa

Sound speed measurements in shock compressed solids have long been valuable for the development of equations of state at extreme conditions, shock-induced phase transformations, and a comprehensive characterization of the thermophysical response of high-pressure standards. We present results from plate impact experiments to 190 GPa to determine the longitudinal sound speed in copper—an important high-pressure standard. Surprisingly, the sound speeds determined using the two most common experimental techniques—the front surface impact (FSI) approach and the release wave overtake (RWO) approach—diverge significantly for stresses greater than ~100 GPa. Further analyses, including numerical simulations, show that the FSI experiments provide the correct sound speeds and that fundamental assumptions underlying the RWO method are likely violated due to the complex release response of shock compressed copper. The sound speeds determined using the FSI approach provide for a more accurate high-pressure description of copper in dynamic compression experiments. The present findings are in contrast to the results for shock compressed silver [Wallace et al., Phys. Rev. B 104, 214106 (2021)], where both methods provided consistent sound speed results. Thus, the findings presented here demonstrate the need to experimentally verify the validity of the RWO method on a case-by-case basis. Finally, we note that even at the high stresses in the present work, the copper unloading response shows a time-dependent, quasielastic response often observed in metals at lower stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hugoniot properties of porous stainless steel: Insights from molecular dynamics simulations

Here, we present a systematic study of Hugoniot properties of porous 316L stainless steel using both a simple interpolation scheme and direct shock simulations in order to probe pore collapse kinetics as well as final thermodynamic states. Both methods indicate that equilibrated Hugoniot properties depend on pore density only and not on the pore distribution or size. We then create a simple porous equation of state model that is shown to be accurate for a range of validation data. This allows us to extend our simulations to make direct comparison to experimental data that have generally significantly larger system sizes and durations. In addition, our direct shock simulations indicate that the relaxation time after hotspot formation is system size dependent and can reach nanosecond timescales for the largest pores investigated in our study, thereby possibly having a measurable effect on fast dynamic loading experiments.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Effect of Grain Boundary Misorientation on Spall Strength in Ta via Shock-Free Simulations with Relatively Few Atoms

A suite of 37 molecular dynamics simulations is conducted at two system sizes to systematically characterize the role of grain boundary (GB) misorientation on spall strength in pure BCC tantalum (Ta). The systems studied consist of bicrystals with a single [110] symmetric tilt grain boundary. Two loading conditions are compared: (i) homogeneous extension under uniaxial strain simulated in this study and (ii) piston/flyer impact of sample, which induces heterogeneous deformation via shockwave propagation along the length of the sample. The piston/flyer impact is taken from the literature and run on the same set of GB misorientation angles using LAMMPS. The major finding here is that both methods result in similar spall strength predictions, but the homogeneous extension method generally requires two to three orders of magnitude fewer atoms and similar reductions in computational costs. Spall strength results systematically overpredict using this method, by about 10% for the dataset three orders of magnitude smaller than piston/flyer simulations, and 5% for the dataset two orders of magnitude smaller. Lastly, the effect of system size and pre-compression magnitude on spall strength is systematically characterized.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2022

Collapse of helium-filled voids in extreme deformation: Dislocation mechanisms

The mechanisms responsible for the collapse of helium-filled bubbles during the passage of shock waves in monocrystalline copper are revealed. Both internal pressure (caused by pre-existing helium atoms) and bubble size are varied in molecular dynamics simulations to understand the atomistic scale deformation as they are subjected to shock compression at pressures of 48, 123, and 170 GPa, corresponding to particle velocities of 1.0, 2.0, and 2.5 km/s. Both empty and helium filled bubbles serve as dislocation sources, generating intense, localized plastic regions. There are distinct differences in the collapse of empty voids compared to He-filled bubbles, the former requiring less stress and generating a greater density of dislocations for a given shock strength. Additionally, a generalized model for dislocation emission is proposed, where the inclusion of shear stress generated by the helium bubble increases the critical stress to generate dislocations at the defect surface, demonstrating the change in plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shock initiation and hot spots in plastic-bonded 1,3,5-triamino-2,4,6-trinitrobenzene (TATB)

TATB (1,3,5-triamino-2,4,6-trinitrobenzene) is a powerful explosive whose dynamical behavior is difficult to study because TATB is so insensitive to initiation by shock waves. We used a tabletop microscope equipped with 0–4.5 km/s laser-launched flyer plates to study shock initiation of TATB, which was fabricated in the form of an array of hundreds of plastic-bonded explosive minicharges (X-TATB = 80% TATB + 20% Sylgard 182 polymer). The 4 ns shocks from the flyer plates were not effective in initiating TATB, but we also developed a two-layer array where flyers first initiated a plastic-bonded PETN (pentaerythritol tetranitrate) charge (X-PETN = 80% PETN + 20% Sylgard), which drove an initiating 25 ns shock into the X-TATB. Thermal emission from shocked X-TATB was used to measure time-dependent temperature profiles with a resolution of 2 ns and to produce high-speed (5 ns) videos. In X-TATB, flyer plates produced 2500–3500 K hot spots and combustion at 2500 K. With X-PETN initiators, X-TATB had 3500–4000 K hot spots and a powerful volume explosion lasting a few nanoseconds. Prospects for producing TATB detonations on a tabletop are discussed.

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

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structural response of α-quartz under plate-impact shock compression

Because of its far-reaching applications in geophysics and materials science, quartz has been one of the most extensively examined materials under dynamic compression. Despite 50 years of active research, questions remain concerning the structure and transformation of SiO 2 under shock compression. Continuum gas-gun studies have established that under shock loading quartz transforms through an assumed mixed-phase region to a dense high-pressure phase. While it has often been assumed that this high-pressure phase corresponds to the stishovite structure observed in static experiments, there have been no crystal structure data confirming this. In this study, we use gas-gun shock compression coupled with in situ synchrotron x-ray diffraction to interrogate the crystal structure of shock-compressed α-quartz up to 65 GPa. Our results reveal that α-quartz undergoes a phase transformation to a disordered metastable phase as opposed to crystalline stishovite or an amorphous structure, challenging long-standing assumptions about the dynamic response of this fundamental material.

58 GEOSCIENCES↗
DOE OSTI2025

Ultrafast x-ray diffraction of high-pressure phases in dynamically compressed TiO 2

Here, we investigate the high-pressure polymorphism of Ti⁢O 2 under laser-shock compression from 54(5) to 137(7) GPa using in situ femtosecond x-ray diffraction. Our results provide experimental evidence of the 𝑃⁢𝑐⁢𝑎⁢2 1 -type distorted fluorite structure formed from polycrystalline Ti⁢O 2 dynamically compressed to 54(5) GPa. Upon higher compression, we observe the direct formation of the ninefold coordinated F⁡e 2 ⁢P-type phase at 68(4) and 78(3) GPa in polycrystalline and [001]-oriented Ti⁢O 2 , respectively. This represents an unprecedented 100 GPa reduction in the shock synthesis pressure of the F⁡e 2⁢ P-type structure relative to quasihydrostatic loading conditions. On pressure release, the F⁡e 2⁢ P-type phase transforms to the α-Pb⁢O 2 structure and, at later times, reverts to rutile. Thus, the rutile →F⁡e 2⁢ P and α-Pb⁢O 2 →rutile transformations are both observed to occur on nanosecond timescales. Our results highlight the unique ability of high-strain-rate uniaxial compression to synthesize novel high-pressure phases and also indicate the importance of in situ atomic-level probes in developing pressure-temperature phase diagrams.

Crystal structures↗
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 OSTI2020

Transcriptomic Profiles of Zymomonas mobilis 8b to Furfural Acute and Long-Term Stress in Both Glucose and Xylose Conditions

Zymomonas mobilis 8b is an ethanologenic bacterium engineered to utilize both glucose and xylose. The impacts of lignocellulosic hydrolyzate inhibitors on the growth of Zymomonas mobilis 8b have been investigated. However, the molecular responses of these inhibitors have not been completely elucidated yet. In this study, molecular responses to furfural were investigated using transcriptomic approaches of both chip-based microarray and a directional mRNA-Seq. Furfural acute shock time-course experiment with 3 g/L furfural supplemented when cells reached exponential phase and stress response experiment in the presence of 2 g/L furfural from the beginning of fermentation were carried out to study the physiological and transcriptional profiles of short-term and long-term effects of furfural on 8b. Furfural negatively affected 8b growth in terms of final biomass and the fermentation time. Transcriptomic studies indicated that the response of 8b to furfural was dynamic and complex, and differences existed between short-term shock and long-term stress responses. However, the gene function categories were similar with most down-regulated genes related to translation and biosynthesis, while the furfural up-regulated genes were mostly related to general stress responses. Several gene candidates have been identified and genetic studies indicated that expression of ZMO0465 and cysteine synthase operon ZMO0003-0006 driven by its native promoter in a shuttle vector enhanced the furfural tolerance of 8b. In addition, the relationship between microarray and mRNA-Seq was compared with good correlations. The directional mRNA-Seq data not only provided the gene expression profiling, but also can be applied for transcriptional architecture improvement to identify and confirm operons, novel transcripts, hypothetical gene functions, transcriptional start sites, and promoters with different strength.

09 BIOMASS FUELS↗
DOE OSTI2021

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Establishing temperature from radiance of dynamically compressed metals

Precise temperature determination is a significant challenge in extreme environments of dynamic compression studies. How can radiance measurements taken in high-pressure shock experiments constrain temperature in a meaningful and physically consistent way? Experiments maintaining sample compression against a transparent window can be tailored to present a uniform measurement area with uncertain spectral emissivity. We compare several methods to analyze radiance collected at multiple wavelengths, applying statistical methods and physical principles to improve temperature inference. With proper radiance collection and analysis, dynamic temperature uncertainties become comparable to thermomechanical ambiguities of the emitting surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A field programmable gate array-based timing and control system for the dynamic compression sector

A field programmable gate array (FPGA)-based timing and trigger control system has been developed for the Dynamic Compression Sector (DCS) user facility located at the Advanced Photon Source (APS) at Argonne National Laboratory. The DCS is a first-of-its-kind capability dedicated to dynamic compression science. All components of the DCS laser shock station—x-ray choppers, single-shot shutter, internal laser triggers, and shot diagnostics—must be synchronized with respect to the arrival of x rays in the hutch. An FPGA synchronized to the APS storage ring radio frequency clock (352 MHz) generates trigger signals for each stage of the laser and x-ray shutter system with low jitter. The developed FPGA-based control system was the first system used to control the laser and the shutter system since its commissioning, and it has been developing since then to improve the timing jitter. The system is composed of a Zynq FPGA, a debug card, line drivers, and a power supply. The delay and offsets of the trigger signals can be adjusted by using a user-friendly graphical user interface with high precision. Here, the details of the system architecture, timing requirements, firmware, and software implementation along with the performance evaluation are presented in this paper. The system offers low timing jitter (15.5 ps rms) with respect to the APS 352 MHz clock, suitable for the 100 ps (FWHM) x-ray bunch duration at the APS.

47 OTHER INSTRUMENTATION↗