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.

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.

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.

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.

An improved convolution perfectly matched layer for elastic second-order wave equation

A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.

THE PARAMETER PERTURBATION METHOD ON ELASTIC WAVE EQUATION IN INHOMOGENEOUS MEDIUM

In this paper, the medium parameters of the elastic wave equation in inhomogeneous medium are rewritten by introducing the referential variables and the perturbational variables, and the wave equation whose sources are the medium parameter perturbational term in homogeneous medium is obtained. By using the Green function theory, the integral equation of the perturbational parameters is obtained. Then the displacement field in homogeneous medium is considered the result of the first iteration, and the displacement field is solved by this integral equation. When the perturbations of medium parameters are about 50 percent, this method can solve the displacement field effectively. from the analysis of the numerical results, the characteristics of wave field in inhomogeneous medium are obtained. The results conform with the local principles of wave function in inhomogeneous medium.

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.

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.

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.

SCATTERING OF HARMONIC ANTI-PLANE SHEAR WAVES BY AN INTERFACE CRACK IN MAGNETO-ELECTRO-ELASTIC COMPOSITES

The dynamic behavior of an interface crack in magneto-electro-elastic composites under harmonic elastic anti-plane shear waves is investigated for the permeable electric boundary conditions. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable was the jump of the displacements across the crack surfaces. To solve the dual integral equations, the jump of the displacements across the crack surface was expanded in a series of Jacobi polynomials. Numerical examples were provided to show the effect of the length of the crack, the wave velocity and the circular frequency of the incident wave on the stress, the electric displacement and the magnetic flux intensity factors of the crack. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for anti-plane shear problem.

DYNAMIC BEHAVIOR OF TWO PARALLEL SYMMETRY CRACKS IN MAGNETO-ELECTRO-ELASTIC COMPOSITES UNDER HARMONIC ANTI-PLANE WAVES

The dynamic behavior of two parallel symmetry cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied by Schmidt method. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable is the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surface were expanded in a series of Jacobi polynomials. The relations among the electric filed, the magnetic flux and the stress field were obtained. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for the dynamic anti-plane shear fracture problem. The shielding effect of two parallel cracks was also discussed.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

RAYLEIGH LAMB WAVES IN MICROPOLAR ISOTROPIC ELASTIC PLATE

The propagation of waves in a homogeneous isotropic micropolar elastic cylindrical plate subjected to stress free conditions is investigated. The secular equations for symmetric and skew symmetric wave mode propagation are derived. At short wave limit, the secular equations for symmetric and skew symmetric waves in a stress free circular plate reduces to Rayleigh surface wave frequency equation. Thin plate results are also obtained. The amplitudes of displacements and microrotation components are obtained and depicted graphically. Some special cases are also deduced from the present investigations. The secular equations for symmetric and skew symmetric modes are also presented graphically.

Theconvolutionaldiferentiatormethodfornu┐mericalmodelingofacousticandelasticwave┐fields

Based on the techniques of forward and inverse Fourier transformation, the authors discussed the design scheme of ordinary differentiator used and applied in the simulation of acoustic and elastic wavefields in isotropic media respectively. To compress Gibbs effects by truncation effectively, Hanning window is introduced in. The model computation shows that, the convolutional differentiator method has the advantages of rapidity, low requirements of computer′s inner storage and high precision, which is a potential method of numerical simulation.

THE SCATTERING OF HARMONIC ELASTIC ANTI-PLANE SHEAR WAVES BY A FINITE CRACK IN INFINITELY LONG STRIP USING THE NON-LOCAL THEORY

In this paper, the scattering of harmonic anti-plane shear wavesby a finite crack in infinitely long strip is studied using thenon-local theory. The Fourier transform is applied and a mixedboundary value problem is formulated. Then a set of dual integralequations is solved using the Schmidt method instead of the first orthe second integral equation method. A one-dimensional non-localkernel is used instead of a two-di- mensional one for the anti-planedynamic problem to obtain the stress occurring at the crack tips.Contrary to the classical elasticity solution, it is found that nostress singularity is present at the crack tip. The non-local dynamicelastic solutions yield a finite hoop stress at the crack tip, thusallowing for a fracture criterion based on the maximum dynamic stresshypothesis. The finite hoop stress at the crack tip depends on thecrack length, the width of the strip and the lattice parameters.

Frequency-domain Elastic wave Simulation Based on the Nonoverlapping Domain Decomposition Method

A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrimination of the elastic wave equation are derived. The algorithm of no overlapping domain decomposition method is given. The method solves the elastic wave equation by iteratively solving sub problems defined on smaller sub domains. Numerical computations both for homogeneous and inhomogeneous media show the effectiveness of the proposed method. This method can be used in the full-waveform inversion.

New Periodic Wave Solutions and Their Interaction for （2＋1）-dimensional KdV Equation

A class of new doubly periodic wave solutions for （2＋1）-dimensional KdV equation are obtained by introducing appropriate Jacobi elliptic functions and Weierstrass elliptic functions in the general solution（contains two arbitrary functions） got by means of multilinear variable separation approach for （2＋1）-dimensional KdV equation. Limiting cases are considered and some localized excitations are derived, such as dromion, multidromions, dromion-antidromion, multidromions-antidromions, and so on. Some solutions of the dromion-antidromion and multidromions-antidromions are periodic in one direction but localized in the other direction. The interaction properties of these solutions, which are numerically studied, reveal that some of them are nonelastic and some are completely elastic. Furthermore, these results are visualized.

Using Refined Theory to Studied Elastic Wave Scattering and Dynamic Stress Concentrations in Plates with Two Cutouts

In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.

SCATTERING OF ANTI-PLANE SHEAR WAVES BY A SINGLE CRACK IN AN UNBOUNDED TRANSVERSELY ISOTROPIC ELECTRO-MAGNETO-ELASTIC MEDIUM

A theoretical treatment of the scattering of anti-plane shear(SH) waves is provided by a single crack in an unbounded transversely isotropic electro-magneto-elastic medium. Based on the differential equations of equilibrium, electric displacement and magnetic induction intensity differential equations, the governing equations for SH waves were obtained. By means of a linear transform, the governing equations were reduced to one Helmholtz and two Laplace equations. The Cauchy singular integral equations were gained by making use of Fourier transform and adopting electro-magneto impermeable boundary conditions. The closed form expression for the resulting stress intensity factor at the crack was achieved by solving the appropriate singular integral equations using Chebyshev polynomial. Typical examples are provided to show the loading frequency upon the local stress fields around the crack tips. The study reveals the importance of the electro-magneto-mechanical coupling terms upon the resulting dynamic stress intensity factor.

The Quasi-Static Approximation of Heat Waves in Anisotropic Thermo-Elastic Media

The equilibrium equations of anisotropic media, coupled to the heat conduction equations, are studied here based on the standard spaces of the physical presentation, in which an new thermo-elastic model based on the second law of thermodynamics is induced. The uncoupled heat wave equation for anisotropic media is deduced. The results show that the equation of heat wave is of the properties of dissipative waves. In final part of this paper, we discuss the propagation behaviour of heat waves for transversely isotropic media.