Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock waves”

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 217 records · Page 12

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↗

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↗

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↗

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↗

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↗

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↗

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Acoustic Observations of the OSIRIS-REx Sample Return Capsule Re-Entry from Wendover Airport

The Origins, Spectral Interpretation, Resource Identification, and Security‐Regolith Explorer sample return capsule (SRC) re‐entered the Earth’s atmosphere at hypersonic speeds from interplanetary space on 24 September 2023. The current work reports on 18 ground‐based acoustic sensors deployed at Wendover Airport, the same location that the Genesis and Stardust SRC re‐entries were recorded. Four different sensors (Chaparral Physics, Gem, Wilson Engineering Research and Development [WERD], and RedVox) were deployed in close proximity to compare their performance. All the sensors captured an N‐wave signal associated with the SRC re‐entry shock wave followed by a broadband coda. The Chaparral Physics array served as the high‐fidelity reference measurement. The N‐wave signal had a peak‐to‐peak amplitude of 4.07 Pa with a fundamental frequency of 4.98 Hz from 167.5° measured clockwise from north, nearly perpendicular to the SRC trajectory. In addition, high coherence in the coda was shown to be associated reflections from the surrounding mountains. In general, the more economical sensors (Gem, WERD, and RedVox) produced results that were consistent with these observations and sensor specifications. Beamforming with these single sensors arranged as an array showed agreement with the high‐fidelity array to within a couple of degrees. Furthermore, the current high‐fidelity results were compared with the measurements during the Genesis and Stardust SRC re‐entries. All three entries produced a broadband fundamental peak at a frequency that was inversely related to the SRC diameter as well as evidence of reflections from the surrounding topography.

KC, Real J. [Oklahoma State University, Stillwater↗

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite Element Simulation of the Acoustic Pressure Inside a Beverage Container for Non-Thermal, Ultrasound-based Pasteurization

The purpose of this effort is to investigate whether large acoustic pressure waves can be transmitted inside beverage containers to enable pasteurization. Acoustic waves are known to induce large nonlinear compressive forces and shock waves in fluids, suggesting that compression waves may be capable of damaging bacteria inside beverage containers without appreciably increasingly the temperature or altering the freshness and flavor of the beverage contents. Although a combined process such as thermosonication (e.g., sonication with heating) is likely more efficient, it is instructive to compute the acoustic pressure field distribution inside the beverage container. The COMSOL simulations used two and three-dimensional models of beverage containers placed in a water bath to compute the acoustic pressure field. A limitation of these COMSOL models is that they cannot determine the bacterial lysis efficiency, rather the models provide an indirect metric of bacterial lysis based on the magnitude of the pressure field and its distribution.

42 ENGINEERING↗

Magnesium oxide at extreme temperatures and pressures studied with first-principles simulations

We combine two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, to determine the equation of state of magnesium oxide in the regime of warm dense matter, with densities ranging from 0.35 to 71 g cm -3 and temperatures ranging from 10 000 K to 5 × 10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock wave compression measurements and inertial confinement fusion experiments. We study the electronic structure of MgO and the ionization mechanisms as a function of density and temperature. We show that the L-shell orbitals of magnesium and oxygen hybridize at high density. This results in a gradual ionization of the L-shell with increasing density and temperature. In this regard, MgO behaves differently from pure oxygen, which is reflected in the shape of the MgO principal shock Hugoniot curve. The curve of oxygen shows two compression maxima, while that of MgO shows only one. We predict a maximum compression ratio of 4.66 to occur for a temperature of 6.73 × 10 7 K. Lastly, we research how multiple shocks and ramp waves can be used to cover a large range of densities and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗

Laser Shock Modeling Archival Discussion

The purpose of this discussion is to provide archival information related to the Laser Shock finite element modeling effort. As discussed in Laser Shock System, Assessing bond strength in layered materials (see Section 8.0), the Laser Shock System creates a high-amplitude shockwave on the frontside of a structure (i.e. an aluminum 6061-T6 plate for this discussion) via a high-energy pulsed laser. The shock wave is monitored on the back surface of the structure as a velocity time history. It is a compressive wave as it comes to the back surface but reflects as a tensile wave. If a bond exists in the structure and the reflected tensile wave exceeds its interface threshold stress, then bond rupture occurs. The desire of the Laser Shock effort is to establish (with multiple tests) an ultimate bond strength which can be used for fuel plate design calculations. In this process, the finite element modeling effort is the catalyst to mimic the Laser Shock tests and provide the damage information at the bond that is useful for fuel plate design calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A finite scale model for shock structure

In this paper we explore the use of a finite scale model of fluid dynamics to predict the finite structure of a shock wave in a perfect gas. We begin by documenting the history and issues that have arisen when Navier–Stokes theory is applied to the shock structure problem, and continue by motivating the improvement that finite scale theory might provide from its representation of inviscid (anomalous) dissipation. Furthermore, our primary results include the formulation of a traveling wave equation from finite scale theory, an analysis of the solutions of that equation as regards shock width and monotonicity properties, and an estimation of the critical parameter of the theory from experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗