Asymmetric elastic wave scattering in elastic plates
Structural Health Monitoring presentation for ASA
SEARCH · Engineering Papers
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.
Structural Health Monitoring presentation for ASA
Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments
Flexomagnetism is the coupling between magnetism and strain gradients and is a technologically relevant phenomenon. We present a theory of elastic wave propagation in a linear elastic flexomagnetic material with microstructure and strain gradient elastic interactions. The expressions of frequency, phase velocity, and group velocity of longitudinal and transverse waves are derived and are shown to depend on the flexomagnetic coefficient and microstructure. We also show that the effect of flexomagnetism and microstructure can lead to some interesting phenomena in wave propagation, which are not observed in classical linear elasticity theory of waves. Specifically, in contrast to classical linear elastic materials, where wave propagation is non-dispersive, flexomagnetic materials with microstructure can exhibit both normal and abnormal dispersion. It is also noteworthy that, in flexomagnetic materials with gradient elasticity, the phase velocities of transverse waves can exceed those of longitudinal waves, which is atypical in classical elasticity. Furthermore, waves can also attenuate for a certain range of wavenumbers that depend on the flexomagnetic coefficient and microstructural parameters. Finally, we explore the possibility of waves exhibiting zero group velocity modes, where waves are non-propagating but have strong local energy confinement, negative group velocity modes, where the wave packet moves in the opposite direction to that of wave propagation, and the phenomenon of wave freezing, where a propagating wave stops in space without diffusing or spreading.
Overcoming the low signal-to-noise ratio (SNR) in laser ultrasonic testing of composite materials remains a significant challenge. Current efforts focus on enhancing SNR by inserting more energy into the material through temporal and/or spatial modulation of the laser beam. However, potential SNR improvements through wave convergence and wave energy manipulation have been overlooked. This paper addresses this gap by demonstrating the convergence of different wave types to a designated point and by showing the feasibility of directing absorbed laser energy into a specific wave type through spatial modulation of the laser beam. To achieve this, mathematical expressions for the convergent laser beams are derived. Various laser beam profiles are then introduced to the thermoelastic equations and solved using the finite element method. The sample under investigation is a transversely isotropic unidirectional carbon fiber reinforced plastic, characterized by anisotropic thermal expansion coefficients and thermal conductivities. Results reveal pronounced convergence of the intended wave type at the center due to laser beam shaping. This study showcases the ability to direct absorbed laser energy toward a specific wave type through spatial modulation of the laser beam and highlights the role of material anisotropy in energy focusing.
Guided elastic waves are a truly cross-disciplinary key enabling technology. For more than five decades, surface acoustic wave (SAW) and bulk acoustic wave devices find widespread applications. Nowadays, different types of guided elastic waves cover the wide spectrum of applications spanning from quantum technologies to the life sciences, from controlling single excitations to macroscopic collective states in condensed matter. Six years after the first 2019 SAW roadmap, we believe it is time to make a step back and take a fresh look at the status of the field and its future challenges. Since the first roadmap in 2019, the spectrum clearly expanded and this new edition presents a current snapshot of the status of this vibrant field and prospects for potential future developments.
We report on the scaling and units of the Rayleigh wave elastic response function A R (ω) that is output from the widely used Computer Programs in Seismology (CPS) to aid in the modeling of ground motion sourced by atmospheric explosions. The program uses mixed units (km, second, km/s, gm/cc) to keep A R (ω) near 10 0 and prevent any numerical underflow or overflow. We compare two models for the response of an elastic half space to the output from CPS. Our application inputs the recommended, mixed unit geological models to determine how researchers must scale this output to obtain physical units for A R (ω) that represents the amplitude scaling for the minimum group velocity (Airy phase) contribution to Rayleigh waves. We determine that a CPS user must scale the output for A R (ω) by 10 -12 (m/km) 2 (g/cc/m 3 /kg) to obtain MKS (meter, kg, second) units and then must multiply this result by the vertical component eigenfunction squared, that is evaluated at the free surface.
We introduce pulse-echo ultrasound as a method for characterizing the impact of powder bed fusion parameters on the properties of additively manufactured stainless-steel components, their material anisotropy, and location-dependent heterogeneity. Our results indicate that accurate characterization requires careful selection of ultrasonic propagation paths, which must consider the direction of additive layering, variations in processing parameters, and the component's geometry. We employed two distinct methods to estimate material properties from ultrasonic data: One assumes isotropy, while the other accounts for anisotropic interactions during the propagation of elastic waves. When applied to samples fabricated with laser energy densities ranging from 24 to 42 J/mm³ , these methods revealed transverse isotropy and weak anisotropy (quantified by small Thomsen parameters, ε = 0.0651 and γ = 0.0092) and less than a ∼6 % change in acoustic impedance. The assumption of isotropy, in this case, leads to small errors (less than 4 % or 1 % for Young's modulus in the build or transverse directions) when estimating orthotropic material properties using ultrasonic data measured along just two orthogonal directions, one of which must align with the build direction. By comparing ultrasonic measurements — which aggregate the spatial variability in material properties along the length of elastic wave propagation into a single value — with localized measurements obtained from surface nanoindentation, we uncovered and spatially profiled significant differences between the surface and interior properties. Specifically, the surface Young's modulus decreased from approximately 210 GPa to 180 GPa within a depth of about 3 mm. We attribute this surface-localized heterogeneity in PBF-fabricated components to distinct thermal histories experienced by the surface and interior regions. Collectively, the results of this study establish a framework for the ultrasonic characterization of material heterogeneity and anisotropy in material properties and demonstrate its application in additively manufactured metal components.
Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.
Electronic systems contained within sealed metallic enclosures, such as nuclear waste containers where wire penetration is not an option, require a through-wall data transfer technology to communicate. One such technology utilizes piezoelectric ultrasonic transducers to communicate through the barrier using elastic waves. Currently demonstrated ultrasonic through-wall data communication techniques have relied on having a direct interface with the ultrasonic transducers on either side. Here, this work studies the coupling of an antenna to such an ultrasonic communication system. A computationally efficient model of a piezobarrier ultrasonic communication system and antenna using a simplified antenna model and a transfer matrix piezobarrier model is utilized to perform a representative parameter space search. The system efficiency and bandwidth trends are shown for the system parameters as well as operating frequency. The tradeoff between efficiency and bandwidth is also presented. A system configuration within the parameter space is fabricated and demonstrated to achieve a data rate of up to 100 kb/s.
Abrupt, local frictional fault failure comprises a displacement that is normally accompanied by acoustic emission (AE)—an impulsive elastic wave broadcast with an amplitude proportional to particle velocity. The aggregate of these displacements is the basic fault motion. In laboratory shear experiments, the examination of a sequence of laboratory earthquakes includes continuous measurements of fault motion and the associated AE that is broadcast. From these measurements, connections between the fault motion and cumulative sum of the AE amplitude can be identified. The composition of the AE broadcasts reveals inhomogeneity in the fault mechanical structure from which they arise. This inhomogeneity can be decomposed into a time invariant AE component and an articulated AE component. The articulated AE component serves as a “state of the fault diagnostic” that follows a distinctive pattern to fault failure. Thus, the articulated AE component can be used directly to monitor the state of the fault.
Laser ultrasound (LU) is a technique that uses a pump laser and a probe laser to optically generate and detect elastic waves in a material. Despite its advantages over traditional contact transducer-based ultrasound, industrial adoption has been limited by complex optical setups and the inability of multi-mode fibers to deliver a stable Gaussian profile for the pump laser. Here, we report a fully fiber-coupled thermoelastic LU system that uses an anti-resonant hollow-core single-mode fiber to deliver 1 mJ nanosecond pulses of 1064 nm light, while preserving the fundamental Gaussian mode (pump laser). When combined with a fiber-coupled interferometer (probe laser), a small, flexible, and environmentally robust sensor capable of optically generating and detecting high frequency broadband ultrasound is realized. We demonstrate such an LU system implemented in situ on a four-axis precision lathe. High-resolution thickness gauging is performed, before and after precision cutting, by exciting and measuring a zero-group velocity guided wave mode. The measurements are verified with ex-situ traceable coordinate measuring machine data. Mean absolute deviations of 0.1%, of nominal thickness, before cutting, and 0.2% and 0.3%, after stepped and tapered cuts, respectively, are reported. A theoretical background for thermoelastic ultrasound generation in an elastic waveguide is also presented. Attention is given to the effect of the pump laser profile on wave generation to elucidate the importance of using single-mode laser light. The fiber-coupled system demonstrated is well-suited for use in scientific and engineering sensing applications and facilitates the adoption of LU for industrial non-destructive testing.
SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.
This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.
The recent development of neural operator (NeurOp) learning for solutions to the elastic wave equation shows promising results and provides the basis for fast large-scale simulations for different seismological applications. In this article, we use the Fourier neural operator (FNO) model to directly solve the 3D Helmholtz wave equation for fast seismic ground-motion simulations on different frequencies and show the frequency bias of the FNO model, that is, it learns the lower frequencies better comparing to the higher frequencies. To reduce the frequency bias, we adopt the multistage FNO training, that is, after training a stage 1 FNO model for estimating the ground motion, we use a second FNO model as the stage 2 to learn from the residual, which greatly reduced the errors on the higher frequencies. By adopting this multistage training, the FNO models show reduced biases on higher frequencies, which enhanced the overall results of the ground-motion simulations. Thus the multistage training FNO improves the accuracy and realism of the ground-motion simulations.
Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.
Low-phase-noise microwave-frequency integrated oscillators provide compact solutions for various applications in signal processing, communications, and sensing. Surface acoustic waves (SAWs), featuring orders-of-magnitude shorter wavelength than electromagnetic waves at the same frequency, enable integrated microwave-frequency systems with much smaller footprint on chip. SAW devices also allow higher-quality (Q) factors than electronic components at room temperature. Here, we demonstrate a low-phase-noise gigahertz-frequency SAW oscillator on 128°Y-cut lithium niobate, where the SAW resonator occupies a footprint of 0.05 mm 2 . Leveraging phononic crystal band-gap-edge modes to balance between Q factors and insertion losses, our 1-GHz SAW oscillator features a low phase noise of -132.5 dBc/Hz at a 10-kHz offset frequency and an overlapping Hadamard deviation of 6.5 × 10 -10 at an analysis time of 64 ms. The SAW resonator-based oscillator holds high potential in developing low-noise sensors and acousto-optic integrated circuits.
Temperature variation is often considered an undesirable factor in ultrasonic testing, as wave velocity is highly sensitive to thermal fluctuations. However, completely eliminating the temperature effect is difficult, particularly in tests requiring precise velocity measurements. Recently, a method called Thermal Modulation of Nonlinear Ultrasonics (TMNL) has been developed. Instead of eliminating the thermal effect, the TMNL method leverages the temperature variation as a driving force to stimulate the nonlinear response of the medium and modulate the ultrasonic waves propagating in it. These modulated waves can then be used to evaluate the nonlinear behaviours of the test medium. This paper presents a focused review of the TMNL technique, including its theoretical foundations, particularly the conceptual challenges in integrating thermal effects into classical acoustoelastic theory, and recent applications in non-destructive evaluation (NDE). Three case studies are presented to demonstrate TMNL’s application in detecting microcracking in concrete, assessing ageing in polymer materials, and enabling temperature compensation in acoustoelastic tests. In conclusion, the review also summarises related studies, including photothermal crack modulation, and discusses current limitations and future directions of TMNL in elastic media.
Deep learning has emerged as a transformative tool for enhancing the extraction of reliable information from seismograms, addressing the increasing demand for precise and efficient seismic data analysis. We introduce an innovative encoder–decoder deep learning model, named WaveDenoiser, designed for noise reduction in the time domain, thereby eliminating the need for spectrogram computations that have been used for existing deep learning tools and significantly improving processing speed. Utilizing the benchmark dataset that is Stanford Earthquake Dataset, we developed three models of varying sizes: base, medium, and large. Notably, the large (referred to as WaveDenoiser) model demonstrated superior performance, achieving a median signal‐to‐noise ratio improvement of 8.8 dB on in‐distribution unseen data (in the same geographic region) and 7.7 dB on out‐distribution unseen data (in a new geographic region), outpacing both the base and medium models. Further evaluation of the WaveDenoiser model revealed a reduction in median arrival‐time errors by 0.02 s for P waves and 0.01 s for S waves when processing waveforms prior to phase picking using PhaseNet on in‐distribution unseen data. When tested on out‐distribution unseen data, the model also effectively reduced the P‐wave median arrival‐time error by 0.02 and 0.01 s in median arrival‐time error for S waves. Importantly, the application of WaveDenoiser resulted in a significant reduction of phase picking outliers by 1.1% to 3.6% for both P and S waves. In addition, we achieved over five times acceleration in processing speed compared with the seisBench implementation of DeepDenoiser. Our findings underscore the potential of WaveDenoiser as a powerful tool for improving seismic data analysis and processing efficiency.