A review of time-domain impedance modelling and applications

Time-domain impedance boundary conditions are reviewed and similarities between the models are shown. The extended Helmholtz resonator model is used with a time-domain impedance boundary condition for impedance eduction and simulation of the noise radiation from a lined aeroengine nozzle. Model parameters obtained by impedance eduction from flow duct measurements and then used to predict the noise radiation for realistic aeroengine nozzle flow and geometry. The prediction of sound radiation from a coaxial jet nozzle at approach conditions with the educed model parameters is performed in 2D and 3D. Various features of the numerical method that enable it to be used for realistic applications are presented. This includes the suppression of flow instabilities in the linear inviscid solution and analysis of the solution to prove the global conservation of acoustic energy in each individual result. It is found that the current approach is too time consuming for optimisation, even for a 2D simulation using parallel processing. Thus, finally the potential of porting a CAA method to the graphics processing unit (GPU) is shown in one example. The GPU based computation is about 100 times faster. The results are encouraging, even though the GPU variant needs further validation.     

An investigation into noise radiation from flow control valves with particular reference to flow rate measurement

Flow control valve noise has been investigated, both analytically and experimentally, for two cases of variable orifice area and variable mass flow rate. The analysis is based on Lighthill's acoustic analogy, together with Curle's theory for noise generation from fluid-solid boundary interactions. Results indicate that monitoring sound levels of valve noise could be useful as a feasible means for flow rate measurements and on-line automatic flow control in pipes.     

Simulation of pressure pulsations in the breathing circuit of a deep-sea diving system by using the four-pole method

A method for calculating pressure pulsations in compressor and engine manifolds known as the “four-pole method” is applied to the breathing circuit of a deep-sea diving system. The response of the diver's helmet to the input of two banks of compressors is formulated. A computer is used to solve the response equation and generate the sound pressure level spectrum in the helmet. Results are shown to be in good agreement with experimental data. The computer simulation is used to examine the effects of such parameters as diver's depth, breathing gas flow rate, cylinder bank phasing and system geometry.     

A High-Efficiency Spectral Method for Two-Dimensional Ocean Acoustic Propagation Calculations

The accuracy and efficiency of sound field calculations highly concern issues of hydroacoustics. Recently, one-dimensional spectral methods have shown high-precision characteristics when solving the sound field but can solve only simplified models of underwater acoustic propagation, thus their application range is small. Therefore, it is necessary to directly calculate the two-dimensional Helmholtz equation of ocean acoustic propagation. Here, we use the Chebyshev–Galerkin and Chebyshev collocation methods to solve the two-dimensional Helmholtz model equation. Then, the Chebyshev collocation method is used to model ocean acoustic propagation because, unlike the Galerkin method, the collocation method does not need stringent boundary conditions. Compared with the mature Kraken program, the Chebyshev collocation method exhibits a higher numerical accuracy. However, the shortcoming of the collocation method is that the computational efficiency cannot satisfy the requirements of real-time applications due to the large number of calculations. Then, we implemented the parallel code of the collocation method, which could effectively improve calculation effectiveness.     

Stability and three-wave interaction of disturbances in supersonic boundary layer with mass exchange on the wall

The interaction of disturbances in a boundary layer of the compressible gas is considered in the linear and nonlinear approximation (the weakly nonlinear theory of stability) in the presence of mass exchange (gas blowing or suction) on the surface. The regimes of moderate (the Mach number M = 2) and high (M = 5.35) supersonic velocities of the flow are considered. The suction from the surface is shown to lead to a considerable variation of the linear evolution of disturbances: the vortex disturbances of the first mode and the acoustic disturbances of the second mode are stabilized, the rate of variation is determined by suction intensity. The nonlinear interactions in three-wave systems between the vortex waves in asymmetric triplets at M = 2 and between the waves of different nature (acoustic and vortex waves) ?? in the symmetric ones at M = 5.35 are considered. The planar acoustic wave is the excitation wave in the latter, which excites the three-dimensional subharmonic components of the vortex nature. It is shown that one can delay considerably the transition region with the aid of suction, thereby one can reduce the skin-friction drag. In the gas blowing regime, strong deformations of the mean fields of boundary layers occur, which lead to the destabilization of the vortex and acoustic waves in the linear region and activate the nonlinear processes in transition region. One can expect that this will lead to the acceleration of tripping in supersonic flow.     

小尺度阻尼边界封闭空间声传递函数的有限元素法计算

如何求解阻尼边界封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键技术，能处理任意形状及复杂边界条件的有限元素法可作为求解该问题的适合方法，以室内声声有源Helmholtz方程及其相应边界方程为基础，本文推导出了用于小尺度阻尼边界封闭空间声传递函数的有限元素求解方法，并编制了相应的计算机程序，在算例中，首先通过与模态叠加法计算结果进行比较，验证了该方法的正确性。最后计算了某型车体内腔中任意两点间声传递函数。     

Determination of the exact solutions to the inhomogeneous burgers equation with the use of the darboux transformation

A method of determining the exact solutions to the Burgers equation on the basis of the Darboux transformation is described. It is shown that a single application of the Darboux transformation to the homogeneous Burgers equation transforms the latter into the inhomogeneous equation describing acoustic wave propagation against transonic flow in the de Laval nozzle. In this case, the contraction ratio of the nozzle is fixed and determined by the viscosity coefficient of the medium. Based on the exact solution of the homogeneous Burgers equation, for the aforementioned problem of the flow in the nozzle, all the possible regular steady-state solutions are presented and the evolution of nonstationary solutions is investigated. The algorithm of a multiple Darboux transformation, which allows an increase in the strength of inhomogeneity, i.e., in the contraction ratio of the nozzle, is determined. This approach leads to a discrete set of possible contraction ratios at which exact solutions can be obtained. The Crum’s theorem is used to derive a formula that allows determination of the exact solutions to the inhomogeneous Burgers equation from the solutions to the homogeneous heat transfer equation. It is noted that, in fact, the proposed algorithm of the multiple Darboux transformation makes it possible to decrease the viscosity coefficient of the medium in a discrete way.     

声波方程数值反演的FR法

本文提出了声波方程反演的FR方法,FR法是一个应用范围很广的波动方程反演方法,它不受空间维数与边界条件的限制,适应于各类波动反问题。文中推导了算法的基本公式,并给出了一维与二维声波反问题的数值算例。     

Computing particle velocity and intensity in diffracting systems

In a previous paper the authors have shown how it is possible to compute the acoustic pressure due to an axi-symmetric radiator with different boundary conditions. The technique involves the conversion of a double integral to a single integral by means of a geometric transform, and may be applied with the Green's function appropriate to each boundary condition.The present work shows how the form of the single integrals makes it possible to derive and compute both the particle velocity and the wave intensity. It has become apparent that the particle velocity must be described as an elliptic vector. A new concept, the ‘diffraction front’, has been defined as a surface in which there is no net energy transfer, and this has been applied to the evaluation of the magnitude and direction of energy flow. Plots have been made of flow lines and diffraction fronts, and also three-dimensional representations of energy intensity have been made for those acoustic fields in which pressure has already been evaluated.     

FINITE ELEMENT AND BOUNDARY ELEMENT MODELLING FOR THE ACOUSTIC WAVE TRANSMISSION IN MEAN FLOW MEDIUM

The frequency characteristics of the acoustic wave transmission in a medium with mean flow are considered. One approach is to solve the Helmholtz equation with mean flow medium in original co-ordinates, which is directly discretized for the one-dimensional and the axisymmetric FEM. Another approach is to transform the equation into the standard Helmholtz equation, which is discretized for the axisymmetric FEM and the three-dimensional BEM. The numerical models are examined first for a straight circular duct. The solutions by the numerical approaches are compared with the analytical solution. The examination is then extended to the case when the mean flow is locally present in the muffler with expansion chamber. To model the spatial mean flow in the BEM model, the partitioned domain approach is also developed. No shear effect between the two regions are included.     

An extended car-following model considering the appearing probability of truck and driver's characteristics

In this paper, the appearing probability of truck is introduced and an extended car-following model is presented to analyze the traffic flow based on the consideration of driver's characteristics, under honk environment. The stability condition of this proposed model is obtained through linear stability analysis. In order to study the evolution properties of traffic wave near the critical point, the mKdV equation is derived by the reductive perturbation method. The results show that the traffic flow will become more disorder for the larger appearing probability of truck. Besides, the appearance of leading truck affects not only the stability of traffic flow, but also the effect of other aspects on traffic flow, such as: driver's reaction and honk effect. The effects of them on traffic flow are closely correlated with the appearing probability of truck. Finally, the numerical simulations under the periodic boundary condition are carried out to verify the proposed model. And they are consistent with the theoretical findings.     

Investigation on Rijke pipe＇s acoustic characteristics by numerical simulation： modeling the pulsing flow field coupled the inner of pipe with its outer space

Based on the results of fluid dynamics, heat transfer and acoustics, a Computational Fluid Dynamics （CFD） method was utilized to study the acoustic characteristics and self-excited pulsation mechanism inside a Rijke pipe. To avoid settling the irrational boundary conditions of the finite-amplitude standing wave in the Rijke thermo-acoustic system, the simulation modeling in the flow field, which coupled the inner of pipe with its outer space, was carried out to replace the traditional way in form of internal flow field numerical investigations. A hypothesis for heat source in energy equation including the relationship on unsteady heat of air around heat source, oscillation pressure and oscillation velocity was presented. To reflect the essence of Rijke pipe, simulation on self-excited oscillation was conducted by means of its own pulsation of pressure, velocity and temperature. This method can make the convergence process steady and effectively avoid divergence. The physical phenomenon of the self-excited Rijke pipe was analyzed. Moreover, the mechanisms on the Rijke pipe＇s self-excited oscillation were explained. Based on this method, comparative researches on the acoustic characteristic of the Rijke pipe with different size and different shape of nozzle were performed. The simulation results agreed with the experimental data satisfactorily. The results show that this numerical simulation can be used to study the sound pressure of nozzle for the engineering application of Rijke pipes.     

三维海洋波导中多声散射的数值解

研究三维海洋波导中多障碍物的声散射问题.基于Galerkin变分原理和DtN映射方法,发展一种自然边界元与有限元耦合的求解方法,提出有效的数值实现方法,避免常规有限元中需要采用吸收边界的缺点.对多种形状多障碍物的数值计算结果表明,所建立的理论方法是有效的.     

Vibro-acoustic behaviors of an elastically restrained double-panel structure with an acoustic cavity of arbitrary boundary impedance

An analytical study on the vibro-acoustic behaviors of a double-panel structure with an acoustic cavity is presented. Unlike the existing studies, a structural–acoustic coupling model of an elastically restrained double-panel structure with an acoustic cavity having arbitrary impedance on sidewalls around the cavity is developed in which the two dimensional (2D) and three dimensional (3D) modified Fourier series are used to represent the displacement of the panels and the sound pressure inside the cavity, respectively. The unknown expansions coefficients are treated as the generalized coordinates and the Rayleigh–Ritz method is employed to determine displacement and sound pressure solutions based on the energy expressions for the coupled structural–acoustic system. The effectiveness and accuracy of the present model is validated by numerical example and comparison with finite element method (FEM) and existing analytical method, with good agreement achieved. The influence of key parameters on the vibro-acoustic behaviors and sound transmission of the double-panel structure is investigated, including: cavity thickness, boundary conditions, sidewall impedance, and the acoustic medium in the cavity.     

A prediction of the acoustic field of a vibrating box by FEM-BEM

I.Mathematicalmodelofacousticfie1dByapplyingtheSommerfeldradiationcondihonThus,anacousticfic1disdescribcdbyHelmholtzformula.Byconsideringaficldpointinthearea,Helmholtzformulaisgivenbythefollowingthreeintegra1equations:(1)Theexternalintcgra1fOrmu1a'(2)Thcintcrnalintegralformu1af(3)Thesurfaceintegralformu1afwhcreP.isthesoundpressureonthcSsurface;rthedistancefrompointqtothee1ementsurfacedSiqtheficldpoint,tuthefrcqucncyoftheacousticradiation;pthedensityofthesurroundingmedium,vsthenormalvibrati…     

Perturbation theory of scattering from irregular bodies

A perturbation solution is derived for the following problem: A time harmonic wave of amplitude ψ, propagating in a medium with wave number k, is incident on an irregular volume V, inside of which the propagation constant k′(r) can be an arbitrary function of | r |, where r is a position vector with origin inside V. The boundary conditions are that both ψ and its normal derivative ∂ψ/∂n may be discontinuous across the surface of V. Special cases of these conditions correspond to acoustic scattering, to B-wave scattering from a dielectric cylinder, or to the classical Dirichlet (ψ = 0) or Neumann (∂ψ/∂n = 0) surface conditions. An integral equation is derived that satisfies the appropriate differential equations both outside and inside the body, and satisfies the boundary conditions as well. This equation is reduced to a set of linear algebraic equations by expansion in a certain basis set and these linear equations are then solved in a perturbation approximation for the case that the surface of the body differs from a sphere or cylinder by a small parameter λ. Comparison is made with formulae in the literature, and except for some minor discrepancies, which are here corrected, there is general agreement.     

Effects of a multi-layered poro-elastic ground on attenuation of acoustic waves and ground vibration

Ground conditions affect the propagation of outdoor sound and vibration. This paper focuses on the interaction between air pressure and porous ground at low frequencies- where mechanisms other than the rigid porous effects used in locally reacting models may be important. A 2-D analytical model has been developed in this study for the calculation of acoustic and acousto-seismic admittances in a multi-layered poro-elastic ground. The model can be used as a prediction tool both for the ground effect on the sound and the generation of ground vibration by the sound. The modelled acoustic admittance is validated successfully against established rigid frame admittance models over a frequency range of 1 Hz-3 kHz. Moreover, the acousto-seismic impedance is verified against full scale airblast field test data measured during the Norwegian Trials full scale field test program. For certain ground types, the predicted acoustic admittance illustrates a different behaviour compared with predictions from the traditional rigid frame acoustic impedance models. This emphasises the importance of including a deformable frame in the model to obtain realistic results for these conditions.     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.