有限长号筒中非线性声传播的数值模拟研究

研究了大振幅活塞声源经有限长号筒向外辐射声波的非线性声学问题。采用具有频散保持特性的高精度计算格式求解了适用于变截面管道的一维非线性声场模型,并考虑到声波的非线性畸变及管口处的声反射,加入了宽频时域声阻抗边界条件。宽频阻抗模型的共轭复数系数采用优化拟合方法近似求解,并采用递推卷积算法快速求解出时域声阻抗。在弱非线性条件下模拟指数形号筒中的声传播取得了与已有实验相一致的结果,表明模型能够描述声波非线性畸变带来的宽频特性。进而采用本模型数值模拟了大振幅活塞声源在双曲形、锥形、指数形和正弦形号筒中的非线性声传播问题,结果表明号筒出口声压级受活塞振动速度、频率以及号筒形状的影响,并分析讨论了波形畸变与号筒几何形状之间的关系。      

Model and estimation method for predicting the sound radiated by a horn loudspeaker - With application to a car horn

This paper deals with a new car horn device made of a sound synthesizer and an electrodynamic horn loudspeaker. It presents an one-dimensional model allowing to predict the loudspeaker efficiency and a specific method to estimate experimentally the model parameters. First, this model aims at reducing the time spent in the design process. Second it aims at correcting the sound emitted by the sound synthesizer in order that the listener hears the sound designed for creating the warning message. The study gives a survey of the vast loudspeaker literature. It is based on the conventional electroacoustic approach used for electrodynamic loudspeakers and on wave propagation models used for characterizing acoustic horns. The estimation of the model parameter values is performed using measurements of the electrical impedance of the loudspeaker and of the acoustic impedance of the horn. The model is assessed by comparing the calculated and measured electrical impedances and horn efficiencies. Results show that the model predicts well the horn efficiency up to 2500 Hz, the limitation being due to the horn radiation impedance modelization.     

Weak nonlinear propagation of sound in a finite exponential horn.

This article presents an approximate solution for weak nonlinear standing waves in the interior of an exponential acoustic horn. An analytical approach is chosen assuming one-dimensional plane-wave propagation in a lossless fluid within an exponential horn. The model developed for the propagation of finite-amplitude waves includes linear reflections at the throat and at the mouth of the horn, and neglects boundary layer effects. Starting from the one-dimensional continuity and momentum equations and an isentropic pressure-density relation in Eulerian coordinates, a perturbation analysis is used to obtain a hierarchy of wave equations with nonlinear source terms. Green's theorem is used to obtain a formal solution of the inhomogeneous equation which takes into account linear reflections at the ends of the horn, and the solution is applied to the nonlinear horn problem to yield the acoustic pressure for each order, first in the frequency and then in the time domain. In order to validate the model, an experimental setup for measuring fundamental and second harmonic pressures inside the horn has been developed. For an imposed throat fundamental level, good agreement is obtained between predicted and measured levels (fundamental and second harmonic) at the mouth of the horn.     

AN ACOUSTIC TRANSFORMER COMPOSED OF HORNS AND ITS APPLICATION

A new acoustic transformer was developed by connecting three horns to improve radiation performance in the frequency region below 500 Hz. The proposed acoustic transformer was evaluated by numerical analysis using the commercial computer program SYSNOISE and by experiment. The acoustic transformer is composed of three horns, one of which was used in an inverted form. A design model was developed by use of Webster's horn equation and showed that the transformer can improve radiation efficiency. This was confirmed by numerical calculation using SYSNOISE. An acoustic projector was designed by use of the developed transformer and a piezoelectric unimorph-type actuator. The sound pressure measured at the mouth of the constructed acoustic projector was compared with the sound pressures evaluated at the same location by numerical calculation to investigate the differences between the numerical simulation model and the actual acoustic projector. Sound pressures generated by several acoustic radiators were calculated numerically and compared with the measured and calculated sound pressures of the developed acoustic projector to evaluate the effects of use of the proposed acoustic transformer. The comparative evaluation shows that the proposed acoustic transformer can provide up to a 10 dB gain over use of a horn in the narrow band low-frequency region from 100 to 200 Hz.     

A study on the multi-freedom broadband impedance model for time-domain simulations

The purpose of this paper is to construct a general broadband impedance model, which is suited for predicting acoustic propagation problems in time domain.A multi-freedom broadband impedance model for sound propagation over impedance surfaces is proposed and the corresponding time domain impedance boundary condition is presented.Basing on the extended Helmholtz resonator,the multi-freedom impedance model is constructed through combing with a sum of rational functions in the form of general complex-conjugate pole-residue pairs and it is proved that the impedance model is well posed.The impedance boundary condition can be implemented into a computational aeroacoustics solver by a recursive convolution technique, which results in a fast and computationally efficient algorithm.The two dimensional and three dimensional benchmark problems are selected to validate the accuracy of the proposed impedance model and time domain simulations.The numerical results are in good agreement with the reference solutions.It is demonstrated that the proposed impedance model can be used to describe the broadband characteristics of acoustic liners,and the corresponding time domain impedance boundary condition is viable and accurate for the prediction of sound propagation over broadband impedance surfaces.     

含气泡液体中的非线性声传播

考虑到分布在液体中的气泡是声波在含气泡液体中传播时引起非线性的一个很重要的因素,本文研究了声波在含气泡液体中的非线性传播.将气体含量的影响引入到声波在液体中传播的方程中,从而得到声波在气液混合物中传播的数学模型.通过对该模型进行数值模拟发现,气体含量、驱动声场声压幅值及驱动声场作用时间均会影响到气液混合物中的声场分布及声压幅值大小.液体中的气泡会"阻滞"液体中声场的传播并将能量"聚集"在声源附近.对于连续大功率的驱动声场来说,液体中的气泡会"阻滞"气液混合物中声场及其能量的传播.     

Finite element simulation of non-linear acoustic generation in a horn loudspeaker

The loudspeaker is an electro-acoustic device for sound reproduction which requires the distortion as small as possible. The distortion may arise from the magnetic non-linearity of the york, the uneven magnetic field distribution, the mechanical non-linearity at the diaphragm suspension and the acoustic non-linearity due to the high sound pressure and velocity in the duct-radiation system. A horn is sometimes provided in front of the vibrating diaphragm radiator, which plays an important role to increase the efficiency by matching the acoustic impedance between the radiator and the ambient medium. The horn is in many cases folded twice or three times to shorten the length, which further degrades the reproduction quality. The sound intensity and velocity are apt to attain very high in the small cross-sectional area in the throat and in the folded regions, which may cause the distortion due to the non-linear effect of the medium. The present paper is to investigate the frequency characteristics of the loudspeaker numerically evaluating the generation of the harmonics and sub-harmonics. An axisymmetric folded horn is considered for which the wave equation with the non-linear term retained is solved by the finite element method. The solution is made in time domain in which the sound pressure calculated at the opening end of the horn is Fourier-transformed to the frequency domain to evaluate the distortion, while the wave marching in the horn is visualized.     

0-3型压电复合材料覆盖层水下 吸声性能的理论研究

目前压电分流阻尼技术在振动和噪声领域的应用得到了广泛的关注. 本文尝试将压电分流阻尼技术应用于水下吸声领域, 以提高覆盖层的吸声性能. 将压电覆盖层厚度模态的机电方程和声波传播的传递矩阵相结合, 建立一维电声模型. 该模型可以用于分析多层压电和非压电水下吸声覆盖层的吸声性能. 采用该模型分析了0-3型压电复合材料覆盖层的水下吸声性能. 压电复合材料的参数是采用Furukawa的模型计算的. 研究结果表明, 采用合适的分流电阻, 负电容分流电路可以在较宽的频率范围显著提高覆盖层的吸声性能. 其原理可以从阻抗匹配的角度解释, 负电容分流电路可以调整压电覆盖层的表面声阻抗, 使之与水的特性声阻抗相匹配.     

Reflectance of acoustic horns and solution of the inverse problem

A method is described for solving the inverse problem of determining the profile of an acoustic horn when time-domain reflectance (TDR) is known only at the entrance. The method involves recasting Webster's horn equation in terms of forward and backward propagating wave variables. An essential feature of this method is a requirement that the backward propagating wave be continuous at the wave-front at all locations beyond the entrance. Derivation of the inverse solution raises questions about the meaning of causality in the context of wave propagation in non-uniform tubes. Exact reflectance expressions are presented for infinite exponential, conical and parabolic horns based on exact solutions of the horn equation. Diameter functions obtained with the inverse solution are a good match to all three horn profiles.     

The use of a hybrid model to compute the nonlinear acoustic performance of silencers for the finite amplitude acoustic wave

In the present study, a hybrid method is proposed for predicting the acoustic performance of a silencer for a nonlinear wave. This method is developed by combining two models: (i) a frequency-domain model for the computation of sound attenuation due to a silencer in a linear regime and (ii) a wavenumber space model for the prediction of the nonlinear time-evolution of finite amplitudes of the acoustic wave in a uniform duct of the same length as the silencer. The present method is proposed under the observation that the physical process of the nonlinear sound attenuation phenomenon of a silencer may be decoupled into two distinct mechanisms: (a) a linear acoustic energy loss that owes to the mismatch in the acoustic impedance between reactive elements and/or the sound absorption of acoustic liners in a silencer; (b) a nonlinear acoustic energy loss that is due to the energy-cascade phenomenon that arises from the nonlinear interaction between components of different frequencies. To establish the validity of the present model for predicting the acoustic performance of silencers, two model problems are considered. First, the performance of simple expansion mufflers with nonlinear incident waves has been predicted. Second, proposed method is applied for computing nonlinear acoustic wave propagation in the NASA Langley impedance duct configuration with ceramic tubular liner (CT57). Both results obtained from the hybrid models are compared with those from computational aero-acoustic techniques in a time-space domain that utilize a high-order finite-difference method. Through these comparisons, it is shown that there are good agreements between the two predictions. The main advantage of the present method is that it can effectively compute the nonlinear acoustic performance of silencers in nonlinear regimes without time-space domain calculations that generally entail a greater computational burden.     

圆孔非线性声阻抗三维时域仿真与特性分析

为研究有限振幅声波作用下圆孔的非线性声学特性,提出了基于三维时域计算流体动力学(CFD)仿真的圆孔非线性声阻抗提取方法,通过求解层流方程来模拟声信号在圆孔及上下游的传播,以及采用横向周期性边界条件来考虑高穿孔率时圆孔之间相互作用的影响。研究了不同幅值声波作用下孔径、厚度和穿孔率对声阻抗的影响规律,通过对质点振速幅值、频率和板厚等组成的无量纲参量进行非线性回归分析,得到了圆孔非线性声阻抗的拟合公式,并将其转换为可考虑多频声波影响的时域模型。最后结合声阻抗时域模型和有限差分方法计算了直通穿孔管消声器在小振幅和有限振幅声波作用下的传递损失,通过与实验测量结果的比较,验证了拟合公式的准确性和实用性。     

Acoustic tomography of the nonlinear parameter by a small number of transducers

A wave approach to solving the problem of reconstructing the acoustic nonlinear parameter distribution by a simple transmitting-receiving system with broadband modulation of primary waves is proposed. The approach uses the effect of sound scattering by sound, which makes it possible to employ a complicated time processing of signals instead of spatial processing and, consequently, to reduce the number of transmitting and receiving transducers. The possibility of practical implementation of a two-dimensional tomographic scheme based on coherent correlation processing of broadband quasi-random signals at combination frequencies is analyzed. Results of model numerical experiments on reconstructing complicated spatial distributions of the nonlinear parameter are presented.     

Counterpropagation of waves with shock fronts in a nonlinear tissue-like medium

A numerical model for describing the counterpropagation of one-dimensional waves in a nonlinear medium with an arbitrary power law absorption and corresponding dispersion is developed. The model is based on general one-dimensional Navier-Stokes equations with absorption in the form of a time-domain convolution operator in the equation of state. The developed algorithm makes it possible to describe wave interactions in the presence of shock fronts in media like biological tissue. Numerical modeling is conducted by the finite difference method on a staggered grid; absorption and sound speed dispersion are taken into account using the causal memory function. The developed model is used for numerical calculations, which demonstrate the absorption and dispersion effects on nonlinear propagation of differently shaped pulses, as well as their reflection from impedance acoustic boundaries.     

Finite element simulation of an axisymmetric acoustic transmission system with a sound absorbing wall

A finite element formulation for the axisymmetric acoustic field problem with an arbitrary sound absorbing wall is developed on the true adjoint system approach. A triangular ring element with a second order polynomial trial function is used for discretization. As a first numerical example, the transmission characteristic of acoustic filters of expansion type is considered, in which the wall is partly treated with sound absorbing material. As a second example the effect of a temperature gradient in the chamber on the transmission characteristic is calculated. The third example is the calculation of the throat impedance of conical and exponential acoustic horns, for which the half sphere into which the radiation travels is assumed to be terminated at an appropriate radius by a hypothetical wall having the characteristic impedance of the acoustic medium. The calculated results are verified by comparison with measured results.     

Numerical simulation of the nonlinear ultrasonic pressure wave propagation in a cavitating bubbly liquid inside a sonochemical reactor

We investigate the acoustic wave propagation in bubbly liquid inside a pilot sonochemical reactor which aims to produce antibacterial medical textile fabrics by coating the textile with ZnO or CuO nanoparticles. Computational models on acoustic propagation are developed in order to aid the design procedures. The acoustic pressure wave propagation in the sonoreactor is simulated by solving the Helmholtz equation using a meshless numerical method. The paper implements both the state-of-the-art linear model and a nonlinear wave propagation model recently introduced by Louisnard (2012), and presents a novel iterative solution procedure for the nonlinear propagation model which can be implemented using any numerical method and/or programming tool. Comparative results regarding both the linear and the nonlinear wave propagation are shown. Effects of bubble size distribution and bubble volume fraction on the acoustic wave propagation are discussed in detail. The simulations demonstrate that the nonlinear model successfully captures the realistic spatial distribution of the cavitation zones and the associated acoustic pressure amplitudes.     

Surface impedance design with ground corrugation for mitigation of large-calibre gun blast noise

The surface impedance design approach is proposed for mitigating large-calibre gun blast noise. Surrounding the blast noise, we employ a group of concentric trenches with critical depths to dampen the propagation of the acoustic wave. These trenches behave like quarter-wavelength resonators and produce acoustic soft surfaces at their openings. The sound pressure is then mitigated over these soft surfaces by destructive interference and the wave attenuates rapidly along the ground surface. To evaluate the overall acoustic performance of such a design, we develop an efficient numerical solver by treating the geometry as a body of revolution (BOR). The symmetry of the structure in the revolution direction allows the 3D boundary integral equation (BIE) for acoustic wave scattering to be reduced to a 2D integral equation by the use of Fourier series expansions. Numerical experiments show that this model can effectively suppress the acoustic wave propagation horizontally and the reduction can reach about 15 dB for large-calibre gun noise with very low-frequency components.     

液体材料超声处理过程中声场和流场的分布规律研究

针对合金熔体等液体材料的超声处理过程,选取水作为透明模型材料,采用数值模拟计算和示踪粒子实验方法,研究了20和490 kHz两种频率超声作用下水中的声场和流场分布.结果表明,增大变幅杆半径能够提高水中声压水平,扩大空化效应的发生区域.当超声频率为20 kHz时,水中声压最大值出现在超声变幅杆下端面处,且声压沿传播距离的增大而显著减小.如果超声频率增加至490 kHz,水中的声压级相比于20 kHz时明显提高,且声压沿着超声传播方向呈现出周期性振荡特征.两种频率超声作用下水中的流场呈现相似的分布特征,且平均流速均随着变幅杆半径增大表现出先升高后降低的趋势.变幅杆半径相同时,20 kHz频率超声作用下水中的平均流速高于490 kHz频率超声.采用示踪粒子图像测速技术实时观察和测定了水中的流速分布,发现其与计算结果基本一致.     

Acoustic particle velocity horns

The paper considers receiving acoustic horns designed for particle velocity amplification and suitable for use in vector sensing applications. Unlike conventional horns, designed for acoustic pressure amplification, acoustic velocity horns (AVHs) deliver significant velocity amplification even when the overall size of the horn is much less than an acoustic wavelength. An AVH requires an open-ended configuration, as compared to pressure horns which are terminated at the throat. The appropriate formulation, based on Webster's one-dimensional horn equation, is derived and analyzed for single conical and exponential horns as well as for double-horn configurations. Predicted horn amplification factors (ratio of mouth-to-throat radii) were verified using numerical modeling. It is shown that three independent geometrical parameters principally control a horn's performance: length l, throat radius R(1), and flare rate. Below a predicted resonance region, velocity amplification is practically independent of frequency. Acoustic velocity horns are naturally directional, providing maximum velocity amplification along the boresight.     

谐振管形状对热声发动机内非线性声场的影响

研究了谐振管一端受活塞声源激励,另一端刚性封闭条件下,管道形状对热声发动机谐振管内部非线性声场的影响。基于流体力学基本方程建立了渐变截面谐振管内一维非线性声场的模型,考虑了黏性耗散及非线性效应的影响。利用伽辽金法数值求解了该模型的速度势方程,分析了谐振管形状、活塞振动速度及激励频率对管内声场的影响。将双曲形、指数形、锥形、正弦形等四种变截面谐振管内的非线性声场与圆柱形直管的情况进行了比较。结果反映了谐振管内声场的压力波动受活塞振动速度及谐振管形状的影响;显示了当活塞振动幅度较大时,谐振管内出现的波形畸变、频率曲线偏移、共振频率滞后等非线性现象;揭示了变截面谐振管在抑制管内的高阶谐波及提高压比等方面的优越性。      

Influence of resonator shape on nonlinear acoustic field in a thermoacoustic engine

The influence of the resonator shape on nonlinear acoustic field in a thermoacoustic engine is studied.The resonator of themoacoustic engine is boundary driving by a piston at one end,and the other end of it is rigid closed.A one-dimensional wave equation that accounts for gas dynamic nonlinearities and viscous dissipation in the resonator is established based on the governing equations of viscous hydromechanics.The nonlinear wave equation is solved using approximate Galerkin method.The nonlinear acoustic field in four different types of shaped resonators including hyperbolical,exponential,conical and sinusoidal are obtained and compared with that of a cylindrical resonator.It is found that the amplitude and waveform of the pressure are strongly affected by the resonator shape,the driving amplitude and the oscillation frequency of the piston.Waveform distortion,resonance frequency shift and hysteresis are observed,when the piston oscillation amplitude is large enough.The advantages of shaped resonator for thermoacoustic engine lie in inhibition of higher order harmonics and improvement of pressure ratio,etc.