Search NASA⌕ Search

SEARCH · Search NASA

Results for “elastic shear 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.

82 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic Strength of Iron at High Pressures and Strain Rates

Accurate modeling of meteorite impacts, and deformation of planetary cores require characterization of the flow strength and in-elasticity of iron in its different phases. In this letter, we investigate the flow strength of both the ambient α-phase and high-pressure ϵ-phase of iron at strain-rates of 1 × 10 5 s –1 and pressures up to 42 GPa using high pressure - pressure shear plate impact experiments. We report the strength of the ϵ-iron to be significantly higher than α-phase but consequently one order smaller than the previously reported dynamic strength at high pressures. As a result, the complete stress-strain response of the ϵ-phase is reported for the first time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Temperature-induced densification in compressed basaltic glass revealed by in-situ ultrasonic measurements

Abstract Acoustic velocities of a model basalt glass (64 mol% CaMgSi2O6 + 36 mol% CaAl2Si2O8) were measured along different pressure-temperature (P-T) paths. One set of experiments involved isothermal compression-decompression cycles, performed at temperatures of 300, 641, 823, and 1006 K and pressures up to 12.2 GPa. The other set of experiments involved constant-load heating-cooling cycles at temperatures up to 823 K and pressures up to 7.5 GPa. Both sets of experiments were performed in a multi-anvil apparatus using a synchrotron-based ultrasonic technique. Our results show that the glass compressed isothermally at 300 K (cold-compression) displays anomalously decreasing compressional (VP) and shear (VS) wave velocities with increasing pressure until ~8 GPa. Beyond 8 GPa, both VP and VS start to increase sharply with pressure and irreversible densification of the glass occurred, producing large hysteresis loops of velocities upon decompression. However, for the glass compressed isothermally at increasingly higher temperatures (hot-compression), the velocity minima gradually shift to lower pressures. At temperature close to the glass transition temperature Tg, the velocity minima disappear completely, displaying a monotonic increase of velocities during compression and higher VP and VS during decompression. In addition, constant-load heating-cooling experiments show that velocities generally decrease slightly with increasing temperature, but start to increase once heated above a threshold temperature (~650 K). During cooling the velocities increase almost linearly with decreasing temperature, resulting in higher velocities (~1.5–2.5% higher) when returned to 300 K. This implies that a temperature-induced densification may have occurred in the glass at high pressures. Raman spectra on recovered samples show that the hot-compressed and high-P heated glasses contain distinctly densified and depolymerized structural signatures compared to the initial glass and the cold-compressed glass below the velocity transition pressure PT (~8 GPa). Such densification may be attributed to the breaking of bridging oxygen bonds and compaction in the intermediate-range structure. Our results demonstrate that temperature can facilitate glass densification at high pressures and point out the importance of P-T history in understanding the elastic properties of silicate glasses. Comparison with melt velocity suggests that hot-compressed glasses may better resemble the pressure dependence of velocity of silicate melts than cold-compressed glasses, but still show significantly higher velocities than melts. If the abnormal acoustic behaviors of cold-compressed glasses were used to constrain melt fractions in the mantle low-velocity regions, the melt fractions needed to explain a given velocity reduction would be significantly underestimated at high pressures.

Geochemistry & Geophysics↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2025

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

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↗
DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗