Search NASA⌕ Search

SEARCH · Search NASA

Results for “elastic wave equations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

In Situ Determination of Thermoelastic Properties of Magnesite at High Pressure and Temperature With Implications to Seismic Detectability of Moderately Carbonated Lithologies in the Earth's Mantle

Abstract Subduction of carbonate‐bearing oceanic plates into Earth's interior recycles carbon from the surface to the deep mantle. The subducted carbonates can significantly affect mantle properties and dynamics. Magnesite is recognized as one of the major potential carbon hosts in the deep mantle because of its stability up to deep lower mantle conditions. However, despite many previous studies on the equation of state and elastic properties of magnesite, large discrepancies still exist for its elastic moduli and their pressure and temperature derivatives. Here we report in situ density and elastic wave velocity measurements on a natural magnesite at simultaneous high pressure‐temperature conditions up to ∼8 GPa–1073 K in a multi‐anvil apparatus using ultrasonic and synchrotron X‐ray techniques. Global fitting of the data set to finite strain equations yieldsK S0 = 114.0 ± 1.2 GPa, = 4.9 ± 0.3, = (−0.019 ± 0.002) GPa/K,G 0 = 68.6 ± 0.4 GPa,G′ = 1.6 ± 0.1, and = (−0.018 ± 0.001) GPa/K. Compared to other major upper mantle and transition zone minerals, magnesite has intermediateV P , the lowestV S , and the lowest density. Thus, magnesite possesses higherV P /V S ratio than other mantle minerals in normal mantle regions, whereas this feature is less pronounced in subduction zone environments. Modeling of the velocity profiles of carbonated lithologies along different geotherms suggests that moderately‐enriched magnesite domains are unlikely to be detected seismically in the Earth's mantle.

Geochemistry & Geophysics↗
DOE OSTI2023

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
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 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 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

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elastic constants in monocrystalline tungsten under quasi-hydrostatic pressures to 11.3 GPa

Compressional (P) and shear (S) wave velocities of tungsten single crystals along the [100] and [110] directions were measured using ultrasonic interferometry at room temperature up to 11.3 GPa. Least-squares fitting of $V$$^{[100]}_{P}$, $V$$^{[100]}_{S}$, $V$$^{[110]}_{P}$, and pressure to finite strain (FS) equations yields the elastic constants: C 11 = 523.5(5) GPa, C 12 = 205.1(32) GPa, and C 44 = 160.8(4) GPa, along with their respective pressure derivatives: $C$$^{′}_{11}$ = 6.20(2), $C$$^{′}_{12}$ = 3.35(1), and $C$$^{′}_{44}$ = 1.65(6). Using the Voigt–Reuss–Hill approximation, the elastic moduli were derived as K S0 = 311.2(22) GPa and G 0 = 160.2(7)GPa, along with their respective pressure derivatives: $K$$^{′}_{S0}$ = 4.30(4) and $G$$^{′}_{0}$ = 1.56(1). The Debye temperature of tungsten was determined to be 380.7(8) K, showing good agreement with previous calorimetric measurements. The elastic anisotropy increases slightly from 1.01(1) at ambient pressure to 1.03(1) at 11.3 GPa, suggesting that tungsten remains nearly isotropic under compression. Poisson's ratio slightly increased from 0.281(3) to 0.288(3) with pressure. Additionally, Pugh's ratio decreased from 0.512(4) to 0.494(4), while Pettifor's ratio increased from 0.148(10) to 0.175(10) as pressure reached 11.3 GPa. These results suggest that tungsten is weakly ductile compared to the critical thresholds for ductile behavior (i.e., Pugh's ratio <0.6 and Pettifor's ratio >0, indicating ductile behavior), but that its ductility increases marginally under pressure.

Wang, Ran [Stony Brook University, NY (United Stat↗
DOE OSTI2022

Shock Hugoniot of Forged and Additively Manufactured 304L Stainless Steel

The purpose of this research was to measure the equation of state for additively manufactured (AM) and forged 304L stainless steel using a novel experimental technique. An understanding of the dynamic behavior of AM metals is integral to their timely adoption into various applications. The Hugoniot of the AM 304L was compared to that of the forged 304L at particle velocities where the material retains a two-wave structure. This comparison enabled us to determine the sensitivity of the equation of state to microstructure as varied due to processing. Our results showed that there was a measurable difference in the measured shock velocity between the AM and forged 304L. The shock wave velocities for the AM 304L were found to be ~3% slower than those for the forged 304L at similar particle velocities. To understand these differences, properties such as densities, sound speeds, and texture were measured and compared between the forged and AM materials. Our results showed that no measurable difference was found in these properties. Additionally, it is possible that differing elastic wave amplitudes may influence shock velocity.

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 OSTI2023

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
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 OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Single-crystal elasticity of (Al,Fe)-bearing bridgmanite up to 82 GPa

Abstract Thermoelastic properties of mantle candidate minerals are essential to our understanding of geophysical phenomena, geochemistry, and geodynamic evolutions of the silicate Earth. However, the lower-mantle mineralogy remains much debated due to the lack of single-crystal elastic moduli (Cij) and aggregate sound velocities of (Al,Fe)-bearing bridgmanite, the most abundant mineral of the planet, at the lower mantle pressure-temperature (P-T) conditions. Here we report single-crystal Cij of (Al,Fe)-bearing bridgmanite, Mg0.88Fe0.1Al0.14Si0.90O3 (Fe10-Al14-Bgm) with Fe3+/ΣFe = ~0.65, up to ~82 GPa using X-ray diffraction (XRD), Brillouin light scattering (BLS), and impulsive stimulated light scattering (ISLS) measurements in diamond-anvil cells (DACs). Two crystal platelets with orientations of (–0.50, 0.05, –0.86) and (0.65, –0.59, 0.48), that are sensitive to deriving all nine Cij, are used for compressional and shear wave velocity (νP and νS) measurements as a function of azimuthal angles over 200° at each experimental pressure. Our results show that all Cij of singe-crystal Fe10-Al14-Bgm increase monotonically with pressure with small uncertainties of 1–2% (±1σ), except C55 and C23, which have uncertainties of 3–4%. Using the third-order Eulerian finite-strain equations to model the elasticity data yields the aggregate adiabatic bulk and shear moduli and respective pressure derivatives at the reference pressure of 25 GPa: KS = 326 ± 4 GPa, µ = 211 ± 2 GPa, KS′ = 3.32 ± 0.04, and µ′ = 1.66 ± 0.02 GPa. The high-pressure aggregate νS and νP of Fe10-Al14-Bgm are 2.6–3.5% and 3.1–4.7% lower than those of MgSiO3 bridgmanite end-member, respectively. These data are used with literature reports on bridgmanite with different Fe and Al contents to quantitatively evaluate pressure and compositional effects on their elastic properties. Comparing with one-dimensional seismic profiles, our modeled velocity profiles of major lower-mantle mineral assemblages at relevant P-T suggest that the lower mantle could likely consist of about 89 vol% (Al,Fe)-bearing bridgmanite. After considering uncertainties, our best-fit model is still indistinguishable from pyrolitic or chondritic models.

Geochemistry & Geophysics↗
DOE OSTI2022

Physics-informed neural networks for identification of material properties using standing waves

A metallic structure in its initial stage of failure involves plastic deformation or environmental degradation that changes the elastic modulus and density. This work presents the detection of change in wave velocity (a function of elastic modulus and density) as a system identification problem. A physics-informed neural network (PINN) is proposed to solve the system identification problem. The PINN takes the spatial coordinates of scanning locations and time as inputs and provides the displacement and wave velocity as outputs. The governing partial differential equation of standing waves in a rod is incorporated into the neural network as physics in the form of a loss function. The wave velocity vector is randomly initiated. During the training of the network, physics is used to determine and update the wave velocity target vector from the network’s displacement predictions. The measured data, comprising sparse displacement response on the rod structure, are used to train the PINN. The wave velocity at the sparse locations on the rod is learned from the predicted displacements during the training. Using the predictions of the trained network, the response of free vibration or material property variation can be reconstructed at unscanned locations on the structure to obtain high-resolution maps for full-field imaging to detect and localize the changes caused by plastic deformation. The PINN’s sparse scanning and simultaneous prediction capability during training can lead to high scanning and data-processing speeds. This capability yields a nondestructive evaluation system that can predict the presence of degraded material locations as the structural vibrations are scanned and processed in real time.

Rathod, Vivek↗
DOE OSTI2020

Impact of nuclear reactor radiation on the performance of AlN/sapphire surface acoustic wave devices

The performance of an AlN/sapphire surface acoustic wave (SAW) delay line device was characterized in real time under irradiation inside a nuclear reactor. Both its resonant frequency and transmission efficiency were observed to respond to a change in reactor power. The response follows an exponentially saturating behavior after a step power increase, followed by an exponentially decaying recovery after reactor shutdown. A sensitivity analysis based on the governing electro-mechanical equations shows that the frequency shift can be attributed to the softening of sapphire’s elastic constants under neutron radiation. Additionally, a kinetic rate equation is adopted to interpret device response and describe its microstructural evolution. These results suggest that the AlN/sapphire SAW device remains functional under irradiation, is sensitive to neutron and gamma ray fluxes, and offers an opportunity for remote sensing and in-situ measurement of material properties when exposed to nuclear reactor environment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗