Wave Dispersion Research Articles

Study of P wave dispersion in patients with paroxysmal atrial fibrillation and its role in prediction of atrial fibrillation recurrence

BackgroundIt has been known that increased P wave duration and P wave dispersion reflect prolongation of intra-atrial and interatrial conduction time and the inhomogeneous propagation of sinus impulses, which are well-known electrophysiologic characteristics in patients with atrial arrhythmias and especially paroxysmal atrial fibrillation. The objective of this study was assessment of P wave dispersion value in cases with paroxysmal atrial fibrillation and its role in predicting recurrence.ResultsForty-eight patients with documented paroxysmal AF were subjected to clinical evaluation, electrocardiogram and routine Doppler echocardiogram. We found that a statistically significant association was detected between P wave dispersion and older age, diabetic and hypertensive cases with positive correlation also detected with left atrial dimension (LAD), left ventricle size and diastolic dysfunction grade. Mean corrected P wave dispersion and corrected QT interval were higher among cases using sotalol, ca channel blockers, among cases using nitrates and among cases with Morris index > 0.04. Higher mean value of corrected QT was associated with biphasic P v1 shape. Old age, female sex, P wave dispersion and QT wave dispersion are statistically significant predictors of PAF recurrence.ConclusionP wave dispersion in patients with paroxysmal atrial fibrillation was strongly correlated to older age, diabetic and hypertensive patients and also with left atrial dimension (LAD), left ventricle size and diastolic dysfunction grade. Also, mean corrected P wave dispersion can predict atrial fibrillation recurrence in patients with Morris index > 0.04, old age, female sex, and QT wave dispersion.

Effect of nonlocality on the dispersion relations of mechanical metamaterials

The advent of mechanical metamaterials, with their unconventional and programmable properties, has revolutionised the field of engineering across length scales. Their design remains an area of interest to researchers requiring careful analysis of unit cells so as to achieve desired mechanical and wave dispersion behaviour. While metamaterials are known to extract their unique properties from the internal geometry of the unit cells, the contributions of the constituent materials are typically accounted for via their local material properties in the form of Young’s modulus and Poisson’s ratio. In this article, we demonstrate that material nonlocality of the constituents can also play a crucial role in modulating the nontrivial properties of metamaterials. Such nonlocality can arise from various sources, such as finer scale inhomogeneity like defects, inclusions etc. Towards this, we propose a nonlocal continuum model by upscaling a discrete potential through a stochastic gradient estimator (SGE) interpreted herein as a generalised Cauchy–Born rule. The nonlocal parameters can also be estimated from experimental data by solving an inverse problem, similar to the estimation of constitutive parameters in classical continuum mechanics. The proposed model accurately captures the dispersion behaviour at the continuum level as compared to its discrete counterpart. We further incorporate the model within the transfer matrix method to derive the dispersion relations of a metabar as a function of the constituent material nonlocality. We envisage that the present study introduces a new handle in the form of constituent material nonlocality for designing architected metamaterials with unprecedented mechanical functionalities.

A numerical investigation of longitudinal wave for thermo-acoustoelastic effects on fluid-embedded in two viscoelastic layer plates at higher temperatures

Ultrasound has become indispensable for inspecting complex structures due to its efficient propagation in both liquids and solid materials. Extensive research has been conducted across various domains, encompassing wave behavior, structural interactions, material dispersion due to heterogeneity, propagation characteristics near boundaries, defect detection, and signal processing. In this context, longitudinal waves are particularly favored for examining thick components exhibiting surface defects or thermal gradients across their depth. In our current work, we employed a numerical approach based on the Transfer Matrix Method (TMM) to model ultrasound signals backscattered under normal incidence by a multilayer structure immersed in water. This method involves representing each layer using a quadrupole formalism that combines stresses and velocities. Simultaneously, we harnessed the thermodynamic equation for aqueous solutions and Fourier’s law to derive the expression for the thermo-acoustic coefficient of longitudinal waves within a fluid embedded between two parallel immersed plates, enabling precise control of the fluid temperature over a range from 0 °C to 700 °C. In the time domain, we evaluated the ultrasonic parameters of a fluid exposed to a temperature range from 0 to 700 °C. This evaluation was supported by a variety of multilayer structures designed to ensure the accurate processing of backscattered signals. These structures were excited by longitudinal waves, with a central frequency fixed at 5 MHz. The results obtained emphasize the sensitivity of longitudinal waves to temperature variations and reveal the presence of wave dispersion induced by temperature fluctuations. In order to investigate this dispersion, we conducted an extensive study of ultrasonic longitudinal wave propagation across a range of excitation frequencies from 2 to 12 MHz. This study encompassed the calculation of phase and group velocities, as well as the corresponding attenuation. The results confirm the relevance of our research to various applications and highlight the complex interaction between longitudinal wave dispersion and temperature variations.

Investigation of phonons and magnons in [Ni80Fe20/Au/Co/Au] N multilayers

The interaction between phonons and magnons is a rapidly developing area of research, particularly in the field of acoustic spintronics. To discuss this interaction, it is necessary to observe two different waves (acoustic and spin waves) with the same frequency and wavelength. In the Ni80Fe20/Au/Co/Au system deposited on a silicon substrate, we observe the interaction between spin waves and surface acoustic waves using Brillouin light scattering spectroscopy. As a result, we can selectively control (activate or deactivate) the magnetoelastic interaction between the fundamental spin wave mode and surface acoustic waves. This is achieved by adjusting the magnetostrictive layer thickness in the multilayer. We demonstrate that by adjusting the number of layers in a multilayer structure, it is possible to precisely control the dispersion of surface acoustic waves while having minimal impact on the fundamental spin wave mode.

Two Effective Degrees of Freedom Can Represent the Dominant Features of Global Rayleigh Wave Dispersion Maps

AbstractObjectively exploring global variations in crust and upper mantle structure helps constrain fundamental aspects of Earth's plate tectonic and convective processes. Here we adopted a Variational Auto‐Encoder to explore the degrees of freedom of global Rayleigh wave dispersion maps at 4–40 mHz. We found that two latent variables sufficiently represent the global variations, suggesting inherent coupling between crustal and mantle seismic properties. We propose that the two extracted latent variables mostly correspond with crustal thickness and upper mantle thermal structure. The first variable shows low values for continental mountain belts and ocean spreading ridges, contrasted by high values for abyssal plains. The second variable shows low values for most oceanic lithosphere and Phanerozoic continental areas, contrasted by high values for Archean cratons. Latent space correlations indicate that continental lithosphere has more strongly coupled depth features than beneath the oceans, which might be a consequence of its longevity.

A generalized peridynamic material correspondence formulation using non-spherical influence functions

Peridynamics theory is a nonlocal continuum mechanics theory. The peridynamic material correspondence formulation provides passage for direct incorporation of any material constitutive models from the classical local continuum mechanics theory. The performance of the material correspondence formulation heavily relies on the influence function whose currently used form is spherical and depends only on the bond length. Together with the way how the nonlocal deformation gradient is constructed based on the material point horizon, this results in the well-known issue of material instability in the conventional peridynamic material correspondence formulation. In this paper, a generalized peridynamic material correspondence formulation that enables use both spherical and non-spherical influence functions is developed. A recently proposed class of parameterized non-spherical influence functions that depend on both the bond length and bond direction is used in the generalized formulation. Numerical examples including wave dispersion, small and finite deformation and linear elastic fracture analyses are studied to check the material stability, convergence characteristic, prediction accuracy and applicability to fracture problems without singularity issue of the proposed formulation. It is found that the proposed formulation is inherently stable, possesses linear convergence rate, yields highly accurate predictions for both small and finite deformation problems, and small horizon can be used to improve computational efficiency.

Forward modeling of quake’s infrasound recorded in the stratosphere on board balloon platforms

Acoustic waves generated by seismic waves contain information on the internal structure of planets, and can be sensed by pressure sensors onboard high-altitude balloons. To identify the various contributions (infrasound signal, noise, balloon response, etc.) in such pressure records, a full waveform modeling is implemented and completed by infrasound ray tracing and additional data analysis. Here, we analyze the Stratéole-2 pressure data associated with two earthquakes (Garcia et al. Geophys. Res. Lett. 49(15):e98844, 2022) and compared these to full waveform simulations by SPECFEM2D-DG-LNS software. Even if our simulations do not precisely reproduce the waveform observed in the frequency range [0.05, 0.3] Hz, we show that the waveform presents more sensitivity to quake and internal structure parameters than to atmospheric structure, and that seismic surface wave dispersion is observed in balloon pressure records. The long-duration pressure oscillations observed after the main infrasonic signal cannot be fully reproduced by our one-dimensional input model even when source time function complexity and aftershocks are considered. These features are ascribed mainly to the complex vertical ground movements below the balloon and partly to late secondary infrasound arrivals excited by the interactions of seismic waves with the topography. These results enhance the advantages and limitations of quake-related infrasound observations on board terrestrial and planetary balloon platforms.Graphical

The single-phase solution and Whitham modulation equations of the defocusing self-induced transparency system

Under investigation in this paper is the defocusing self-induced transparency system which can describe propagation characters of short pulses in Erbium doped fibers with a two-level medium. The single-phase periodic solution and Whitham modulation equations will be derived via finite-gap integration method and Whitham modulation theory. In addition, the oscillatory structure of dispersive shock waves will be constructed and analyzed.

Controllable multi-soliton spectral tunneling by airy pulses in photonic crystal fiber

A straightforward method to induce tunable multi-soliton spectral tunneling (SST) using a finite energy Airy pulse in a photonic crystal fiber (PCF) featuring three zero-dispersion wavelengths is proposed numerically in the paper. The study reveals that the distinctive multi-envelope structure of the Airy pulse enables the control of tunneling distances, properties, and soliton numbers through modulation of its attenuation factor and peak intensity. Notably, the SST process involving the Airy pulse generates multiple blue-shifted and red-shifted dispersion waves, along with associated trapping waves, contributing to the formation of an ultra-flat and ultra-wide supercontinuum. Furthermore, the velocity and efficiency of multi-SST can be manipulated to enhance the overall tunability of the phenomenon. The findings underscore the potential application of Airy pulse-driven SST in secure multi-frequency signal transmission and controlled supercontinuum generation, with the tunability of these processes being governed by the specific attenuation factor and peak power settings of the input Airy pulse. This work advances the understanding of multi-SST dynamics in PCFs and provides a theoretical framework for exploiting these dynamics in practical information security and optical spectrum management contexts.

Computation and analysis of surface wave dispersion and attenuation in layered viscoelastic-vertical transversely isotropic media by the generalized R/T coefficient method

SUMMARY In this study, we propose a systematic and effective method, that is, an extended version of the generalized reflection/transmission (R/T) coefficient method, for computing the phase-velocity (${c}_r$) dispersion curves, attenuation coefficient ($\alpha $) curves, and eigenfunctions of both Rayleigh and Love waves as well as the ellipticity of Rayleigh waves in layered viscoelastic-vertical transversely isotropic (VTI) media. The numerical scheme of combining the root-searching method with the local optimization method is designed for determining the complex-valued modal solutions (i.e. complex wavenumber $k = {\omega {/ {\vphantom {\omega {{c}_r - i\alpha }}}} {{c}_r - i\alpha }}$) of surface waves. The near-surface sedimentary geological environment is taken as the model example because it is typical viscoelastic-VTI media. Besides the anisotropic-viscoelastic (AV) media, our algorithm can also compute surface waves in isotropic-elastic (IE), isotropic-viscoelastic (IV) and anisotropic-elastic (AE) media by resetting the corresponding parameters. Using the six-layer half-space models and in these four media, we verify the correctness of our algorithm by benchmarking the modal solutions against those from other methods. In the four-layer half-space model, by comparing the results of IE, IV, AE and AV media, we analyse the effects of velocity anisotropy, viscoelasticity and attenuation anisotropy on the dispersion and attenuation characteristics of both Rayleigh and Love waves in detail. Our study can provide a theoretical basis and useful tool for surface wave imaging considering the anisotropy and/or viscoelasticity of the medium, which has the potential to better investigate the solid Earth's internal structure.

Theoretical study on the effects of Mn ion doping and applied magnetic field in (In,Mn)As

Diluted magnetic semiconductors are a recent research area due to their ability to enhance ferromagnetic properties and facilitate the electrical detection of magnetoresistance and polarization. (In, Mn)As dilute magnetic semiconductor has potential application in the field of spintronic devices, such as spin field-effective transistors, spin laser light-emitted diodes, modern technology, multi-functional devices, green technology, and nanotechnology. For this study, we have considered the RKKY interaction between Mn2+ spins via delocalized carriers. The effect of manganese ion concentration and applied magnetic field on ferromagnetic diluted magnetic semiconductor properties such as dispersion relation, Curie temperature, and reduced magnetization of (In, Mn)As are studied. We have developed a spin-wave model using the Holstein-Primakaff transformation. Based on the developed model, the number of ferromagnetic magnons, dispersion relations, and Curie temperatures were calculated with or without an applied magnetic field. The reduced magnetization is also calculated. The graph of Curie temperature and magnetization of (In,Mn)As versus temperature with applied field up to 6 T and manganese ion concentration from 0.01 to 0.1 are plotted. The graph of spin wave dispersion of (In,Mn)As versus a wave vector with varying manganese ion concentration with and without an applied magnetic field up to 6 T. In this study, an InMnAs Curie temperature of 290.68 K is found without an applied magnetic field with a 0.1 manganese ion concentration, which is near room temperature. Moreover, with an applied magnetic field of 6 T at 0.1 manganese ion concentration, a Curie temperature of 342.466 K is found, which is above room temperature. Hence, these temperatures are suitable for the field of next-generation spintronic new technology.

Some gradient theories in linear visco-elastodynamics towards dispersion and attenuation of waves in relation to large-strain models

Some gradient theories in linear visco-elastodynamics towards dispersion and attenuation of waves in relation to large-strain models

Integrable approximations of dispersive shock waves of the granular chain

In the present work we revisit the shock wave dynamics in a granular chain with precompression. By approximating the model by an α-Fermi–Pasta–Ulam–Tsingou chain, we leverage the connection of the latter in the strain variable formulation to two separate integrable models, one continuum, namely the KdV equation, and one discrete, namely the Toda lattice. We bring to bear the Whitham modulation theory analysis of such integrable systems and the analytical approximation of their dispersive shock waves in order to provide, through the lens of the reductive connection to the granular crystal, an approximation to the shock wave of the granular problem. A detailed numerical comparison of the original granular chain and its approximate integrable-system-based dispersive shocks proves very favorable in a wide parametric range. The gradual deviations between (approximate) theory and numerical computation, as amplitude parameters of the solution increase are quantified and discussed.

Inversion of shear wave slowness and attenuation from dipole acoustic logging data based on an equivalent tool theory

The presence of a sonic logging tool affects the propagation characteristics of dipole-flexural waves in a fluid-filled borehole. As a result, the tool effect should be considered during the processing and inversion of the dipole-flexural waves for the acoustic properties of a shear wave (S wave). We show that the presence of the tool not only changes the slowness dispersion of the flexural wave but also results in a shift of its frequency-dependent attenuation toward lower frequencies. The equivalent tool theory (ETT) can realistically model the effect of an acoustic tool on guided wave propagation. With the ETT framework, we develop a unified workflow for estimating the S-wave slowness and attenuation from dipole array acoustic waveforms. The workflow is an extension of the model-based dispersive processing of dipole logging data. Our workflow involves two inversion procedures. By matching the actual dispersion data with the theoretical dispersion predicted by the ETT, we can invert the S-wave slowness and tool parameters. With the estimated S-wave slowness and other borehole parameters, we further calculate the partition coefficients of flexural waves using the ETT. Based on the concept of energy partitioning, we formulate a linear optimization problem and invert for the S-wave attenuation from the frequency-dependent dipole attenuation. We use synthetic and field data to demonstrate the effectiveness and applicability of the developed method. The results show that it is important to incorporate the tool effect when processing dipole acoustic logging data for S-wave slowness and attenuation, especially in a fast formation.

Research and application of Rayleigh wave imaging based on the Born–Jordan time‐frequency distribution

AbstractCurrently, the horizontal resolution of Rayleigh wave exploration is low. In this study, we propose the Born–Jordan time‐frequency distribution to analyse Rayleigh waves. The seismic signal was filtered with a wavelet transform for denoising, and the Rayleigh wave was separated in the time domain. Using the Born–Jordan time‐frequency distribution, the time waveform of each frequency comprising the Rayleigh wave from every seismic channel was obtained, and the time difference of the Rayleigh wave with the same frequency was calculated, based on which the dispersion curve between the two channels was obtained. Combined with the multichannel Rayleigh wave dispersion curve, phase velocity and frequency imaging under the seismic arrangement were obtained. Applying this method to detect abnormal geological bodies in engineering investigations showed that hard geologic bodies, such as comcrete rocks, have high velocity and frequency, whereas weak ones have low velocity and frequency. This strategy facilitated the detection of fractured zones, underground goafs and obstacles during pipe‐jacking construction near the surface.

Reservoir Mud Releasing May Suboptimize Fluvial Sand Supply to Coastal Sediment Budget: Modeling the Impact of Shihmen Reservoir Case on Tamsui River Estuary

AbstractRegular release of sediment from reservoir has been increasingly adopted as a strategy for sustainable management. Here, we use a process‐based morphodynamic model to simulate the estuarine sediment dynamics impacted by turbidity current venting implemented by the Shihmen Reservoir during three typhoon events in 2008. Upon validation with the post‐event bathymetries, the model hindcasts reveal that mud releasing can be effective in mitigating reservoir siltation, yet may be a suboptimal strategy for alleviating coastal sediment deficit. A vast majority of the released muds were delivered through the estuary and exported to offshore by flood advection, wave dispersion, and tidal flushing. The flood‐driven sands, sourced mainly from downstream tributaries, were instead the major contributor to coastal sediment budget. However, mud mantling (covering and immobilizing sand deposits by the reservoir‐released muds) reduced sand availability and thus sand delivery to the coast. For the present case, 25% of the released muds were deposited along the way, presence of these mud covers reduced sand delivery by 15%, compared to a hypothetical scenario of clear‐water flood releases. The relative sand transport deficit is found to increase linearly with the degree of bed mud saturation, 1–D/R, with D/R the ratio of single‐event mud deposit to release. Given broad relevance to global reservoirs encountering the problems of siltation and coastal sediment deficit, our findings highlight that sustainable management needs to look beyond just a bulk amount of sediment, but it is critical to consider how different sediment fractions are interacting and impacted by human activities.

Adjoint Synthesis for Trans‐Oceanic Tsunami Waveforms and Simultaneous Inversion of Fault Geometry and Slip Distribution

AbstractTsunamis propagate over long distances and can cause widespread damage even after crossing ocean basins. Prediction of tsunamis in distant areas based on observations near their sources is critical to mitigating damage. In recent years, the accuracy of numerical models of trans‐oceanic tsunami propagation has improved significantly due to the incorporation of effects such as the solid earth response to tsunami loading and wave dispersion. However, these models are computationally expensive and have not been fully utilized for real‐time prediction. Here, we derive the adjoint operator for the linear set of equations describing deep‐ocean tsunami propagation and show how a pre‐computed database of adjoint states can achieve rapid synthesis of tsunami waveforms at target sites from nonpoint arbitrary tsunami sources. The adjoint synthesis method allows for an exhaustive parameter search for tsunami source estimation. A method for simultaneous inversion of fault geometry and slip distribution using adjoint synthesis with Sequential Monte Carlo method was proposed and applied to the 2012 Haida Gwaii earthquake tsunami. The influence of model accuracy and the amount of observed data on the estimation of tsunami sources and waveforms was examined. It was found that with a highly accurate propagation model, using only a limited amount of observed data produced source and waveform estimates very similar to the final models obtained with much larger data sets. The final inferred fault model involved megathrust slip distributed between the Haida Gwaii trench and the Queen Charlotte fault. The proposed method can also quantify the uncertainty of the waveform forecasts.

Lithospheric Evolution of the South‐Central United States Constrained by Joint Inversion of Receiver Functions and Surface Wave Dispersion

AbstractIn the present study, we use broadband seismic data recorded by 190 stations of the EarthScope program's Transportable Array to construct a 3‐D shear wave velocity model for the upper 180 km using a non‐linear Bayesian Monte‐Carlo joint inversion of receiver functions (RFs) and Rayleigh wave dispersion curves. Ambient noise and teleseismic data are used for obtaining Rayleigh wave phase velocity dispersion curves. A resonance removal filtering technique is applied to the RFs contaminated by reverberations from the thick sedimentary layers that cover most of the region. Our observations of the higher crustal shear velocities (∼3.40 km/s) beneath the Sabine Block (SB), along with the estimated relatively thicker crust (∼34.0 km) and lower crustal Vp/Vs estimates (∼1.80) in comparison with the rest of the Gulf Coastal Plain (GCP) (∼3.10 km/s for crustal shear velocities, ∼29.0 km for crustal thickness, and ∼1.90 for crustal Vp/Vs estimates), indicating that this crustal block has different crustal properties from the surrounding coastal plain regions. The southern Ouachita Mountains have a thin crust (∼30.0 km) and low mean crustal Vp/Vs value (∼1.73), suggesting that lower crustal delamination has occurred in this region. Low velocities in the upper mantle beneath most of the GCP are interpreted as a combined result of thin lithosphere, higher‐than‐normal temperatures, and possibly compositional variations.

Numerical Investigation of Evaluating the Degree of Damage in Concrete Plate Cold Joints through Surface Wave Dispersion Velocity Profiles

This study aims to develop a technology for assessing the interfacial condition of cold joints in concrete plates, focusing on the dispersion velocity profile of surface waves. The methodology employs a small impact hammer embedded with a piezoelectric element and a displacement receiver strategically positioned on either side of the seam. In this initial investigation, numerical models were constructed using ANSYS LS-DYNA. The cold joint zone was simulated by randomly removing elements with a predetermined void ratio. Displacement waveforms at the sensor node underwent analysis using a short-time Fourier Transform and the reassigned method to derive the dispersive velocity profile of surface waves. The study delves into the impacts of varying test-line lengths, impact durations, and void distributions. Numerical findings indicate that cold joints are detectable with a void ratio of 40% or greater. For the 60% and 80% void ratios, surface wave velocities register notably lower for wavelengths below 0.1 m and 0.2 m, respectively. Under the same porosity ratio for surface-observable cold joints, if their extension depth is smaller, the wave velocity on the dispersion curve will be lower than that of a fully penetrated cold joint. The current study reveals that a contact duration of 60 μs and an impactor-receiver distance of 0.82 m yield the most distinguishable results for detecting the extension and porosity ratio of a cold joint. Accompanied by tests with a hammer with a contact duration of 20 μs, additional information about the joint can be unveiled. Finally, various void distributions in the cold joint predominantly influence the wave velocity at shorter wavelengths.

Disparate Far‐Field Responses of Deep Pacific‐Plate Subduction Beneath Northeastern Asia: Implication for the North–South Segmented Crustal Modification of the Great Khingan Range

AbstractThe far limit of a plate subduction process and its related far‐field dynamic process are fundamental topics for plate tectonics. The Great Khingan Range (GKAR) in the western flank of NE China is currently under the far‐field influence of the Pacific plate subduction. Benefiting from the newly deployed seismic arrays in the northern GKAR, we take full advantage of seismic data from both temporary and permanent networks and employ an improved scheme of the H‐κ stacking method by introducing surface wave dispersion to obtain the integrated maps of Moho depth and crustal bulk Vp/Vs ratio in this region. Our results show that the Moho depth has a giant step near the North–South Gravity Lineament (NSGL). Meanwhile, the distribution of bulk Vp/Vs ratio presents a north–south difference roughly by 50°N, where the south subregion consists of a collage of high and low Vp/Vs ratio; by contrast, the north subregion is characterized by unified high values. The east‐west structural discrepancies across the NSGL from the Earth's surface to the mantle transition zone indicate the difference in the strength of modifications between the near and far‐fields from the Pacific plate subduction. The complicated distribution of the Vp/Vs ratio in the south subregion supports that secondary small‐scale upwellings underneath the Cenozoic volcanic groups dominate the tectonic reworking in this area. The lack of Cenozoic volcanism and a more straightforward distribution of the Vp/Vs ratio in the north subregion might allude to a tectonically inactive part of NE China.