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.     

Broadband optimization of a pentamode-layered spherical acoustic waveguide

Transformational acoustics offers the theoretical possibility of cloaking obstacles within fluids, provided metamaterials having continuously varying bulk moduli and densities can be found or constructed. Realistically, materials with the proper, continuously varying anisotropies do not presently exist. Discretely layered cloaks having constant material parameters within each layer are a viable alternative, but due to their discrete nature, may become ineffective outside of narrow frequency ranges. Because of such limitations, there is interest in finding discretely layered systems that can be effective in as wide as possible bandwidth without the need for unrealizable material properties within each layer. The present work introduces a novel methodology for finding optimal material parameters for use in such layered cloaks. In principle, the technique could be applied to any acoustic or electromagnetic scattering problem, but for purposes of demonstration, this paper considers a fluid-loaded acoustically hard sphere with a cloak that comprised layered pentamodes, whose material properties are constrained to lie within reasonable ranges relative to the density and bulk modulus of water.     

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.     

Elastodynamic behavior of mechanical cloaks designed by direct lattice transformations

Steering waves in elastic solids is more demanding than steering waves in electromagnetism or acoustics. As a result, designing material distributions which are the counterpart of optical invisibility cloaks in elasticity poses a major challenge. Waves of all polarizations should be guided around an obstacle to emerge on the downstream side as though no obstacle were there. Recently, we have introduced the direct-lattice-transformation approach. This simple and explicit construction procedure led to extremely good cloaking results in the static case. Here, we transfer this approach to the dynamic case, i.e., to elastic waves or phonons. We demonstrate broadband reduction of scattering, with best suppressions exceeding a factor of five when using cubic coordinate transformations instead of linear ones. To reliably and quantitatively test these cloaks efficiency, we use an effective-medium approach.     

Antiplane elastic wave cloaking using metamaterials,homogenization and hyperelasticity

We consider the problem of how to cloak objects from antiplane elastic waves using two alternative techniques. The first is the use of a layered metamaterial in the spirit of the work of Torrent and Sanchez-Dehesa (2008) who considered acoustic cloaks, motivated by homogenization theories, whilst the second is the use of a hyperelastic cloak in the spirit of the work of Parnell et al. (2012). We extend the hyperelastic cloaking theory to the case of a Mooney–Rivlin material since this is often considered to be a more realistic constitutive model of rubber-like media than the neo-Hookean case studied by Parnell et al. (2012), certainly at the deformations required to produce a significant cloaking effect. Although not perfect, the Mooney–Rivlin material appears to be a reasonable hyperelastic cloak. This is clearly encouraging for applications. We quantify the effectiveness of the various cloaks considered by plotting the scattering cross section as a function of frequency, noting that this would be zero for a perfect cloak.     

大延伸非均匀介质中地震波全弹性散射理论II-弹性波多次散射理论

推出大延伸非均匀连续介质的弹性波能量传输表达式,并以此为基础建立并发展了相应的弹性动力学能量传递理论.构造了一个基于非均匀薄层或非均匀相屏单次散射迭代法的多次散射模型.该模型既适用于弱散射,也适用于强散驰 既适用于普通散射,也适用于转换散射;尽管高频情况下只考虑普通散射及前向散射.应用该模型计算了弹性多次散射的能通量,处理了散射衰减问题。数值实验的结果表明,短周期地震图上的尾波主要来源于S波散射。     

Cloaking of a circular cylindrical elastic inclusion from antiplane elastic waves and resonance effects

Cloaking of a circular cylindrical elastic inclusion embedded in a homogeneous linear isotropic elastic medium from antiplane elastic waves is studied. The transformation or change-of-variables method is used to determine the material properties of the cloak and the homogenization theory of composites is used to construct a multilayered cloak consisting of many bi-material cells. The large system of algebraic equations associated with this problem is solved by using the concept of multiple scattering with wave expansion coefficient matrices. Numerical results for cloaking of an elastic inclusion and a rigid inclusion are compared with the case of a cavity. It is found that while the cloaking patterns for the three cases are similar, the major difference is that standing waves are generated in the elastic inclusion and the multilayered cloak cannot prevent the motion inside the elastic inclusion, even though the cloak seems nearly perfect. Waves can penetrate into and cause vibrations inside the elastic inclusion, where the amplitude of standing waves depend on the material properties of the inclusion but are very much reduced when compared to the case when there is no cloak. For a prescribed mass density, the displacements inside the elastic cylinder decrease as the shear modulus increases. Moreover, the cloaking of the elastic inclusion over a range of wavenumbers is also investigated. There is significant low frequency scattering even if the cloak consists of a large number of layers. When the wavenumber increases, the multilayered cloak is not effective if the cloak consists of an insufficient number of layers. Resonance effects that occur in cloaking of elastic inclusions are also discussed.     

Disturbance evolution in the leading-edge region of a blunt flat plate in supersonic flow

The spatio-temporal dynamics of small disturbances in viscous supersonic flow over a blunt flat plate at freestream Mach number M_∞=2.5 is numerically simulated using a spectral approximation to the Navier–Stokes equations. The unsteady solutions are computed by imposing weak acoustic waves onto the steady base flow. In addition, the unsteady response of the flow to velocity perturbations introduced by local suction and blowing through a slot in the body surface is investigated. The results indicate distinct disturbance/shock-wave interactions in the subsonic region around the leading edge for both types of forcing. While the disturbance amplitudes on the wall retain a constant level for the acoustic perturbation, those generated by local suction and blowing experience a strong decay downstream of the slot. Furthermore, the results prove the importance of the shock in the distribution of perturbations, which have their origin in the leading-edge region. These disturbance waves may enter the boundary layer further downstream to excite instability modes.     

Transport equations for elastic and other waves in random media

We derive and analyze transport equations for the energy density of waves of any kind in a random medium. The equations take account of nonuniformities of the background medium, scattering by random inhomogeneities, polarization effects, coupling of different types of waves, etc. We also show that diffusive behavior occurs on long time and distance scales and we determine the diffusion coefficients. The results are specialized to acoustic, electromagnetic, and elastic waves. The analysis is based on the governing equations of motion and uses the Wigner distribution.     

Analytical extension of Finite Element solution for computing the nonlinear far field of ultrasonic waves scattered by a closed crack

Nonlinear scattering of ultrasonic waves by closed cracks subject to contact acoustic nonlinearity (CAN) is determined using a 2D Finite Element (FE) coupled with an analytical approach. The FE model, which includes unilateral contact with Coulomb friction to account for contact between crack faces, provides the near-field solution for the interaction between in-plane elastic waves and a crack of different orientations. The numerical solution is then analytically extended in the far-field based on a frequency domain near-to-far field transformation technique, yielding directivity patterns for all linear and nonlinear components of the scattered waves. The proposed method is demonstrated by application to two nonlinear acoustic problems in the case of tone-burst excitations: first, the scattering of higher harmonics resulting from the interaction with a closed crack of various orientations, and second, the scattering of the longitudinal wave resulting from the nonlinear interaction between two shear waves and a closed crack. The analysis of the directivity patterns enables us to identify the characteristics of the nonlinear scattering from a closed crack, which provides essential understanding in order to optimize and apply nonlinear acoustic NDT methods.     

Design and analysis of nonlinear-transformation-based broadband cloaking for acoustic wave propagation

A coordinate-transformation method can be used to design invisibility cloaks for many types of waves, including acoustic waves. The traditional method for designing a cloak depends on a transformation from a virtual space to a physical space. Previous acoustic cloaks that are mainly designed with linear-transformation-based acoustics have drawbacks that acoustic wave trajectories in the cloaks cannot be controlled and tuned. This work uses a nonlinear mapping from a ray trajectory perspective to construct acoustic cloaks with tunable non-singular material properties. Use of a ray trajectory equation is a straightforward and alternate way to study propagation characteristics of different types of waves, which allows more flexibility in controlling the waves. A broadband cylindrical cloak for acoustic waves in an inviscid fluid is realized with layered non-singular, homogeneous, and isotropic materials based on a nonlinear transformation. Some advantages and improvements of the invisibility nonlinear-transformation cloak over a traditional linear-transformation cloak are analyzed. The invisibility capability of the nonlinear-transformation cloak can be tuned by adjusting a design parameter that is shown to have influence on the acoustic wave energy flowing into the region inside the cloak. Numerical examples show that the nonlinear-transformation cloak is more effective for making a domain undetectable by acoustic waves in an inviscid fluid and shielding acoustic waves from outside the cloak than the linear-transformation cloak in a broad frequency range. The methodology developed here can be used to design nonlinear-transformation cloaks for other types of waves.     

圆弧形凹陷地形表面覆盖层对入射平面P波的影响

利用Fourier-Bessel级数展开法给出了表面具有覆盖层的圆弧形凹陷地形对入射平面P波散射问题的一个解析解，并利用该解分析了不同形状凹陷地形表面覆盖层刚度和厚度对入射P波的影响。结果表明，凹陷地形表面覆盖层的存在，即使厚度很薄，对入射P波的散射也具有很大覆盖，覆盖层刚度和厚度的变化可显著改善凹陷地形场地的动力特性。     

Normal incidence of a plane acoustic wave on a rigid wall covered with a thin compressible layer

The one-dimensional unsteady problem of the variation of the pressure on a rigid wall covered with a thin compressible layer upon which a plane acoustic wave impinges is investigated. The investigation is carried out from two standpoints: without allowance for wave processes in the layer (in this case the layer is modeled by means of a special boundary condition [1] and the pressure on the wall is a continuous function of time) and with allowance for the waves transporting the pressure perturbation from the outer edge of the layer to the wall and back (in this case the pressure on the wall is a piecewise-continuous function of time). A criterion of the proximity of the results of the two approaches is the smallness of the acoustic impedance ratio before interaction begins. This holds true even when the high intensities of the incident waves lead to considerable compression of the layer and an increase in its acoustic impedance.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 139–148, July–August, 1988.     

Theory of scattering from multilayered bodies of arbitrary shape

Green 's functions and boundary integral equation methods are used to derive a matrix set of equations for scattering from a multilayered homogeneous elastic body embedded in an infinite elastic material. The surfaces separating the layers have arbitrary shape. The formalism for the three-layer material is derived in detail and generalized to N-layers. A matrix factorization method (MFM) is shown to considerably simplify the computational problem. The relation to the problems of acoustic waves in fluids and electromagnetic waves in a dielectric material is briefly indicated.     

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.     

Exterior cloaking with active sources in two dimensional acoustics

We cloak a region from a known incident wave by surrounding the region with three or more devices that cancel out the field in the cloaked region without significantly radiating waves. Since very little waves reach scatterers within the cloaked region, the scattered field is small and the scatterers are for all practical purposes undetectable. The devices are multipolar point sources that can be determined from Green's formula and an addition theorem for Hankel functions. The cloaking devices are exterior to the cloaked region.     

微槽-吸气组合控制平板边界层多模态稳定性研究

边界层转捩会使高超声速飞行器壁面摩阻和热流显著增加,因此在高超声速飞行器设计过程中往往占据重要地位.针对高超声速飞行器多模态转捩控制问题,提出了微槽道(1 mm)与边界层吸气的组合控制方法,并通过直接数值模拟和线性稳定性理论研究了Ma=4.5平板边界层的稳定性及组合控制效果.边界层在无控状态时,同时存在失稳的第一、二模态波,且二维第二模态波最不稳定;单纯施加微槽道控制时,边界层第二模态波会被抑制但第一模态波会被略微激发.对比而言,采用“微槽-吸气”组合控制后,不仅增强了对第二模态波的抑制效果,而且减弱了第一模态波的激发程度;同时随着吸气强度的增加,第二模态波不稳定区域明显收缩、频率显著增高,而第一模态波则变化不明显.相较于单纯的微槽道,吸气增强了“微槽吸收”与“声波散射”作用,因此中等吸气强度下该组合控制方法对第一和第二模态波的增长率分别实现了12.63%和28.02%的抑制效果.以上结果表明“微槽-吸气”组合控制手段具有适用宽频、布置区域灵活的优点,展现出了一定的多模态控制效果.     

功能梯度磁电弹耦合结构中的间隙波

研究了由均匀磁电弹半空间和功能梯度磁电弹层组成的耦合结构中间隙波的传播特性.假定功能梯度层的材料性能沿厚度方向呈指数变化,且其表面机械自由,但承受两种电磁边界条件.首先推导了频散方程,然后结合数值算例分析了功能梯度层材料性能的梯度变化、厚度及电磁边界条件对相速度的影响,结果对功能梯度磁电弹材料在声波器件中的应用具有参考价值.      

GAS WAVE IN FUNCTIONALLY GRADED MAGNETO-ELECTRO-ELASTIC COUPLED STRUCTURE

研究了由均匀磁电弹半空间和功能梯度磁电弹层组成的耦合结构中间隙波的传播特性.假定功能梯度层的材料性能沿厚度方向呈指数变化,且其表面机械自由,但承受两种电磁边界条件.首先推导了频散方程,然后结合数值算例分析了功能梯度层材料性能的梯度变化、厚度及电磁边界条件对相速度的影响,结果对功能梯度磁电弹材料在声波器件中的应用具有参考价值.     

Resonance of Feshbach-type and explosive instability of magnetoelastic waves in solids

In this paper, we theoretically describe a magnetoelastic system with an explosive instability in which the magnetic energy can be converted into the mechanical one via an extremely efficient channel. The principle is based on the amplification of acoustic waves traveling in an antiferromagnetic crystal under the action of transversal electromagnetic pumping. It is known that such parametric interaction can produce an exponential acoustic wave amplification. Our finding consists in the addition of another pumping channel by means of a discrete resonant acoustic mode whose impact is similar to the Feshbach resonance observed in another physical system (ultra-cold gazes). The addition of the second pumping mechanism results in the explosive instability having the dynamics of a singularity at a finite time instead of the exponential growth. In our example the traveling waves are Lamb waves in an antiferromagnetic plate, and the additional pumping is a shear resonance. We establish equations governing the considered three-phonon parametric instability, theoretically analyze the instability conditions, and give a relevant numerical example illustrating the explosive magneto-acoustic dynamics. It is shown that the explosive scenario can occur with a very low signal level i.e. Lamb waves amplitudes comparable to spontaneous noise in the system.