Exact solution of three-dimensional acoustic field in a wedge with perfectly reflecting boundaries

An exact approach is presented to compute the three-dimensional(3D) acoustic field in a homogeneous wedge-shaped ocean with perfectly reflecting boundaries. This approach applies the Fourier synthesis technique, which reduces a 3D point-source ideal wedge problem into a sequence of two-dimensional(2D) line-source ideal wedge problems, whose analytical solution is well established. A comparison of numerical efficiency is provided between this solution and the solution proposed by Buckingham,which is obtained by a sequence of integral transforms. The details of numerical implementation of these two solutions are also given. To validate the present approach and at the same time compare numerical efficiency between this approach and Buckingham's analytical solution, two numerical examples are considered. One is the Acoustical Society of America(ASA) benchmark wedge problem and the other is a wide-angle wedge problem. Numerical results indicate that the present approach is efficient and capable of providing accurate 3D acoustic field results for arbitrary receiver locations, and hence can serve as a benchmark model for sound propagation in a homogeneous wedge-shaped ocean.     

Benchmark solutions for sound propagation in an ideal wedge

Sound propagation in a wedge-shaped waveguide with perfectly reflecting boundaries is one of the few range-dependent problems with an analytical solution, and hence provides an ideal benchmark for a full two-way solution to the wave equation. An analytical solution for the sound propagation in an ideal wedge with a pressure-release bottom was presented by Buckingham and Tolstoy [Buckingham and Tolstoy 1990 J. Acoust. Soc. Am. 87 1511]. The ideal wedge problem with a rigid bottom is also of great importance in underwater acoustics. We present an analytical solution to the ideal wedge problem with a perfectly reflecting bottom, either rigid or pressure-release, which may be used to provide a means for investigating the sound field in depth-varying channels, and to establish the accuracy of numerical propagation models. Closed-form expressions for coupling matrices are also provided for the ideal waveguides characterized by a homogeneous water column bounded by perfectly reflecting boundaries. A comparison between the analytical solution and the numerical solution recently proposed by Luo et al. [Luo W Y, Yang C M and Zhang R H 2012 Chin. Phys. Lett. 29 014302] is also presented, through which the accuracy of this numerical model is illustrated.     

On the Method of Source Images for the Wedge Problem Solution in Ocean Acoustics: Some Corrections and Appendices

In this study the method of source images for the problem of sound propagation in a penetrable wedge [G. Deane and M. Buckingham, J. Acoust. Soc. Am. 93 (1993) 1319–1328] is revisited. This solution is very important three-dimensional (3D) benchmark in computational underwater acoustics, since a wedge bounded from above by the sea surface and overlying a sloping penetrable bottom is the simplest model of a shallow-sea waveguide near the coastline. The corrected formulae for the positions of the source images and bottom images are presented together with the explanation of their derivation. The problem of branch choice in the reflection coefficient is thoroughly discussed, and the corresponding explicit formulae are given. In addition, numerical validation of the proposed branch choice schemes and the resulting wedge problem solutions are presented. Finally, source images solution is computed for a series of examples with different ratios of shear and bulk moduli in the bottom. The interplay between the acoustic-elastic waves coupling and the horizontal refraction in the wedge is demonstrated.     

Line source and site characterizations for defining the sound transmission loss of building facades

An analytical model is presented for defining the sound transmission loss of building facades exposed to noise from line sources. The model describes the non-diffuse sound field incident upon the facade in terms of both source and site parameters. The effects of facade orientation relative to the line source and the sound propagation with distance are introduced as a single term in the definition of the facade sound transmission loss. This term defines a mean angle of incidence for the exterior sound field that is equivalent to a point source location relative to a point on the facade. Numerical results are presented illustrating the magnitudes of these effects and it is shown that alternative methods for conducting field measurements of building facade sound transmission loss may be related by using this model.     

水平变化双层波导声传播问题的耦合简正波解

为了验证现有模型的精度,导出了全反射下边界双层波导中简正波耦合矩阵的解析表达式,并将其应用到全局矩阵耦合简正波模型(Direct Global Matrix Coupled-Mode)中,使得该模型可以提供水平变化双层波导问题的标准解。文中首先利用COUPLE的简正波及耦合矩阵数值解验证了该简正波及耦合矩阵解析表达式的正确性;其次,采用改进的DGMCM模型求解了双层波导海山声传播损失,结果表明,改进后的DGMCM模型可以非常精确地求解水平变化双层波导问题,可作为求解此类问题的标准模型使用。     

水平变化波导中的全局矩阵耦合简正波方法

提出了一种新的水平变化波导中声场的耦合简正波求解方法,该方法能够处理二维点源和线源问题,提供声场的双向解。该方法利用全局矩阵(DGM)一次性求解耦合模式的系数,消除了传播矩阵递推求解中存在的误差累积问题;此外,改善了现有模型中对距离函数的归一化方法,从而避免了泄露模式指数增长导致的数值溢出问题。本文还给出了绝对软海底理想波导中耦合矩阵的闭合表达式,并分析了单个阶梯下简正波耦合现象。此外,本文还计算了理想楔形波导中的声传播问题(ASA标准问题),并与解析解及COUPLE07计算结果进行了比较,结果表明该方法是一种稳定、精确的水平变化波导中的声场计算方法。      

一种水平变化可穿透波导中声传播问题的耦合简正波方法

通过利用标准简正波程序KRAKEN计算本地简正波解及耦合矩阵, 进一步发展了求解水平变化波导中声场的全局矩阵耦合简正波方法(Luo et al., "A numerically stable coupled-mode formulation for acoustic propagation in range-dependent waveguides," Sci. China-Phys. Mech. Astron. 55, 572 (2012)), 使得该方法可以处理具有可穿透海底及随深度变化声速剖面等实际问题, 并提供声场的完全双向解. 本文还给出了双层波导中耦合矩阵的解析表达式, 并利用其验证了本方法中耦合矩阵数值算法的精度. 最后, 利用改善后的全局矩阵耦合简正波模型(DGMCM)计算了美国声学学会(ASA)提出的可穿透楔形波导标准问题, 将所得数值解与参考解比较, 结果表明DGMCM方法可以精确处理水平变化波导中声传播实际问题. 关键词： 耦合简正波理论 全局矩阵方法 可穿透楔形波导     

Diffraction of sound from a dipole source near to a barrier or an impedance discontinuity

Pierce's formulation for the diffraction of spherical waves by a hard wedge has been extended to the case of the sound field due to a dipole source. The same approach is also used to extend a semiempirical model for sound propagation above an impedance discontinuity due to a dipole source. The resulting formulas have been validated by comparing their numerical solutions with that computed by summing the sound fields due to two closely spaced monopole sources of equal magnitude but opposite in phase. These new formulations are then used to develop a simple model for calculating the dipole sound field diffracted by a barrier above an impedance ground. Applications of these models relate to transportation noise prediction, particularly railway noise abatement, for which dipole sources are commonly used. The numerical predictions have been found to compare reasonably well with indoor measurements using piezoceramic transducers as dipole sources.     

Analytical solution based on the wavenumber integration method for the acoustic field in a Pekeris waveguide

An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.     

Non-line-of-sight sound source localization using matched-field processing

Acoustic diffraction allows sound to travel around opaque objects and therefore may allow beyond-line-of-sight sensing of remote sound sources. This paper reports simulated and experimental results for localizing sound sources based on fully shadowed microphone array measurements. The generic geometry includes a point source, a solid 90° wedge, and a receiving array that lies entirely in the shadow defined by the source location and the wedge. Source localization performance is assessed via matched-field (MF) ambiguity surfaces as a function of receiving array configuration, and received signal-to-noise ratio for the Bartlett and minimum variance distortionless (MVD) MF processors. Here, the sound propagation model is developed from a Green's function integral treatment. A simple 16 element line array of microphones is tested in three mutually orthogonal orientations. The experiments were conducted using an approximate 50-to-1-scaled tabletop model of a blind city-street intersection and produced ambiguity surfaces from source frequencies between 17.5 and 19 kHz that were incoherently summed. The experimental results suggest that a sound source may be localized by the MVD processor when using fully shadowed arrays that have significant aperture parallel to the edge of the wedge. However, this performance is reduced significantly for signal-to-noise ratios below 40 dB.     

Observationally constrained modeling of sound in curved ocean internal waves: examination of deep ducting and surface ducting at short range

A study of 400 Hz sound focusing and ducting effects in a packet of curved nonlinear internal waves in shallow water is presented. Sound propagation roughly along the crests of the waves is simulated with a three-dimensional parabolic equation computational code, and the results are compared to measured propagation along fixed 3 and 6 km source/receiver paths. The measurements were made on the shelf of the South China Sea northeast of Tung-Sha Island. Construction of the time-varying three-dimensional sound-speed fields used in the modeling simulations was guided by environmental data collected concurrently with the acoustic data. Computed three-dimensional propagation results compare well with field observations. The simulations allow identification of time-dependent sound forward scattering and ducting processes within the curved internal gravity waves. Strong acoustic intensity enhancement was observed during passage of high-amplitude nonlinear waves over the source/receiver paths, and is replicated in the model. The waves were typical of the region (35 m vertical displacement). Two types of ducting are found in the model, which occur asynchronously. One type is three-dimensional modal trapping in deep ducts within the wave crests (shallow thermocline zones). The second type is surface ducting within the wave troughs (deep thermocline zones).     

Noise attenuation by a hard wedge-shaped barrier

This paper is concerned with the problem of sound screening by a wedge-like barrier. The sound source is assumed to be point like, and the receiver is located in the shadow of the source sound field, so that according to geometrical optics only the field diffracted by the edge of the barrier is considered. First, for the hard wedge in space, three models are used for calculating the amplitude of the edge-diffracted field. These are the uniform theory of diffraction (UTD), the Hadden-Pierce model, both in the frequency domain, and the Biot-Tolstoy theory of diffraction which is a time domain formulation. It is first shown that even at relatively low frequencies, the frequency domain models perform quite satisfactorily as compared to the exact time domain theory. Hence, and due to its relative simplicity the UTD is proposed as an accurate calculation scheme for solving problems with edge diffraction by hard wedges. It is also proved from theoretical calculations that the amplitude of the edge-diffracted field increases for an increasing angle of the wedge, and consequently the hard half-plane gives the lowest field amplitude in the shadow zone. Some applications are then considered for evaluating the performance of a barrier on a flat ground, either completely hard or with mixed homogeneous boundary conditions. An improvement of the scheme for calculating the sound field in the all-hard case is achieved through considering the multiple diffraction, in this case only to the second order, between the top of the wedge barrier and its base. The results show that for usually occurring situations, increasing the angle of the hard wedge barrier affects negatively its efficiency through diminishing its insertion loss. These conclusions are also supported by the results of some experimental measurements conducted at a scale-model level.     

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.     

The diffraction of sound by an impedance sphere in the vicinity of a ground surface

The problem of sound diffraction by an absorbing sphere due to a monopole point source was investigated. The theoretical models were extended to consider the case of sound diffraction by an absorbing sphere with a locally reacting boundary or an extended reaction boundary placed above an outdoor ground surface of finite impedance. The separation of variables techniques and appropriate wave field expansions were used to derive the analytical solutions. By adopting an image method, the solutions could be formulated to account for the multiple scattering of sound between the sphere and its image near a flat acoustically hard or an impedance ground. The effect of ground on the reflected sound fields was incorporated in the theoretical model by employing an approximate analytical solution known as the Weyl-van der Pol formula. An approximation solution was suggested to determine the scattering coefficients from a set of linearly coupled complex equations for an absorbing sphere not too close to the ground. The approximate method substantially reduced the computational time for calculating the sound field. Preliminary measurements were conducted to characterize the acoustical properties of an absorbing sphere made of open cell polyurethane foam. Subsequent experiments were carried out to demonstrate the validity of the proposed theoretical models for various source/receiver configurations around the sphere above an acoustically hard ground and an impedance ground. Satisfactory comparative results were obtained between the theoretical predictions and experimental data. It was found that the theoretical predictions derived from the approximate solution agreed well with the results obtained by using the exact solutions.     

楔形海域中三维效应对声源距离估计的影响分析*

楔形海域中声传播会出现显著的三维效应,导致水平方向上特定接收位置的声信号出现模态多次到达或者模态影区。利用射线模型,首先对楔形波导中的传播损失进行计算,分析三维效应的存在对于声波传播的影响, 结果表明,随着距离的增加,水平折射效应愈加明显,所得的三维结果与利用二维模型计算所得相差逐渐增大。然后对比二维声场和三维声场中的接收波形,并对其进行模态分离,证明在某些接收位置将会出现明显的模态多次到达或模态影区。最后通过匹配全波解以及分离出的各号简正波所对应的脉冲波形进行声源距离估计研究,并分析三维效应对于声源距离估计结果的影响。结果表明,由于三维效应的存在,使用三维模型能够得到准确估计声源距离,而利用二维模型进行估计得到的结果与实际距离相差较多,总体上偏大。     

Three-dimensional steady state Green function for a layered isotropic plate

The elastodynamic response of a layered isotropic plate to a source point load having an arbitrary direction is studied in this paper. A decomposition technique is developed within each homogeneous isotropic lamina to simplify the general three-dimensional plane-wave propagation problem as a separate plane-strain problem and an anti-plane-wave propagation problem. The accuracy of computation is assured by cross-checking the numerical results by different methods. Results are checked numerically for a vertical point load acting on a homogeneous and a layered plate by using a hybrid method. On the other hand, results are checked for a horizontal point load by using dynamic reciprocal identities. Results are presented for both a homogeneous as well as a layered plate.     

Measurement and modeling of diffraction from an edge of a thin panel

The diffraction of sound from an edge of a thin chipboard panel was measured in an anechoic chamber, and compared to simulations based on the diffraction formulation developed by Svensson et al. [Svensson UP, Fred RI, Vanderkooy J. An analytic secondary source model of edge diffraction impulse responses. J Acoust Soc Am 1999;106(5):2331-44]. The measurements and simulations were performed for a line of receiver positions below the panel to include cases for which the direct sound had an unobstructed propagation path to the receivers, as well as cases for which the direct sound was occluded by the panel. Comparison of the measured and simulated responses is provided in both the time and frequency-domains, and shows that the differences between them are small over the entire audible frequency range. This case study verifies that the applied diffraction-modeling method gives accurate results, and that the assumptions of ideal source and wedge characteristics inherent in the method do not preclude its use in simulations of realistic scenarios.     

海洋声学中三维抛物方程非均匀网格模型

在海洋声学中,三维抛物方程模型可以有效考虑三维空间的声传播效应。然而,采用三维抛物方程模型分析三维空间内的声传播问题时,计算时间较长,并且需要消耗较大的计算机内存,因此给远距离声场的快速精确计算带来了很大困难。为此,将非均匀网格Galerkin离散化方法用于三维直角坐标系下的水声抛物方程模型中,深度算子和水平算子Galerkin离散方式由均匀网格变为非均匀网格。仿真结果表明,三维直角坐标系下非均匀网格离散的抛物方程模型,在保持计算精度、提高计算速度的同时,可以实现远距离声场的快速预报。另外,针对远距离局部海底地形与距离有关的三维声传播问题,给出了声场快速计算方法;在海底保持水平的区域,采用经典Kraken模型,重构抛物方程算法的初始场,随后依次递推求解地形与距离有关海底下的三维声场。采用改进模型,证明了远距离楔形波导声强增强效应。      

Least-squares Legendre spectral element solutions to sound propagation problems

This paper presents a novel algorithm and numerical results of sound wave propagation. The method is based on a least-squares Legendre spectral element approach for spatial discretization and the Crank-Nicolson [Proc. Cambridge Philos. Soc. 43, 50-67 (1947)] and Adams-Bashforth [D. Gottlieb and S. A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications (CBMS-NSF Monograph, Siam 1977)] schemes for temporal discretization to solve the linearized acoustic field equations for sound propagation. Two types of NASA Computational Aeroacoustics (CAA) Workshop benchmark problems [ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics, edited by J. C. Hardin, J. R. Ristorcelli, and C. K. W. Tam, NASA Conference Publication 3300, 1995a] are considered: a narrow Gaussian sound wave propagating in a one-dimensional space without flows, and the reflection of a two-dimensional acoustic pulse off a rigid wall in the presence of a uniform flow of Mach 0.5 in a semi-infinite space. The first problem was used to examine the numerical dispersion and dissipation characteristics of the proposed algorithm. The second problem was to demonstrate the capability of the algorithm in treating sound propagation in a flow. Comparisons were made of the computed results with analytical results and results obtained by other methods. It is shown that all results computed by the present method are in good agreement with the analytical solutions and results of the first problem agree very well with those predicted by other schemes.     

矩形表面波探头声场的高斯声束叠加法

利用矩形压电晶片和有机玻璃楔块折射可激励出超声表面波,广泛用于固体近表面缺陷检测和材料特性测量.由于描述表面波三维声场的理论方法还鲜有报道,因而主要采用简化的表面波二维声场模型来定量分析这类问题.高斯声束模型近些年被广泛应用于解决超声体波传播的各种复杂问题,然而,目前还没有将其扩展应用到超声表面波的声场的计算中.通过结合表面波格林方程和矩形换能器的高斯声束模型,推导出基于高斯声束叠加的表面波三维声场解析解.进一步,将该方法与点源叠加的数值解进行了分析比较,计算结果表明表面波声场的高斯声束叠加方法在具有较好计算精度的同时,还具有更快的计算效率.