Mean wave propagation in a slab of one-dimensional discrete random medium

A study has been made of the propagation of time harmonic waves through a one-dimensional medium of discrete scatterers randomly positioned over a finite interval L. The random medium is modeled by a Poisson impulse process with density λ. The invariant imbedding procedure is employed to obtain a set of initial value stochastic differential equations for the field inside the medium and the reflection coefficient of the layer. By using the Markov properties of the Poisson impulse process. exact integro-differential equations of the Kolmogorov-Feller type are derived for the probability density function of the reflection coefficient and the field. When the concentration of the scatterers is low, a two variable perturbation method in small λ is used to obtain an approximate solution for the mean field. It is shown that this solution, which varies exponentially with respect to λL, agrees exactly with the mean field obtained by Feldy's approximate method.     

A variational principle for waves in discrete random media

A new variational principle is derived for the Green's function of the linear harmonic response of a scalar wave field in a discretely heterogeneous medium. The variational principle is derived directly from exact multiple scattering integral equations and is stated in terms of a certain functional of trial configuration dependent fields. The functional is found to be stationary with respect to small variations in the fields when those fields have their correct configuration dependence. A certain trial dependence in the fields is shown to generate the Lax quasicrystalline multiple scattering equations. It is furthermore clear that these equations are optimal in the sense that any more realistic form for the trial fields unavoidably generates approximate equations entailing the generally unknown high order correlation functions. At its stationary point the proposed functional takes on a physically significant value. It becomes the change in the medium admittance due to the introduction of the scatterers. In a nonlossy medium this is related to the spectral density of modes. As errors in the trial fields cancel to first order, the final evaluation of the functional at its stationary point is especially accurate. Other related functionals are proposed and discussed as well.     

Discriminating linear from nonlinear elastic damage using a nonlinear time reversal DORT method

The DORT method is a selective detection and focusing technique originally developed to detect defects and damages which induce linear changes of the elastic moduli. It is based on the time reversal (TR) where a signal collected from an array of transducers is time reversed and then back-propagated into the medium to obtain focusing on selected targets. TR is based on the principle of spatial reciprocity. Attenuation, dispersion, multiple scattering, mode conversion, etc. do not break spatial reciprocity. The presence of defects or damage, may cause materials to show nonlinear elastic wave propagation behavior that will break spacial reciprocity. Therefore the DORT method will not allow focusing on nonlinear elastic scatterers. This paper presents a new method for the detection and identification of multiple linear and nonlinear scatterers by combining nonlinear elastic wave spectroscopy, time reversal and DORT method. In the presence of nonlinear hysteretic elastic scatterers, forcing the solid with a harmonic excitation, the time reversal operator can be obtained not only at the fundamental frequency of excitation, but also at the odd harmonics. At the fundamental harmonic, either inhomogeneities and linear damages can be individually selected but only at odd harmonics nonlinear hysteretic elastic damages exist. A procedure was developed where by decomposing the operator at the odd harmonics, it was possible to focus on nonlinear scatterers and to differentiate them from the linear inhomogeneities. A complete mathematical nonlinear DORT formulation for 1 and 2D structures is presented. To model the presence of nonlinear elastic hysteretic scatterers a Preisach–Mayergoyz (PM) material constitutive model was used. Results relative to 1 and 2 dimensional structures are reported showing the capability of the method to focus and discern selectively linear and nonlinear scatterers. Furthermore, an analysis was conducted to study the influence of the number of sources and their location on the imaging process showing that using a higher numbers of sensors does not automatically bring to a minor uncoupled behaviour between the nonlinear targets.     

Theory of scattering from a nearly transparent anomaly

Summary Recent advances in the design of uhf and microwave communication equipment, coupled with the increasing need for expanding the usable frequency spectrum, have generated considerable interest in making use of scattered signals which are propagated beyond the radio horizon. Although a large number of experimental measurements have been reported in the literature, much theoretical work in interpreting the basic scattering phenomena remains to be done. In order to gain a more detailed insight into the scattering mechanism, an approximate equation of propagation of electromagnetic energy in a nearly transparent medium is applied to the case in which the medium contains a dielectric inhomogeneity in the form of an isolated gaussian-shaped perturbation in the refractive index. Equations for the scattering caused by the perturbing blob are illustrated graphically. The energy extracted from the incident wave by the blob is illustrated graphically by a plot of the total scattering cross-section as a function of blob size.     

Multichannel resonant scattering theory applied to the acoustic scattering by an eccentric elastic cylindrical shell immersed in a fluid

The multichannel resonant scattering theory is developed in the field of the acoustic scattering. Because it takes into account mode conversions, this theory deals with a nondiagonal scattering S matrix. It is used here for the study of the acoustic scattering by an elastic eccentric shell. In this case the mode conversions are due to the fact that the Lamb type waves propagating around the shell are submitted to a reflection-refraction phenomenon when passing through the thinner part of the shell. All informations on the resonances are provided by the study of Argand's diagrams. In particular, the partition of the resonant energy over all the vibration modes of the scatterer is obtained. A numerical validation of the multichannel resonant scattering theory is then given. We focus our attention on the bifurcation of resonances, and the existence of angular diagrams with an odd numbers of lobes; both are due to mode conversions. Experiments have been performed which are in agreement with the theoretical results.     

Space–time focusing of acoustic waves on unknown scatterers

Consider a propagative medium, possibly inhomogeneous, containing some scatterers whose positions are unknown. Using an array of transmit–receive transducers, how can one generate a wave that would focus in space and time near one of the scatterers, that is, a wave whose energy would confine near the scatterer during a short time? The answer proposed in the present paper is based on the so-called DORT method (French acronym for: decomposition of the time reversal operator) which has led to numerous applications owing to the related space-focusing properties in the frequency domain, i.e., for time-harmonic waves. This method essentially consists in a singular value decomposition (SVD) of the scattering operator, that is, the operator which maps the input signals sent to the transducers to the measure of the scattered wave. By introducing a particular SVD related to the symmetry of the scattering operator, we show how to synchronize the time-harmonic signals derived from the DORT method to achieve space–time focusing. We consider the case of the scalar wave equation and we make use of an asymptotic model for small sound-soft scatterers, usually called the Foldy–Lax model. In this context, several mathematical and numerical arguments that support our idea are explored.     

Scattering of obliquely incident waves by straight features in a plate

The scattering of obliquely incident waves by straight features in a plate is solved analytically. The reflection matrix of a free straight edge and the scattering matrix of a straight thickness step are obtained respectively by modal decomposition based on a real orthogonal relation. The formulas are illustrated through numerical examples. The matrices are found to be Hermitian for the propagating modes; thus, the mode conversions are reciprocal in terms of energy. The matrices can be used to determine the scattering from more complicated straight features if the cross sections are approximated as sequences of stairs and steps.     

A numerical study of super-resolution through fast 3D wideband algorithm for scattering in highly-heterogeneous media

We present a wideband fast algorithm capable of accurately computing the full numerical solution of the problem of acoustic scattering of waves by multiple finite-sized bodies such as spherical scatterers in three dimensions. By full solution, we mean that no assumption (e.g. Rayleigh scattering, geometrical optics, weak scattering, Born single scattering, etc.) is necessary regarding the properties of the scatterers, their distribution or the background medium. The algorithm is also fast in the sense that it scales linearly with the number of unknowns. We use this algorithm to study the phenomenon of super-resolution in time-reversal refocusing in highly-scattering media recently observed experimentally (Lemoult et al., 2011), and provide numerical arguments towards the fact that such a phenomenon can be explained through a homogenization theory.     

Reconstruction of multiply-scattered arrivals from the cross-correlation of waves excited by random noise sources in a heterogeneous dissipative medium

Imagining a medium composed of an arbitrary distribution of point-like heterogeneities, we study the reconstruction of scattered waves in Green's function derived from the cross-correlation function of waves excited by random noise sources of which the distribution is stationary and homogeneous. We show that the reconstruction process is intimately related to generalized forms of the optical theorem. The role of absorption in the formulation of the theorem is discussed. The reconstruction of multiply-scattered arrivals from the cross-correlation of two random wavefields is demonstrated to all orders of scattering for the simple case of two point scatterers, through application of the optical theorem for a single scatterer. In the case of N point scatterers, the cross-correlation of two Green's functions is expressed in the form of Feynman-like diagrams. The wavepaths that contribute to the reconstruction of an arbitrary multiply-scattered arrival of Green's function are identified. Repeated application of the generalized optical theorem, formulated as a diagrammatic rule, demonstrates the destructive interference between all spurious multiply-scattered arrivals.     

双介质耦合刚性基弹性层平面应变型导波模式及界面散射能量分配

过渡段动力稳定性问题已成为制约400 km/h及以上高铁路基设计的关键难题,亟需从波动和能量的角度探究由基础非均匀引发的线路系统动力响应放大机理.文章将轨下基础简化为上表面自由、底端固定的刚性基弹性层,将高铁过渡段车致弹性波传播问题提炼为非均匀介质刚性基弹性层中波的散射问题,建立双介质耦合刚性基弹性层平面应变模型,优化该类波导结构频散方程在复平面求根方法,并结合岩土类介质特征展开刚性基弹性层频散分析,以明确其多模式导波特性及散射能量分配,最后,围绕弹性层厚度、刚度比等影响因素开展对比分析.结果表明:刚性基弹性层各模式导波均具有截止频率,弹性层厚度越小,杨氏模量越大,各阶导波模式的截止频率越高;入射波在双介质刚性基弹性层发生散射后,透射场基阶模式导波会占据主体能量,随着高阶导波模式被逐一激发,反射场及透射场高阶模式能量占比会在全频率范围呈现“此消彼长”状态;交换两侧弹性层材料,改变弹性层厚度及两弹性层刚度比不会显著改变能量分布规律,但总体来看,能量更易集中在较软侧弹性层中,各模式导波在激发初始频段会更为活跃,可分配到更多能量.     

The scattering of electroelastic waves by an ellipsoidal inclusion in piezoelectric medium

The scattering problem for a single ellipsoidal piezoelectric inclusion embedded in piezoelectric medium is investigated. Based on the polarization method, the extended displacements are expressed in terms of integral equations, whose kernels are obtained from the Green’s functions of homogenous matrix. In this paper, the 3D dynamic Green’s functions are derived by means of the Radon transform technique. To illustrate the use of the equations, scattering by a piezoelectric, ellipsoidal inhomogeneity in a piezoelectric medium is considered in the low frequency and an asymptotic formula for this scattering cross-section is obtained. Numerical results of the scattering cross-sections are carried out for a spheroidal BaTiO₃-inclusion in a PZT-5H-matrix.     

缺陷岩体中弹性波传播的多次散射效应分析

弹性波在岩体中传播时与岩体缺陷相互作用形成复杂的传播图案。为研究缺陷对弹性波多次散射作用的影响,建立了双椭圆缺陷模型,基于Green函数基本解,采用边界积分的计算方法,得到了反映缺陷界面条件的刚度矩阵,分析了弹性波在双椭圆缺陷间的多次散射效应。结果表明：与单椭圆缺陷模型相比,双缺陷的相互作用使得弹性波频散和衰减效应增强,定量给出了缺陷的影响区域,从而明确了多次散射效应的尺度界限。进一步探讨了弹性波传播的多尺度效应,结果表明频散的Rayleigh峰、Mie峰和衰减的峰值频率同椭圆长轴和入射波波长两个尺度密切相关,存在明确的定量关系。相应的数值模拟结果表明,弹性波和缺陷相互作用在缺陷界面上诱发界面波,该界面波也存在频率相关性,影响了弹性波宏观传播的频散和衰减特征。     

Multiple scattering of elastic waves by pinned dislocation segments in a continuum

The coherent propagation of elastic waves in a solid filled with a random distribution of pinned dislocation segments is studied to all orders in perturbation theory. It is shown that, within the independent scattering approximation, the perturbation series that generates the mass operator is a geometric series that can thus be formally summed. A divergent quantity is shown to be renormalizable to zero at low frequencies. At higher frequencies said quantity can be expressed in terms of a cut-off with dimensions of length, related to the dislocation length, and physical quantities can be computed in terms of two parameters, to be determined by experiment. The approach used in this problem is compared and contrasted with the scattering of de Broglie waves by delta-function potentials as described by the Schrödinger equation.     

Asymptotic inverse scattering

Several inverse techniques are developed for determining the shape of an unknown scattering surface by analyzing backscattered acoustic or electromagnetic waves. These techniques are based on asymptotic high frequency representations of the fields and may be divided into three categories. The first one is the geometrical imaging method where the surface is reconstructed by means of a travel-time analysis which is here specified to the far field by utilizing Minkowski's support function. Furthermore, a geometrical method is developed for localizing edges from mid field data measured along a curve. The second category is called quasigeometrical imaging and uses geometric optics or higher order amplitudes for the reconstruction. It is shown that cross-polarized electromagnetic far field amplitude measurements permit one to deduce the complete quadratic approximation of the surface at the specular points from which the surface can be reconstructed pointwise. The third category may be subsumed under ‘asymptotic inverse scattering identities’. Here, asymptotic relations between scattered fields and distributions associated with the geometry of the scatterer are established. It is shown that the physical optics far field inverse scattering identity is only a leading order asymptotic relation but as such is also valid for non-convex scatterers. Furthermore, asymptotic inverse scattering identities are deduced which relate the singular function of a closed surface to the backscattered field data measured on a sphere enclosing the scatterer. This generalizes far field results of Cohen and Bleistein (Wave Motion 1 (1979), p. 153) to the mid field.     

Flexural wave scattering in a quarter-infinite thin plate with circular scatterers

Scattering of flexural waves by circular scatterers in a quarter-infinite thin plate is formulated using the wave expansion method together with the method of images. The scattered waves are expressed as a summation series of wave functions and the unknown scattering coefficients are determined by enforcing boundary conditions at the scatterers. Both holes and rigid scatterers are studied. Simply-supported and roller-supported boundary conditions on the quarter-infinite thin plate are also considered. The analysis can be used to determine the stress concentration caused by circular scatterers in quarter-infinite thin plates.     

Bremmer series for the multimodal sound propagation in inhomogeneous waveguides

The method of Bremmer series is presented and implemented for the multimodal sound propagation in a waveguide with varying cross-section. The solution is constructed iteratively by summing terms of increasing order, with each order representing the number of scattering events. This formulation generalizes to higher dimensions the series decomposition proposed by Bremmer in one-dimensional inhomogeneous media, and it also gives an insight on the complex wave scattering of the coupled guided modes. The accuracy and convergence of the solution are inspected and a comparison is made with the one-way coupled mode equation derived assuming no reflected waves. It is notably shown that the first-order Bremmer series, simply obtained by quadrature, is a relevant alternative to classical WKB or one-way approximations.     

Two dimensional wave problems in rotating elastic media

The rotation of an elastic medium makes it act anisotropically and dispersively. The eigenvectors for plane wave propagation are in general complex and thus the waves are elliptically polarized. In general the waves are neither pure shear nor pure compressional waves, and their speeds depend on the ratio of rotational frequency of the medium and the angular frequency of the wave.The class of problems discussed here involves waves propagating perpendicularly to the axis of rotation and in particular we discuss plane strain modes. The reflection and refraction of plane waves is considered.The plane waves are used to construct a general solution in cylindrical coordinates. The solution is given in terms of Bessel functions. The cylindrical solution is applied to scattering by circular cylinders. The problem of free oscillations is mentioned briefly.     

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.     

Scattering of water waves by vertical porous barriers: An analytical approach

An analytical approach is proposed here to study scattering of deep water waves by a submerged or a surface piercing vertical porous barrier. It involves a connection between two wave potentials of which one is the solution of a boundary value problem associated with wave scattering by the porous barrier and the other is the solution of a complementary type problem where barrier and gap positions are interchanged and solid barrier takes the position of the porous barrier. The connection also involves an auxiliary or a connection wave potential. The potential for the solid barrier problem involves incident wave forcing while the auxiliary potential describes a solid barrier type problem that involves a non-physical forcing. The solution procedure of Ursell (Ursell, 1947) is chosen to solve these boundary value problems explicitly in the case of normal wave incidence as it also determines necessary exact behavior of the potential at the barrier edge. The reflection coefficients are also connected and the reflection amplitudes of the normally incident wave against the vertical porous barriers are obtained analytically. Numerical results for reflection and transmission coefficients are presented.