An analogical study of wave equations,physical quantities,conservation and reciprocity equations between electromagnetic and elastic waves

This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differential operations.Using the electromagnetic-to-elastic substitution,the analogous relations of the conservation laws of energy and momentum are investigated between these two physical fields.Moreover,the energy-based and momentum-based reciprocity theorems for an elastic wave are also derived in the time-harmonic state,which describe the interaction between two elastic wave systems from the perspectives of energy and momentum,respectively.The theoretical results obtained in this analysis can not only improve our understanding of the similarities of these two linear systems,but also find potential applications in relevant fields such as medical imaging,non-destructive evaluation,acoustic microscopy,seismology and exploratory geophysics.

Reconfigurable mechanism-based metamaterials for ternary-coded elastic wave polarizers and programmable refraction control

Elastic metamaterials with unusual elastic properties offer unprecedented ways to modulate the polarization and propagation of elastic waves.However,most of them rely on the resonant structural components,and thus are frequency-dependent and unchangeable.Here,we present a reconfigurable 2D mechanism-based metamaterial which possesses transformable and frequency-independent elastic properties.Based on the proposed mechanism-based metamaterial,interesting functionalities,such as ternarycoded elastic wave polarizer and programmable refraction,are demonstrated.Particularly,unique ternary-coded polarizers,with 1-trit polarization filtering and 2-trit polarization separating of longitudinal and transverse waves,are first achieved.Then,the strong anisotropy of the proposed metamaterial is harnessed to realize positive-negative bi-refraction,only-positive refraction,and only-negative refraction.Finally,the wave functions with detailed microstructures are numerically verified.

On Klein tunneling of low-frequency elastic waves in hexagonal topological plates

Incident particles in the Klein tunnel phenomenon in quantum mechanics can pass a very high potential barrier.Introducing the concept of tunneling into the analysis of phononic crystals can broaden the application prospects.In this study,the structure of the unit cell is designed,and the low frequency(<1 k Hz)valley locked waveguide is realized through the creation of a phononic crystal plate with a topological phase transition interface.The defect immunity of the topological waveguide is verified,that is,the wave can propagate along the original path in the cases of impurities and disorder.Then,the tunneling phenomenon is introduced into the topological valley-locked waveguide to analyze the wave propagation,and its potential applications(such as signal separators and logic gates)are further explored by designing phononic crystal plates.This research has broad application prospects in information processing and vibration control,and potential applications in other directions are also worth exploring.

Mathematical framework of nonlinear elastic waves propagating in pre-stressed media

Acoustic nonlinearity holds potential as a method for assessing material stress.Analogous to the acoustoelastic effect,where the velocity of elastic waves is influenced by third-order elastic constants,the propagation of nonlinear acoustic waves in pre-stressed materials would be influenced by higher-order elastic constants.Despite this,there has been a notable absence of research exploring this phenomenon.Consequently,this paper aims to establish a theoretical framework for governing the propagation of nonlinear acoustic waves in pre-stressed materials.It delves into the impact of pre-stress on higher-order material parameters,and specifically examines the propagation of one-dimensional acoustic waves within the contexts of the uniaxial stress and the biaxial stress.This paper establishes a theoretical foundation for exploring the application of nonlinear ultrasonic techniques to measure pre-stress in materials.

Evaluation of soil fabric using elastic waves during load-unload

It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using three granular materials in the laboratory.The variations of compression and shear wave velocities(Vp and Vs)are investigated during load-unload cycles under dry and drained conditions.Supplementary discrete element method(DEM)simulations are performed to understand the evolution of soil fabric during the equivalent load-unload cycles using spherical particles.Vp and Vs are not always reversible even though the stress state regains its isotropic condition after unload,indicating that Vp and Vs are governed by not only the stress state but also the fabric change.The variations of Vp/Vs are density-and stress-dependent;a higher level of stress ratio(s01/s03)threshold is observed for denser packings to trigger a significant change in wave velocity ratio(Vp/Vs)for experimental results using spherical glass beads and simulation data using spherical particles.Considering the particle shape,a higher s01/s03 threshold is found for more angular particles than rounded particles.The DEM result reveals that Vp/Vs of spherical particles can be correlated linearly with the evolution of fabric ratio(Fver/Fhor)during loadunload in a pre-peak range under dry and drained conditions.

The mechanism to reform dynamic performance of an elastic wave-front in a piezoelectric semiconductor by the wave-carrier interaction induced from static biasing fields

The propagation of an elastic wave(EW)in a piezoelectric semiconductor(PSC)subjected to static biasing fields is investigated.It is found that there exist two coupling waves between electric field and charge carriers.One is stimulated by the action of the polarized electric field in the EW-front on charge carriers(EFC),and the other is stimulated by the action of initial electric field in biasing fields on dynamic carriers(IEC).Obviously,the latter is a man-made and tunable wave-carrier interaction.A careful study shows that IEC can play a leading role in remaking dynamic performance of the wave-front and an inter-medium role in transferring energy from biasing fields to EW-fronts.Hence,a method is proposed to reform the EW performance by biasing-fields:reforming the dispersivity of EW-fronts by promoting competition between IEC and EFC and inverting the dissipation by the IEC to transfer energy from biasing fields to EWfronts.The corresponding tuning laws on the phase-frequency characteristics of an EW show that the wave velocity can be regulated smaller than the pure EW velocity at a lowfrequency and larger than the pure piezoelectric wave velocity at a high-frequency.As for regulating the amplitude-frequency characteristics of the EW by the IEC,analyses show that EWs can obtain amplification only for those with relatively high vibration frequencies(small wave lengths).The studies will provide guidance for theoretical analysis of waves propagating in PSCs and practical application and design of piezotronic devices.

Investigation of long-wavelength elastic wave propagation through wet bentonite-filled rock joints

The saturation of the compacted bentonite buffer in the deep geological repository can cause bentonite swelling,intrusion into rock fractures,and erosion.Inevitably,erosion and subsequent bentonite mass loss due to groundwater inflow can aggravate the overall integrity of the engineered barrier system.Therefore,the coupled hydro-mechanical interaction between the buffer and rock during groundwater inflow and bentonite intrusion should be evaluated to guarantee the long-term safety of deep geological disposal.This study investigated the effect of bentonite erosion and intrusion on the elastic wave propagation characteristics in jointed rocks using a quasi-static resonant column test.Jointed rock specimens with different joint conditions(i.e.joint surface saturation and bentonite filling)were prepared using granite rock discs sampled from the Korea Underground Research Tunnel(KURT)and Gyeongju bentonite.The long-wavelength longitudinal and shear wave velocities were measured under different normal stress levels.A Hertzian-type power model was used to fit the wave velocities,and the relationship between the two fitted parameters provided the trend of joint conditions.Numerical simulations using three-dimensional distinct element code(3DEC)were conducted to better understand how the long-wavelength wave propagates through wet bentonite-filled rock joints.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Tunable three-dimensional nonreciprocal transmission in a layered nonlinear elastic wave metamaterial by initial stresses

In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.

3D elastic wave equation forward modeling based on the precise integration method

The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid, Finite Difference Scheme

This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

PML and CFS-PML boundary conditions for a mesh-free finite difference solution of the elastic wave equation

Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.

Harten-Lax-van Leer-contact(HLLC)approximation Riemann solver with elastic waves for one-dimensional elastic-plastic problems

A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

Influence of bedding structure on stress-induced elastic wave anisotropy in tight sandstones

To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum principal stress.P-wave velocities along 64 different directions on each specimen were monitored frequently to understand the anisotropy change at various stress levels by fitting Thomsen’s anisotropy equation.The results show that the elastic wave anisotropy is very sensitive to mechanical loading.Under hydrostatic loading,the magnitude of anisotropy is reduced in all three specimens.However,under deviatoric stress loading,the evolution of anisotropic characteristics(magnitude and orientation of the symmetry axis)is bedding orientation dependent.Anisotropy reversal occurs in specimens with bedding normal/oblique to the maximum principal stress.P-wave anisotropyε0 is linearly related to volumetric strain Sv and dilatancy,indicating that stress-induced redistribution of microcracks has a significant effect on P-wave velocity anisotropy.The closure of initial cracks and pores aligned in the bedding direction contributes to the decrease of the anisotropy.However,opening of new cracks,aligned in the maximum principal direction,accounts for the increase of the anisotropy.The experimental results provide some insights into the microstructural behavior under loading and provide an experimental basis for seismic data interpretation and parameter selection in engineering applications.

Propagation behavior of elastic waves in one-dimensional piezoelectric/piezomagnetic phononic crystal with line defects

Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating energy reflection/transmittion coefficients of quasi-pressure and quasi-shear waves. Line defects are created by the re- placement of PE or PM constituent layer. The defect modes existing in the first gap are considered and the influences on defect modes of the material properties and volume fraction of the defect layers, the type of incident waves, the location of defect layer and the number of structural layers are discussed in detail. Numerical results indicate that defect modes are the most obvious when the defect layers are inserted in the middle of the perfect PCs; the types of incidence wave and material properties of the defect layers have important effects on the numbers, the location of frequencies and the peaks of defect modes, and the defect modes are strongly de- pendent on volume fraction of the defect layers. We hope this paper will be found useful for the design of PE/PM acoustic filters or acoustic transducer with PCs structures.

ELASTIC WAVEFIELD CALCULATION FOR HETEROGENEOUS ANISOTROPIC POROUS MEDIA USING THE 3-D IRREGULAR-GRID FINITE-DIFFERENCE

Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous media.The ‘slow’P wave in porous media wave simulation is highly dispersive.Finer grids are needed to get a precise wavefield calculation for models with curved interface and complex geometric structure.Fine grids in a global model greatly increase computation costs of regular grids scheme.Irregular fine or coarse grids in local fields not only cost less computing time than the conventional velocity-stress FDM,but also give a more accu- rate wavefield description.A dispersion analysis of the irregular-grid finite difference operator has confirmed the stability and high efficiency.The absorbing boundary condition is used to elimi- nate artificial reflections.Numerical examples show that this new irregular-grid finite difference method is of higher performance than conventional methods using regular rectangular grids in simulating elastic wave propagation in heterogeneous anisotropic porous media.

Theoretical and numerical investigation of HF elastic wave propagation in two-dimensional periodic beam lattices

The elastic wave propagation phenomena in two-dimensional periodic beam lattices are studied by using the Bloch wave transform. The numerical modeling is applied to the hexagonal and the rectangular beam lattices, in which, both the in-plane （with respect to the lattice plane） and out-of-plane waves are considered. The dispersion relations are obtained by calculating the Bloch eigenfrequencies and eigenmodes. The frequency bandgaps are observed and the influence of the elastic and geometric properties of the primitive cell on the bandgaps is studied. By analyzing the phase and the group velocities of the Bloch wave modes, the anisotropic behaviors and the dispersive characteristics of the hexagonal beam lattice with respect to the wave prop- agation are highlighted in high frequency domains. One im- portant result presented herein is the comparison between the first Bloch wave modes to the membrane and bend- ing/transverse shear wave modes of the classical equivalent homogenized orthotropic plate model of the hexagonal beam lattice. It is shown that, in low frequency ranges, the homog- enized plate model can correctly represent both the in-plane and out-of-plane dynamic behaviors of the beam lattice, its frequency validity domain can be precisely evaluated thanks to the Bloch modal analysis. As another important and original result, we have highlighted the existence of the retro- propagating Bloch wave modes with a negative group veloc- ity, and of the corresponding ＂retro-propagating＂ frequency bands.

Research on propagation properties of elastic waves in two-phase anisotropic media

With the development of seismic engineering and seismic exploration of energy, the underground media that westudy are more and more complicated. Conventional anisotropy theory or two-phase isotropy theory is difficult todescribe anisotropic media containing fluid, such as fractures containing gas, shales containing water Based onBlot theory about two-phase anisotropy, with the use of elastic plane wave equations, we get Christoffel equations.We calculate and analyze the effects of frequency on phase velocity, attenuation, amplitude ratio and polarizationdirection of elastic waves of two-phase, transversely isotropic media. Results show that frequency affects slow Pwave the greatest among the four kinds of waves, i.e., fast P wave, slow P wave, fast S wave and slow S wave.Fluid phase amplitude to solid phase amplitude ratio of fast P wave, fast S wave and slow S wave approaches unitfor large dissipation coefficients. Polarization analysis shows that polarization direction of fluid phase displacement is different from, not parallel to or reverse to, that of solid phase displacement in two-phase anisotropic media.