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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2022

Sound speed and refractive index of amorphous CaSiO3 upon pressure cycling to 40 GPa

Abstract Brillouin spectroscopy at room temperature and pressures up to 40 GPa documents nearly identical elasticity and refractive index of amorphous CaSiO3 created by two different methods: temperature-quenching the melt at ambient pressure and pressure-amorphizing crystalline wollastonite at room temperature. We find reproducible hysteresis of 0 to 8% on pressure cycling that is small relative to the 30 to 60% changes in shear and longitudinal wave velocities over this pressure range. Together with observed changes in refractive index and previous results from Raman spectroscopy, these measurements reveal a continuous and reversible change in atomic packing induced by pressure. Unlike many other silicate glasses, amorphous CaSiO3 exhibits highly reproducible properties, behaving smoothly and reversibly under pressure cycling and possessing similar structure and elasticity regardless of synthesis paths for the starting material, which suggests that the amorphous solid may mimic the liquid over the pressure range investigated.

Geochemistry & Geophysics↗
DOE OSTI2024

Effects of irradiation damage on the hardness and elastic properties of quaternary and high entropy transition metal diborides

Multi-principal component transition metal (TM) diborides represent a class of high-entropy ceramics (HECs) that have received considerable interest in recent years owing to their promising properties for extreme environment applications that include thermal/ environmental barriers, hypersonic vehicles, turbine engines, and next-generation nuclear reactors. While the addition of chemical disorder through the random distribution of TM elements on the cation sublattice has offered opportunities to tailor elastic stiffness and hardness, the effects of irradiation-induced structural damage on the physical properties of these complex materials have remained largely unexplored. To this end, changes in the hardness and elastic moduli of a high-entropy TM diboride (Hf 0.2 Nb 0.2 Ta 0.2 Ti 0.2 Zr 0.2 )B 2 and three of its quaternary subsets following irradiation with 10 MeV gold (Au) ions to fluences of up to 6 × 10 15 Au cm -2 are investigated at the micrometer and sub-micrometer length-scales via the dispersion of laser-generated surface acoustic waves (SAW) and nanoindentation, respectively. The nanoindentation measurements show that the TM diborides exhibit an initial increase in hardness following irradiation with energetic Au ions, with a subsequent decrease in hardness following further irradiation. One quaternary composition, (Hf 1/3 Ta 1/3 Ti 1/3 )B 2 , exhibits a notable exception to the trend and continues to exhibit an increase in hardness with ion irradiation fluence. Although differences in the absolute values of the effective elastic moduli obtained from the measured SAW dispersion and nanoindentation are observed (and attributed to microstructural variations at the measurement length-scale), both techniques yield similar trends in the form of an initial reduction and subsequent saturation in the elastic modulus with increasing ion irradiation fluence. The quaternary TM diboride (Hf 1/3 Ta 1/3 Ti 1/3 )B 2 again exhibits a departure from this trend. The high-entropy TM diboride (Hf 0.2 Nb 0.2 Ta 0.2 Ti 0.2 Zr 0.2 )B 2 exhibits the greatest recovery in hardness and modulus when irradiated to high ion fluences following initial changes at low fluence, indicating superior resistance to radiation-induced damage over its quaternary counterparts. Opportunities for designing HECs with superior hardness and modulus for enhanced radiation resistance (compared to their single constituent counterparts) by tailoring chemical disorder and bond character in the lattice are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sound Velocities of Stishovite at Simultaneous High Pressure and High Temperature Suggest an Eclogite‐Rich Layer Beneath the Hawaii Hotspot

Compressional and shear wave velocities of polycrystalline stishovite (SiO 2 ) have been measured at simultaneous high pressures and temperatures up to 14.5 GPa and 800°C. By fitting velocities to the finite strain equations, the elastic moduli and density were determined to be K S0 = 306.6(46) GPa, $K$$^{′}_{S}$ = 4.92(10), ∂K S /∂T = −0.024(1) GPa/K, G 0 = 229.0(34) GPa, G′ = 1.07(10), ∂G/∂T = −0.017(1) GPa/K, ρ 0 = 4.287(2) g/cm 3 . Our modeling suggested that, in the eclogite, coesite-stishovite transition can increase P and S wave velocities by 2.4% and 3.5%, respectively. A comparison between geophysical observations and our model shows that the coesite-stishovite phase transition in the eclogite can potentially be responsible for the occurrence of the X discontinuity beneath Hawaii. In addition, our current results suggest an eclogite-rich layer between 340 and 450 km depth beneath Hawaii. The eclogite concentration at the top and bottom of the layer is 41–55 vol% and >77 vol%, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Clarifying the definition of ‘transonic’ screw dislocations

We report a number of recent Molecular Dynamics (MD) simulations have demonstrated that screw dislocations in face centred cubic (fcc) metals can achieve stable steady state motion above the lowest shear wave speed (v shear ) which is parallel to their direction of motion (often referred to as transonic motion). This is in direct contrast to classical continuum analyses which predict a divergence in the elastic energy of the host material at a crystal geometry dependent ‘critical’ velocity v crit . Within this work, we first demonstrate through analytic analyses that the elastic energy of the host material diverges at a dislocation velocity (v crit ) which is greater than v shear , i.e. v crit >v shear . We argue that it is this latter derived velocity (v crit ) which separates ‘subsonic’ and ‘supersonic’ regimes of dislocation motion in the analytic solution. In addition to our analyses, we also present a comprehensive suite of MD simulation results of steady state screw dislocation motion for a range of stresses and several cubic metals at both cryogenic and room temperatures. At room temperature, both our independent MD simulations and the earlier works find stable screw dislocation motion only below our derived v crit . Nonetheless, in real-world polycrystalline materials v crit cannot be interpreted as a hard limit for subsonic dislocation motion. In fact, at very low temperatures our MD simulations of Cu at 10 Kelvin confirm a recent claim in the literature that true ‘supersonic’ screw dislocations with dislocation velocities v>v crit are possible at very low temperatures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Sound velocities of iron-nickel (Fe90Ni10) alloy up to 8 GPa and 773 K: The effect of nickel on the elastic properties of bcc-iron at high P-T

Sound velocities of iron and iron-based alloys at high pressure and high temperature are crucial for understanding the composition and structure of Earth’s and other telluric planetary cores. In this study, we performed ultrasonic interferometric measurements of both compressional (ν P ) and shear (ν S ) velocities on a polycrystalline body-centered-cubic (bcc)-Fe 90 Ni 10 up to 8 GPa and 773 K. The elastic moduli and their pressure and temperature derivatives are derived from least-square fits to third-order finite strain equations, yielding K S0 = 154.2(8) GPa, G 0 = 73.2(2) GPa, K' S0 = 4.6(2), G' 0 = 1.5(1), ∂K S /∂T = –0.028(1) GPa/K, and ∂G/∂T = –0.023(1) GPa/K. Furthermore, a comparison with literature data on bcc-Fe suggests that nickel not only decreases both P and S wave velocities but also weakens the temperature effects on the elastic moduli of Fe-Ni alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermo-acoustoelastic effect of Rayleigh wave: Theory and experimental verification

Previous studies showed that the thermally-induced ultrasonic bulk wave velocity change could be used to measure acoustoelastic coefficients and third-order elastic constants of elastic materials. This method is naturally immune from the ambient temperature effect and has improved sensitivity and a simpler test setup than the conventional acoustoelastic test. However, Rayleigh wave is preferred for thick components or structures with only one accessible surface. In this work, the thermo-hyperelastic constitutive equation, along with acoustoelastic theory, is used to derive the expression of the thermo-acoustoelastic coefficient (TAEC) of Rayleigh wave. The numerical relationship between the TAEC of Rayleigh wave and Murnaghan constants ( l, m and n ) are given for common metals. Further, the TAEC expressions for Rayleigh wave and shear wave are similar, and both are dominated by the constant m . The TAEC of Rayleigh wave was measured on an aluminum 6061 specimen using the thermal modulation experiment in a temperature range of 22 ~ 35°C. The measured TAEC shows good agreement with the theoretical calculation. Then the third-order elastic constants were calculated based on TAECs of bulk waves and Rayleigh wave.

42 ENGINEERING↗
DOE OSTI2022

Seismic strain energy partitioning: estimating the strain energy of seismic body waves.

This report details a method to estimate the energy content of various types of seismic body waves. The method is based on the strain energy of an elastic wavefield and Hooke’s Law. We present a detailed derivation of a set of equations that explicitly partition the seismic strain energy into two parts: one for compressional (P) waves and one for shear (S) waves. We posit that the ratio of these two quantities can be used to determine the relative contribution of seismic P and S waves, possibly as a method to discriminate between earthquakes and buried explosions. We demonstrate the efficacy of our method by using it to compute the strain energy of synthetic seismograms with differing source characteristics. Specifically, we find that explosion-generated seismograms contain a preponderance of P wave strain energy when compared to earthquake-generated synthetic seismograms. Conversely, earthquake-generated synthetic seismograms contain a much greater degree of S wave strain energy when compared to explosion-generated seismograms.

58 GEOSCIENCES↗
DOE OSTI2022

Implications of Sound Velocities of Natural Topaz on the Seismic L–Discontinuity

Compressional and shear wave velocities of topaz have been measured by ultrasonic interferometry in a multi-anvil apparatus up to 13 GPa at room temperature. By fitting the measured velocities under pressure to finite strain equations, the elastic moduli of topaz and their pressure derivatives were obtained, yielding K S0 = 165.5 (25) GPa, K S ' = 2.9 (1), G 0 = 116.5 (17) GPa, and G' = 1.0 (1). Modeling of the velocity profiles of subducted sediments suggests that at 8 GPa (~240 km depth) the P- and S- wave velocities of the subducted sediments exhibit a first order increase of 10% and 11%, respectively. Furthermore, the comparison between seismic observations and velocity contrasts of sediment + MORB model at 240 km shows dehydration of a 1.5–3.9 km thick sediment layer can potentially contribute to the seismically observed Lehmann discontinuity in Tonga, Japan, Philippine Sea, South America, Indonesia, and South Shetland subduction zones.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Field-angle dependence of sound velocity in the Weyl semimetal TaAs

The elastic modulus c 44 of a single crystal of the Weyl semimetal TaAs was investigated by measuring relative changes in the sound velocity under application of a magnetic field up to 10 T. Using an ultrasonic pulsed-echo technique, we studied the shear response of the crystal when the angle between the sound wave propagation and the magnetic field is changed. We observe a broken tetragonal symmetry at fields above 6 T, an anisotropy that is likely related to a longitudinal negative magnetoresistance and therefore might provide evidence of the chiral anomaly, one of the main topological signatures of this class of materials. Furthermore, we also observe quantum oscillations in the sound velocity whose frequencies vary with magnetic field orientation. A fan diagram of Landau level indices reveals topological and trivial Berry phases, depending on the field orientation, indicating a sensitivity to different Fermi surface pockets that do or do not enclose Weyl nodes, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elasticity of davemaoite as a primary contributor to lower-mantle heterogeneities

Geophysical detection of subducted mid–ocean ridge basalt (MORB) in the lower mantle is hindered by uncertainties in the elasticity of Fe,Al,Mg,Ti–bearing davemaoite, a key MORB component. Using Brillouin spectroscopy and x-ray diffraction, we determined the elasticity of a Ca 0.906(1) Fe 2+ 0.027(1) Fe 3+ 0.042(1) Mg 0.033(1) Al 0.072(1) Ti 0.020(1) Si 0.912(1) O 3 davemaoite up to 113 gigapascals and 2294 K. Here, we found that it exhibited a shear wave velocity 10 to 20% slower than end-member davemaoite, making it the slowest phase among major lower-mantle minerals. Our models show that MORB, containing 20 to 25 volume percent davemaoite, potentially contributes to large low-shear-velocity provinces (LLSVPs), whereas a cumulate layer enriched in davemaoite crystallized from basal magma ocean may comprise ultralow-velocity zones (ULVZs). Davemaoite’s ability to host incompatible and heat-producing elements possibly links LLSVPs and ULVZs to mantle plume initiation and geochemical signatures of ocean island basalts.

Zhou, Wen-Yi [Texas A & M Univ., College Station, ↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2021

Elastic moduli of hexagonal diamond and cubic diamond formed under shock compression

The relative stiffnesses and strengths of hexagonal diamond and cubic diamond constitute long standing unresolved issues, because relevant experimental results are lacking. Laser interferometry was used to experimentally determine the longitudinal sound speeds and moduli in hexagonal diamond and cubic diamond formed during the shock compression of graphite. The hexagonal diamond longitudinal moduli are significantly larger than the cubic diamond longitudinal moduli, and even exceed averaged cubic diamond single crystal values. Here, the measured hexagonal diamond longitudinal moduli, combined with high-pressure bulk moduli for cubic diamond single crystals, show that shock-formed hexagonal diamond shear moduli are larger than the shear moduli for cubic diamond single crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗
DOE OSTI2023

Atomistic insight into the ferroelastic post-stishovite transition by high-pressure single-crystal X-ray diffraction

Abstract The post-stishovite transition is a classic pseudo-proper typed ferroelastic transition with a symmetry-breaking spontaneous strain. This transition has been studied using high-pressure spontaneous strains, optic modes, and elastic moduli (Cij) based on the Landau modeling, but its atomistic information and structural distortion remain poorly understood. Here we have conducted synchrotron single-crystal X-ray diffraction measurements on stishovite crystals up to 75.3 GPa in a diamond-anvil cell. Analysis of the data reveals atomic positions, bond lengths, bond angles, and variations of SiO6 octahedra across the transition at high pressure. Our results show that the O coordinates split at ~51.4 GPa, where the apical and equatorial Si-O bond lengths cross over, the SiO6 octahedral distortion vanishes, and the SiO6 octahedra start to rotate about the c axis. Moreover, distortion mode analysis shows that an in-plane stretching distortion (GM1+ mode) occurs in the stishovite structure at high pressure while a rotational distortion (GM2+ mode) becomes dominant in the post-stishovite structure. These results are used to correlate with elastic moduli and Landau parameters (symmetry-breaking strain e1–e2 and order parameter Q) to provide atomistic insight into the ferroelastic transition. When the bond lengths of two Si-O bonds are equal due to the contribution from the GM1+ stretching mode, C11 converges with C12, and the shear wave VS1[110] polarizing along [110] and propagating along [110] vanishes. Values of e1–e2 and Q are proportional to the SiO6 rotation angle from the occurrence of the GM1+ rotational mode in the post-stishovite structure. Our results on the pseudo-proper type transition are also compared with that for the proper type in albite and improper type in CaSiO3 perovskite. The symmetry-breaking strain, in all these types of transitions, arises as the primary effect from the structural angle (such as SiO6 rotation or lattice constant angle) and its relevant distortion mode in the low-symmetry ferroelastic phase.

Geochemistry & Geophysics↗