Numerical simulation of parametric sound generation and its application to length-limited sound beam

In this study we propose a simulation model for predicting the nonlinear sound propagation of ultrasound beams over a distance of a few hundred wavelengths, and we estimate the beam profile of a parametric array. Using the finite-difference time-domain method based on the Yee algorithm with operator splitting, axisymmetric nonlinear propagation was simulated on the basis of equations for a compressible viscous fluid. The simulation of harmonic generation agreed with the solutions of the Khokhlov–Zabolotskaya–Kuznetsov equation around the sound axis except near the sound source. As an application of the model, we estimated the profiles of length-limited parametric sound beams, which are generated by a pair of parametric sound sources with controlled amplitudes and phases. The simulation indicated a sound beam with a narrow truncated array length and a width of about one-quarter to half that of regular a parametric beam. This result confirms that the control of sound source conditions changes the shape of a parametric beam and can be used to form a torch like low-frequency sound beam.     

A beam tracing method for interactive architectural acoustics

A difficult challenge in geometrical acoustic modeling is computing propagation paths from sound sources to receivers fast enough for interactive applications. This paper describes a beam tracing method that enables interactive updates of propagation paths from a stationary source to a moving receiver in large building interiors. During a precomputation phase, convex polyhedral beams traced from the location of each sound source are stored in a "beam tree" representing the regions of space reachable by potential sequences of transmissions, diffractions, and specular reflections at surfaces of a 3D polygonal model. Then, during an interactive phase, the precomputed beam tree(s) are used to generate propagation paths from the source(s) to any receiver location at interactive rates. The key features of this beam tracing method are (1) it scales to support large building environments, (2) it models propagation due to edge diffraction, (3) it finds all propagation paths up to a given termination criterion without exhaustive search or risk of under-sampling, and (4) it updates propagation paths at interactive rates. The method has been demonstrated to work effectively in interactive acoustic design and virtual walkthrough applications.     

一种预测扩散声场在非规则形状刚性壁面附近干涉图样的方法

扩散声场会在反射边界附近形成干涉图样,研究方法包括平面波模型、简正模态分析、渐进模态分析等,但仅适用于尺度远大于声波波长的矩形声腔。提出一种预测扩散声场在非规则刚性壁面结构附近形成的干涉图样的数值方法,表明结构附近“受挡”声压的互谱矩阵取决于:(1)假定该结构在自由空间中振动辐射声音时其表面法向振速到表面及场点声压的边界元系数矩阵;(2)假定结构置于自由空间中且表面刚性时,点声源辐射声波入射到结构表面上产生的散射声场的边界元系数矩阵;(3)扩散声场均方声压。仿真表明,该途径预测的干涉图样与理论值完全吻合。该预测方法还可用于混响环境下声源附近直达声压均方值的空间分布估计,为混响环境下设备的声源定位提供帮助。      

Calculation of sound radiant efficiency and sound radiant modes of arbitrary shape structures by BEM and general eigenvalue decomposition

In this paper, a numerical method is presented to calculate sound radiant efficiency and radiant modes of arbitrary shape structures. Some methods have been proposed to compute sound radiant efficiencies and sound radiant modes of plates and beams. However, there is not a valid method to calculate for arbitrary shape structures except for measurement at the present time. The method proposed can predicate the sound radiant efficiencies and the sound radiant modes for arbitrary shape structures by boundary element method (BEM) and general eigenvalue decomposition. The validity of this method is demonstrated by two numerical examples of pulsating sphere and radiation cube.     

Acousto-optic interaction of diverging Gaussian beams in anisotropic media

A strong acousto-optic interaction of bounded light and sound beams of a Gaussian shape is considered for different geometries permitting long-term interaction. The anisotropic spreading of an acoustic beam in the course of its propagation is taken into account. The spectra of light beams are described by a Gaussian-polynomial expression, and a set of differential equations is derived which allows one to describe various geometries of acousto-optic diffraction. It is demonstrated that the transmission function of an acousto-optic cell essentially depends on the ratio between the dimensions of the light and sound beams and on the angle between the propagation direction of the diffracted beam and that of the incident beam. The degree of spreading of the acoustic beam noticeably affects the suppression of the transmission side lobes.     

基于光传播Monte Carlo模拟的三维模型特征提取方法

在现代工业设计、人工智能、软件设计等领域,三维模型正展现强劲需求与活力,传统三维模型特征提取方法仅能提取模型表面特征,难以满足复杂模型特征提取需求,为提高三维模型的特征提取精度,基于光谱分析以及可见光传播特性,提出一种具有高区分度的三维模型特征提取方法。首先利用光散射系数、吸收系数、各向异性等光学特性参数,量化分析光在不同介质中透射、散射、反射概率,并确定最佳光谱模拟波段;其次采用Monte-Carlo法模拟光子束在三维模型中的传播历程,获得光子束传播轨迹的角度、距离、能量等多种统计量,计算不同统计量权重,经过统计分析后完成特征提取;然后在ESB国际通用三维模型库中,测试不同光子束数量、约束空间形状对特征提取效果影响,从而确定最佳模拟参数;最后使用多种特征提取方法与光传播模拟法进行特征提取效果比较,并采用查准率、查全率以及E测度评价指标对特征提取效果进行定量测试。实验结果表明,三维模型特征提取准确度对光传播约束空间形态较为敏感,光子传播的最佳约束空间为球体;三维特征提取效率随着光子束数量升高而降低,在保证特征提取精度的前提下,10000~25000区间是光子束模拟数量的最佳取值范围;基于光传播模拟的特征提取准确度高于小波变换、距离夹角以及D2分布方法,能够满足三维模型检索需求,更适合复杂三维模型的离线特征提取与应用。结合光谱分析与可见光传播特性的模型特征提取法拓宽了光谱分析的应用范围,能够提取出融合三维模型表面特性与内部形态的有效特征,为高精度特征提取技术研究注入新动力。     

Incorporating source directionality into outdoor sound propagation calculations

Many outdoor sound sources, such as aircraft or ground vehicles, exhibit directional radiation patterns. However, long-range sound propagation algorithms are usually formulated for omnidirectional point sources. This paper describes two methods for incorporating directional sources into long-range sound propagation algorithms. The first is the equivalent source method (ESM), which determines a compact distribution of omnidirectional point sources reproducing a given directivity pattern in the far field. This method can be used with any propagation algorithm because it explicitly reconstructs a source function as a set of point sources with certain amplitudes and positions. The second is a directional starter method (DSM), which is developed specifically for the parabolic equation (PE) algorithms. This method derives narrow- or wide-angle directional starter fields, corresponding to a given source directivity pattern, without reconstructing the equivalent source distribution. Although the ESM can also be used for the PE, the DSM is simpler and can be more convenient, especially if the sound propagation is calculated only for one or a few azimuthal directions. While these two methods are found to produce generally distinct starter fields, they nonetheless yield identical directivity patterns.     

深海海底斜坡环境下的声传播

海底地形变化对声传播具有很大影响,在南海深海区域海底斜坡环境下进行了一次声传播实验,实验显示倾斜海底环境下声传播损失出现了一些不同于平坦海底环境下的现象,分析并解释了海底地形变化对产生声传播差异的原因.结果表明,海底斜坡对声波的反射增强作用可使斜坡上方的声传播损失减少约5 d B.当声波第一次入射到达的海底位置有较小幅度的山丘(凸起高度小于1/10海深)时,海底小山丘即可对声波有反射遮挡作用,导致在其反射区特定传播距离和深度上出现倒三角声影区,比平坦海底环境下相同影区位置处的传播损失增大约8 d B,影响深度可达海面以下1500 m.而海底斜坡对声波的反射阻挡作用使得从海面反射及水体向下折射的会聚区结构消失,只剩下从水体向上折射的会聚结构.因此,海底地形对深海声传播影响较大,在水下目标探测和性能评估等应用中应予以重视.     

Noise attenuation in factories

The cost of reducing noise levels in factories by covering large surfaces with sound absorbents is high. It is therefore important to be able to calculate in advance the effectiveness of absorbents and to determine how absorbents may be chosen and distributed for maximal noise reduction for the invested capital. For this purpose a mathematical model of sound propagation and attenuation in factories has been developed on certain simplifying assumptions. The interrelationship between the different parameters is found to be rather intricate and the mathematical model must be evaluated using a computer program. The influence on the noise levels of sound scattering objects, the shape and size of the factory and absorbents on the room surfaces have been studied. Especially if the ceiling is high, acoustic baffles are found to be more efficient than absorbents on the ceiling in reducing noise levels.     

Model equation for strongly focused finite-amplitude sound beams

A model equation that describes the propagation of sound beams in a fluid is developed using the oblate spheroidal coordinate system. This spheroidal beam equation (SBE) is a parabolic equation and has a specific application to a theoretical prediction on focused, high-frequency beams from a circular aperture. The aperture angle does not have to be small. The theoretical background is basically along the same analytical lines as the composite method (CM) reported previously [B. Ystad and J. Berntsen, Acustica 82, 698-706 (1996)]. Numerical examples are displayed for the amplitudes of sound pressure along and across the beam axis when sinusoidal waves are radiated from the source with uniform amplitude distribution. The primitive approach to linear field analysis is readily extended to the case where harmonic generation in finite-amplitude sound beams becomes significant due to the inherent nonlinearity of the medium. The theory provides the propagation and beam pattern profiles that differ from the CM solution for each harmonic component.     

Propagation-independent fields

Nondiffracting beams are solutions of the wave equation which keep their shape during free propagation. Because of this property, they are widely studied, both as mechanical and electromagnetic waves. Another scheme, however, can be devised in order to find beams that, even if do not maintain their shape upon propagation, are independent on the propagation distance, apart from a phase factor. In this sense it can be said that these beams are also nondiffracting. In this work, the scheme leading to this different kind of nondiffracting fields, which we term propagation-independent beams, is presented. Closed-form expressions for a few of these fields, for 1-D sources, are found. A method, which in principle is able to produce beams of this type, is given. This also clarifies some of their properties. Finally, the item of finite realizations of these beams, which in principle require infinite-dimension sources, is analyzed through numerical simulations. It is also shown that the practical feasibility of these fields poses some limitations on the bandwidth they can present.     

Modeling the propagation of nonlinear three-dimensional acoustic beams in inhomogeneous media

A three-dimensional model of the forward propagation of nonlinear sound beams in inhomogeneous media, a generalized Khokhlov-Zabolotskaya-Kuznetsov equation, is described. The Texas time-domain code (which accounts for paraxial diffraction, nonlinearity, thermoviscous absorption, and absorption and dispersion associated with multiple relaxation processes) was extended to solve for the propagation of nonlinear beams for the case where all medium properties vary in space. The code was validated with measurements of the nonlinear acoustic field generated by a phased array transducer operating at 2.5 MHz in water. A nonuniform layer of gel was employed to create an inhomogeneous medium. There was good agreement between the code and measurements in capturing the shift in the pressure distribution of both the fundamental and second harmonic due to the gel layer. The results indicate that the numerical tool described here is appropriate for propagation of nonlinear sound beams through weakly inhomogeneous media.     

非傍轴矢量异常空心光束的传输特性

 将Cai提出的异常空心光束的理论模型推广到非傍轴范畴,推导出非傍轴矢量异常空心光束传输的解析表达式,用以研究它在自由空间中的传输特性。研究表明,异常空心光束在传输过程中光束形状会发生变化。与高斯光束不同,非傍轴异常空心光束傍轴近似成立条件依赖于传输距离,这与异常空心光束光强分布随传输距离的变化有关。非傍轴异常空心光束远场的光束质量可用桶中功率来描述,桶中功率随f参数(波长与束腰宽度的比值)的增大而减小。     

Efficient calculation of two-dimensional periodic and waveguide acoustic Green's functions

New representations and efficient calculation methods are derived for the problem of propagation from an infinite regularly spaced array of coherent line sources above a homogeneous impedance plane, and for the Green's function for sound propagation in the canyon formed by two infinitely high, parallel rigid or sound soft walls and an impedance ground surface. The infinite sum of source contributions is replaced by a finite sum and the remainder is expressed as a Laplace-type integral. A pole subtraction technique is used to remove poles in the integrand which lie near the path of integration, obtaining a smooth integrand, more suitable for numerical integration, and a specific numerical integration method is proposed. Numerical experiments show highly accurate results across the frequency spectrum for a range of ground surface types. It is expected that the methods proposed will prove useful in boundary element modeling of noise propagation in canyon streets and in ducts, and for problems of scattering by periodic surfaces.     

Development of a hybrid computer model for simulating the complicated virtual sound field in enclosures

A hybrid computer model, SOFIS, has been developed for the simulation of an enclosed virtual sound field in an arbitrary shaped enclosure. It can be used to calculate the impulse response and acoustical parameters of different positions in a virtual enclosure. This paper describes the way in which SOFIS models the sound source, the receiver and the sound propagation throughout the enclosed space with curved surfaces or barriers. A phase tracing method and the calculation of acoustic indexes are also discussed in this paper.     

Comparison of two approximate methods for hard-edged diffracted flat-topped light beams

The integral resulted in an infinite series of Bessel functions and expanding a hard aperture into a complex-Gaussians shape are proposed as two methods for studying the propagation properties of the hard-edged diffraction flat-topped light beam. Using the two methods, the corresponding analytical propagation equations of flat-topped light beams through a circular apertured ABCD optical system are obtained. Some numerical calculations and comparative analyses by using the two methods and the diffraction integral formulae are made. It is shown that the first method of an infinite series of Bessel functions is superior to the second of expanding a hard aperture function into a complex-Gaussians shape at the aspect of calculation accuracy, but the second method is superior to the first method at the aspect of the improvement in the calculation efficiency.     

The extended Fourier pseudospectral time-domain method for atmospheric sound propagation

An extended Fourier pseudospectral time-domain (PSTD) method is presented to model atmospheric sound propagation by solving the linearized Euler equations. In this method, evaluation of spatial derivatives is based on an eigenfunction expansion. Evaluation on a spatial grid requires only two spatial points per wavelength. Time iteration is done using a low-storage optimized six-stage Runge-Kutta method. This method is applied to two-dimensional non-moving media models, one with screens and one for an urban canyon, with generally high accuracy in both amplitude and phase. For a moving atmosphere, accurate results have been obtained in models with both a uniform and a logarithmic wind velocity profile over a rigid ground surface and in the presence of a screen. The method has also been validated for three-dimensional sound propagation over a screen. For that application, the developed method is in the order of 100 times faster than the second-order-accurate FDTD solution to the linearized Euler equations. The method is found to be well suited for atmospheric sound propagation simulations where effects of complex meteorology and straight rigid boundary surfaces are to be investigated.     

Numerical modeling of the exterior-to-interior transmission of impulsive sound through three-dimensional,thin-walled elastic structures

Exterior propagation of impulsive sound and its transmission through three-dimensional, thin-walled elastic structures, into enclosed cavities, are investigated numerically in the framework of linear dynamics. A model was developed in the time domain by combining two numerical tools: (i) exterior sound propagation and induced structural loading are computed using the image-source method for the reflected field (specular reflections) combined with an extension of the Biot–Tolstoy–Medwin method for the diffracted field, (ii) the fully coupled vibro-acoustic response of the interior fluid–structure system is computed using a truncated modal-decomposition approach. In the model for exterior sound propagation, it is assumed that all surfaces are acoustically rigid. Since coupling between the structure and the exterior fluid is not enforced, the model is applicable to the case of a light exterior fluid and arbitrary interior fluid(s). The structural modes are computed with the finite-element method using shell elements. Acoustic modes are computed analytically assuming acoustically rigid boundaries and rectangular geometries of the enclosed cavities. This model is verified against finite-element solutions for the cases of rectangular structures containing one and two cavities, respectively.     

浅海涌浪对表面声道声传播的影响

受海面强风和海-气相互作用影响,表面声道普遍存在于冬季海洋环境中,是一种天然有利于声传播的波导.但是海面波浪使得海表形成粗糙界面,会严重破坏这种优良性能.本文利用南海北部海区的一次冬季声传播实验数据,研究表面声道声传播特性.研究表明,海底底质对表面声道内声传播的影响较弱,当海面风较小时,涌浪造成的影响为主要原因.实验数据显示,考虑涌浪后的粗糙海面给70km远处带来了10dB的传播损失增长.因此在考察南海北部海区冬季声场特性时,不仅要考虑海面风浪的影响,更需要考虑周围海域传来的涌浪的影响.研究涌浪存在时的声传播特性对提升声纳设备在海况较差时的使用性能具有重要意义.