Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock structure”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measuring the structure and equation of state of polyethylene terephthalate at megabar pressures

Abstract We present structure and equation of state (EOS) measurements of biaxially orientated polyethylene terephthalate (PET, $$({\hbox {C}}_{10} {\hbox {H}}_8 {\hbox {O}}_4)_n$$ ( C 10 H 8 O 4 ) n , also called mylar) shock-compressed to ( $$155 \pm 20$$ 155 ± 20 ) GPa and ( $$6000 \pm 1000$$ 6000 ± 1000 ) K using in situ X-ray diffraction, Doppler velocimetry, and optical pyrometry. Comparing to density functional theory molecular dynamics (DFT-MD) simulations, we find a highly correlated liquid at conditions differing from predictions by some equations of state tables, which underlines the influence of complex chemical interactions in this regime. EOS calculations from ab initio DFT-MD simulations and shock Hugoniot measurements of density, pressure and temperature confirm the discrepancy to these tables and present an experimentally benchmarked correction to the description of PET as an exemplary material to represent the mixture of light elements at planetary interior conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2025

High-Entropy Alloys for Accelerator Beam Window Applications

Development of novel high-entropy alloys (HEAs) is currently underway for potential use as beam windows in future multi-megawatt target systems at Fermilab. HEAs encompass a new class of materials with a vast design space allowing for material properties to be tailored for particular applications and to potentially offer improved resistance to beam-induced radiation damage and thermal shock effects. The alloy systems being studied consist of several compositions of AlCoCrMnTiV with 4 6 component elements. These alloys are all predicted by CALPHAD simulation to have a single-phase BCC crystal structure and low density, with some compositions displaying ordered, nanoscale precipitates. This presentation will briefly discuss alloy design and synthesis before giving a detailed description of the characterization studies of these HEAs in both the pristine state and post-irradiation by low-energy heavy ions to high damage levels. Electron microscopy techniques to quantify elemental homogeneity and composition, determine grain size, shape, and orientation, and quantify lattice parameters, defect structures and precipitate phases are all being used to study alloy microstructures. Mechanical properties of the alloys at the microscale will be reported. The evolution of these properties as a function of radiation damage will also be described. Bulk thermal characteristics of these HEAs have been tested to measure specific heat capacity and coefficient of thermal expansion as a function of temperature. To determine bulk tensile properties a miniature tensile testing apparatus is under development; it s commissioning will be covered briefly. The talk will conclude with our plans for alloy down-selection.

Burleigh, A. [Fermilab]↗
DOE OSTI2024

The Physics of Core-Collapse Supernovae: Explosion Mechanism and Explosive Nucleosynthesis

Recent developments in multi-dimensional simulations of core-collapse supernovae have considerably improved our understanding of this complex phenomenon. In addition to that, one-dimensional (1D) studies have been employed to study the explosion mechanism and its causal connection to the pre-collapse structure of the star, as well as to explore the vast parameter space of supernovae. Nonetheless, many uncertainties still affect the late stages of the evolution of massive stars, their collapse, and the subsequent shock propagation. In this review, we will briefly summarize the state-of-the-art of both 1D and 3D simulations and how they can be employed to study the evolution of massive stars, supernova explosions, and shock propagation, focusing on the uncertainties that affect each of these phases. Finally, we will illustrate the typical nucleosynthesis products that emerge from the explosion.

Boccioli, Luca (ORCID:000000024819310X)↗
DOE OSTI2022

Samarium: from a distorted-fcc phase to melting under dynamic compression using in-situ x-ray diffraction

Abstract Lattice and electronic structure interactions for f -electrons are fundamental challenges for lanthanide equation of state development. Difficulties in first-principles calculations, such as density functional theory (DFT), emphasize the need for well-characterized experimental data. Here, we measure in-situ x-ray diffraction of shocked samarium (Sm) and temperature along the Hugoniot for the first time, providing direct evidence for phase transitions. We report direct evidence of a distorted fcc (dfcc) phase at 23 GPa. Shocked samarium melts from the dfcc phase starting at 33 GPa (1333 K), with complete melt at 40 GPa (1468 K). Previous work indicated shock melt at 27 GPa (1200 K), underscoring the significance of x-ray measurements for detecting phase transitions. Interestingly, our observed melting is in sharp contrast with the melting reported by a diamond anvil cell study. These experimental data can tightly constrain first principles calculations and serve as key touchstones for equation of state modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evidence for flare-accelerated particles in large scale loops in the behind-the-limb gamma-ray solar flare of September 29, 2022

We report on the detection of the gamma-ray emission above 100 MeV from the solar flare of September 29, 2022, by Fermi LAT with simultaneous coverage in HXR by Solar Orbiter STIX. The Solar Orbiter-Earth separation was 178° at the time of the flare as seen from Earth, with Solar Orbiter observing the east limb. Based on STIX imaging, the flare was located 16° behind the eastern limb as seen from Earth. The STIX and GBM non-thermal emission and the LAT emission above 100 MeV all show similarly shaped time profiles, and the Fermi profiles peaked only 20 s after the STIX signal from the main flare site, setting this flare apart from all the other occulted flares observed by Fermi LAT. The radio spectral imaging based on the Nançay Radioheliograph and ORFEES spectrograph reveal geometries consistent with a magnetic structure that connects the parent active region behind the limb to the visible disk. We studied the basic characteristics of the gamma-ray time profile, in particular, the rise and decay times and the time delay between the gamma-ray and HXR peak fluxes. We compared the characteristics of this event with those of four Fermi LAT behind-the-limb flares and with an on-disk event and found that this event is strikingly similar to the impulsive on-disk flare. Based on multiwavelength observations, we find that the gamma-ray emission above 100 MeV originated from ions accelerated in the parent active region behind the limb and was transported to the visible disk via a large magnetic structure connected to the parent active region behind the limb. Our results strongly suggest that the source of the emission above 100 MeV from the September 29, 2022 flare cannot be the CME-driven shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The 2024 July 16 solar event: a challenge to the coronal mass ejection origin of long-duration gamma-ray flares

We present a multi-spacecraft analysis of the 2024 July 16 long-duration gamma-ray flare (LDGRF) detected by the Large Area Telescope on the Fermi satellite. The measured > 100 MeV γ-ray emission persisted for over seven hours after the flare impulsive phase, and was characterized by photon energies exceeding 1 GeV and a remarkably hard parent-proton spectrum. In contrast, the phenomena related to the coronal mass ejection (CME)-driven shock linked to this eruption were modest, suggesting an inefficient proton acceleration unlikely to achieve energies well above the 300 MeV pion-production threshold to account for the observed γ-ray emission. Specifically, the CME was relatively slow (∼600 km/s) and the accompanying interplanetary type-II/III radio bursts were faint and short-lived, unlike those typically detected during large events. In particular, the type-II emission did not extend to kilohertz frequencies and disappeared ∼5.5 hours prior to the LDGRF end time. Furthermore, the associated solar energetic particle (SEP) event was very weak, short-duration, and limited to a few tens of MeV, even at magnetically well-connected spacecraft. These findings demonstrate that a very fast CME resulting in a high-energy SEP event is not a necessary condition for the occurrence of LDGRFs, challenging the idea that the high-energy γ-ray emission is produced by the back-precipitation of shock-accelerated ions into the solar surface. The alternative origin scenario based on local particle trapping and acceleration in large-scale coronal loops is instead favored by the observation of giant arch-like structures of hot plasma over the source region that persisted for the entire duration of this LDGRF.

Sun: UV radiation↗
DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Recovery of a high-pressure phase formed under laser-driven compression

The recovery of metastable structures formed at high pressure has been a long-standing goal in the field of condensed matter physics. While laser-driven compression has been used as a method to generate novel structures at high pressure, to date no high-pressure phases have been quenched to ambient conditions. Here we demonstrate, using in situ x-ray diffraction and recovery methods, the successful quench of a high-pressure phase which was formed under laser-driven shock compression. We show that tailoring the pressure release path from a shock-compressed state to eliminate sample spall, and therefore excess heating, increases the recovery yield of the high-pressure ω phase of zirconium from 0% to 48%. Our results have important implications for the quenchability of novel phases of matter demonstrated to occur at extreme pressures using nanosecond laser-driven compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

X-Ray Diffraction of Ramp-Compressed Silicon to 390 GPa

Silicon (Si) exhibits a rich collection of phase transitions under ambient-temperature isothermal and shock compression. This report describes in-situ diffraction measurements of ramp compressed Si between 40 and 389 GPa. Angle-dispersive x-ray scattering reveals that Si assumes an hexagonal close-packed (hcp) structure between 40 and 93 GPa, and, at higher pressure, a face-centered cubic (fcc) structure which persists to at least 389 GPa, the highest pressure for which the crystal structure of Si has been investigated. Furthermore, the range of hcp stability extends to higher pressures and temperatures than predicted by theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Shock wave formation from head-on collision of two subsonic vortex rings

Vortex ring collisions have attracted intense interest in both water and air studies (Baird in Proc R Soc Lond Ser Math Phys Sci 409:59–65, 1987, Poudel et al. in Phys Fluids 33:096105, 2021, Lim and Nickels in Nature 357:225, 1992, New et al. in Exp Fluids 57:109, 2016, Suzuki et al. in Geophys Res Lett 34, 2007, Yan et al. in J Fluids Eng 140:054502, 2018, New et al. in J Fluid Mech 899, 2020, Cheng et al. in Phys Fluids 31:067107, 2019, Hernández and Reyes in 29:103604, 2017, Mishra et al. in Phys Rev Fluids, 2021, Zednikova et al. in Chem Eng Technol 42:843–850, 2019, Kwon et al. in Nature 600:64–69, 2021). These toroidal structures spin around a central axis and travel in the original direction of impulse while spinning around the core until inertial forces become predominant causing the vortex flow to spontaneously decay to turbulence (Vortex Rings, https://projects.iq.harvard.edu/smrlab/vortex-rings). Previous studies have shown the collision of subsonic vortex rings resulting in reconnected vortex rings, but the production of a shock wave from the collision has not been demonstrated visibly (Lim and Nickels in Nature 357:225, 1992, Cheng et al. in Phys Fluids 31:067107, 2019). Here we present the formation of a shock wave due to the collision of explosively formed subsonic vortex rings. As the vortex rings travel at Mach 0.66 toward the collision point, they begin to trap high pressure air between them. Upon collision, high pressure air was imploded and released radially away from the axis of the collision, generating a visible shock wave traveling through and away from the colliding vortices at Mach 1.22. Our results demonstrate a pressure gradient with high pressure release creating a shock wave. We anticipate our study to be a starting point for more explosively formed vortex collisions. For example, explosives with different velocities of detonation could be tested to produce vortex rings of varying velocities.

42 ENGINEERING↗
DOE OSTI2023

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

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

Early structure formation constraints on the ultralight axion in the postinflation scenario

Many works have concentrated on the observable signatures of the dark matter being an ultralight axion-like particle (ALP). We concentrate on a particularly dramatic signature in the late-time cosmological matter power spectrum that occurs if the symmetry breaking that establishes the ALP happens after inflation -- white-noise density fluctuations that dominate at small scales over the adiabatic fluctuations from inflation. Here, these fluctuations alter the early history of nonlinear structure formation. We find that for symmetry breaking scales of f A ~ 10 13 – 10 15 GeV, which requires a high effective maximum temperature after inflation, ALP dark matter with particle mass of m A ~ 10 –13 – 10 –20 eV could significantly change the number of high-redshift dwarf galaxies, the reionization history, and the Lyα forest. We consider all three observables. We find that the Lyα forest is the most constraining of current observables, excluding f A ≳ 10 15 GeV (m A ≲ 10 –17 eV) in the simplest model for the ALP and considerably lower values in models coupled to a hidden strongly interacting sector (f A ≳ 10 13 GeV and m A ≲ 10 –13 eV). Observations that constrain the extremely high-redshift tail of reionization may disfavor similar levels of isocurvature fluctuations as the forest. Future z ~ 20–30 21cm observations have the potential to improve these constraints further using that the supersonic motions of the isocurvature-enhanced abundance of ~10 4 M ⊙ halos would shock heat the baryons, sourcing large BAO features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2025

A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations (2006–2011)

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense.

astrometry↗
DOE OSTI2024

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Selected 3D Flash-X Checkpoints for the Long-Time Evolution of a 9.6 Solar-Mass Core-Collapse Supernova Model

This dataset contains selected 3D Flash-X checkpoint files from the long-time evolution of a low-energy core-collapse supernova explosion of a 9.6 Msun zero-metallicity, low-mass iron-core progenitor. The checkpoints span the shock-breakout phase through the young-remnant phase, ending at approximately 3 yr after core bounce. The dataset is intended for follow-up analysis and post-processing, especially radiation-transport calculations using the hydrodynamic and compositional structure of the ejecta. For a full description of the numerical setup, physical assumptions, limitations, and interpretation of the simulation, users should refer to the associated paper.

79 ASTRONOMY AND ASTROPHYSICS↗