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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Micro-strains, local stresses, and coherently diffracting domain size in shock compressed Al(100) single crystals

In situ x-ray diffraction (XRD) measurements and their analysis in Al single crystals shock compressed along the [100]-direction were utilized to examine shock wave induced microstructural heterogeneities. High-resolution XRD line profiles for the 200, 400, and 600 Al peaks were measured in uniaxial strain compression states to either 5.6 or 11.7 GPa and partial stress release to 3.5 or 6.6 GPa, respectively. Broadening of the XRD line profiles was analyzed to determine the magnitude of the longitudinal micro-strain distribution (0.195% and 0.28% full width at half maximum for 3.5 and 6.6 GPa stresses, respectively) and the size of coherently diffracting domains (CDDs) (0.125 and 0.07 μm for 3.5 and 6.6 GPa stresses, respectively). From the longitudinal micro-strain distributions, the distribution of local stress differences (or stress deviators) was obtained in the shocked state. The full width at half maximum of this distribution, a measure of the local stress inhomogeneities, is greater than half of the macroscopic stress difference for both 3.5 and 6.6 GPa peak stresses, suggesting considerable variation in local stress deviators. The CDD sizes determined here are comparable to characteristic length scales for defect-free regions determined from defect density measurements in post-shock recovery experiments. As a result, the present work represents an important step in understanding material microstructure and inhomogeneities in the shock-compressed state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Design and optimization of supersonic turbines for detonation combustors

Detonation-based engines offer a potential surge in efficiency for compact thermal power systems. However, these cycles require ad-hoc components adapted to the high outlet velocity from the detonation combustors. This paper presents the design methodology of turbine stages suitable for supersonic inlet conditions and provides a detailed analysis of optimized turbine geometries. A reduced-order solver examines the supersonic blade rows’ functional design space, quantifies the turbine’s non-isentropic performance, and budgets the turbine loss for different optimized leading-edge designs and chord to pitch ratios. The shock-wave interactions were identified as the predominant contributor to turbine losses, and optimal pitch-chord ratios were determined for various inlet Mach numbers. Finally, with this tool, the specific-power output for a wide range of design configurations was computed; and the metal angle that ensures flow starting and maximizes power extraction was calculated. The detailed numerical study describes the flow interactions in a supersonic turbine and offers new correlations to guide the design of future supersonic turbines.

42 ENGINEERING↗
DOE OSTI2020

Calibrating the Galactic Cepheid Period–Luminosity Relation from the Maximum-likelihood Technique

We present a realization of the maximum-likelihood technique, which is one of the latest modifications of the Baade–Becker–Wesselink method. Our approach is based on nonlinear calibrations of the effective temperature and bolometric correction, which take into account metallicity and surface gravity. It allows one to estimate the key Cepheid parameters, the distance modulus, and the interstellar reddening, combining photometric and spectroscopic data (including the effective temperature data). This method is applied to a sample of 44 Galactic Cepheids for which multiphase temperature measurements are available. The additional data correction is performed to subtract the impact of the component in binary/multiple systems. We also study the effect of shock waves, whose presence in the stellar atmosphere distorts the observational data and leads to systematic errors in the obtained parameters. We determine the optimal restriction on the input data to eliminate this effect. This restriction provides accurate period–radius and period–luminosity relations that are consistent with the results in previous studies. We found the following relations: log R = (0.68 ± 0.03) · log P + (1.14 ± 0.03), M {sub v} = − (2.67 ± 0.16) · (log P − 1) − (4.14 ± 0.05).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Atmospheric clearing of solid-particle debris using femtosecond filaments

Through nonlinear self-focusing, femtosecond pulses can propagate several kilometers beyond diffraction limits, forming an ionization channel in air known as a laser filaments. It has been demonstrated that in the wake of the filament, aerosols can be effectively cleared to improve the transmission of subsequent laser pulses or secondary light sources, pertinent to applications in atmospheric sensing. However, the current understanding of aerosol clearing is founded on interactions with droplets to simulate fogs and clouds and thus does not extend to solid particles or atmospheric debris. Using optical trapping, we isolate both graphite and silica microparticles and directly measure the subsequent displacement caused by the filament using time-resolved shadowgraphy. The shock wave from the filament is demonstrated to propel particles away from the filament, directly contributing to atmospheric debris clearing. Particles exposed to the laser light in either the intense filament core or the surrounding energy reservoir are axially displaced along the beam path. Furthermore, it is found that the optomechanical properties of the particle largely influence the axial displacement induced by laser exposure through mechanisms such as radiation pressure, mass ejection from ablation or optical damage, and particle deagglomeration.

42 ENGINEERING↗
DOE OSTI2021

Peak states of molybdenum single crystals shock compressed to high stresses

In this work, to determine crystal anisotropy effects at high stresses, peak states behind the plastic shock waves were examined in BCC single crystals. Using plate impact experiments, molybdenum (Mo) single crystals were shock compressed up to 190 GPa elastic impact stress along [100], [110], and [111] orientations. Laser interferometry was used to measure wave velocities and particle velocity profiles at the Mo–LiF window interface. These data were analyzed to obtain in-material quantities in the peak states. The Hugoniots for [100] and [110] orientations were comparable, but the Hugoniot for the [111] orientation was different from the other two orientations. Also, these Mo single crystal Hugoniots display differences from the polycrystalline Mo Hugoniots. Although none of the differences can be considered large, the present results demonstrate that, unlike FCC metal single crystals (Cu, Al), some anisotropy is preserved in Mo single crystal Hugoniots even at high stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2023

Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Dissertation]

This thesis presents the computational design, modeling, and analysis of three experiments in high-energy density physics (HEDP), all of them concerning fundamental radiation flows. The first experiment is a laboratory astrophysics experiment to investigate the role of the Kelvin-Helmholtz instability (KHI) in the process of galactic filaments supplying gas to galactic halos. The achieved goal was to provide a first study in which the role of the instability is maximal and predict behavior in future iterations of an experiment accessing a more radiative regime where the role of the instability is stifled. This experiment would help answer how certain galaxies are able to grow so rapidly and produce many stars as the KHI limits this process. The second experiment, COAX, is a radiation flow experiment with a novel spectroscopy diagnostic configuration, designed to spatially measure the temperature of a radiation wave as it travels down a doped foam. A key result of this work was the development of a synthetic spectroscopy application and application of modern spectroscopy comparison techniques to provide our first temperature reconstructions from the experimental data. This experimental platform serves as the launching ground for a number of new experiments that vary the basic premise and thus is foundational to our ongoing research. The final experiment is a full integration of modeling, design, and theoretical development for the Radishock experiment. This experimental platform studies the head-on collision of a radiation wave with a counter-propagating shock, and like COAX, uses spectroscopy to diagnose and detect the interaction. My research analyzes the successful shots, indicating aspects of successful detections and suggests improvements to future iterations of the design.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Potential energy surfaces for high-energy N + O 2 collisions

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.

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

Simulation and analysis of the 3D double shock-induced recompression of spall using a hybrid PAGOSA with FLIP+MPM

In this work, an innovative hybrid algorithm PAGOSA with FLIP + MPM is first presented by coupling PAGOSA with the particle FLIP + MPM, and applied to systematically exploring the 3D recompression of spall induced by double-shock waves in ductile materials as the few existing studies are currently limited to the pure longitudinal recompression of spall, and ignore the 3D effects. This algorithm solves the difficulties the grid-based methods encounter when applied to capturing the fracture in material. We first validate the capabilities of PAGOSA with FLIP + MPM to predict different fracture/fragmentation cases by solving two benchmark problems and comparing with the analytical solutions or the experiment results. The convergences and the infinity-norm errors are also investigated. Subsequently, the 3D effects are demonstrated by simulating the spallation in material using PAGOSA with FLIP + MPM and comparing with those in the longitudinal cases. Additionally, some factors that affect the spallation with 3D effects are also discussed. Finally, the recompressions of spall in material with/without 3D effects are systematically simulated and analyzed. The recompression of spall driven by a high explosive is also simulated to show the ability of PAGOSA with FLIP + MPM to handle this kind of problem. The conditions that can lead to a recompression of spall are explored. Also, numerical results show that PAGOSA with FLIP + MPM can accurately predict different fracture cases, that the 3D effects play an important role in fracture/fragmentation, and the recompression of spall is very sensitive to the occurrence conditions, shows great differences when considering 3D effects in real applications. Moreover, the present hybrid method can be easily extended to other grid-based techniques employed for fracture in materials.

3D effects↗
DOE OSTI2026

Improving Chirped Fiber Bragg Grating Resolution for Position-Sensitive Sensors in Shock- and Detonation-Driven Experiments

Chirped fiber Bragg gratings (CFBGs) are robust diagnostic sensors that are widely used to track detonation-driven and shock wave propagation. CFBGs are inscribed with a linearly chirped periodic index of refraction changes that alter the Bragg wavelength along the length of the probe. The light return of each individual Bragg element is captured by a detector at a unique time to map the full reflected spectrum. The CFBG spectrum is measured with a dispersive Fourier transform of the reflected light that temporally stretches the spectrum to increase spatial resolution and make a one-to-one map of the wavelength on a time axis. Here, we propose an improvement of CFBG temporal resolution by incorporating two co-linear laser pulses with orthogonal polarization states and a 5 ns time offset. The two separate signals were split and tracked by two separate detectors. An oscilloscope captured good separation in the signals, and two separate spectrograms were generated and interleaved in the post-processing of the data. This novel technique doubled the CFBG temporal resolution and led to a doubled location resolution. As a proof-of-concept of this technique, the resolution improvement was compared between standard CFBG measurements and the two polarization states method on a position-sensitive CFBG sensor. CFBG resolution doubling will advance sensor capabilities and will have a direct impact on improving capture and analysis in dynamic, high-explosive experiments.

42 ENGINEERING↗
DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Hugoniot states and optical response of soda lime glass shock compressed to 120 GPa

In contrast to relatively pure silica glass (fused silica-FS), commercial silica-rich glasses contain significant fractions of additional oxide components. In particular, soda-lime glass (SLG) consists of approximately 71% SiO 2 by weight, which raises the question: what is the effect of additional cations on the shock compression response of silica-rich glasses? To address this concern, plate impact experiments were conducted to determine the high-pressure Hugoniot states for shocked SLG (37 to 120 GPa) and compared with recently reported results on FS. Using laser interferometry, particle velocity profiles were measured at the impact surface and the SLG/LiF window interface. In all experiments, the transmitted profiles show a single shock wave, with no features indicative of a phase transformation. The Hugoniot states determined from the wave profiles are described well using a linear shock velocity–particle velocity relation. Interferometry measurements (using 1550-nm wavelength laser) at the impact surface show that SLG is transparent when shocked to 55 GPa, transparent for tens of nanoseconds after impact between 55 and 81 GPa, and opaque beyond 81 GPa. From impact surface measurements, a linear relationship between the apparent and true particle velocity was observed, resulting in a linear relationship between the refractive index (at 1550 nm) and density. At 120 GPa, the SLG density is nearly twice its ambient value, indicating that SLG can achieve highly dense amorphous states. When compared to FS, shocked SLG is much less compressible and likely does not transform to a crystalline phase. A plausible explanation for this difference is indicated.

58 GEOSCIENCES↗