Elastic wave propagation characteristics under anisotropic squirt-flow condition

On the generalized self-consistent model for elastic wave propagation in composite materials

Wave propagation in glass fibre-reinforced polymer (GFRP) bars subjected to progressive damage—Experimental and numerical results

Shell type structures are used in a wide variety of industrial components, especially in automotive or aeronautic industries. A full understanding and predictive modeling of elastic waves propagation through those structures can be of great importance at their design stage for several reasons, such as the noise reduction in a vehicle or the control of high energy shock waves triggered by pyrotechnic cut in a space launcher. This paper deals with the theoretical and numerical modeling of elastic wave propagation in shells. The aim is to discuss the relevance of theoretical models and to develop appropriate and efficient numerical tools.

This paper considers the propagation of elastic waves in periodic two-dimensional mass–spring structures with diagonal springs. The second-neighbor interactions in non-diagonal directions are included to account for the nonlocal effect. The influences of the spring stiffness in the diagonal directions and the nonlocal effect on the propagation characteristics of elastic waves are then scrutinized. Through the dispersion relation curve and the equi-frequency contours, it is seen that when the diagonal spring stiffness increases, the slope of the second curve in the –M direction will not always be positive, meaning that the negative group velocity occurs. Therefore, an incident wavevector with a chosen angle to the negative group velocity can lead to the negative refraction phenomenon in the two-dimensional mass–spring structure. Another interesting phenomenon called directional radiation of elastic waves can also be achieved by adjusting the nonlocal effect. Within a certain range, the stronger the nonlocal effect in a specific direction is, the more obviously the elastic waves propagate along this direction. In this paper, we theoretically analyze and numerically simulate the phenomena of negative refraction and directional wave propagation by choosing a proper set of parameters of the two-dimensional mass–spring structure.

Seismic exploration of deep oil/gas reservoirs involves the propagation of seismic waves in high-pressure media. Traditional elastic wave equations are not suitable for describing such media. The theory of acoustoelasticity establishes the dynamic equation of wave propagating in prestressed media through constitutive relation using third-order elastic constants. Many studies have been carried out on numerical simulations for acoustoelastic waves, but are mainly limited to 2D cases. A standard staggered-grid (SSG) finite-difference (FD) approach and the perfectly matched layer (PML) absorbing boundary are combined to solve 3D first-order velocity-stress equations of acoustoelasticity to simulate wave propagating in 3D prestressed solid medium. Our numerical results are partially validated by plane-wave analytical solution through the comparison of calculated and theoretical P-/S-wave velocities as a function of confining prestress. We perform numerical simulations of acoustoelastic waves under confining, uniaxial and pure-shear prestressed conditions. The results show the stress-induced velocity anisotropy in acoustoelastic media, which is closely related to the direction of prestresses. Comparisons to seismic simulations based on the theory of elasticity illustrate the limitation of conventional elastic simulations for prestressed media. Numerical simulations prove the significant effect of prestressed conditions on seismic responses.

Continuous wave acoustic method for determination of moisture content in agricultural soil

A simple model is proposed to study wave propagation in hard sphere-reinforced elastic random composites. Compared to existing related models, the proposed model is featured by a modified form of classical elastodynamic equations in which the inertia term is substituted by the acceleration field of the mass centre of a representative unit cell, supplied with a derived simple differential relation between the displacement field of the composite and the displacement field of the mass centre of a representative unit cell. The present model enjoys conceptual and mathematical simplicity although it is restricted to hard sphere-filled elastic composites in which the elastic moduli of embedded spheres are much (at least 4-5 times) stiffer than those of a softer matrix. Explicit formulas are derived for the attenuation coefficient and the effective phase velocity of plane longitudinal P-waves and transverse S-waves. The efficiency and reasonable accuracy of the present model are demonstrated by reasonably good agreement between the predicted results and some established known data. The proposed model could offer a potential general method to study various three-dimensional dynamic problems of hard sphere-filled elastic random composites.

In this study, we use the interferometric ghosts of the scattered surface waves for estimating the location of a near-surface scatterer. The scatterer is embedded in a homogeneous halfspace. We perform finite difference modeling of elastic wave propagation to calculate the seismograms. We choose different locations for virtual sources and obtain interferometric estimates of the scattered surface waves for these locations. We solve the inverse problem by using the interferometric traveltimes of the modeled scattered surface waves. The results show that the location of the scatterer is reasonably well estimated.

This paper is concerned with the numerical modeling of elastic wave propagation in layered media. It considers two isotropic homogeneous elastic solids in perfect contact. The interface is parallel to the free surface. Two finite difference methods are developed. The usefulness of the methods are investigated for long time simulations and the accuracy of the results are compared with the response from an approximate model.

Three‐dimensional (3D) elastic wave propagation modeling in the velocity‐stress formulation using finite differences (FDs) have been investigated for a homogeneous medium covered by a representative, relatively steep surface topography consisting of a 1D square root function. This scenario using various numerical implementations is explored. The behavior with regard to stability of long simulations is expected to be indicative of each numerical implementation's robustness for other types of topographies/media. Employing various combinations of the FD order was found only to change the time of the first incidence of instability. On the other hand, nonequidistant grids in the horizontal and vertical directions are found to be extremely useful for long‐term stability of 3D wave propagation modeling with our free‐surface boundary condition for single‐valued topographies. In particular dz ≥ (3/2) dx is found completely stable for all tested vP/vSratios. Such relationships of using dz > dx are also favorable for more accurate Rayleigh‐wave modeling. Setting the density equal to 1/10 of its interior value at one layer only at the surface is another simple means of achieving stability.

Effect of negative permeability on elastic wave propagation in magnetoelastic multilayered composites

To improve the performance of stencil computations tend to be challenging due to their inherent memory access pattern. A condition that remains valid for wave propagation simulation engines. Moving towards using elastic waves instead of acoustic ones (e.g., used in medical imaging), results in computationally more expensive processes along with an increment on memory usage. Despite the computational demand, the elevated cost of exploration joined to a low success rate drove the Oil & Gas industry to rely on elastic anisotropic wave propagation models as the core of many geophysical imaging mechanisms to explore subsurface without an excessive investment. Moreover, exploration codes are expected to run fast and efficiently on modern architectures. The Intel® Xeon Phi™ processors emerge as an energy efficient solution that provides a good trade-off between market price and computing capability. In this paper, we study the effect of several optimization techniques using the YASK framework to implement and evaluate a 25-points stencil of an elastic wave propagation engine for Intel Xeon Phi processors. The results showed improvements of up to 7× in computations and 8× in memory bandwidth with respect to the non-tuned version, reaching up to 75% of the attainable performance. We collected performance metrics for a set of the most representative optimizations and revealed the relation between each strategy and fundamental characteristics of both code and hardware.

We consider the propagation of elastic waves in gas-filled porous media at small but non-zero values of Knudsen numbers $$ {\text{Kn}} $$ , where $$ {\text{Kn}} = \lambda /l $$ , $$ \lambda $$ is the mean free path of gas molecules; $$ l $$ is the characteristic size of inclusion (the so-called slip regime). In this case, it is possible to apply the classic equations of hydrodynamics with modified boundary conditions at solid walls. We have assumed that the gas molecules distribution function is satisfied at the modified Maxwell boundary conditions (Struchtrup 2013; Mohammadzadeh and Struchtrup 2015). We have obtained the expressions for drag and added mass coefficients for the Biot equations of poroelasticity for a system of randomly oriented gas-filled cylindrical capillaries. Our calculations have shown that the drag and added mass coefficients depend considerably on the Knudsen number and the properties of the surface. The influence of the interfacial slip effect on the velocities of the compressional wave of the first kind and shear wave is small, but the velocity and attenuation of the compressional wave of the second kind are considerably influenced by this effect. The results obtained show the fundamental possibility of the determination of the accommodation coefficient by measuring the velocity of the compressional wave of the second kind for different values of the Knudsen number.

A coupled three-dimensional fluid-elastic wave propagation mathematical model has been developed to handle environmental interactions in the ocean between a fluid medium and an elastic bottom. The existing model combines three-dimensional fluid and elastic wave propagation models with the incorporation of a set of horizontal fluid-elastic interface conditions. This paper extends the above model to consider an irregular fluid-elastic interface. The theoretical development of the irregular fluid-elastic interface equations is presented. Comparison of the fluid-elastic irregular interface to the horizontal case is one of the major objectives of this paper. Another major objective of this paper is the construction of a complete set of three-dimensional fluid-elastic wave equations, including the irregular interface, in a form suitable for stable numerical solution.

Ordered Ca3TaGa3Si2O14 crystals from langasite family are grown. Full sets of elastic, piezoelectric and dielectric constant are determined. Some characteristics of elastic wave propagation are calculated.