水下掩埋目标的散射声场计算与实验

水下掩埋目标声散射问题是识别和探测掩埋目标的理论基础, 是声散射研究领域的热点问题. 本文基于射线声学推导了掩埋情况下目标声散射计算的格林函数近似式, 并在此基础上进一步给出了相应的远场积分公式. 在有限元方法的基础上, 将推导得到的公式写入有限元仿真软件, 对软件功能进行拓展, 构建二维轴对称目标的声散射模型, 并计算掩埋情况下弹性实心球在不同条件下的目标强度, 获得了其散射声场随频率、掩埋深度、沙层吸收系数等参数的变化规律. 开展实心球的自由空间和浅掩埋条件下水池声散射实验, 利用共振隔离技术处理实验数据, 提取目标声散射的纯弹性共振特征进行分析, 结果表明可将其用于掩埋目标识别和探测. 最后利用总散射声场与理论计算结果进行对比, 验证了理论仿真的正确性.     

Realization of a Gravity-Capillary Wave Surface Using Sub-Nyquist Spectral Samples

A method for numerical realization of a full gravity-capillary wave surface that is specified by a wave-spectral model and a dispersive relationship of the surface wave is developed on the basis of the angular spectral representation of a random water surface. A significant aspect of the method is that it requires a smaller number of spectral samples than that required by the Whitteker-Shannon sampling theorem for the complete generation of a full gravity-capillary wave surface, without resulting in any appreciable errors in auto-correlation functions for the surface displacement or surface slopes. The method enables the unified treatment of gravity and capillary waves in numerical studies of higher order characteristics of the thermal radiation emitted from, and the light scattered by, the wave surface.     

Exact wave field simulation for finite-volume scattering problems

An exact boundary condition is presented for scattering problems involving spatially limited perturbations of arbitrary magnitude to a background model in generally inhomogeneous acoustic media. The boundary condition decouples the wave propagation on a perturbed domain while maintaining all interactions with the background model, thus eliminating the need to regenerate the wave field response on the full model. The method, which is explicit, relies on a Kirchhoff-type integral extrapolation to update the boundary condition at every time step of the simulation. The Green's functions required for extrapolation through the background model are computed efficiently using wave field interferometry.     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.     

Gauss声束对离轴椭圆柱的声辐射力矩

利用部分波展开法求解得到了Gauss声束入射下刚性和非刚性椭圆柱的声散射系数,推导了一般情况下的声辐射力矩表达式.在此基础上,通过一系列数值仿真详细分析了离轴距离、入射角度和束腰半径对声辐射力矩的影响.结果表明:正向与负向声辐射力矩均可以在一定条件下存在;低频情况下刚性椭圆柱比非刚性椭圆柱更容易产生较强的声辐射力矩;特定频率的入射声场可以激发出非刚性椭圆柱不同阶的共振散射模式,因而非刚性椭圆柱的声辐射力矩峰值与频率的关系更密切;增加束腰半径有利于扩大散射截面,进而增加椭圆柱的声辐射力矩.该研究结果预期可以为利用声辐射力矩实现粒子的可控旋转和流体黏度的反演提供一定的理论指导.     

Modeling of Coulomb coupling and acoustic wave propagation in LiNbO3

The angular spectrum propagation technique is applied to modeling of wave propagation through piezoelectric media. The calculations of the angular spectrum propagator are based on the relevant equation for the slowness surface resulting from the solution of the Christoffel equation with piezoelectrically stiffened elastic constants. A two-dimensional FFT algorithm is applied in the final field superposition. We concentrate on the case of Coulomb coupling through local electrical point contacts on both the excitation and detection side. To model that case we superpose solutions for acoustic Green's functions of different propagation modes convoluted with equivalent distributed effective sources. Calculated results are in good agreement with the measured ones.     

雪层覆盖地面电磁散射的矩量法

为实现分层介质粗糙面电磁散射的矩量法研究,给出一种分层介质粗糙面电磁积分方程的区域分解方法.将格林定理应用于粗糙面所分的各子空间,结合波动方程和格林函数推导分层粗糙面的电磁积分方程,利用矩量法对其进行离散,数值计算得到雪层覆盖地面散射系数的角分布曲线,其中,粗糙表面由一维带限Weierstrass分形谱和Monte Carlo方法模拟.通过与时域有限差分法数值结果的比对,验证该方法的准确性,并分析散射系数随雪和地面参数、介质参数以及入射波参数的变化,获得了较完整的散射特征.     

Modeling of wave dispersion along cylindrical structures using the spectral method

Algorithm and code are presented that solve dispersion equations for cylindrically layered media consisting of an arbitrary number of elastic and fluid layers. The algorithm is based on the spectral method which discretizes the underlying wave equations with the help of spectral differentiation matrices and solves the corresponding equations as a generalized eigenvalue problem. For a given frequency the eigenvalues correspond to the wave numbers of different modes. The advantage of this technique is that it is easy to implement, especially for cases where traditional root-finding methods are strongly limited or hard to realize, i.e., for attenuative, anisotropic, and poroelastic media. The application of the new approach is illustrated using models of an elastic cylinder and a fluid-filled tube. The dispersion curves so produced are in good agreement with analytical results, which confirms the accuracy of the method. Particle displacement profiles of the fundamental mode in a free solid cylinder are computed for a range of frequencies.     

分层吸声结构的声学设计与性能分析

本文给出了声波在多层介质中的传播方程，通过数值计算，分析了各层参数对分层吸声结构声学性能的影响规律。用参数匹配方法，对多层吸声结构各层材料的参数进行设计。结果表明，通过合理分布各层材料，多层吸声结构可以在较宽频段上达到满意的吸声性能。     

A Green's function method for surface acoustic waves in functionally graded materials

Acoustic wave propagation in anisotropic media with one-dimensional inhomogeneity is discussed. Using a Green's function approach, the wave equation with inhomogeneous variation of elastic property and mass density is transformed into an integral equation, which is then solved numerically. The method is applied to find the dispersion relation of surface acoustic waves for a medium with continuously or discontinuously varying elastic property and mass density profiles.     

利用矢量海洋环境噪声提取声场格林函数

考虑到矢量水听器在垂直方向上具有8字形指向性,能够有效抑制远方非平稳噪声源的干扰,提出了一种矢量环境噪声相关函数(NCF)提取声场时域格林函数(TDGF)的方法。基于简正波理论建立了声压和垂直振速垂直相关性模型。在此基础上,给出了声压和垂直振速相关函数提取声场纵向格林函数的过程.数值仿真对比和实验数据分析表明,相对于声压提取方法,垂直振速提取方法能够有效消除直达波前出现的亮纹与亮区干扰。此外,对于同等时间长度噪声序列,声压提取方法只提取到直达波路径,而垂直振速提取方法还提取到了我们更为关心的海底反射路径。利用直达波与海底反射波到达时延差估计的海深与实测海深吻合较好。      

驻波声场中悬浮临界密度及稳定性研究

本文以声场中物体为研究对象,理论上得到行波和驻波场中的声辐射压力方程.在驻波声场中引入临界悬浮密度概念,可作为物体能否在非线性声场中悬浮的判据,同时给出谐振腔移动速度的最大范围.更进一步,以实验参数作为数值计算的输入来指导实验,并结合实验结果讨论了驻波声场中样品密度和大小、发射面和反射面形状以及两者之间的距离、反射面的尺寸等因素对物体悬浮稳定性的影响,发现当物体尺寸和密度确定时,调控好谐振腔的长度,增加波腹处的声压是提升声悬浮稳定性的有效手段.     

Enhanced scattering of acoustic waves at interfaces

We propose a general method to realize a total scattering of an incident acoustic wave at interfaces between different media while allowing the flow of air, fluids and/or particles. This originates from the enlargement of the equivalent acoustic scattering cross section of an embedded object coated with acoustic metamaterials, which causes the coated object to behave as a scatterer bigger than its physical size. We theoretically design a model circular cylindrical object coated with such metamaterials whose properties are determined according to two different, but identical, methods. The desired function is confirmed for both far-field and near-field cases with full wave simulations based on the finite element method. This work reveals a promising way to achieve noise shielding and naval camouflage.     

矢量拖线阵水听器流噪声响应特性

针对传统拖线阵流噪声理论的局限性, 建立了完善的矢量拖线阵流噪声理论分析方法, 可全面准确地揭示矢量拖线阵流噪声响应特性. 基于细长圆柱的湍流边界层压力起伏Carpenter模型, 采用波数-频率谱分析方法对矢量拖线阵流噪声响应特性进行了理论研究, 导出了圆柱形矢量水听器流噪声响应的声压和振速自功率谱及其互功率谱的解析表达式, 定量分析了流噪声响应功率谱与拖曳速度、水听器尺寸、套管尺寸和材料等参数之间变化规律; 另外, 还讨论了圆柱形矢量水听器偏离护套轴线时矢量拖线阵流噪声响应, 导出了流噪声响应的声压、径向和轴向振速自功率谱及其互功率谱的解析表达式, 数值计算结果表明: 轴线偏移距离对声压和轴向振速的高频噪声的影响要大于对低频噪声的影响, 而对径向振速的全频段噪声都有明显影响, 且对振速分量影响要远大于对声压影响.     

In-plane vibration of thin elliptic plates submitted to uniform pulsed microwave irradiations

The technique of acoustic generation by microwave excitation in structures is applied here to study the in-plane vibration of full or hollowed elliptic plates. The absorption of pulsed microwave irradiations by a material causes a sudden rise of its temperature and the generation of an acoustic wave by thermoelastic effect. A semi-analytic theoretical model is developed to predict the in-plane displacement fields in elliptic thin plates submitted to a uniform temperature rise. It is assumed that the isotropic and viscoelastic plate constitutive material is submitted to a thermoelastic excitation under a plane stress state. The wave equations that govern the Helmholtz displacement potentials are resolved in an elliptic cylindrical system of coordinates by means of infinite angular and radial Mathieu functions series. The displacement field is finally obtained by taking into account the zero stress conditions on the boundaries of the plates. The comparison between the theoretical and the experimental responses of full and hollowed elliptic plates shows a good agreement that permits the validation of the developed model.     

Elastic wave field computation in multilayered nonplanar solid structures: a mesh-free semianalytical approach

Multilayered solid structures made of isotropic, transversely isotropic, or general anisotropic materials are frequently used in aerospace, mechanical, and civil structures. Ultrasonic fields developed in such structures by finite size transducers simulating actual experiments in laboratories or in the field have not been rigorously studied. Several attempts to compute the ultrasonic field inside solid media have been made based on approximate paraxial methods like the classical ray tracing and multi-Gaussian beam models. These approximate methods have several limitations. A new semianalytical method is adopted in this article to model elastic wave field in multilayered solid structures with planar or nonplanar interfaces generated by finite size transducers. A general formulation good for both isotropic and anisotropic solids is presented in this article. A variety of conditions have been incorporated in the formulation including irregularities at the interfaces. The method presented here requires frequency domain displacement and stress Green's functions. Due to the presence of different materials in the problem geometry various elastodynamic Green's functions for different materials are used in the formulation. Expressions of displacement and stress Green's functions for isotropic and anisotropic solids as well as for the fluid media are presented. Computed results are verified by checking the stress and displacement continuity conditions across the interface of two different solids of a bimetal plate and investigating if the results for a corrugated plate with very small corrugation match with the flat plate results.     

Surface response of a viscoelastic medium to subsurface acoustic sources with application to medical diagnosis

The response at the surface of an isotropic viscoelastic medium to buried fundamental acoustic sources is studied theoretically, computationally and experimentally. Finite and infinitesimal monopole and dipole sources within the low audible frequency range (40-400 Hz) are considered. Analytical and numerical integral solutions that account for compression, shear and surface wave response to the buried sources are formulated and compared with numerical finite element simulations and experimental studies on finite dimension phantom models. It is found that at low audible frequencies, compression and shear wave propagation from point sources can both be significant, with shear wave effects becoming less significant as frequency increases. Additionally, it is shown that simple closed-form analytical approximations based on an infinite medium model agree well with numerically obtained "exact" half-space solutions for the frequency range and material of interest in this study. The focus here is on developing a better understanding of how biological soft tissue affects the transmission of vibro-acoustic energy from biological acoustic sources below the skin surface, whose typical spectral content is in the low audible frequency range. Examples include sound radiated from pulmonary, gastro-intestinal and cardiovascular system functions, such as breath sounds, bowel sounds and vascular bruits, respectively.     

Spectral Green's Functions for Layered Gyrotropic Structures Through a Transmission-Line Approach

In this work, a straightforward method is shown to derive spectral Green's functions in classes of planar structures including gyrotropic media, based on a recently proposed transmission-line formalism. In particular, the analysis is expressly carried out for a current line source in layered configurations with transversely-magnetized ferrites, useful for microwave and millimeter-wave waveguiding and radiating devices.     

Modal vibrations of a cylindrical radiator over an impedance plane

The problem of acoustic radiation from an infinite cylinder undergoing harmonic modal surface vibrations near a locally reacting planar boundary is considered. The formulation utilizes the appropriate wave field expansions, the classical method of images, and the translational addition theorem for cylindrical wave functions, along with a simple local surface reaction model involving a complex amplitude wave reflection coefficient applied to simulate the relevant boundary conditions for the given configuration. The analytical results are illustrated with a numerical example in which the cylindrical surface is immersed near a layer of fibrous material set on an impervious rigid wall. The numerical results reveal the important effects of interface local surface reaction and source position on the computed modal impedance component values and the radiated on-axis far-field pressure. The benchmark solution presented can lead to a better understanding of acoustic radiation from near-interface two-dimensional sources, which are commonly encountered problems in outdoor acoustics and noise control engineering. Eventually, it could be used to validate those found by numerical approximation techniques.     

Scholte Wave at Air-Metal Interface Generated by a Pulsed Disc-Like Source

The normal displacement and pressure of Scholte and Leaky Rayleigh waves at air-metal interface generated by a pulsed disc-like source are simulated theoretically by the Cagniard-de Hoop method and studied by laser ultrasound technique experimentally. It is found that the Scholte wave detected by a photorefractive interferometer is mainly contributed by the surface pressure and the Leaky Rayleigh wave is dominated by the surface displacement. It is also proven that the pulse width of these interface waves is mainly determined by the acoustic time delay on the generating source size under our experimental conditions.