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.     

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.     

New results for isotropic point scatterers: Foldy revisited

Since 1945, Foldy’s method has been used to predict velocity and attenuation for various types of scatterers. In this paper, it is shown that Foldy’s method also yields predictions of reflection and transmission of scalar waves by a random distribution of point or line scatterers contained in a slab. Results are given in two and three dimensions, and for normal and oblique incidences. Formulae are also obtained for the reflection and transmission of longitudinal waves by point or line scatterers distributed in an elastic (non-viscous) fluid. Energy equations are derived, and expressions are obtained for the energy dissipated in the slab on average over one period. Curves for the reflection and transmission coefficients are presented in the case of solid cylindrical bars immersed in a fluid. The results obtained in this paper are expected to be valid for a low density of scatterers. Potential applications of this work occur in ultrasonic evaluation of materials, seismic exploration and medical ultrasonics, where reflected (or backscattered) data are used to construct maps or images of the materials (metals, composites, earth subsurface, tissue). The formulae of this work are expected to provide useful tools for better and more efficient mapping or imaging.     

Elastic wave scattering by a three-dimensional inhomogeneity in an elastic half space

In this paper we will consider scattering of elastic waves in a half space. The half space is an isotropic, linear and homogeneous medium except for a finite inhomogeneity. The T-matrix method (also called the “extended boundary condition method” or “null field approach”) is extended to derive expressions for the elastic field inside the half space and the surface field on the interface. The assumptions on the source that excites the half space are fairly weak. In the numerical applications found in this paper we assume a Rayleigh surface wave to be the incoming field, and we only compute the surface displacements. We make illustrations on some simple types of scatterers (spheres and spheroids; the latter ones can be arbitrarily oriented).     

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.     

Asymptotic approximations for Bloch waves and topological mode steering in a planar array of Neumann scatterers

We study the canonical problem of wave scattering by periodic arrays, either of infinite or finite extent, of Neumann scatterers in the plane; the characteristic lengthscale of the scatterers is considered small relative to the lattice period. We utilise the method of matched asymptotic expansions, together with Fourier series representations, to create an efficient and accurate numerical approach for finding the dispersion curves associated with Floquet–Bloch waves through an infinite array of scatterers. The approach lends itself to direct scattering problems for finite arrays and we illustrate the flexibility of these asymptotic representations on topical examples from topological wave physics.     

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.     

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.     

A novel design scheme for acoustic cloaking of complex shape based on region partitioning and multi-origin coordinate transformation

Acoustic cloaking is an important application of acoustic metamaterials. This article proposes a novel design scheme for acoustic cloaking based on the region partitioning and multi-origin coordinate transformation. The cloaked region is partitioned into multiple narrow strips. For each strip, a local coordinate system is established with the local origin located at the strip center, and a coordinate transformation in the local coordinate system is conducted to squeeze the material along the strip length direction to form the cloaked region. To facilitate the implementation of the acoustic cloak, the multilayer effective medium is used to approximate the non-uniform anisotropic material parameters. The effectiveness of the proposed coordinate transformation method is verified by comparing the results from our method with those in the literature. Firstly, the results of a circular acoustic cloak in the literature are reproduced by using our finite element (FE) simulations for validation. Then, a comparison is made between the traditional coordinate transformation scheme and our new scheme for simulating an elliptical acoustic cloak. The results indicate that the proposed multi-origin coordinate transformation method has a better cloaking effect on the incident wave along the ellipse minor axis direction than the traditional method. This means that for the same object, an appropriate transformation scheme can be selected for different incident wave directions to achieve the optimal control effect. The validated scheme is further used to design an arch-shaped cloak composed of an upper semicircular area and a lower rectangular area, by combining the traditional single-centered coordinate transformation method for the semicircular area and the proposed multi-origin method for the rectangular area. The results show that the designed cloak can effectively control the wave propagation with significantly reduced acoustic pressure level. This work provides a flexible acoustic cloak design method applicable for arbitrary shapes and different wave incident directions, enriching the theory of acoustic cloaking based on coordinate transformation.     

Transformation elastodynamics and cloaking for flexural waves

The paper addresses an important issue of cloaking transformations for fourth-order partial differential equations representing flexural waves in thin elastic plates. It is shown that, in contrast with the Helmholtz equation, the general form of the partial differential equation is not invariant with respect to the cloaking transformation. The significant result of this paper is the analysis of the transformed equation and its interpretation in the framework of the linear theory of pre-stressed plates. The paper provides a formal framework for transformation elastodynamics as applied to elastic plates. Furthermore, an algorithm is proposed for designing a broadband square cloak for flexural waves, which employs a regularised push-out transformation. Illustrative numerical examples show high accuracy and efficiency of the proposed cloaking algorithm. In particular, a physical configuration involving a perturbation of an interference pattern generated by two coherent sources is presented. It is demonstrated that the perturbation produced by a cloaked defect is negligibly small even for such a delicate interference pattern.     

Analytical solutions for crack-like scatterers and sources in isotropic elastic plates

A common theoretical framework can be employed for modelling both active and passive approaches to structural health monitoring based on guided waves. The same fundamental solutions can be used to represent the scattered field due to crack-like damage in an active approach, or the acoustic emission due to various failure mechanisms in a passive approach. It is shown that the reciprocal theorem can be used to derive compact analytical formulae for those fundamental solutions expressed as sums over the propagating modes. The only detailed calculations that are required are those involved in determining the relevant dispersion curves and through-thickness mode functions. For clarity, the approach is first described for two-dimensional configurations. The subsequent extension to three-dimensional requires the stationary phase approximation. It is shown that the modal coefficients for crack-like sources (or scatterers) are proportional to the conjugate stress component of the corresponding mode. Furthermore, the approach can also be used for an edge crack or delamination in a semi-infinite plate, and the theoretical formulae are compared with experimentally validated computational results for that case. Further extensions are briefly noted.     

Sphere scattering of nonlinearly interacting acoustic waves

Sphere scattering of the field of nonlinearly interacting plane acoustic waves when the sphere is located in the region of nonlinear interaction between the primary pumping waves of a parametric antenna is considered. An analytic expression for the secondary field pressure at the difference frequency is obtained. This expression describes the process of nonlinear interaction of the incident and scattered waves. The secondary-field total pressure components, which characterize the interaction between the incident plane waves and scattered spherical waves are analyzed. The numerical results and experimental data are given.Taganrog. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 4–12, March–April, 1995.     

Turbulence in broken lee waves

In this paper, the waves' breaking in the lee waves in successfully simulated by the atmospheric mesoscale numerical model with a second-order turbulent closure. It is further proved that the turbulence in the wave-breaking region plays the role of intense mixing for the average field, which leads to the trapping of upward propagating waves and thus promotes the development of the downslope wind. The turbulent structure in the wave-breaking region is discussed and the following conclusions are obtained: (1) In the wave-breaking region, the turbulent heat fluxes transfer from inside to outside and the turbulent momentum fluxes transfer from outside to inside. (2) In the wave-breaking region, the turbulent energy mainly comes from the wind shear and the buoyancy promotes the turbulent development only in part of the region. (3) In the upper part of the wave-breaking region, the turbulent momentum fluxes behave as a counter-gradient transfer. (4) The turbulent mixing in the wave-breaking region is non-local. recommended by Professor Li Jiachun.     

Analysis of nonlinear acoustoelastic effect of surface acoustic waves in laminated structures by transfer matrix method

The propagation of surface acoustic waves in a layered half-space is investigated in this paper, where a thin cubic Ge film is perfectly bonded to an isotropic elastic Si half-space. Application of the transfer matrix and by solving the coupled field equations, solutions to the mechanical displacements are obtained for the film and elastic substrate, respectively. The phase velocity equations for surface acoustic waves are obtained. Effects of the homogeneous initial stresses induced by the mismatch of the film and substrate are discussed in detail. The results are useful for the design of acoustic surface wave devices.     

Effect of viscosity on the propagation of nonlinear Alfvén waves along a tangential magnetohydrodynamic discontinuity in an incompressible fluid

It is proposed to consider the propagation of surface waves along a tangential magnetohydrodynamic discontinuity in the particular case where the fluid velocities on both sides of the interface are equal to zero. In [1] it was shown that waves called surface Alfvén waves may be propagated along the surface separating a semi-infinite region without a field from a region with a uniform magnetic field. The linear theory of surface Alfvén waves in a compressible medium was considered in [2]. In [3] the damping of surface Alfvén waves as a result of viscosity and heat conduction was investigated. The propagation of low-amplitude nonlinear surface Alfvén waves in an incompressible fluid in the absence of dissipative processes is described by the integrodifferential equation obtained in [4]. By means of a numerical solution of this equation it was shown that a perturbation initially in the form of a sinusoidal wave will break. The breaking time was determined. In this paper the equation derived in [4] is extended to the case of a viscous fluid. It is shown that the equation obtained does not have steady-state solutions. The propagation of periodic disturbances is investigated numerically. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 94–104, November–December, 1986. The author wishes to thank L. S. Fedorov for assisting with the calculations.     

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.     

初应力对压电层状结构声表面波传播性能的影响

研究了压电层状结构中初应力对广义Ｒａｙｌｅｉｇｈ波传播相速度和机电耦合性能的影响,通过求解含初应力的运动微分方程,对自由界面电学开路和短路两种情况得到了相应的相速度方程。给出了具体的数值算例,所得结果对于提高和改善声表面波器件性能有参考意义。     

Coupled magnetoacoustic waves in a conducting paramagnetic fluid

We consider the propagation of small disturbances in a paramagnetic conducting fluid in a uniform constant magnetic field. Because of coupling of the mechanical and magnetic effects, coupled magnetoacoustic oscillations of a wave nature develop in a certain (resonant) frequency region. The usual MHD waves and uniform magnetization oscillations occur far from resonance. Dissipative processes are accounted for.The equations of motion for a conducting paramagnetic fluid in which interaction of the hydrodynamic velocity with the magnetization and the magnetic field was taken into account phenomenologically were obtained in [1], One of the consequences of this interaction is the propagation of coupled magnetoelastic waves in the fluid; this phenomenon is examined in the present paper.     

Wilton ripples between two uniform streaming magnetic fluids

An analysis is made of the small-amplitude capillary-gravity waves which occur on the interface of two incompressible inviscid magnetic fluids of different densities. The waves arise as a result of second harmonic resonance. The fluids moving with uniform velocities parallel to their interface are stressed by an oblique magnetic field. The linear relations between the oblique magnetic field and the instability criteria of the linear waves are analyzed. At the stability region (away from the neutral curve) of the linear theory, a pair of coupled non-linear partial differential equations are presented. On the neutral curve, a pair of coupled non-linear partial differential equations are introduced. The last pair of equations may be regarded as the counterparts of the single Klein-Gordon equation which occurs in the non-resonant case. In all cases, the wave profile and its stability conditions are obtained. These conditions are discussed analytically and graphically.     

Reduction of the sonic boom level by heating the flow in front of the body

A numerical study of the possibility of reducing the sonic boom level in the case of local heat release to a supersonic gas flow at Mach number equal to 2 ahead of a body is described. The computations are performed for a spherical heat supply zone located on the flight trajectory ahead of the tip of an axisymmetric thin body. For the numerical study the combined method of “phantom bodies” is used. Different magnitudes of heat supply to the incoming flow are tested. These calculations are performed with allowance for interaction of shock waves emanating from the heated gas region and from the body in the far field. The computational results show that the local heat supply to a supersonic gas flow ahead of a body can reduce the sonic boom level by more than 20 %. The reduction of the sonic boom level is ensured by changing the free-stream parameters ahead of the body and by preventing the coalescence of shock waves from the heat supply zone and from the body in the far field.