轴对称体声振耦合的边界子波谱与有限元耦合方法

探讨了子波在Helmholtz积分方程及声振耦合中的应用,在建立了求解轴对称Helmholtz积分方程的子波谱方法的基础上,构造了轴对称子波谱与轴对称有限元的耦合方法,该方法可以处理轴对称问题的任意边界条件.进行了声振耦合问题的模态分析.     

Use of principle velocity patterns in the analysis of structural acoustic optimization

This work presents an application of principle velocity patterns in the analysis of the structural acoustic design optimization of an eight ply composite cylindrical shell. The approach consists of performing structural acoustic optimizations of a composite cylindrical shell subject to external harmonic monopole excitation. The ply angles are used as the design variables in the optimization. The results of the ply angle design variable formulation are interpreted using the singular value decomposition of the interior acoustic potential energy. The decomposition of the acoustic potential energy provides surface velocity patterns associated with lower levels of interior noise. These surface velocity patterns are shown to correspond to those from the structural acoustic optimization results. Thus, it is demonstrated that the capacity to design multi-ply composite cylinders for quiet interiors is determined by how well the cylinder be can designed to exhibit particular surface velocity patterns associated with lower noise levels.     

Automotive panel noise contribution modeling based on finite element and measured structural-acoustic spectra

The radiated noise contributions of automotive body panels to the interior sound pressure levels are modeled using an approximate spectral formulation that applies the theoretical interior acoustic sensitivity terms derived from a finite element model and measured spatial-averaged structural-acoustic spectra. The finite element calculation is validated by comparison to a set of experimental acoustic transfer functions. A measurement set-up for the sound intensity and structural-acoustic response is applied to acquire the cross and auto power spectra needed to predict the relative mean-squared velocity term of each control plane near the panel surface, and to obtain the individual panel contribution function. The proposed approach also computes the noise spectra in 1/12 octave band form at selected positions in the passenger compartment, which matches well with the overall experimental results. Through an actual passenger car application, the approximate computational scheme is proven to be generally quite robust and effective for analyzing higher frequency interior noise problems.     

A boundary integral equation method using auxiliary interior surface approach for acoustic radiation and scattering in two dimensions

This paper presents an effective solution method for predicting acoustic radiation and scattering fields in two dimensions. The difficulty of the fictitious characteristic frequency is overcome by incorporating an auxiliary interior surface that satisfies certain boundary condition into the body surface. This process gives rise to a set of uniquely solvable boundary integral equations. Distributing monopoles with unknown strengths over the body and interior surfaces yields the simple source formulation. The modified boundary integral equations are further transformed to ordinary ones that contain nonsingular kernels only. This implementation allows direct application of standard quadrature formulas over the entire integration domain; that is, the collocation points are exactly the positions at which the integration points are located. Selecting the interior surface is an easy task. Moreover, only a few corresponding interior nodal points are sufficient for the computation. Numerical calculations consist of the acoustic radiation and scattering by acoustically hard elliptic and rectangular cylinders. Comparisons with analytical solutions are made. Numerical results demonstrate the efficiency and accuracy of the current solution method.     

Axial time-averaged acoustic radiation force on a cylinder in a nonviscous fluid revisited

Objective

The present research examines the acoustic radiation force of axisymmetric waves incident upon a cylinder of circular surface immersed in a nonviscous fluid. The attempt here is to unify the various treatments of radiation force on a cylinder with arbitrary radius and provide a formulation suitable for any axisymmetric incident wave.

Method and results

Analytical equations are derived for the acoustic scattering field and the axial acoustic radiation force. A general formulation for the radiation force function, which is the radiation force per unit energy density per unit cross-sectional surface, is derived. Specialized forms of the radiation force function are provided for several types of incident waves including plane progressive, plane standing, plane quasi-standing, cylindrical progressive diverging, cylindrical progressive converging and cylindrical standing and quasi-standing diverging waves (with an extension to the case of spherical standing and quasi-standing diverging waves incident upon a sphere).

Significance and some potential applications

This study may be helpful essentially due to its inherent value as a canonical problem in physical acoustics. Potential applications include particle manipulation of cylindrical shaped structures in biomedicine, micro-gravity environments, fluid dynamics properties of cylindrical capillary bridges, and the micro-fabrication of new cylindrical crystals to better control light beams.     

Numerical difficulties with boundary element solutions of interior acoustic problems

Although boundary element methods have been applied to interior problems for many years, the numerical difficulties that can occur have not been thoroughly explored. Various authors have reported low-frequency breakdowns and artificial damping due to discretization errors. In this paper, it is shown through a simple example problem that the numerical difficulties depend on the solution formulation. When the boundary conditions are imposed directly, the solution suffers from artificial damping, which may potentially lead to erroneous predictions when boundary element methods are used to evaluate the performance of damping materials. This difficulty can be alleviated by first computing an impedance or admittance matrix, and then using its reactive component to derive the solution for the acoustic field. Numerical computations are used to demonstrate that this technique eliminates artificial damping, but does not correct errors in the reactive components of the impedance or admittance matrices, which then causes nonexistence and nonuniqueness difficulties at the interior resonance frequencies for hard-wall and pressure release boundary conditions, respectively. It is shown that the admittance formulation is better suited to boundary element computations for interior problems because the resonance frequencies for pressure release boundary conditions do not begin until the smallest dimension of the boundary surface is at least one half the acoustic wavelength. Aside from producing much more accurate predictions, the admittance matrix is also much easier to interpolate at low frequencies due to the absence of interior resonances. For the example problem considered, only the formulation using the reactive component of the admittance matrix produces accurate solutions as long as the surface element discretization satisfies the standard six-element-per-wavelength rule.     

A new integral equation formulation for an axisymmetric structure

A new method is developed to compute the acoustic field outside an axisymmetric structure from the normal velocity values on the surface. Surface pressure and normal velocity are expanded in a series of functions that are orthonormal on the surface of the structure and have a constant ratio of pressure to normal derivative of pressure at vanishing frequency. The Helmholtz integral equation is next used to compute the field everywhere outside the structure. The method is tested by applying it to scattering from a rigid cylinder with hemispherical endcaps. The series is shown to converge very rapidly.     

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.     

Data completion method for the characterization of sound sources

This paper presents a different approach to solve the inverse acoustic problem. This problem is an "ill-posed" problem since the solution is very sensitive to measurement precision. A classical way to solve this problem consists in inversing a propagation operator which relates structure quantities (acoustic pressures or gradients) to near-field quantities (acoustic pressures or gradients). This can be achieved by using near-field acoustical holography (NAH) in separable coordinate systems. In order to overcome this limitation, the inverse boundary element method (IBEM) can be implemented to recover all acoustic quantities in a three-dimensional space and on an arbitrary three-dimensional source surface. In this paper, the data completion method (DCM) is developed: the acoustic gradients and pressures are known on a surface surrounding the source, but are unknown on its structure. The solution is given by the resolution of the Helmholtz formulation applied on the empty domain between the two boundaries made by the measurements quantities and the structure of the source. The conventional method applies directly the integral formulation for the empty domain. Another way of solving this Helmholtz formulation can be achieved by splitting it in two well-posed subproblems in a Steklov-Poincare?'s formulation. The data completion method allows one to solve the problem with acoustic perturbations due to sources on the exterior domain, or due to a confined domain, without altering the results.     

On the effect of viscosity and thermal conductivity on sound propagation in ducts: A re-visit to the classical theory with extensions for higher order modes and presence of mean flow

The dispersion equation for the axisymmetric modes of viscothermal acoustic wave propagation in uniform hard-walled circular ducts containing a quiescent perfect gas is classical. This has been extended to cover the non-axisymmetric modes and real fluids in contemporary studies. The fundamental axisymmetric mode has been the subject of a large number of studies proposing approximate solutions and the characteristics of the propagation constants for narrow and wide ducts with or without mean flow is well understood. In contrast, there are only few publications on the higher order modes and the current knowledge about their propagation characteristics is rather poor. On the other hand, there is a void of papers in the literature on the effect of the mean flow on the quiescent modes of propagation. The present paper aims to contribute to the filling of these gaps to some extent. The classical theory is re-considered with a view to cover all modes of acoustic propagation in circular ducts carrying a real fluid moving axially with a uniform subsonic velocity. The analysis reveals a new branch of propagation constants for the axisymmetric modes, which appears to have escaped attention hitherto. The solution of the governing wave equation is expressed in a modal transfer matrix form in frequency domain and numerical results are presented to show the effects over wide ranges of frequency, viscosity and mean flow parameters on the propagation constants. The theoretical formulation allows for the duct walls to have finite impedance, but no numerical results are presented for lined ducts or ducts carrying a sheared mean flow.     

Transient response of an acoustic medium by an excited submerged spherical shell.

An exact closed-form solution is obtained for the transient response of an acoustic fluid due to an excited submerged spherical shell. Step axisymmetric stress acting on the inner surface of the shell is expanded into modal Fourier components in which for each modal stress systematic exact expressions for modal radial displacement of the shell, fluid pressure, fluid particle velocity and displacement are presented. The superiority of the formulations herein is its applicability to any time and spatial distance in the fluid as well as on the shell. The formulation also does not require additional numerical computations associated with other methods.     

Measured wavenumber: frequency spectrum associated with acoustic and aerodynamic wall pressure fluctuations

Direct measurements of the wavenumber-frequency spectrum of wall pressure fluctuations beneath a turbulent plane channel flow have been performed in an anechoic wind tunnel. A rotative array has been designed that allows the measurement of a complete map, 63×63 measuring points, of cross-power spectral densities over a large area. An original post-processing has been developed to separate the acoustic and the aerodynamic exciting loadings by transforming space-frequency data into wavenumber-frequency spectra. The acoustic part has also been estimated from a simple Corcos-like model including the contribution of a diffuse sound field. The measured acoustic contribution to the surface pressure fluctuations is 5% of the measured aerodynamic surface pressure fluctuations for a velocity and boundary layer thickness relevant for automotive interior noise applications. This shows that for aerodynamically induced car interior noise, both contributions to the surface pressure fluctuations on car windows have to be taken into account.     

The forward and backward projection of acoustic fields from axisymmetric ultrasonic radiators using impulse response and Hankel transform techniques

A generalized impulse response formulation to evaluate the harmonic pressure field of ultrasonic planar vibrators having axisymmetric nonuniform surface velocity distributions is presented. The harmonic pressure is expressed as a Fourier transform of a generalized impulse response which is a function of the spatially nonuniform velocity of the vibrator. A backward projection method is then developed to reconstruct the normal surface velocity of axisymmetric vibrators from harmonic field pressures using an angular spectrum or Hankel transform formulation. The numerical accuracy of the backward projection technique is evaluated using the impulse response formulation to evaluate the pressure fields for several velocity distributions on disk vibrators. Experiments were performed to reconstruct the velocity distributions over the surface of a uniformly driven piezoelectric ceramic disk and ceramic ring using farfield measurements of the complex pressure. The experimental results were in good agreement with theoretical results based on the electrode patterns of the transducers.     

A finite element approach for predicting nozzle admittances

A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and irrotational, and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation by using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. The boundary conditions are described for both straight and curved sonic lines. A two dimensional FEM with linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.     

An axisymmetric poroelastic finite element formulation.

In the past, various two- and three-dimensional Cartesian, poroelastic finite element formulations have been proposed and demonstrated. Here an axisymmetric formulation of a poroelastic finite element is presented. The intention of this work was to develop a finite element formulation that could easily and efficiently model axisymmetric sound propagation in circular structures having arbitrary, axially dependent radii, and that are lined or filled with elastic porous sound absorbing materials such as foams. The formulation starts from the Biot equations for an elastic porous material expressed explicitly in axisymmetric form. By following a standard finite element development, a u-U formulation results. Procedures for coupling the axisymmetric elements to an adjacent acoustical domain are described, as are the boundary conditions appropriate for unfaced foams. Calculations described here show that the present formulation yields predictions as accurate as a Cartesian, three-dimensional model in much reduced time. Predictions made using the present model are also compared with measurements of sound transmission through cylindrical foam plugs, and the predicted results are shown to agree well with the measurements. Good agreement was also found in the case of sound transmission through a conical foam plug.     

Analytic formulation and numerical implementation of an acoustic pressure gradient prediction

Two new analytical formulations of the acoustic pressure gradient have been developed and implemented in the PSU-WOPWOP rotor noise prediction code. The pressure gradient can be used to solve the boundary condition for scattering problems and it is a key aspect to solve acoustic scattering problems. The first formulation is derived from the gradient of the Ffowcs Williams–Hawkings (FW–H) equation. This formulation has a form involving the observer time differentiation outside the integrals. In the second formulation, the time differentiation is taken inside the integrals analytically. This formulation avoids the numerical time differentiation with respect to the observer time, which is computationally more efficient. The acoustic pressure gradient predicted by these new formulations is validated through comparison with available exact solutions for a stationary and moving monopole sources. The agreement between the predictions and exact solutions is excellent. The formulations are applied to the rotor noise problems for two model rotors. A purely numerical approach is compared with the analytical formulations. The agreement between the analytical formulations and the numerical method is excellent for both stationary and moving observer cases.     

Hybrid near-field acoustic holography

Hybrid near-field acoustical holography (NAH) is developed for reconstructing acoustic radiation from an arbitrary object in a cost-effective manner. This hybrid NAH is derived from a modified Helmholtz equation least squares (HELS) formula that expands the acoustic pressure in terms of outgoing and incoming waves. The expansion coefficients are determined by solving an overdetermined linear system of equations obtained by matching the assumed-form solution to measured acoustic pressures through the least squares. Measurements are taken over a conformal surface around a source at close range so that the evanescent waves can be captured. Next, the modified HELS is utilized to regenerate as much acoustic pressures on the conformal surface as necessary and take them as input to the Helmholtz integral formulation implemented numerically by boundary element method (BEM). The acoustic pressures and normal velocities on the source surface are reconstructed by using a modified Tikhnov regularization (TR) with its regularization parameter determined by generalized cross validation (GCV) method. Results demonstrate that this hybrid NAH combines the advantages of HELS and inverse BEM. This is because a majority of the input data are regenerated but not measured, thus the efficiency of reconstruction is greatly enhanced. Meanwhile, the accuracy of reconstruction is ensured by the Helmholtz integral theory and modified TR together with GCV method, provided that HELS converges fast enough on the measurement surface. Numerical examples of reconstructing acoustic quantities on the surface of a simplified engine block are demonstrated. [Work supported by NSF.]     

Acoustical cavity excitation.

A tutorial exposition is given of cavity excitations and it is shown that relatively simple analytical solutions for the excitation of acoustic fields in cavities with lossy walls can be obtained for a broad class of such cavities. The standard analytic procedure whereby one expresses the acoustic field as a sum over the eigenmodes of the cavity is reviewed and its limitations are pointed out. An alternate procedure is to seek to make use of Green's functions depending on a single coordinate and to not express such Green's functions as sums over eigenfunctions. Instead, one expresses the Green's function as a closed form function that has a discontinuous slope at the value of the coordinate corresponding to the source location. Although the class of cavities for which this approach is possible is limited, it is sufficiently broad to be of intrinsic interest and the results may be helpful to those seeking benchmark solutions and analytic insight. Displayed examples are for cavities whose shapes are spheres and cylinders, with interior monopole, dipole, line, and distributed sources. The formulation yields simple expressions for the total power absorbed by the cavity walls.     

封闭腔体噪声控制中亥姆霍兹共振器的优化设计方法

提出用亥姆霍兹共振器控制声腔内噪声时计算共振器最优阻尼比和最优工作带宽的理论公式,并进行实验验证。首先,建立共振器与待控腔体的声学耦合方程,以最小化腔体内目标声压幅值为参考,对共振器的阻尼比和工作带宽进行理论分析,求出最优阻尼比和最优工作带宽的计算公式。接着,提出在声腔噪声控制中使用最优亥姆霍兹共振器的实施步骤。最后,以一维声学腔体内的噪声为控制对象,通过对比控制前后的理论结果与实测数据,验证最优阻尼比和最优工作带宽的理论公式。结果表明,本文开发的亥姆霍兹共振器优化设计方法能准确地预报共振器的最优阻尼比与最优工作带宽,在声腔中低频噪声控制中有广泛的应用前景。