非均匀流体介质内部散射声场重建的逐层计算方法*

入射声波激励下非均匀流体介质内部散射声场的重建方法对超声层析成像具有重要意义。以往采用矩量法求解,但该方法全域离散形成的复数满秩矩阵规模随着分辨率与计算精度的提高而急剧增大,对算力具有很高的要求,一定程度上限制了其在实际中的应用。为克服上述缺陷,本文以逐层离散、逐层计算为核心思想,以声散射基本公式与近场声全息理论为基础,推导出逐层计算非均匀流体介质内部散射声场的理论公式并给出对应的几何离散模型。为验证该方法的可行性,以矩量法为参照,对同样的介质模型进行介质内部声场重构仿真。结果表明,逐层算法不仅可以有效地重建非均匀流体介质内部散射声场,且大幅度减小了求解规模。     

An effective quiescent medium for sound propagating through an inhomogeneous,moving fluid

The idea of similarity between acoustic fields in a moving fluid and in a certain "effective" quiescent medium, first put forward by Lord Rayleigh, proved very helpful in understanding and modeling sound propagation in an atmosphere with winds and in an ocean with currents, as well as in other applications involving flows with small velocity compared to sound speed. Known as effective sound speed approximation, the idea is routinely utilized in the contexts of the ray theory, normal mode representation of the sound field, and the parabolic approximation. Despite the wide use of the concept of effective sound speed in acoustics of moving media, no theoretical justification of Rayleigh's idea was published that would be independent of the chosen representation of the sound field and uniformly apply to distinct propagation regimes. In this paper, we present such a justification by reducing boundary conditions and a wave equation governing sound fields in the inhomogeneous moving fluid with a slow flow to boundary conditions and a wave equation in a quiescent fluid with effective sound speed and density. The derivation provides insight into validity conditions of the concept of effective quiescent fluid. Introduction of effective density in conjunction with effective sound speed is essential to ensure accurate reproduction of acoustic pressure amplitude in the effective medium. Effective parameters depend on sound speed, flow velocity, and density of the moving fluid as well as on sound propagation direction. Conditions are discussed under which the dependence on the propagation direction can be avoided or relaxed.     

Computer models in room acoustics: The ray tracing method and the auralization algorithms

Computer algorithms are described for constructing virtual acoustic models of various rooms that should satisfy some specific sound quality criteria. The algorithms are based on the ray tracing method, which, in the general case, allows calculation of the amplitude of an acoustic ray that survived multiple reflections from arbitrary curved surfaces. As a result, calculations of room acoustics are reduced to tracing the trajectories of all the acoustic rays in the course of their propagation with multiple reflections from reflecting surfaces to the point of their complete decay. For this approach to be used, the following physical properties of a room should be known: the geometry of the reflecting surfaces, the absorption and diffusion coefficients on each of these surfaces, and the decay law for rays propagating in air. The proposed models allow for the solution of the important problem of architectural acoustics called the auralization problem, i.e., to predict how any given audio segment will sound in any given hall on the basis of computer simulation alone, without any full-scale testing in specific halls.     

海洋声学环境下水中声源的时延估计法定位精度分析

水中声源的定位精度受到海洋声学环境的重要影响。结合海上试验的实际应用,分析了水下观测平台采用时延估计法对声源的定位精度问题。根据理论分析,计算了时延估计误差、海洋中声速不均匀、平台非稳性、及声传播起伏等因素引起的俯仰角和方位角误差。利用误差传递公式,获得了上述因素引起的不同平台深度下,不同距离声源的定位误差。比较了采用平面阵与立体阵、是否补偿声线弯曲效应等条件下定位误差的变化,并通过海上试验结果进行了部分验证。研究结果表明,海洋声速不均匀对定位误差的贡献最大。采用立体阵代替平面阵、测量海洋声速剖面并补偿声线弯曲引起的定位误差,在1000m距离上可使定位相对误差从最大30%降低到约10%,有效提高了较远距离上的定位精度。研究结果对于采取措施提高水中声源的定位精度有指导意义。      

Efficiency of acoustic heating produced in the thermoviscous flow of a fluid with relaxation

Instantaneous acoustic heating of a fluid with thermodynamic relaxation is the subject of investigation. Among others, viscoelastic biological media described by the Maxwell model of the viscous stress tensor, belong to this type of fluid. The governing equation of acoustic heating is derived by means of the special linear combination of conservation equations in differential form, allowing the reduction of all acoustic terms in the linear part of the final equation, but preserving terms belonging to the thermal mode responsible for heating. The procedure of decomposition is valid for weakly nonlinear flows, resulting in the nonlinear terms responsible for the modes interaction. Nonlinear acoustic terms form a source of acoustic heating in the case of dominative sound, which reflects the thermoviscous and dispersive properties of a fluid. The method of deriving the governing equations does not need averaging over the sound period, and the final governing dynamic equation of the thermal mode is instantaneous. Some examples of acoustic heating are illustrated and discussed, conclusions about efficiency of heating caused by different sound impulses are made.     

大尺寸弹性材料中声速的实验室测量方法研究

针对大尺寸弹性材料中声速的测量问题,将传统穿透测量法与分层介质中的射线声学理论相结合,设计了一种在水池中反演测量材料中纵、横波声速的方法.该方法使用水听器接收透过待测材料后到达的高频水声信号,依据实测信号传播时间与理论值比较代价函数,最终利用寻优算法计算材料中的两类声速。本方法既克服了传统方法不能同时获得材料内纵、横波声速的缺点,又可保证待测材料的完整性。利用开发的测量系统,水池实际测量结果的不确定度小于3.5%      

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.     

Research on sound velocity measurement in a laboratory for large elastic materials

A measurement scheme carried out in a tank is designed to obtain the compressionaland shear-wave velocities of a large elastic material.A hydrophone is used to receive the high frequency acoustic signals which penetrate the tested material,in order to determine the transmission time from the source to the hydrophone,the transmission time is also calculated according to the ray acoustic theory in layered media.A cost function is built based on the measured and the calculated transmission time,then the compressional- and shear-wave velocities can be obtained using the optimization algorithm.Compared with the traditional measurement scheme,this approach can not only get the 2 kinds of sound velocities in the tested material at the same time,but also keep the integrality of the tested material.With the proposed measurement system,the uncertainty of measurement results is less than 3.5%.     

Equations for finite-difference, time-domain simulation of sound propagation in moving inhomogeneous media and numerical implementation

Finite-difference, time-domain (FDTD) calculations are typically performed with partial differential equations that are first order in time. Equation sets appropriate for FDTD calculations in a moving inhomogeneous medium (with an emphasis on the atmosphere) are derived and discussed in this paper. Two candidate equation sets, both derived from linearized equations of fluid dynamics, are proposed. The first, which contains three coupled equations for the sound pressure, vector acoustic velocity, and acoustic density, is obtained without any approximations. The second, which contains two coupled equations for the sound pressure and vector acoustic velocity, is derived by ignoring terms proportional to the divergence of the medium velocity and the gradient of the ambient pressure. It is shown that the second set has the same or a wider range of applicability than equations for the sound pressure that have been previously used for analytical and numerical studies of sound propagation in a moving atmosphere. Practical FDTD implementation of the second set of equations is discussed. Results show good agreement with theoretical predictions of the sound pressure due to a point monochromatic source in a uniform, high Mach number flow and with Fast Field Program calculations of sound propagation in a stratified moving atmosphere.     

Ultrasonic flow measurement and wall acoustic impedance effects

An examination of the influence of wall acoustic impedance effects on sound propagation in flowing liquids confined by cylindrical walls is presented. Special focus is given to the importance of the wall acoustic impedance value for ultrasonic flow meter performance. The mathematical model presented allows any radially-dependent axial flow profile to be examined in the linear flow acoustics regime where fluid flow speed is much smaller than the fluid sound speed everywhere in the fluid medium.     

Generation of the vorticity mode by sound in a Bingham plastic

This study investigates interaction between acoustic and non-acoustic modes, such as vorticity mode, in some class of a non-newtonian fluid called Bingham plastic. The instantaneous equations describing interaction between different modes are derived. The attention is paid to the nonlinear effects in the field of intense sound. The resulting equations which describe dynamics of both sound and the vorticity mode apply to both periodic and aperiodic sound of any waveform. They use only instantaneous quantities and do not imply averaging over the sound period. The theory is illustrated by an example of acoustic force of vorticity induced in the field of a Gaussian sound beam. Some unusual peculiarities in both sound and the vorticity induced in its field as compared to a newtonian fluid, are discovered.     

浅海周期起伏海底环境下的声传播

海底粗糙对水下声传播及水声探测等应用具有重要影响.利用黄海夏季典型海洋环境,分析了同时存在海底周期起伏和强温跃层条件下的声传播特性,结果表明:由于海底周期起伏的存在,对于低频(<1 kHz)、近程(10 km)的声信号,传播损失可增大5—30 dB.总结了声传播损失及脉冲到达结构随声源深度、海底起伏周期及起伏高度等因素变化的规律.当海底起伏周期不变时,起伏高度越大引起的异常声传播的影响随之变大;当起伏高度不变时,随着起伏周期变大,其对声传播的影响逐渐变小.用射线理论分析了其影响机理,由于海底周期起伏改变了声波与海底的入射和反射角度,使得原本小掠射角入射到海底的声线变为大掠射角,导致海底的反射损失增大;另一方面,声线反射角度的改变会使得原本可以到达接收点的声能量,由于与海底作用次数增加或变为反向传播而大幅度衰减.在浅海负跃层环境下,声源位于跃层上比位于跃层下对声传播影响更大.周期起伏海底对脉冲声传播的影响表现在引起不同角度的声线(或简正波号数)之间的能量发生转化,一些大角度声线能量衰减加大,多途结构变少.多途结构到达时间及相对幅度的变化进而影响声场的频谱,会使得基于匹配场定位的方法性能受到影响.所以,声呐在实际浅海环境中应用时,应对起伏海底的影响予以重视.此外,研究结果对海底地形测绘空间精度的提高也具有重要参考意义.     

Passive acoustic tomography of the ocean using arrays of unknown shape

It was recently established that ocean acoustic tomography based on an inversion of ray travel times can be implemented without use of any dedicated sound sources by cross-correlating the ambient noise recorded on two line arrays, the shapes of which are known. In contrast to active tomography, the amount of useful information from noise interferometry is proportional to the product of the numbers of receivers in the two arrays. In our study based on the 2D and 3D numerical experiments, we examine a hypothesis concerning the feasibility of simultaneous performance of a passive ray tomography and passive positioning of arrays through cross-correlation of ambient or shipping noise. The numerical experiments are conducted under conditions close to those of a field experiment on passive ocean tomography. It is demonstrated that, when using arrays of 20–40 hydrophones, the sound velocity profile and the array shape can be found from noise correlation to an accuracy adequate for oceanological and acoustic applications.     

Ray dynamics of the propagation of sound in nonuniform moving media

This paper presents a new ray theory for the propagation of sound waves in nonuniformly moving media. It is found that the ray equations in weakly inhomogeneous and slowly moving media are analogous to the equations of motion of charged particles in nonuniform electric and magnetic fields. The adiabatic approximation is used to study the problem of the propagation of sound rays in a model of near-ocean-bottom waveguide with horizontal flow and slowly varying parameters along the direction of propagation of the wave. A general formula is derived that describes the transverse displacement of the trajectory of the ray relative to the direction of propagation of the wave.     

声波作用下球形颗粒外声流分布的数值模拟

综合考虑声学边界层内的热损失和黏性损失,建立处于平面驻波声压波节位置二维球形颗粒外声流计算模型,利用分离时间尺度的数值方法对颗粒外声流流场特征进行模拟.将模拟结果与相应的解析解和实验结果对比,验证了数值模拟的可靠性.在此基础上,研究了雷诺数Re和斯特劳哈尔数Sr对球形颗粒声学边界层内二阶声流流场结构、涡流强度及范围的影...     

双层界面下非近轴近似多元高斯模型的相控阵声场模拟

非近轴近似多元高斯模型克服了近轴近似条件的限制,能保证较大偏转角度下模拟声场的精确度和效率。根据双层介质的瑞利积分模型和单层介质的非近轴近似多元高斯模型,推导出双层介质中基于非近轴近似多元高斯模型的单阵元辐射声场计算模型。并在得到各阵元声束偏转聚焦的延迟时间基础上,累加得到基于非近轴近似高斯模型的超声相控阵横波检测辐射声场计算模型。模拟计算带丙烯酸树脂楔块的超声相控阵探头在钢中横波检测的声束偏转与偏转聚焦辐射声场,并与近轴近似多元高斯模型的横波检测声束的偏转聚焦辐射声场进行比较分析。对比分析结果表明非近轴近似多元高斯模型计算速度更快,用时约为近轴近似多元高斯模型的1/13。此外,非近轴近似高斯模型方法横波检测的聚焦声束覆盖区域更大,能量更集中,更适于远场区域检测。      

A theoretical analysis and experimental study of the behaviour of sound in corridors

The acoustic ray image theory is developed to predict relative sound levels along a corridor or in adjacent rooms caused by a sound source in the corridor. From experimental results in full-size corridors and in an acoustic scale model, general comments about the behaviour of sound in corridors—and the effect of small intruding walls and the position of absorbent material on sound levels in the corridor—are made.     

Variations of mode amplitudes in a range-dependent waveguide

The ray method of calculating the mode amplitudes is used to analyze the sound fields in deep-water acoustic waveguides with two types of inhomogeneities of the refractive index: (i) weak inhomogeneities that cause small perturbations of ray trajectories and (ii) strong inhomogeneities with large spatial scales. Simple analytical relations are derived for describing the variations of the mode structure of the sound field in the presence of the aforementioned inhomogeneities. A new criterion defining the validity of the adiabatic approximation is formulated. To illustrate and test the results obtained, a numerical simulation of the sound fields is performed on the basis of the parabolic equation method.

     

采用小生境遗传算法反演浅海声速剖面研究

针对浅海声速剖面反演问题,采用小生境遗传算法,结合声线搜索的最快本征声线匹配反演,实现浅海负跃层条件下的声速剖面估计.利用经验正交函数对声速剖面的多参数不确定性降维,依据声场计算模型获取的最快特征声线传播时延与观测声传播时长进行匹配,采用小生境遗传优化算法,获取最优经验正交函数估计,实现声速剖面反演.按上述方法反演处理...     

Wave analysis of ray chaos in underwater acoustics

The dispersion of a wave packet in an acoustic medium is considered in the paraxial wave approximation, where the effective potential, due to variation of the speed of propagation, varies both with depth and propagation distance. The analysis of the resulting parabolic equation, similar to the Schrodinger equation, clearly demonstrates the role of ray chaos in enhancing the dispersion of the initial packet. However, wave coherence effects are also seen that suppress the effects of the ray chaos in a manner analogous to the effects of quantum chaos. (c) 1999 American Institute of Physics.