Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock loading”

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 181 records · Page 10

Compression rate effects on the Bi-III stability field in dynamic diamond anvil cell XFEL studies of bismuth

This study aims to resolve conflicting observations of bismuth’s high-pressure phase transitions across static, intermediate, and shock compression regimes. We probed its high-pressure structural sequence using the dynamic diamond anvil cell with and time-resolved X-ray diffraction with microsecond resolution at an X-ray free-electron laser. At room temperature and pressures up to 20 gigapascals, bismuth evolves through the same structural sequence previously identified under static compression. However, the transformation behavior differs in important ways under dynamic loading. In particular, the transition from the incommensurate intermediate-pressure phase to the high symmetry high pressure phase begins at pressures about 2 to 4 gigapascals lower than in static experiments, and the onset of this transformation depends on the compression rate. As a result, the stability field of the intermediate-pressure phase is reduced under rapid compression. Our results suggest that at sufficiently fast loading rates the intermediate phase may be bypassed entirely, consistent with previous shock-compression observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

In situ measurement of damage evolution in shocked magnesium as a function of microstructure

Accurate modeling and prediction of damage induced by dynamic loading in materials have long proved to be a difficult task. Examination of postmortem recovered samples cannot capture the time-dependent evolution of void nucleation and growth, and attempts at analytical models are hindered by the necessity to make simplifying assumptions, because of the lack of high-resolution, in situ, time-resolved experimental data. We use absorption contrast imaging to directly image the time evolution of spall damage in metals at ~1.6-μm spatial resolution. We observe a dependence of void distribution and size on time and microstructure. The insights gained from these data can be used to validate and improve dynamic damage prediction models, which have the potential to lead to the design of superior damage-resistant materials.

36 MATERIALS SCIENCE↗

Time-resolved velocity and ion sound speed measurements from simultaneous bow shock imaging and inductive probe measurements

We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (“b-dot”) probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large (R M > 10), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach–Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the sound speed to infer the value of Z̄T e , where Z̄ is the average ionization and T e is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic (M S ~ 8), super-Alfvénic (M A ~ 2) aluminum plasma generated during the ablation stage of an exploding wire array on the Magpie generator (1.4 MA, 250 ns). The velocity and Z̄T e measurements agree well with the optical Thompson scattering measurements reported in the literature and with 3D resistive magnetohydrodynamic simulations in GORGON.

47 OTHER INSTRUMENTATION↗

Hands-on Introduction to Shock Physics (2020 Report, RR2020-IPD-1)

Los Alamos National Laboratory (LANL) researchers have played key roles in the development of high-pressure science and shock physics since the 1950s. Central to the Laboratory’s experimental capability are multiple impact systems and explosive loading capabilities that span impact velocities in excess of 8 km/s. More recently, our efforts have taken us to 3rd and 4th generation light sources where these loading platforms are coupled to X-ray diagnostics that are now providing new insights into matter-at-extremes. As the availability of these new diagnostics and loading techniques increases, there is an increased need within the high pressure community for researchers with the knowledge base and expertise to use them. Few institutions offer any level of theoretical or experimental coursework designed to train students for work in the field of shock physics or dynamic compression science, and those that do cannot meet the current demand of the national laboratories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic Field Amplification by a Plasma Cavitation Instability in Relativistic Shock Precursors

Abstract Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

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

Crystal structures↗

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Design and Characterization of a Lens-Coupled System for Dynamic X-Ray Diffraction

X-ray diffraction (XRD) is a necessary technique for understanding states of materials under static and dynamic loading conditions. The higher-pressure Equation of State (EOS) of many materials can only be explored via shock or ramp compression at temperatures and pressures of interest. While static XRD work has yielded EOS measurements in the 100 - 200 GPa regime, dynamic X-ray diffraction (DXRD) can explore EOS phases in the TPa regime, which closely resembles inner-core planetary conditions. DXRD hinges on the ability to measure the exact phase or phase change of a material while under dynamic loading conditions. Macroscopic diagnostic systems (e.g. velocimetry and pyrometry) can infer a phase change but not identify the specific phase entered by a material. While microscopic (atomic-level) diagnostic systems (e.g. DXRD) have been designed and implemented in Department of Energy’s (DOE) National Laboratories complex, the unique nature of Sandia National Laboratories’ Pulsed Power Facility (Z Machine) prohibits the use of such devices. The destructive nature of Z experiments presents a challenge to data capture and retrieval. Furthermore there are electromagnetic interference, X-ray background, and mechanical constraints to consider. Thus, a multi-part X-ray diagnostic for use on the Z Machine and Z-Beamlet Laser system has been designed and analyzed. Portions of this new DYnamic SCintillator Optic (DYSCO) have been built, tested and fielded. A data analysis software has been written. Finally, the radiance profile of the DYSCO’s scintillator has been characterized through experiments performed at the University of Arizona.

36 MATERIALS SCIENCE↗

Introduction to the Demonstration of Pyroshock Response Prediction Methods Using NASA ShockSat Data Virtual Collection

Mechanical shocks, particularly pyroshocks, are environments that all spacecraft and launch vehicles experience. Pyroshocks can be challenging design environments for space hardware because these events are often high frequency, high intensity transients. Furthermore, since design decisions must be made relatively early in the design phase, understanding, and accurately defining pyroshock excitation and the loads that it imparts on components is important to ensure the high reliability required of spacecraft and rockets. Here, general guidelines about shock transmission through spacecraft structures are provided in NASA and MIL handbooks and standards, but these are relatively broad guidelines, and recent studies have shown that more accurate and less-conservative predictions are needed. However, predicting shock response accurately is still a challenge because of the high frequencies in a pyroshock excitation.

Babuska, Vit [Sandia National Laboratories (SNL-NM↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

Stress field measurements using quantitative schlieren

Quantitative schlieren analysis is extended here to optically transparent solids in quasi-static and dynamic experiments to measure stress distributions. The quasi-static experiments in polymethyl methacrylate (PMMA) compared refraction angles and stress gradients calculated from schlieren images to the analytical Flamant solution of a line load on a half-space. The quantitative schlieren measurements of the stress field in the thin sample with a load compared well to the analytical solution. The analysis method was then extended to explosive induced shock waves in PMMA. The explosive induced response of PMMA was experimentally studied using high-speed schlieren to visualize the shock propagation in conjunction with Photon Doppler Velocimetry (PDV) to record surface velocity histories. The stress state estimated from the schlieren images was compared to the stress calculated from the PDV measurements. High-speed imaging limitations caused the shock wave to not be fully resolved in the images, but was resolved in the PDV measurement. The stress state behind the shock calculated from the high-speed images followed a similar trend to the stress calculated from the PDV measurements.

Physics↗

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium

The “genomic shock” hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession “Bhyb26.” We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species’ similarity in transposable element load may account for the subtlety of the observed genome dominance.

59 BASIC BIOLOGICAL SCIENCES↗

Material effects on dynamics in triple-nozzle gas-puff Z pinches

Here, the gas-puff Z pinch has a long history with myriad applications as an efficient neutron or x-ray source. Its simplicity as a load configuration makes it suitable for studying fundamental plasma physics phenomena such as instabilities and energy transport. For example, the implosion of cylindrical shells onto a fusion fuel are inherently susceptible to instability growth on their external surfaces; if such instabilities are unmitigated, then the consequences in terms of degraded performance can be substantial. Similarly, mitigating heat transport from a hot fuel to its colder surrounding container can make fusion conditions more easily achievable. Here we have conducted a systematic study of triple-nozzle (outer liner, inner liner, fuel) gas puffs using two-dimensional (2D) magnetohydrodynamic simulations to investigate the effect of load material on the relevant dynamics. Analogous to past studies on spherical blast waves and converging shock waves, a trend emerges linking increased radiative cooling, lower adiabatic index, and increased magneto-Rayleigh-Taylor instability growth. Notably, our results suggest that, for the present configuration, Ar radiates less than both Ne and Kr during the early stages of the implosion while mass is being swept up and perturbations begin to seed instability growth. Consequently, pinches with Ar on the outer surface exhibit more stable 2D behavior. Here we also present a parameter scan of thermonuclear neutron yield, Y, as a function of peak current, I pk and dopant concentration with Ne or Ar, depending on the inner liner material. Above 6 MA, our results suggest Y ∝ $I$$^{5}_{pk}$ and even substantial mixing (10% by volume) of Ne into the fuel does not drastically reduce yield, suggesting an Ar/Ne/fuel configuration may reliably achieve DD thermonuclear yields of 10 13 – 10 14 /cm in the 10–20 MA range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oxidation of an iso-paraffinic alcohol-to-jet fuel and n-heptane mixture: an experimental and modelling study

Oxidation of ATJ/n-heptane blends was studied over a wide range of test conditions, using single pulse shock tubes. Test conditions were designed to study the effect of pressure (4 and 50 bar), and fuel loading (~100-1400 ppm) on the oxidation of the blends across a wide range of temperatures (800-1300 K). These effects were observed by measuring concentrations of intermediate species formed using in-line Gas Chromatography (GC) and GCxGC TOF-MS. Results showed that increasing the initial fuel load does not have a significant impact on the results. However, increasing the pressure shifts the mixture reactivity to a lower temperature by about 150 K and causes the fuel to be oxidized instead of decomposing pyrolytically as at lower pressures. Additional experiments for pure ATJ and pure n-heptane were performed at conditions matching the ATJ/n-heptane 50 bar experiments to analyze the differences between the pure and mixed fuels. Speciation measurements were compared against predictions from a detailed kinetic model. The ability of the kinetic model to capture the effects of varying experimental test conditions on the evolution of intermediate species is discussed and kinetic analyses have been conducted to identify the important reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗

A compact x-ray diffraction system for dynamic compression experiments on pulsed-power generators

Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.

47 OTHER INSTRUMENTATION↗