Scattering by penetrable acoustic targets

An acoustic target of constant density ?_t and variable index of refraction is imbedded in a surrounding acoustic fluid of constant density ?_a. A time harmonic wave propagating in the surrounding fluid is incident on the target. We consider two limiting cases of the target where the parameter ε ≡ ?_a/?_t → 0 (the nearly rigid target) or ε → ∞ (the nearly soft target). Wh en the frequency of the incident wave is bounded away from the ‘in-vacuo’ resonant frequencies of the target, the resulting scattered field is essentially the field scattered by the rigid target for ε = 0 or the soft target if ε → ∞. However, when the frequency of the incident wave is near a resonant frequency,the target oscillates and its interaction with the surrounding fluid produces peaks in the scattered field amplitude. In this paper we obtain asymptotic expansions of the solutions of the scattering problems for the nearly rigid and the nearly soft targets as ε → 0 or ε → ∞, respectively, that are uniformly valid in the incident frequency. The method of matched asymptotic expansions is used in the analysis. The outer and inner expansions correspond to the incident frequencies being far or near to the resonant frequencies, respectively. We have applied the method only to simple resonant frequencies, but it can be extended to multiple resonant frequencies. The method is applied to the incidence of a plane wave on a nearly rigid sphere of constant index of refraction. The far field expressions for the scattered fields, including the total scattering cross-sections, that are obtained from the asymptotic method and from the partial wave expansion of the solution are in close agreement for sufficiently small values of ε.     

Scattering of elastic waves by an arbitrary small imperfection in the surface of a half-space

T , the first of two articles, is concerned with the scattering of elastic waves by arbitrary surface-breaking or near surface defects in an isotropic half-plane. We present an analytical solution, by the method of matched asymptotic expansions, when the parameter , which is the ratio of a typical length scale of the imperfection to the incident radiation's wavelength, is small. The problem is formulated for a general class of small defects, including cracks, surface bumps and inclusions, and for arbitrary incident waves. As a straightforward example of the asymptotic scheme we specialize the defect to a two-dimensional circular void or protrusion, which breaks the free surface, and assume Rayleigh wave excitation ; this inner problem is exactly solvable by conformal mapping methods. The displacement field is found uniformly to leading order in , and the Rayleigh waves which are scattered by the crack are explicitly determined. In the second article we use the method given here to tackle the important problem of an inclined edge-crack. In that work we show that the scattered field can be found to any asymptotic order in a straightforward manner, and in particular the Rayleigh wave coefficients are given to O(²).     

Scattering of flexural wave in a thin plate with multiple circular holes by using the multipole Trefftz method

The scattering of flexural wave by multiple circular holes in an infinite thin plate is analytically solved by using the multipole Trefftz method. The dynamic moment concentration factor (DMCF) along the edge of circular holes is determined. Based on the addition theorem, the solution of the field represented by multiple coordinate systems centered at each circle can be transformed into one coordinate system centered at one circle, where the boundary conditions are given. In this way, a coupled infinite system of simultaneous linear algebraic equations is derived as an analytical model for the scattering of flexural wave by multiple holes in an infinite plate subject to the incident flexural wave. The formulation is general and is easily applicable to dealing with the problem containing multiple circular holes. Although the number of hole is not limited in our proposed method, the numerical results of an infinite plate with three circular holes are presented in the truncated finite system. The effects of both incident wave number and the central distance among circular holes on the DMCF are investigated. Numerical results show that the DMCF of three holes is larger than that of one, when the space among holes is small and meanwhile the specified direction of incident wave is subjected to the plate.     

Scattering attenuation of elastic waves due to low-contrast inclusions

Seismic scattering attenuation due to random lithospheric heterogeneity has been theoretically modeled using two approaches. One approach is the Born approximation theory (BAT), which is primarily used to treat weak continuous heterogeneity, and the other approach is the Foldy approximation theory (FAT), which deals with sparsely distributed discrete inclusions. We apply the BAT to elastic wave scattering due to inclusions having low contrast with the matrix, and compare the results with those predicted by the FAT. We thus investigate the valid wavenumber range of the BAT based on a reasonable assumption that the inclusions are distributed so sparsely that the FAT is effectively correct for any wavenumber. For simplicity, we consider a specific type of round inclusion, which is either two- or three-dimensional and has a two-valued wave velocity and/or mass density. Both theories are confirmed to yield essentially equivalent results below a certain wavenumber limit, depending on the contrast. This is known as the Rayleigh-Gans scattering regime. Beyond the wavenumber limit, the BAT overestimates the attenuation for common-mode scattering due to wave-velocity contrast, but remains valid with respect to the attenuation for scattering due to mass-density contrast and/or conversion scattering. These conclusions are independent of the spatial dimensions of the media as well as the modes of the elastic waves (P or S). Some advantages of the BAT over the FAT for application to low-contrast inclusions are discussed.     

SH波绕界面孔的散射

用波函数展开方法研究了SH波绕界面孔的散射问题。由入射、反射和透射波组成的自由波场与孔的散射场叠加成总波场。按照一定方式将两个半平面散射波场延拓于全平面，通过Hankel-Fourier展开方法求得了任意形状孔散射场的级数解。以椭圆形孔为例计算了孔边缘的动应力集中系数。     

Scattering and depolarization of electromagnetic waves by a horizontally weakly inhomogeneous medium

A general theory of wave scattering from a weakly inhomogeneous medium is developed for the case where the inhomogeneity varies parallel to the boundary plane. The method of small perturbation is used and terms are carried up to and including the second order. It is found that the scattered waves are depolarized and present in all directions. In the special case of forward- or backscattering the depolarized fields are of the second order and are seen to result from a multiple scattering process; while in other directions, these fields could be of the first order, and result from a single scattering process.     

A Refined Asymptotic Model of Fluid-Structure Interaction in Scattering by Elastic Shells

A refined asymptotic model of fluid-structure interaction in scattering by elastic shells is proposed. The model takes into consideration transverse compression of a shell by a fluid and some other phenomena. As an illustration, scattering of a plane acoustic wave by a circular cylindrical shell is considered. Comparison of numerical data corresponding various approximate approaches is provided. This revised version was published online in July 2006 with corrections to the Cover Date.     

Asymptotic modeling of Helmholtz resonators including thermoviscous effects

We systematically employ the method of matched asymptotic expansions to model Helmholtz resonators, with thermoviscous effects incorporated starting from first principles and with the lumped parameters characterizing the neck and cavity geometries precisely defined and provided explicitly for a wide range of geometries. With an eye towards modeling acoustic metasurfaces, we consider resonators embedded in a rigid surface, each resonator consisting of an arbitrarily shaped cavity connected to the external half-space by a small cylindrical neck. The bulk of the analysis is devoted to the problem where a single resonator is subjected to a normally incident plane wave; the model is then extended using “Foldy’s method” to the case of multiple resonators subjected to an arbitrary incident field. As an illustration, we derive critical-coupling conditions for optimal and perfect absorption by a single resonator and a model metasurface, respectively.     

薄膜涂层材料中币形界面裂纹的弹性波散射

研究了薄膜涂层材料中币形界面裂纹的弹性波散射问题,建立了含有币形界面裂纹的覆层半空间模型,采用Hankel积分变换,将裂纹对弹性波散射的问题转化为求解矩阵形式的奇异积分方程。结合渐近分析和围道积分技术得到积分方程的解,进一步推导了散射波的应力场和位移场,以及动应力强度因子的理论计算公式。在数值算例中,分析了不同材料组合和裂纹尺寸情况下动应力强度因子与入射波频率的关系,并给出了裂纹张开位移的结果。为薄膜涂层材料的动态破坏分析提供了一定的理论基础。     

Scattering far field solution of SH-wave by movable rigid cylindrical interface inclusion

The Green's function is used to solve the scattering far field solution of SH-wave by a movable rigid cylindrical interface inclusion in a linear elastic body. First, a suitable Green's function is developed, which is the fundamental displacement solution of an elastic half space with a movable rigid half-cylindrical inclusion impacted by out-of-plane harmonic line source loaded at any point of its horizontal surface. By using the Green's function, a series of Fredholm integral equations of the first kind which determine the scattering far field can be set up. Then the paper gives the expressions on the far field including the displacement mode of scattering wave and the solution of scattering cross-section. Finally, some examples and numerical results are discussed to analyze the influence of the combination of different media parameters on the answer of far field.     

Scattering of a longitudinal wave by a circular crack in a fluid-saturated porous medium

Physical properties of many natural and man-made materials can be modelled using the concept of poroelasticity. Some porous materials, in addition to the network of pores, contain larger inhomogeneities such as inclusions, cavities, fractures or cracks. A common method of detecting such inhomogeneities is based on the use of elastic wave scattering. We consider interaction of a normally incident time-harmonic longitudinal plane wave with a circular crack imbedded in a porous medium governed by Biot’s equations of dynamic poroelasticity. The problem is formulated in cylindrical co-ordinates as a system of dual integral equations for the Hankel transform of the wave field, which is then reduced to a single Fredholm integral equation of the second kind. It is found that the scattering that takes place is predominantly due to wave induced fluid flow between the pores and the crack. The scattering magnitude depends on the size of the crack relative to the slow wave wavelength and has it’s maximum value when they are of the same order.     

Scattering of elastic waves by a small inclined surface-breaking crack

I we examine the scattering of Rayleigh waves by an inclined two-dimensional plane surface-breaking crack in an isotropic elastic half-plane. We follow the method already introduced by the authors (A and W , 1992a, J. Mech. Phys. Solids 40, 1683) to obtain an analytical solution when the parameter , which is the ratio of a typical length scale of the imperfection to the incident radiation's wavelength, is small. The procedure employed is the method of matched asymptotic expansions, which involves determining the scattered wave field both away from and near the crack. The outer solution is specialized from the general expansion in the first part of this work (A and W , 1992a, J. Mech. Phys. Solids 40, 1683), and the inner problem is exactly solved by the Wiener-Hopf technique. The displacement field and scattered Rayleigh waves are found uniformly to third order in , and concluding remarks are made about the general method as well as the results presented here.     

Stability for the Electromagnetic Scattering from Large Cavities

Consider the scattering of electromagnetic waves from a large rectangular cavity embedded in the infinite ground plane. There are two fundamental polarizations for the scattering problem in two dimensions: TM (transverse magnetic) and TE (transverse electric). In this paper, new stability results for the cavity problems are established for large rectangular shape cavities in both polarizations. For the TM cavity problem, an asymptotic property of the solution and a stability estimate with an improved dependence on the high wavenumber are derived. In the TE case, the first stability result is established with an explicit dependence on the wave number.     

Scattering of guided waves by circumferential cracks in composite cylinders

A somewhat generalized numerical procedure is used in this paper to study the problem of wave scattering by circumferential cracks in composite pipes. The study is motivated by the need to develop a model for the quantitative, ultrasonic non-destructive evaluation of cracks in pipes. For this purpose, a stiffness-based Rayleigh–Ritz type approach is employed first to obtain the approximate wave numbers and wave modes. Using the wave function expansions of the incident and scattered fields in the axial direction and decomposing the problem into separate symmetric and anti-symmetric problems, a three-dimensional wave scattering problem is reduced to two, independent two-dimensional problems over the circular cross-section. Both these problems can be reduced further to quasi-one-dimensions by discretizing the cross-section into finite elements and using a transfer matrix approach in the circumferential direction. This simplification greatly reduces the computational time. A comparison of the results for an isotropic pipe demonstrates the reliability and accuracy of the modified numerical procedure. Numerical results for the reflection and transmission coefficients of different incident wave modes are also presented for a 2-ply composite pipe with a crack. The crack may have an arbitrary circumferential length and radial depth. Simple extrapolations from one wave to another wave, separately incident on a crack, are demonstrated to be impossible due to different mode conversions by the crack.     

Scattering by an anisotropic circle

The scattering by a circle is considered when the outside medium is isotropic and the inside medium is anisotropic (orthotropic). The problem is a scalar one and is phrased as a scattering problem for elastic waves with polarization out of the plane of the circle (SH wave), but the solution is with minor modifications valid also for scattering of electromagnetic waves. The equation inside the circle is first transformed to polar coordinates and it then explicitly contains the azimuthal angle through trigonometric functions. Making an expansion in a trigonometric series in the azimuthal coordinate then gives a coupled system of ordinary differential equations in the radial coordinate that is solved by power series expansions. With the solution inside the circle complete the scattering problem is solved essentially as in the classical case. Some numerical examples are given showing the influence of anisotropy, and it is noted that the effects of anisotropy are generally strong except at low frequencies where the dominating scattering only depends on the mean stiffness and not on the degree of anisotropy.     

流体饱和弹性半空间裂隙对平面P_I波的衍射

采用间接边界元法,求解了饱和半空间裂隙对平面PI波的二维衍射问题。基于单层位势理论,将边界离散并直接在边界单元上施加虚拟荷载(水平作用力、竖向作用力和流量源的叠加)以构造散射波场,并由边界条件确定虚拟荷载密度,总波场由自由波场和散射波场共同组成。通过参数分析研究了入射波频率、入射倾角、埋深、孔隙率、边界渗透条件等因素对饱和半空间中裂隙对平面PI波衍射的影响规律。结果表明:裂隙随埋深增大,地表位移谱振荡加剧,峰值有所降低;随着入射频率增加,孔隙率影响逐渐增大;垂直入射时,水平位移的放大区域主要分布在裂隙两端,斜入射时,主要集中在裂隙正上方地表附近;透水和不透水两种情况下的地表位移幅值和相位差别较小,但干土情况与饱和情况下的位移幅值相差较大。     

Scattering of harmonic waves by a disk-shaped rigid inclusion in a three-dimensional elastic matrix

We consider a three-dimensional problem on the interaction of harmonic waves with a thin rigid movable inclusion in an infinite elastic body. The problem is reduced to solving a system of two-dimensional boundary integral equations of Helmholtz potential type for the stress jump functions on the opposite surfaces of the inclusion. We propose a boundary element method for solving the integral equations on the basis of the regularization of their weakly singular kernels. Using the asymptotic relations between the amplitude-frequency characteristics of the wave farzone field and the obtained boundary stress jump functions, we determine the amplitudes of the shear plane wave scattering by a circular disk-shaped inclusion for various directions of the wave incident on the inclusion and for a broad range of wave numbers.     

Scattering of an incident wave from an interface separating two fluids

We consider the scattering of an incident pulse from an interface separating two fluids. The interface can be either an elastic membrane or a two-fluid interface with surface tension. By considering the limit where the ratio of acoustic wavelength to the surface wavelength is small, we systematically derived a boundary condition relating the scattered wave and the surface deformation. This condition is local and can be used to derive a partial differential equation for the deformation of the interface. This equation includes the contribution of the acoustic waves induced by the motion of the interface and once it is solved it can be used to determine the scattered field. At leading order in our analysis we find the plane wave approximation. The addition of the next order terms results in an on surface condition equivalent to that of Kriegsmann and Scandrett. We present numerical calculations to show that our results are in good agreement with the exact numerical solution as well as that of Kriegsmann and Scandrett. Physical situations where the conditions of our analysis are valid are presented.     

Focusing of an ultrasonic beam by a curved interface

We calculate the wavefields near the caustics, and their cusps, formed when a well collimated, ultrasonic beam scatters from a concave fluid-solid interface. The radius of curvature of the interface is assumed to be sufficiently large and the angle of incidence sufficiently small that only reflection and transmission of the beam need be considered. The incident beam is modeled as a two-dimensional wavefield whose initial profile is rectangular. The aperture is assumed to be sufficiently large (in wavelengths) that a well collimated beam is radiated, and the interface is assumed to lie close enough to the aperture that it is struck by a wavefield that has not yet evolved into a cylindrical wave. The scattered wavefields are represented as multiple integrals and are evaluated using a combination of asymptotic and numerical analysis. Special attention is given to the sometimes competing effects of the shadow boundary of each scattered beam with its corresponding caustic, and cusp if one is formed.