中频声振耦合的波函数统计能量法

为了提高中频声振耦合的计算效率,提出了波函数-统计能量法的结构-声学耦合方法,该方法从波动理论的角度出发,将波函数法(WBM)和统计能量法(SEA)结合,通过在耦合面分别施加声压激励和速度边界条件,推导了耦合面参数理论计算公式。将该方法用到长方体声腔和钢板耦合的模型中,并对100~1000 Hz的计算结果进行了实验验证。WBM-SEA模型与参考FEM-SEA模型以及实验模型的频响曲线对比结果表明,WBM-SEA与FEM-SEA以及实验结果吻合很好,验证了混合WBM-SEA的有效性。通过收敛性分析发现混合WBM-SEA方法计算时间比混合FEM-SEA方法更少。从而可以得出结论:混合波函数-统计能量法方法对中频声振耦合预测是有效的,且比FEM-SEA更加高效。      

A novel hybrid ES-FE-SEA for mid-frequency prediction of Transmission losses in complex acoustic systems

The widely-used numerical modeling approaches such as the finite element method (FEM) and statistical energy analysis (SEA) often have limited applicability to the transmission loss prediction in mid-frequency range. In this paper, a novel hybrid edge-based smoothed FEM coupled with statistical energy analysis (ES-FE-SEA) method is proposed to further improve the accuracy of “mid-frequency” transmission loss predictions. The application of ES-FEM will “soften” the well-known ‘‘overly-stiff’’ behavior in the standard FEM solution and reduce the inherent numerical dispersion error. While the SEA approach deals with the physical uncertainty in the relatively higher frequency range. The plate of interest is appropriately described by an ES-FEM model, due to its relative robustness to perturbations. Its adjacent reverberation cavities are modeled by employing the SEA approach, because of their high model density. The coupling and interaction between SEA subsystems and the FE subsystem is governed by the “reciprocity relationship” theorem. A standard numerical example for benchmarking is examined and excellent agreement was achieved between the prediction and reference results. The proposed ES-FE-SEA is also verified by various numerical examples. The method is finally applied to the modeling a complicated engineering problem–acoustic fields on both sides of the front windshield in a passenger car.     

A fast direct solver for the integral equations of scattering theory on planar curves with corners

We describe an approach to the numerical solution of the integral equations of scattering theory on planar curves with corners. It is rather comprehensive in that it applies to a wide variety of boundary value problems; here, we treat the Neumann and Dirichlet problems as well as the boundary value problem arising from acoustic scattering at the interface of two fluids. It achieves high accuracy, is applicable to large-scale problems and, perhaps most importantly, does not require asymptotic estimates for solutions. Instead, the singularities of solutions are resolved numerically. The approach is efficient, however, only in the low- and mid-frequency regimes. Once the scatterer becomes more than several hundred wavelengths in size, the performance of the algorithm of this paper deteriorates significantly. We illustrate our method with several numerical experiments, including the solution of a Neumann problem for the Helmholtz equation given on a domain with nearly 10000 corner points.     

The DORT solution in acoustic inverse scattering problem of a small elastic scatterer

The DORT (French acronym for Décomposition de l’Opérateur de Retournement Temporel) method is a novel approach for active detection and focusing of acoustic waves on the targets in the scattering medium. This technique involves the determination of the invariant of the time-reversal operator obtained by measurement of the scattering data in a pulse-echo mode. In this paper, a proposed approach based on the DORT method is developed to solve the acoustic inverse scattering problem of a small metallic scatterer. The proposed approach not only estimates the position of the scatterer, but also determines the physical properties of an unknown metallic scatterer such as the shape (cylinder or sphere), the material (density), and the size (radius) in an anisotropic scattering case. Theoretical and numerical simulation results are also studied and investigated to show that the proposed approach can simultaneously characterize all those properties of an unknown metallic scatterer. Moreover, the advantage of the proposed approach is to avoid the complex iterative scheme in solving the direct scattering problem and results in smaller computational load and faster implementation.     

Frequency domain finite-element and spectral-element acoustic wave modeling using absorbing boundaries and perfectly matched layer

Wave propagation modeling as a vital tool in seismology can be done via several different numerical methods among them are finite-difference, finite-element, and spectral-element methods (FDM, FEM and SEM). Some advanced applications in seismic exploration benefit the frequency domain modeling. Regarding flexibility in complex geological models and dealing with the free surface boundary condition, we studied the frequency domain acoustic wave equation using FEM and SEM. The results demonstrated that the frequency domain FEM and SEM have a good accuracy and numerical efficiency with the second order interpolation polynomials. Furthermore, we developed the second order Clayton and Engquist absorbing boundary condition (CE-ABC2) and compared it with the perfectly matched layer (PML) for the frequency domain FEM and SEM. In spite of PML method, CE-ABC2 does not add any additional computational cost to the modeling except assembling boundary matrices. As a result, considering CE-ABC2 is more efficient than PML for the frequency domain acoustic wave propagation modeling especially when computational cost is high and high-level absorbing performance is unnecessary.     

Design and Analysis of Planar Photonic Band Gap Devices

The need for a highly efficient numerical simulation platform for designing photonic band gap (PBG) structures is outlined in the context of various functional device topologies. In this paper we therefore introduce the Method of Auxiliary Sources (MAS) as a semi-analytical, frequency-domain method for computational optics, which has already proven its accuracy and efficiency in various other fields of electrodynamics. The proposed software package provides an easy-to-handle approach to full-wave analysis of two-dimensional (2D) PBG circuits, PBG-based antennas as well as to dense-integrated optics components that contain optical waveguides, scatterers, resonators and other functional elements. Experimental verifications of the numerical results have been conducted along large-scale prototypes in the microwave frequency range for several device topologies.     

Side-branch resonators modelling with Green׳s function methods

This paper deals with strategies for computing efficiently the propagation of sound waves in ducts containing passive components. In many cases of practical interest, these components are acoustic cavities which are connected to the duct. Though standard Finite Element software could be used for the numerical prediction of sound transmission through such a system, the method is known to be extremely demanding, both in terms of data preparation and computation, especially in the mid-frequency range. To alleviate this, a numerical technique that exploits the benefit of the FEM and the BEM approach has been devised. First, a set of eigenmodes is computed in the cavity to produce a numerical impedance matrix connecting the pressure and the acoustic velocity on the duct wall interface. Then an integral representation for the acoustic pressure in the main duct is used. By choosing an appropriate Green?s function for the duct, the integration procedure is limited to the duct–cavity interface only. This allows an accurate computation of the scattering matrix of such an acoustic system with a numerical complexity that grows very mildly with the frequency. Typical applications involving Helmholtz and Herschel–Quincke resonators are presented.     

Vibro-acoustic modelling of a railway bridge crossed by a train

This paper presents a vibro-acoustic modelling of a railway bridge excited by a moving train. The modelling of the bridge-train system is carried out by the modal superposition method taking into account the train mass, the viscoelastic suspension of the vehicles with an unlimited number of trucks and spans of the bridge. The numerical resolution of the coupled equations of motion is carried out by the Newmark’s method with an iterative process. we studied on one hand the influence of the track irregularities on the dynamic behaviour of the bridge-train system on the other hand the noise radiated by the bridge due to the passage of the train. The acoustic pressure is obtained by solving the wave equation which has as excitations source the bridge acceleration, they are considered as acoustic monopoles.     

Computation of acoustic radiation from vibrating structures in motion

The noises from a vibrating structure in motion are often encountered in engineering practice, for example, tire noise and pass-by noise of moving vehicles. Consequently, research on the radiation characteristics of moving acoustic source is of significance. In this paper, a new computational method based on the wave superposition approach is developed for the acoustic field from a vibrating structure in motion. It inherits the advantages of the wave superposition approach in the acoustic computation, and in which a method of moving simple sources is used to eliminate the influence of the Doppler effect. By the proposed method, the acoustic radiation from the moving vibrating structure can be calculated easily with the same implementation process as the conventional wave superposition approach performed in the stationary acoustic field. Finally, the validity and feasibility of the proposed method are verified by the numerical results.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.     

Observer-based beamforming algorithm for acoustic array signal processing

In the field of noise identification with microphone arrays, conventional delay-and-sum (DAS) beamforming is the most popular signal processing technique. However, acoustic imaging results that are generated by DAS beamforming are easily influenced by background noise, particularly for in situ wind tunnel tests. Even when arithmetic averaging is used to statistically remove the interference from the background noise, the results are far from perfect because the interference from the coherent background noise is still present. In addition, DAS beamforming based on arithmetic averaging fails to deliver real-time computational capability. An observer-based approach is introduced in this paper. This so-called observer-based beamforming method has a recursive form similar to the state observer in classical control theory, thus holds a real-time computational capability. In addition, coherent background noise can be gradually rejected in iterations. Theoretical derivations of the observer-based beamforming algorithm are carefully developed in this paper. Two numerical simulations demonstrate the good coherent background noise rejection and real-time computational capability of the observer-based beamforming, which therefore can be regarded as an attractive algorithm for acoustic array signal processing.     

A Consistent Characteristic Boundary Condition for General Fluid Mixture and Its Implementation in a Preconditioning Scheme

Characteristic boundary conditions that are capable of handling general fluid mixtures flow at all flow speeds are developed. The formulation is based on fundamental thermodynamics theories incorporated into an efficient preconditioning scheme in a unified manner. Local one-dimensional inviscid (LODI) relations compatible to the preconditioning system are proposed to obtain information carried by incoming characteristic waves at boundaries accurately. The approach has been validated against a variety of sample problems at a broad range of fluid states and flow speeds. Both acoustic waves and hydrodynamic flow features can pass through the boundaries of computational domain transparently without any unphysical reflection or spurious distortion. The approach can be reliably applied to fluid flows at extensive thermodynamic states and flow speeds in numerical simulations. Moreover, the use of the boundary condition shows to improve the computational efficiency.     

Deterministic-statistical analysis of a structural-acoustic system

The purpose of this paper is to develop an efficient approach for vibro-acoustic analysis. Being simple and representative, an exited plate-acoustic system is selected as a validation case for the vibro-acoustic analysis as the system presents one two-dimensional statistical component (modal dense structure panel—plate) connected to the other component (deterministic acoustic volume—cavity) through the area junction over a surface domain, rather than at a line boundary. Potential industrial applications of the system vibro-acoustic analysis would be in acoustic modelling of vehicle body panels such as the cabin roof panel, and door panels for the boom noise analysis.A new deterministic-statistical analysis approach is proposed from a combination or hybrid of deterministic analysis and statistical energy analysis (SEA) approaches. General theory of the new deterministic-statistical analysis approach is introduced. The main advantage of the new deterministic-statistical analysis approach is its possibility in place of the time consuming Monte Carlo simulation. In order to illustrate and validate the new deterministic-statistical analysis approach, three approaches of the deterministic analysis, the statistical energy analysis and the new deterministic-statistical analysis are then applied to conduct the plate-acoustic system modelling, and their results will be compared. The vibro-acoustic energy coupling characteristic of the plate-acoustic system will be studied. The most suitable frequency range for the new approach will be identified in consideration of computational accuracy, information and speed.     

A hybrid formulation for mid-frequency analysis of assembled structures

A new formulation able to predict the behaviour of structures in the mid-frequency range is presented in this paper. The mid-frequency field is a hybrid domain for which assembled structures exhibit simultaneously low- and high-frequency behaviours, depending on the material and geometrical properties of different subsystems. Thus, dealing with the mid-frequency field requires simulation methods which are able to account the differences in behaviour of different subsystems. The hybrid formulation is based on the coupling of two different formulations, the finite elements for the low-frequency behaving subparts and a probabilistic formulation, the smooth integral formulation, applied to the high-frequency subsystems. The hybrid method enables to correctly predict the deterministic response of the low-frequency parts which is not affected by randomness, and the smooth trend of the contributions of the high-frequency parts. The paper is concluded with several numerical examples computed for coupled one- and two-dimensional structures.     

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.     

复杂结构声辐射模态的计算

结构辐射阻矩阵的特征向量就是其声辐射模态,复杂结构的辐射阻矩阵获取困难阻碍了声辐射模态理论的进一步应用。为此提出一种基于等效源方法的复杂结构辐射阻矩阵计算方法,由此形成了一套完整的计算复杂结构声辐射模态方法。最后通过仿真算例对影响结构辐射阻矩阵计算精度的因素进行了分析,给出了复杂结构辐射阻矩阵构造的一般原则。理论推导和数值仿真算例表明所提出的复杂结构声辐射模态计算方法是有效的。该方法简化了复杂结构声辐射模态的计算,计算效率较高。     

Fast prediction of aerodynamic noise induced by the flow around a cylinder based on deep neural network

Accurate and fast prediction of aerodynamic noise has always been a research hotspot in fluid mechanics and aeroacoustics. The conventional prediction methods based on numerical simulation often demand huge computational resources, which are difficult to balance between accuracy and efficiency. Here, we present a data-driven deep neural network (DNN) method to realize fast aerodynamic noise prediction while maintaining accuracy. The proposed deep learning method can predict the spatial distributions of aerodynamic noise information under different working conditions. Based on the large eddy simulation turbulence model and the Ffowcs Williams-Hawkings acoustic analogy theory, a dataset composed of 1216 samples is established. With reference to the deep learning method, a DNN framework is proposed to map the relationship between spatial coordinates, inlet velocity and overall sound pressure level. The root-mean-square-errors of prediction are below 0.82 dB in the test dataset, and the directivity of aerodynamic noise predicted by the DNN framework are basically consistent with the numerical simulation. This work paves a novel way for fast prediction of aerodynamic noise with high accuracy and has application potential in acoustic field prediction.     

A fast wave superposition spectral method with complex radius vector combined with two-dimensional fast Fourier transform algorithm for acoustic radiation of axisymmetric bodies

Based on the two-dimensional fast Fourier transform (2D FFT) algorithm, a wave superposition spectral method with complex radius vector has been proposed to efficiently analyze the acoustic radiation from an axisymmetric body. First, the complex Fourier series are used along both circumferential and meridian directions, to expand the integral kernel function and unknown source strength density distributed function. Then, by means of the rectangular integral formula, the radiation sound pressure is described in the form of two-dimensional discrete Fourier transform and generalized through 2D FFT algorithm. Finally, several numerical examples are performed to verify the accuracy and efficiency of the present method. Comparing with the other existing analysis ways, the present method has three different characteristics: (1) there is no singular integral in the numerical computation; (2) the unique solution can be given for all eigen wavenumbers owing to the application of the virtual boundary technology with complex radius vector; and (3) the computational efficiency is improved remarkably because all Fourier terms are calculated simultaneously through 2D FFT algorithm.     

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.