SEARCH · Engineering Papers

Results for “Shock wave experiments”

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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Viscoelastic analysis of mussel threads reveals energy dissipative mechanisms

Mussels use byssal threads to secure themselves to rocks and as shock absorbers during cyclic loading from wave motion. Byssal threads combine high strength and toughness with extensibility of nearly 200%. Researchers attribute tensile properties of byssal threads to their elaborate multi-domain collagenous protein cores. Because the elastic properties have been previously scrutinized, we instead examined byssal thread viscoelastic behaviour, which is crucial for withstanding cyclic loading. By targeting protein domains in the collagenous core via chemical treatments, stress relaxation experiments provided insights on domain contributions and were coupled with in situ small-angle X-ray scattering to investigate relaxation-specific molecular reorganizations. Results show that when silk-like domains in the core were disrupted, the stress relaxation of the threads decreased by nearly 50% and lateral molecular spacing also decreased, suggesting that these domains are essential for energy dissipation and assume a compressed molecular rearrangement when disrupted. A generalized Maxwell model was developed to describe the stress relaxation response. The model predicts that maximal damping (energy dissipation) occurs at around 0.1 Hz which closely resembles the wave frequency along the California coast and implies that these materials may be well adapted to the cyclic loading of the ambient conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Shock compression of poly(methyl methacrylate) PMMA in the 1000 GPa regime: Z machine experiments

Hydrocarbon polymers are used in a wide variety of practical applications. In the field of dynamic compression at extreme pressures, these polymers are used at several high energy density (HED) experimental facilities. One of the most common polymers is poly(methyl methacrylate) or PMMA, also called Plexiglass® or Lucite®. Here, we present high-fidelity, hundreds of GPa range experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for PMMA to more than threefold compression up to 650 GPa and re-shock Hugoniot states up to 1020 GPa in an off-Hugoniot regime, where experimental data are even sparser. These data can be used to put additional constraints on tabular equation of state (EOS) models. The present results provide clear evidence for the need to re-examine the existing tabular EOS models for PMMA above ~120 GPa as well as perhaps revisit EOSs of similar hydrocarbon polymers commonly used in HED experiments investigating dynamic compression, hydrodynamics, or inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Collective optical Thomson scattering in pulsed-power driven high energy density physics experiments (invited)

Optical collective Thomson scattering is used to diagnose magnetized high energy density physics experiments at the Magpie pulsed-power generator at Imperial College London. The system uses an amplified pulse from the 2nd harmonic of a Nd:YAG laser (3 J, 8 ns, 532 nm) to probe a wide diversity of hightemperature plasma objects; with densities in the range of 10 17 -10 19 cm -3 and temperatures between 10 eV and a few keV. The scattered light is collected from 100 μm-scale volumes within the plasmas, which are imaged onto optical fiber arrays. Multiple collection systems observe these volumes from different directions, providing simultaneous probing with different scattering K-vectors (and different associated α-parameters, typically in the range 0.5 – 3) allowing independent measurements of separate velocity components of the bulk plasma flow. The fiber arrays are coupled to an imaging spectrometer with a gated ICCD. The spectrometer is configured to view the ion-acoustic waves (IAWs) of the collective Thomson scattered spectrum. Fits to the spectra with the theoretical spectral density function S(K,ω) yield measurements of the local plasma temperatures and velocities. Fitting is constrained by independent measurements of the electron density from laser interferometry, and the corresponding spectra for different scattering vectors. Furthermore, this TS diagnostic has been successfully implemented on a wide range of experiments, revealing temperature and flow velocity transitions across magnetized shocks, inside rotating plasma jets and imploding wire arrays, as well as providing direct measurements of drift velocities inside a magnetic reconnection current sheet.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

A simple hydride model for cerium ejecta particles

Cerium ejecta particles form after a shock wave impacts a cerium plate with a known roughness or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle relative to the gas. A simple model is developed, assuming that the particle remains intact and spherical, to capture the heat and mass transfer that occurs during this process. Model performance is evaluated by comparing with cerium ejecta experiments where the particles are ejected into deuterium gas at initial pressures of 4 and 8 atm. Overall, the model is able to capture the approximately 400 K increase in particle temperature above the surrounding gas temperature.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A simple hydride model for cerium ejecta particles

36 MATERIALS SCIENCE↗

DOE OSTI2025

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

X-ray phase-contrast imaging of strong shocks on OMEGA EP

The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. Furthermore, the interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Electrical conductivity of Sn at high pressure and temperature

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

36 MATERIALS SCIENCE↗

DOE OSTI2019

Design and Scaling of an Omega-EP Experiment to Study Cold Streams Feeding Early Galaxies

Galaxies form in the centers of dark matter halos, and grow by accreting gas from the halos. The gas supply to the galaxy is the bottleneck for star formation and ultimately sets its stellar mass and star formation rate. As gas falls from the cosmic web into the halos, the highly supersonic infalling gas (Mach ~20 for Milky Way-type halos) eventually forms a shock. If cooling is inefficient, this shock quickly stabilizes as a quasi-static slowly expanding accretion shock at the outer edge of the halo. If, however, cooling is important, a recent theory suggests that filaments of gas from the cosmic web will freefall unshocked through the halos, and shock only at the outskirts of the central galaxies. Furthermore, these "cold streams" allow for a steady supply of gas into the galaxies and for more efficient star formation. The cold dense filament flowing into a hot less dense environment is potentially Kelvin–Helmholtz unstable. This instability may hinder the ability of the stream to deliver gas deeply enough into the halo. A design of a well-scaled laser experiment on Omega-EP meant for studying this phenomena is presented in the current work, with relevant theory informed by 2D hydrodynamic simulations of the experiment. We establish the hydrodynamic scaling analysis between the cosmological system and its experimental analog, presenting the experiment as an adiabatic upper limit to informing the role of mixing due to the Kelvin–Helmholtz instability.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Performance Characteristics of a Rotating Detonation Combustor Exiting Into a Pressurized Plenum to Simulate Gas Turbine Inlet

Abstract The present study aims to experimentally characterize the performance of a rotating detonation combustion (RDC) system integrated with a pressurized downstream plenum to simulate the high-pressure inlet conditions of power-generating gas turbines. A thorough understanding of the operational behavior including wave mode behavior, static pressure profile along the combustor length, and dynamic features of pressure fluctuations is crucial for successful integration of RDC with the turbine. In this study, two RDC configurations are investigated, RDC with a constant area annulus and RDC with a converging nozzle. In both cases, the RDC flow exited into a plenum chamber kept at pressures varying from 155 kPa to 330 kPa. RDC was operated on methane and oxygen-enriched air to represent reactants used in land-based power generation. Experiments were conducted for the two RDCs configurations operated at three reactant mass flow rates (0.23, 0.32, and 0.46 kg/s). The RDC performance is characterized by time-averaged static pressure measurements, and wave velocity determined by ionization probes. In addition, dynamic pressure measurements were recorded both inside and near the exit of RDC channel to investigate wave interactions between RDC and downstream plenum. Results show that the RDC with the converging nozzle achieved superior performance while minimizing detrimental interactions with the reflected shock and/or acoustic waves from the downstream plenum.

Engineering↗

DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Insight into the chemistry of TNT during shock compression through ultrafast absorption spectroscopies

Thin films of trinitrotoluene (TNT) were shock compressed using the ultrafast laser shock apparatus at Los Alamos National Laboratory. Visible (VIS) and mid-infrared (MIR) transient absorption spectroscopies were simultaneously performed to probe for electronic and vibrational changes during shock compression of TNT. Additionally, three shock pressures (16 GPa, 33 GPa, and 45 GPa) were selected to observe no reaction, incipient reaction, and strongly developed reactions for TNT within the experimental time scale of <250 ps. Negligible absorption changes in MIR or VIS absorptions were observed at 16 GPa. At 33 GPa, MIR absorptions in the 3000 cm -1 –4000 cm -1 range were observed to increase during the shock and continue to increase during the rarefaction, in contrast to the VIS absorption measurements, which increased during the shock and almost fully recovered during rarefaction. At 45 GPa, both VIS and MIR absorptions were strong and irreversible. The intense and spectrally broad MIR absorptions were attributed to short lived intermediates with strong, spectrally broad absorptions that dominate the spectral response. The MIR and VIS absorption changes observed at 33 GPa and 45 GPa were credited to shock induced chemistry, most likely including the formation of a very broad hydrogenic stretch feature. The results from these experiments are consistent with the chemical mechanisms that include O–H or N–H formation such as CH 3 oxidation or C–N homolysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗