Wave field synthesis of a sound field described by spherical harmonics expansion coefficients

Near-field compensated higher order Ambisonics (NFC-HOA) and wave field synthesis (WFS) constitute the two best-known analytic sound field synthesis methods. While WFS is typically used for the synthesis of virtual sound scenes, NFC-HOA is typically employed in order to synthesize sound fields that have been captured with appropriate microphone arrays. Such recorded sound fields are essentially represented by the coefficients of the underlying surface spherical harmonics expansion. A sound field described by such coefficients cannot be straightforwardly synthesized in WFS. This is a consequence of the fact that, unlike in NFC-HOA, it is critical in WFS to carefully select those loudspeakers that contribute to the synthesis of a given sound source in a sound field under consideration. In order to enable such a secondary source selection, it is proposed to employ the well-known concept of decomposing the sound field under consideration into a continuum of plane waves, for which the secondary source selection is straightforward. The plane wave representation is projected onto the horizontal plane and a closed form expression of the secondary source driving signals for horizontal WFS systems of arbitrary convex shape is derived.     

Scattering of a plane wave by a nonlinear load array on a perfectly conducting shield covered by a dielectric layer

The problem of plane wave scattering by a nonlinear load array on a perfectly conducting shield covered by a homogeneous dielectric layer is solved by the integral equation method. The influence of the coating parameters on the array reflection coefficients at the fundamental and third harmonics of the incident wave are studied numerically. Taganrog Radiotechnical University, Taganrog, Russia Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 43, No. 6, pp. 545–551, June, 2000.     

Double pass wave structure function in weak to strong optical turbulence

The authors have recently published one-way horizontal path wave structure function expressions for plane, spherical, and Gaussian-beam waves. Those results are based on a modification of the Rytov approximation method and, by comparison with experimental data, are believed to be valid under all irradiance fluctuation conditions. In this paper, the earlier results are extended to the double pass case in which the received wave is reflected from a smooth target. The point target case is considered for incident plane, spherical, and Gaussian-beam waves. The finite target case is considered for an incident spherical wave with reflection from a plane mirror or a retroreflector. Comparisons are made between the results derived here, the standard Rytov approximation results, and the asymptotic results given by Banakh for conditions of strong irradiance fluctuations.     

Double pass wave structure function in weak to strong optical turbulence

The authors have recently published one-way horizontal path wave structure function expressions for plane, spherical, and Gaussian-beam waves. Those results are based on a modification of the Rytov approximation method and, by comparison with experimental data, are believed to be valid under all irradiance fluctuation conditions. In this paper, the earlier results are extended to the double pass case in which the received wave is reflected from a smooth target. The point target case is considered for incident plane, spherical, and Gaussian-beam waves. The finite target case is considered for an incident spherical wave with reflection from a plane mirror or a retroreflector. Comparisons are made between the results derived here, the standard Rytov approximation results, and the asymptotic results given by Banakh for conditions of strong irradiance fluctuations.     

Interaction of a plane progressive sound wave with a functionally graded spherical shell

An exact analysis is carried out to study interaction of a time-harmonic plane progressive sound field with a radially inhomogeneous thick-walled elastic isotropic spherical shell suspended in and filled with compressible ideal fluid mediums. Using the laminated approximation method, a modal state equation with variable coefficients is set up in terms of appropriate displacement and stress functions and their spherical harmonics. Taylor’s expansion theorem is then employed to obtain the solution to the modal state equation ultimately leading to calculation of a global transfer matrix. Numerical example is given for a water-submerged/air-filled Aluminum/Zirconia elastic spherical sandwich shell containing a functionally graded interlayer and subjected to an incident progressive plane sound wave. The mechanical properties of the interlayer are assumed to vary smoothly and continuously across the thickness with the change of volume concentration of its constituents. The effect of incident wave frequency, thickness and compositional gradient of the interlayer on the form function amplitude and the average radiation force acting on the composite shell are examined. Limiting cases are considered and fair agreements with well-known solutions are established.     

一种利用激波进行落点测量的新方法

目前用于确定弹刃或导弹落点的声测方法多数是基于这样一个假定：其所发出的声波是一个平面波。但当其下落速度大于声速时,这个假定就不成立了,这时产生了一个锥面激波＾（４）,本文介绍了一种利用正方形基阵对其进行测量的方法。     

单层传声器阵列信号空间重采样的声波分离方法

为了重构非自由声场中目标声源的声场响应,提出单层传声器阵列信号空间重采样的声波分离方法.以球面波函数为基函数,建立由系列球面波函数叠加表达的声场数学模型.基于近场声全息原理,利用单层传声器阵列面上空间重采样形成的两组声压测量信号,求解基函数系数,并重构出传声器阵列两侧声源各自的声场响应,实现声波分离.使用脉动球和振动球共同作用的非自由声场,检验了数学模型以及传声器信号信噪比、传声器阵列形状和面积、声源中心位置、频率等关键参数对声波分离精度的影响,并在全消声室内进行了实验验证.最后,对单层传声器阵列重采样的声波分离方法的实施给出了建议.     

Acoustic centering of sources measured by surrounding spherical microphone arrays

The radiation patterns of acoustic sources have great significance in a wide range of applications, such as measuring the directivity of loudspeakers and investigating the radiation of musical instruments for auralization. Recently, surrounding spherical microphone arrays have been studied for sound field analysis, facilitating measurement of the pressure around a sphere and the computation of the spherical harmonics spectrum of the sound source. However, the sound radiation pattern may be affected by the location of the source inside the microphone array, which is an undesirable property when aiming to characterize source radiation in a unique manner. This paper presents a theoretical analysis of the spherical harmonics spectrum of spatially translated sources and defines four measures for the misalignment of the acoustic center of a radiating source. Optimization is used to promote optimal alignment based on the proposed measures and the errors caused by numerical and array-order limitations are investigated. This methodology is examined using both simulated and experimental data in order to investigate the performance and limitations of the different alignment methods.     

Asymptotic fourier analysis of a “sawtooth like” wave for dual frequency source excitation

An exact weak shock zone solution was found previously for Burgers' equation for plane waves. An approximate solution for spherical waves also was obtained. Both of these solutions held for dual frequency excitation at the source. In the work reported in the present paper, an asymptotic formula was derived which yields harmonic amplitudes at a point beyond the “sawtooth like” zone. The formula's predicted harmonic amplitudes have been found to agree well with numerical results for a specific dual frequency case. Excellent agreement also was obtained with the well-known solution for a monochromatic source. The new formula predicts that for certain frequency ratios of the primary signals the amplitude of the parametrically generated difference signal will exceed that of a directly projected signal at the same frequency and with the same total input power. This result holds for infinite planar and spherical source geometries.The formula should be useful in estimating beam axis values of the harmonic amplitudes for practical source geometries by applying the plane wave version in the near field and the spherical version in the far field. Such information can be of value to those who are constructing mathematical models of parametric array behavior.     

Transmission of a spherical sound wave through a single-leaf wall: Mass law for spherical wave incidence

This paper examines the sound insulation of a single-leaf wall driven by a spherical wave. The transmitted sound field of an infinite elastic plate under a spherical wave incidence is theoretically analyzed and insulation mechanisms are considered. The displacement of the plate is formulated using the Hankel transform in wavenumber space and the transmitted sound pressure in the far-field is obtained by Rayleigh’s formula in an explicit closed form. Moreover, a reduction index is also derived in a closed form by introducing an approximation into the vibration characteristics of the plate. Deterioration of the insulation performance under the spherical wave incidence is caused by an apparent decrease of wall impedance that depends on the directivity of the transmitted sound wave. The mass law for a spherical wave incidence is different from that for a normal plane wave incidence: doubling the weight of the wall or the frequency gives an increase of 3 dB (c.f. 6 dB for a normal plane wave incidence), which is also smaller than the field incidence mass law.     

Two-particle scattering on the lattice: Phase shifts,spin-orbit coupling,and mixing angles

We determine two-particle scattering phase shifts and mixing angles for quantum theories defined with lattice regularization. The method is suitable for any non-relativistic effective theory of point particles on the lattice. In the center-of-mass frame of the two-particle system we impose a hard spherical wall at some fixed large radius. For channels without partial-wave mixing the partial-wave phase shifts are determined from the energies of the nearly spherical standing waves. For channels with partial-wave mixing further information is extracted by decomposing the standing wave at the wall boundary into spherical harmonics, and we solve coupled-channels equations to extract the phase shifts and mixing angles. The method is illustrated and tested by computing phase shifts and mixing angles on the lattice for spin-1/2 particles with an attractive Gaussian potential containing both central and tensor force parts.     

Lateral spatial resolution of widened focused transducer array for optoacoustic tomography

Numerical investigation of a point spread function for a optoacoustic transducer array is described. Analysis of the minimal reconstructed dimension for an optoacoustic image of a point source of spherical waves is performed within the proposed approach. The influence of the array geometrical parameters, number of array transducers, and the frequency band of a single transducer on the lateral resolution in the image plane is investigated in detail. It is demonstrated that the lateral resolution obtained with the help of a transducer array in the image plane is determined unambiguously by the frequency band of a transducer, the flare angle of the array, and the transducer width and does not depend on the number of transducers.     

A spectral method for the wave equation of divergence-free vectors and symmetric tensors inside a sphere

The wave equation for vectors and symmetric tensors in spherical coordinates is studied under the divergence-free constraint. We describe a numerical method, based on the spectral decomposition of vector/tensor components onto spherical harmonics, that allows for the evolution of only those scalar fields which correspond to the divergence-free degrees of freedom of the vector/tensor. The full vector/tensor field is recovered at each time-step from these two (in the vector case), or three (symmetric tensor case) scalar fields, through the solution of a first-order system of ordinary differential equations (ODE) for each spherical harmonic. The correspondence with the poloidal–toroidal decomposition is shown for the vector case. Numerical tests are presented using an explicit Chebyshev-tau method for the radial coordinate.     

Plane wave interpolation in direct collocation boundary element method for radiation and wave scattering: numerical aspects and applications

The classical boundary element formulation for the Helmholtz equation is rehearsed, and its limitations with respect to the number of variables needed to model a wavelength are explained. A new type of interpolation for the potential is then described in which the usual boundary element shape functions are modified by the inclusion of a set of plane waves, propagating in a range of directions. This is termed the plane wave basis boundary element method. The modifications needed to the classical procedures, in terms of integration of the element matrices, and location of collocation points are described. The well-known Singular Value Decomposition solution technique, which is adopted here for the solution of the system matrix equation in its complex form, is briefly outlined. The conditioning of the system matrix is analysed for a simple radiation problem. The corresponding diffraction problem is also analysed and results are compared with analytical and classical boundary element solutions. The CHIEF method is adopted to enhance the quality of the solution, particularly in the vicinity of irregular frequencies. The plane wave basis boundary element method is then applied to two problems: scattering of plane waves by an elliptical cylinder and the multiple circular cylinder plane wave scattering problem. In both cases results are compared with analytical solutions. The results clearly demonstrate that the new method is considerably more efficient than the classical approach. For a given number of degrees of freedom, the frequency for which accurate results can be obtained, using the new technique, can be up to three or four times higher than that of the classical method. This makes the method a powerful new addition to our tools for tackling high-frequency radiation and scattering problems.     

单层传声器阵列重建相干声场的近场声全息方法研究

针对传声器阵列两侧存在相干声源的非自由声场重建问题,提出基于球面谐波函数扩展近场声全息理论的相干声场重建方法。该方法在已知测量面两侧声源几何位置时,使用单层传声器阵列获取测量面处的声压分布,通过最小二乘法获得与目标声源和干扰噪声源响应对应的最优球波函数扩展项数和最优系数向量,结合测点位置的空间坐标进行声波分解,并分别重建出各声源在测量面上的声压分布。为了验证方法的有效性,分别给出了相干噪声源为球形声源和非球形声源的仿真验证,并在全消声室内对双扬声器产生的相干声场的重建进行了实验验证。结果表明:该方法对球形声源和非球形声源干扰下的声场重建都具有较好的效果,球形声源干扰下的重建精度更高。      

Wave envelopes method for description of nonlinear acoustic wave propagation

A novel, free from paraxial approximation and computationally efficient numerical algorithm capable of predicting 4D acoustic fields in lossy and nonlinear media from arbitrary shaped sources (relevant to probes used in medical ultrasonic imaging and therapeutic systems) is described. The new WE (wave envelopes) approach to nonlinear propagation modeling is based on the solution of the second order nonlinear differential wave equation reported in [J. Wójcik, J. Acoust. Soc. Am. 104 (1998) 2654-2663; V.P. Kuznetsov, Akust. Zh. 16 (1970) 548-553]. An incremental stepping scheme allows for forward wave propagation. The operator-splitting method accounts independently for the effects of full diffraction, absorption and nonlinear interactions of harmonics. The WE method represents the propagating pulsed acoustic wave as a superposition of wavelet-like sinusoidal pulses with carrier frequencies being the harmonics of the boundary tone burst disturbance. The model is valid for lossy media, arbitrarily shaped plane and focused sources, accounts for the effects of diffraction and can be applied to continuous as well as to pulsed waves. Depending on the source geometry, level of nonlinearity and frequency bandwidth, in comparison with the conventional approach the Time-Averaged Wave Envelopes (TAWE) method shortens computational time of the full 4D nonlinear field calculation by at least an order of magnitude; thus, predictions of nonlinear beam propagation from complex sources (such as phased arrays) can be available within 30-60 min using only a standard PC. The approximate ratio between the computational time costs obtained by using the TAWE method and the conventional approach in calculations of the nonlinear interactions is proportional to 1/N2, and in memory consumption to 1/N where N is the average bandwidth of the individual wavelets. Numerical computations comparing the spatial field distributions obtained by using both the TAWE method and the conventional approach (based on a Fourier series representation of the propagating wave) are given for circular source geometry, which represents the most challenging case from the computational time point of view. For two cases, short (2 cycle) and long (8 cycle) 2 MHz bursts, the computational times were 10 min and 15 min versus 2 h and 8 h for the TAWE method versus the conventional method, respectively.     

球面电磁波斜入射建筑物内场的数值分析

为了分析窄带电磁脉冲源近、远场辐照大型建筑物内部空间电磁场分布特点,采用时域有限差分方法,对比分析平面波、球面波场源正面斜入射多层多单元建筑物的计算结果,并较为全面地分析了各房间各水平面中心电场强度幅值、各层各水平面电场强度最大值分布情况。建筑物各层相对应高度面上场强分布近似,在电磁波传播方向上窗户房间内部区域场强较强,两种波场源入射结果基本一致;在其余区域场强较弱,特别是在较大空间内的区域,球面波场源入射下场强相对更弱。其中通过与待模拟建筑物外形相紧凑性匹配的球面波场源构建,模拟仿真了窄带球面电磁波在真空空间的辐射传播,计算结果验证所构建球面波近场源准确可行。     

Excitation of a tube wave in a borehole by an external isotropic point source

The excitation of a tube wave in an infinite fluid-filled borehole by an external isotropic point source is considered. The solution to the problem is obtained in the form of a double integral with respect to the ray parameter (slowness) and frequency. The integral with respect to the slowness is transformed to a contour integral in the complex slowness plane and then reduced to the integral over the edges of the cut of the vertical slowness function and the semiresidues at the poles. An asymptotic expression for the wave field in the borehole is obtained with allowance for the radiation condition at infinity. It is shown that, when a longitudinal spherical wave is incident on the borehole, only one low-frequency Stoneley wave is excited and not two, as was assumed earlier [1].     

A theoretical study of ultrasonic wave transmission through a fluid-solid interface

This article develops a model for the study of the transient ultrasonic waves radiated by a transducer in a liquid and transmitted into a solid through a plane interface. The method is an extension to the transient case, of the angular spectrum method previously developed for the monochromatic case. It is based on the decomposition of the ultrasonic field, in impulse plane waves. The radiated waveform is calculated at any point in the field by a simple summation of these impulse plane waves, where the propagation delay and the refraction have been taken into account. These plane waves are, first of all, delayed by an amount of time corresponding to the travel time up to the considered field point. The transmission through the plane interface is taken into account by using Snell refraction laws and transmission coefficients. In the obtained results all the waves previously described by other authors are highlighted: direct wave, edge waves, head waves as well as subsurface waves with a clear resolution between compression and shear waves.     

A low frequency stable plane wave addition theorem

The multilevel fast multipole algorithm (MLFMA) is a well known and very successful method for accelerating the matrix-vector products required for the iterative solution of Helmholtz problems. The MLFMA has an important drawback, namely its inability to handle scattering problems with a lot of subwavelength detail due to the low frequency (LF) breakdown of the MLFMA. There is a need to extend the MLFMA to LF, since alternative methods are less efficient (multipole methods) or hard to implement (spectral methods). In this paper a new addition theorem will be developed that does not suffer from an LF breakdown. Instead it suffers from a high-frequency (HF) breakdown. The new method relies on a novel set of distributions, the so-called pseudospherical harmonics, closely related to the spherical harmonics. These allow the discretization points and translation operators to be calculated in closed form. Hence the method presented in this paper allows the easy implementation of a method that is stable at LF. Furthermore, a combination of the traditional MLFMA and the new method allows for the construction of a broadband MLFMA.