Numerical simulation of scattering wave imaging in a goaf

Goafs are threats to safe mining.Their imaging effects or those of other complex geological bodies are often poor in conventional reflected wave images.Hence,accurate detection of goafs has become an important problem,to be solved with a sense of urgency.Based on scattering theory,we used an equivalent offset method to extract Common Scattering Point gathers,in order to analyze different scattering wave characteristics between Common Scattering Point and Common Mid Point gathers and to compare stack and migration imaging effects.Our research results show that the scattering wave imaging method is more efficient than the conventional imaging method and is therefore a more effective imaging method for detecting goafs and other complex geological bodies.It has important implications for safe mining procedures and infrastructures.

Scattering wave field around a cavity with circular cross-section embedded in saturated soil using boundary element method

Based on Biot’s theory and considering the properties of a cavity,the boundary integral equations for the numerical simulation of wave scattering around a cavity with a circular cross-section embedded in saturated soil are obtained using integral transform methods.The Cauchy type singularity of the boundary integral equation is discussed.The effectiveness of the properties of soil mass and incident field on the dynamic stress concentration and pore pressure concentration around a cavity is analyzed.Our results are in good agreement with the existing solution.The numerical results of this work show that the dynamic stress concentration and pore pressure concentration are influenced by the degree of fluid–solid coupling as well as the pore compressibility and water permeability of saturated soil.With increased degree of fluid–solid coupling,the dynamic stress concentration improves from 1.87 to 3.42 and the scattering becomes more significant.With decreased index of soil mass compressibility,the dynamic stress concentration increases and its maximum reaches 3.67.The dynamic stress concentration increases from 1.64 to 3.49 and pore pressure concentration improves from 0.18 to 0.46 with decreased water permeability of saturated soil.

SCATTERING FAR FIELD SOLUTION OF SH-WAVE BY MOVABLE RIGID CYLINDRICAL INTERFACE INCLUSION

The Green's function is used to solve the scattering far fieldsolution of SH-wave by a mov- able rigid cylindrical interfaceinclusion in a linear elastic body. First, a suitable Green'sfunction is devel- oped, which is the fundamental displacementsolution of an elastic half space with a movable rigid half-cylin-drical inclusion impacted by out-of-plane harmonic line source loadedat any point of its horizontal surface.

scattering and refracting of plane strain wave by a cylindrical inclusion in fluidsaturated soils

The problems involving scattering and refracting of plane strain wave by a cylinder in fluidsaturated soils are studied in this paper. The theoretical solutions for the amplitudes of scattering and refracting are derived on the basis of modified Biot′s theory for twophase medium in engineering, in which the compressibility of soil grains and the viscous coupling effect between the skeletal frame and the pore fluid are taken into account. Moreover, the amplitude equations of potential function are presented when the fluidsaturated soil cylindrical inclusion is degenerated into cylindrical cavity, rigid cylinder, single phase elastic cylinder and single phase fluid cylinder.

Dynamic Stress Intensity Factor and Scattering of SH Wave by Circle Canyon and Crack

The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.

SCATTERING OF FLEXURAL WAVES AND DYNAMIC STRESS CONCENTRATIONS IN MINDLIN'S THICK PLATES WITH A CUTOUT

Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique,the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations.As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular,elliptic cutout are graphically presented in sequence.

SHEAR WAVE SCATTERING FROM A PARTIALLY DEBONDED PIEZOELECTRIC CYLINDRICAL INCLUSION

The scattering of SH wave by a cylindrical piezoelectric inclusion partially debonded from its surrounding piezoelectric material is investigated using the wave function expansion method and singular integral equation technique. The debonding regions are modeled as mul- tiple arc-shaped interface cracks with non-contacting faces. By expressing the scattered ?elds as wave function expansions with unknown coe?cients, the mixed boundary value problem is ?rstly reduced to a set of simultaneous dual series equations. Then dislocation density functions are introduced as unknowns to transform these dual series equations into Cauchy singular integral equations of the ?rst type, which can be numerically solved easily. The solution is valid for arbi- trary number and size of the debonds. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond and two debonds. The e?ects of incidence direc- tion, crack con?guration and various material parameters on the dynamic stress intensity factors are respectively discussed. The solution of this problem is expected to ?nd applications in the investigation of dynamic fracture properties of piezoelectric materials with cracks.

Oblique Water Wave Scattering by Bottom Undulation in a Two-layer Fluid Flowing Through a Channel

The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.

Wave Scattering by Double Slotted Barriers in A Steady Current:Experiments

The adoption of slotted breakwaters can be an ideal option in the protection of very large near-shore floating struc-trees that may extend offshore to a considerable water depth. In this paper, we experimently investigated the behaviour of wave transmission and reflection coefficients of double slotted barriers in the presence of a steady opposing current. The experimental results show that opposing currents have only minor effects on wave reflection, but can significantly reduce the wave transmission through double slotted barriers. The experimental results suggest that coastal currents should be taken into consideration for an economical design of slotted breakwaters.

A model for the scattering of long waves by slotted breakwaters in the presence of currents

Slotted breakwaters have been used to provide economical protection from waves in harbors where surface waves and currents may co-exist. In this paper, the effects of currents on the wave scattering by slotted breakwaters are investigated by using a simple model. The model is based on a long wave approximation. The effects of wave height, barrier geometry and current strength on the reflection and transmission coefficients are examined by the model. The model results are compared with recent experimental data. It is found that both the wave-following and wave-opposing currents can increase the reflection coefficient and reduce the transmission coefficient. The model can be used to study the interaction between long waves and slotted breakwaters in coastal waters.

Wave scattering by undulating bed topography in a two-layer ocean

The problem of wave scattering by undulating bed topography in a two-layer ocean is investigated on the basis of linear theory. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating at both the free surface of the upper layer and the interface between the two layers. Due to a wave train of a particular mode incident on an obstacle which is bottom-standing on the lower layer, reflected and transmitted waves of both modes are created by the obstacle. For small undulations on the bottom of the lower layer, a perturbation method is employed to obtain first-order reflection and transmission coefficients of both modes for incident wave trains of again both modes in terms of integrals involving the bed-shape fimction. For sinusoidal undulations, numerical results are presented graphically to illustrate the energy transfer between the waves of different modes by the undulating bed.

Scattering of seismic waves by three-dimensional large-scale hill topography simulated by a fast parallel IBEM

To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual(GMRES) algorithm and constructing an Open MP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.

Elastic Wave Scattering by Two-Dimensional Periodical Array of Cylinders

We extend the multiple-scattering theory (MST) for elastic wave scattering and propagating in two-dimensional composite. The formalism for the band structure calculation is presented by taking into account the full vector character of the elastic wave. As a demonstration of application of the formalism we calculate the band structure of elastic wave propagating in a two-dimensional periodic arrangement of cylinders. The results manifest that the MST shows great promise in complementing the plane-wave (PW) approach for the study of elastic wave.

SCATTERING OF ELASTIC WAVES IN AN ELASTIC MATRIX CONTAINING AN INCLUSION WITH INTERFACES

Based on the theory of elastic dynamics, the scattering of elasticwaves and dynamic stress concentration of fiber-reinforced compositewith interfaces are studied. Analytical expressions of elastic wavesin different medium areas are presented and an analytic method ofsolving this problem is established. The mode coefficients aredetermined by means of the continuous conditions of displacement andstress on the boundary of the interfaces. The influence of materialproperties and structural size on the dynamic stress con- centrationfactors near the interfaces is analyzed.

Wave Scattering by Porous Bottom Undulation in a Two Layered Channel

The scattering of plane surface waves by bottom undulations in channel flow consisting of two layers is investigated by assuming that the bed of the channel is composed of porous material. The upper surface of the fluid is bounded by a rigid lid and the channel is unbounded in the horizontal directions. There exists only one wave mode corresponding to an internal wave. For small undulations, a simplified perturbation analysis is used to obtain first order reflection and transmission coefficients in terms of integrals involving the shape function describing the bottom. For sinusoidal bottom undulations and exponentially decaying bottom topography, the first order coefficients are computed. In the case of sinusoidal bottom the first order transmission coefficient is found to vanish identically. The numerical results are depicted graphically in a number of figures.

Scattering of Circularly Polarized Terahertz Waves on a Graphene Nanoantenna

We present a surface current method to model the graphene rectangular nanoantenna scattering in the terahertz band with Comsol. Compared with the equivalent thin slab method, the results obtained by the surface current method are more accurate and efficient. Then the electromagnetic scattering of circularly polarized terahertz waves on graphene nanoantennas is numerically analyzed by utilizing the surface current method. The depen- dences of the antenna resonant frequency with the circularly polarized wave on width and length are consistent with those for the linear polarized waves. These results are proved to be useful to design et^cient nanoantennas in terahertz wireless communications.

A Linear Theory for Wave Scattering by Double Slotted Barriers in Weak Steady Currents

Recent experimental results have shown that the presence of a steady current can significantly reduce the energy of transmitted waves. In this paper, a theory is developed to study the wave scattering by single or double vertical slotted barriers in the presence of a weak uniform current. The quasi-linear theory is based on an eigenfunction expansion method. Comparisons between theory and existing experimental results for both single slotted barrier and double slotted barriers show satisfactory agreements. In consideration of wave propagation in a weak current it is found that the friction factor used to characterize the head loss at the slotted barrier depends on both the geometry of the slotted barrier and the strength of the steady current.

Waves scattering induced by an interface crack in a coated material

This paper deals with the two-dimensional problem of elastic wave scattering from a finite crack at the interface between a coated material layer and its substrate. By adopting the Fourier transform method and introducing the crack opening displacement function, the boundary value problem is simplified for numerically solving a system of Cauchy-type singular integral equations by means of Jacobi polynomial expansion. The stress intensity factors and the crack opening displacements are defined in terms of the integral equations solutions. The influence of the dimensionless wave number and the ratio of crack length to layer thickness on the stress intensity factors and crack opening displacements are discussed.

Water Wave Scattering by a Nearly Circular Cylinder Submerged Beneath an Ice-cover

Assuming linear theory, the two-dimensional problem of water wave scattering by a horizontal nearly circular cylinder submerged in infinitely deep water with an ice cover modeled as a thin-elastic plate floating on water, is investigated here. The cross-section of the nearly circular cylinder is taken as r=a（1＋δC（θ））, where a is the radius of the corresponding circular cross-section of the cylinder, δ is a measure of small departure of the cross-section of the cylinder from its circularity and C（θ） is the shape function. Using a simplified perturbation technique the problem is reduced to two independent boundary value problems up to first order in δ. The first one corresponds to water wave scattering by a circular cylinder submerged in water with an ice-cover, while the second problem describes wave radiation by a submerged circular cylinder and involves first order correction to the reflection and transmission coefficients. The corrections are obtained in terms of integrals involving the shape function. Assuming a general Fourier expansion of the shape function, these corrections are evaluated approximately. It is well known that normally incident wave trains experience no reflection by a circular cylinder submerged in infinitely deep water with an ice cover. It is shown here that the reflection coefficient also vanishes up to first order for some particular choice of the shape function representing a nearly circular cylinder. For these cases, full transmission occurs, only change is in its phase which is depicted graphically against the wave number in a number of figures and appropriate conclusions are drawn.

ELASTIC WAVE SCATTERING AND DYNAMIC STRESS IN COMPOSITE WITH FIBER

Based on the theory of elastic dynamics, multiple scattering of elastic waves and dynamic stress concentrations in fiber_reinforced composite were studied. The analytical expressions of elastic waves in different region were presented and an analytic method to solve this problem was established. The mode coefficients of elastic waves were determined in accordance with the continuous conditions of displacement and stress on the boundary of the multi_interfaces. By making use of the addition theorem of Hankel functions, the formulations of scattered wave fields in different local coordinates were transformed into those in one local coordinate to determine the unknown coefficients and dynamic stress concentration factors. The influence of distance between two inclusions, material properties and structural size on the dynamic stress concentration factors near the interfaces was analyzed. It indicates in the analysis that distance between two inclusions, material properties and structural size has great influence on the dynamic properties of fiber_reinforced composite near the interfaces. As examples, the numerical results of dynamic stress concentration factors near the interfaces in a fiber_reinforced composite are presented and discussed.