Acoustic scattering by multiple elliptical cylinders using collocation multipole method

This paper presents the collocation multipole method for the acoustic scattering induced by multiple elliptical cylinders subjected to an incident plane sound wave. To satisfy the Helmholtz equation in the elliptical coordinate system, the scattered acoustic field is formulated in terms of angular and radial Mathieu functions which also satisfy the radiation condition at infinity. The sound-soft or sound-hard boundary condition is satisfied by uniformly collocating points on the boundaries. For the sound-hard or Neumann conditions, the normal derivative of the acoustic pressure is determined by using the appropriate directional derivative without requiring the addition theorem of Mathieu functions. By truncating the multipole expansion, a finite linear algebraic system is derived and the scattered field can then be determined according to the given incident acoustic wave. Once the total field is calculated as the sum of the incident field and the scattered field, the near field acoustic pressure along the scatterers and the far field scattering pattern can be determined. For the acoustic scattering of one elliptical cylinder, the proposed results match well with the analytical solutions. The proposed scattered fields induced by two and three elliptical–cylindrical scatterers are critically compared with those provided by the boundary element method to validate the present method. Finally, the effects of the convexity of an elliptical scatterer, the separation between scatterers and the incident wave number and angle on the acoustic scattering are investigated.     

部分浸没圆柱壳声固耦合计算的半解析法研究

部分浸没圆柱壳-流场耦合系统的声振分析是一种典型的半空间域内声固耦合问题,其振动及声学计算目前主要依赖于数值方法求解,但无论从检验数值法还是从机理上揭示其声固耦合特性,解析或半解析方法的发展都是不可或缺的.本文提出了一种半解析方法,先将声场坐标系建立在自由液面上,采用正弦三角级数来满足自由液面上的声压释放边界条件;接着基于二维Flügge薄壳理论建立了以圆柱圆心为坐标原点的壳-液耦合系统的控制方程;然后再利用Galerkin法处理声固耦合界面的速度连续条件,推导得到声压幅值与壳体位移幅值之间的关系矩阵并求解该耦合系统的振动和水下声辐射.与有限元软件Comsol进行了耦合系统自由、受迫振动和水下辐射噪声计算结的对比分析,表明本文方法准确可靠.本文的研究为解析求解弹性结构与声场部分耦合的声振问题提供了新的思路.     

插管结构对膨胀腔消声器声学性能的影响

将数值模态匹配法(NMM)拓展应用于计算和分析外插管膨胀腔消声器的声学性能,推导了相应的理论公式并编写了计算程序。使用二维有限元法提取横向波数和本征向量,应用模态匹配法计算消声器的传递损失。使用数值模态匹配法和三维有限元法(FEM)研究了插管长度和进出口位置对带有外插进出口管椭圆形非同轴膨胀腔消声器声学性能的影响,两种方法计算结果吻合良好,从而验证了本文数值模态匹配法的正确性。研究结果表明,设置特定的插管长度和进出口位置可以消除消声器的通过频率,进而改善消声器中低频的消声性能。      

Acoustic cavity modes in lens-shaped structures

A method for solving exactly the Helmholtz equation in parabolic rotational coordinates is presented using separability of the eigenfunctions and the Frobenius power series expansion technique. Two examples of interest in acoustics are considered and analyzed quasianalytically: The acoustic pressure in a cavity defined by two paraboloids (forming a lens-shaped structure) with (I) rigid wall boundary conditions and (II) pressure-release boundaries. The rigid-wall (pressure-release) acoustic enclosure problem is a Neumann (Dirichlet) boundary condition problem. In both cases, eigenfunctions and eigenmodes are calculated and the shape dependence of the eigenvalue for the ground state is examined.     

Gauss声束对离轴椭圆柱的声辐射力矩

利用部分波展开法求解得到了Gauss声束入射下刚性和非刚性椭圆柱的声散射系数,推导了一般情况下的声辐射力矩表达式.在此基础上,通过一系列数值仿真详细分析了离轴距离、入射角度和束腰半径对声辐射力矩的影响.结果表明:正向与负向声辐射力矩均可以在一定条件下存在;低频情况下刚性椭圆柱比非刚性椭圆柱更容易产生较强的声辐射力矩;特定频率的入射声场可以激发出非刚性椭圆柱不同阶的共振散射模式,因而非刚性椭圆柱的声辐射力矩峰值与频率的关系更密切;增加束腰半径有利于扩大散射截面,进而增加椭圆柱的声辐射力矩.该研究结果预期可以为利用声辐射力矩实现粒子的可控旋转和流体黏度的反演提供一定的理论指导.     

Acoustic analysis of a rectangular cavity with general impedance boundary conditions

A Fourier series method is proposed for the acoustic analysis of a rectangular cavity with impedance boundary conditions arbitrarily specified on any of the walls. The sound pressure is expressed as the combination of a three-dimensional Fourier cosine series and six supplementary two-dimensional expansions introduced to ensure (accelerate) the uniform and absolute convergence (rate) of the series representation in the cavity including the boundary surfaces. The expansion coefficients are determined using the Rayleigh-Ritz method. Since the pressure field is constructed adequately smooth throughout the entire solution domain, the Rayleigh-Ritz solution is mathematically equivalent to what is obtained from a strong formulation based on directly solving the governing equations and the boundary conditions. To unify the treatments of arbitrary nonuniform impedance boundary conditions, the impedance distribution function on each specified surface is invariantly expressed as a double Fourier series expansion so that all the relevant integrals can be calculated analytically. The modal parameters for the acoustic cavity can be simultaneously obtained from solving a standard matrix eigenvalue problem instead of iteratively solving a nonlinear transcendental equation as in the existing methods. Several numerical examples are presented to demonstrate the effectiveness and reliability of the current method for various impedance boundary conditions, including nonuniform impedance distributions.     

Vibro-acoustic analysis of a rectangular cavity bounded by a flexible panel with elastically restrained edges

A coupled system consisting of an acoustic cavity and an elastic panel is a classical problem in structural acoustics and is typically analyzed using modal approaches based on in vacuo structural modes and the rigidly walled acoustic modes which are pre-determined based on separate component models. Such modeling techniques, however, tend to suffer the following drawbacks or limitations: (a) a panel is only subjected to ideal boundary conditions such as the simply supported, (b) the coupling between the cavity and panel is considered weak, and (c) the particle velocity cannot be correctly predicted from the pressure gradient on the contacting interface, to name a few. Motivated by removing these restrictions, this paper presents a general method for the vibro-acoustic analysis of a three-dimensional (3D) acoustic cavity bounded by a flexible panel with general elastically restrained boundary conditions. The displacement of the plate and the sound pressure in the cavity are constructed in the forms of standard two-dimensional and 3D Fourier cosine series supplemented by several terms introduced to ensure and accelerate the convergence of the series expansions. The unknown expansions coefficients are treated as the generalized coordinates and determined using the Rayleigh-Ritz procedure based on the energy expressions for the coupled structural acoustic system. The accuracy and effectiveness of the proposed method are demonstrated through numerical examples and comparisons with the results available in the literature.     

Calculation of the Critical Distances in a System of Two Colliding Nuclei Beyond the Monopole Approximation

A method for calculating electronic levels in compact superheavy nuclear quasimolecules based on solving the Dirac equation in spherical coordinates using a multipole expansion of a two-center potential is developed. It is shown that, for internuclear distances up to ∼100 fm, such a technique enables fast convergence, which allows one to calculate the electronic energy levels with an accuracy better than 10−6. Moreover, all the multipole moments can be calculated analytically. The critical distances between the similar colliding nuclei have been calculated in the range Z ∼ 87–100 for bottom electronic levels: even 1σg and odd 1σu, respectively.     

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 Chebyshev–Lagrangian method for acoustic analysis of a rectangular cavity with arbitrary impedance walls

A general Chebyshev–Lagrangian method is proposed to obtain the analytical solution for a rectangular acoustic cavity with arbitrary impedance boundary conditions. The originality of the present paper is the successful attempt of applying orthogonal polynomials, such as Chebyshev polynomials of the first kind, to the analysis of a rectangular sound field with general wall impedance. The sound pressure is uniformly expressed as triplicate Chebyshev polynomial series which is independent in each direction. The Chebyshev polynomial series solution is obtained using the Rayleigh–Ritz procedure after considering the influence of boundary impedance on the cavity as the work done by the impedance surfaces in the Lagrangian function. The accuracy and reliability of the proposed method are validated against the analytical solutions and some numerical results available in the literature. Excellent orthogonality and complete properties of the Chebyshev polynomials ensure the rapid convergence, numerical stability, high accuracy of the current solution. The simplicity and low computational cost of the present approach make it preferable to obtain the results of complex models even in the relative high frequency range by choosing enough truncated terms in the sound pressure expression. Numerous cases with various uniform or non-uniform impedance boundary conditions are analyzed numerically and some of the results can be used as benchmark. It is shown that the impedance boundary condition can effectively influence or modify the acoustic characteristics and response of a cavity.     

多重散射方法研究轴对称空腔覆盖层的声学特性

发展了一种多重散射方法研究声学覆盖层的半数值半解析模型,分析了影响轴对称空腔结构声学性能的主要能量耗散机制。在球坐标条件下推导出轴对称空腔结构的位移和应力场基函数,通过对空腔表面基函数的数值积分,得到散射波和入射波之间的传输矩阵方程,结合分层介质声传播理论计算了周期性空腔结构覆盖层的反射、透射和吸声性能。研究结果表明;空腔共振是低频能量耗散的主要形式,边界条件对材料空腔结构的谐振特性影响很大,利用双空腔耦合共振可以拓宽材料的低频吸声频带;背衬对材料的高频吸声影响较小,材料的高频能量损耗取决于空腔的散射和波型转换特性。      

Thickness optimization of a multilayered structure on the coupling surface between a structure and an acoustic cavity

This paper describes a new design method to optimize thickness distribution of a multilayered structure which is located on the coupling surface between a structure and an acoustic cavity. The design method is based on the concept of the density approach in topology optimization incorporating a transfer matrix for a multilayered structure that includes a poroelastic media layer. The one-dimensional transfer matrix adopted here is an approximate representation addressing vibro-acoustic effects inherent in a multilayered structure, and balances calculation resources and desired accuracy. Applying the transfer matrix representation as boundary conditions on the coupling surface between a structure and an acoustic cavity, the modified equilibrium equation of the vibro-acoustic system is derived which is approximately but efficiently solved by the modal approach. In this study, the problem of minimizing the acoustic pressure within the cavity over the prescribed frequency range is formulated under the volume constraint of the poroelastic media layer. The continuous approximation of thickness distribution is assumed, and the thickness of the poroelastic media layer at each nodal point is chosen as design variables. Numerical results show that an acoustic response is significantly reduced by the optimal thickness distribution having a total weight equal to or less than that in the initial uniform thickness. These demonstrate that the proposed method is effective to design the optimal thickness distribution of a multilayered structure.     

Extension of the generalized multipole technique to anisotropic medias

This paper presents an extension of the generalized multipole technique (GMT) for 2D anisotropic scatterers. New expansions similar to the Bessel multipole expansion are derived for arbitrary anisotropic media. Numerical simulations prove the accuracy and the rapid convergence of these expansions. As the results obtained are extremely accurate, this technique is most helpful for the evaluation of reference solutions and for the understanding of the physical interaction of light with arbitrary anisotropic media.     

Vibro-acoustic response of an elliptical plate-cavity coupled system to external shock loads

A general fully coupled three-dimensional vibro-acoustic model is developed to investigate the forced non-stationary acousto-structural response of a thin elastic plate of elliptical planform which is backed by a reverberant, rigid, and finite (closed) elliptic cylindrical acoustic enclosure, while under the action of general external transverse loads of arbitrary temporal and spatial variations. The Laplace transform with respect to the time coordinate is invoked, and the classical method of separation of variables in elliptic coordinates is used to obtain the transformed solutions as a linear combination of even and odd modes in terms of products of radial and angular Mathieu functions. A linear system of coupled algebraic equations is ultimately obtained, which is truncated and then solved numerically by implementing Durbin’s numerical Laplace transform inversion scheme. Detailed numerical simulations are conducted for the temporal histories of plate center-point displacement and on-axis cavity acoustic pressure for air-coupled elliptic aluminum plates of selected aspect ratios when subjected to external loadings of practical interest (i.e., an impulsive point load, a uniformly distributed pulse load, and a blast load). Also, acoustic radiation into the backing enclosure is examined by using appropriate 2D images of the internal sound field for selected cavity depths and plate eccentricities. The presented results confirm that the acousto-elastic characteristics of the coupled plate-cavity system are significantly influenced by the plate aspect ratio, cavity depth and the transverse loading configuration. Validity of the work is established through the computations made by using a commercial finite element package.     

不确定声场分析的二阶区间摄动有限元法

针对一阶区间摄动有限元法在声场参数不确定程度增大时误差过大的缺陷,在二阶Taylor展开的基础上推导了声学二阶区间摄动有限元法,并将其应用于区间不确定声场的声压响应分析。该方法先对声学区间有限元方程的声压响应向量进行二阶Taylor展开,获取声压响应的二阶近似响应向量;再根据二次函数极值定理获得声压响应向量的上下界。二维管道声场与轿车声腔模型的数值分析算例表明,与一阶区间摄动有限元法相比,二阶区间摄动有限元法有效提高了计算精度。因此二阶区间摄动有限元适合不确定度更大的区间不确定声场声压响应分析,具有良好的工程应用前景。      

任意形状金属腔体内电磁脉冲的三维计算

提出了模拟任意形状腔体中的内电磁脉冲的三维直角坐标系时域有限差分(FDTD)算法。该算法采用FDTD共形网格技术模拟任意形状腔体的边界,可以解决腔体内非对称的边界问题。推导了射线斜入射的差分方程,进行了三维数值计算,并采用直角坐标系FDTD算法和柱坐标系FDTD算法计算了射线斜入射圆柱腔体产生的内电磁脉冲,二者吻合很好,验证了直角坐标系FDTD算法正确性。     

水下近水面锥柱组合壳声固耦合效应的半解析求解方法

针对水下近水面锥柱组合壳声固耦合多借助于数值方法求解的现状,本文提出一种半解析方法从机理上分析此类问题。首先基于能量泛函和Sanders壳体理论、虚拟弹簧法以及力与力矩平衡条件建立锥柱组合壳的结构模型;然后采用Legendre谱元法和二维傅里叶变换得到含自由液面的水下声场模型;最后由非线性迭代法和高斯积分求解耦合系统声振控制方程。通过与参考文献和数值方法结果的对比,验证了本文方法的收敛性、正确性和可靠性。研究结果表明,结构参数、浸没深度和激励频率与远场辐射声压密切相关。本文工作可推广到水下含内部结构的复杂旋转组合壳在不同结构边界及声边界下的声固耦合问题。      

A mode matching approach for modeling two dimensional porous grating with infinitely rigid or soft inclusions

This work investigates the acoustical properties of a multilayer porous material in which periodic inclusions are embedded. The material is assumed to be backed by a rigid wall. Most of the studies performed in this field used the multipole method and are limited to circular shape inclusions. Here, a mode matching approach, more convenient for a layered system, is adopted. The inclusions can be in the form of rigid scatterers of an arbitrary shape, in the form of an air-filled cavity or in the form of a porous medium with contrasting properties. The computational approach is validated on simple geometries against other numerical schemes and with experimental results obtained in an anechoic room on a rigid grating embedded in a porous material made of 2 mm glass beads. The method is used to study the acoustic absorption behavior of this class of materials in the low frequency range and at a range of angles of incidence.     

Shape design sensitivity analysis for the radiated noise from the thin-body

Many industrial applications generally use thin-body structures in their design. To calculate the radiated noise from vibrated structure including thin bodies, the conventional boundary element method (BEM) using the Helmholtz integral equation is not an effective resolution. Thus, many researchers have studied to resolve the thin-body problem in various physical fields. No major study in the design sensitivity analysis (DSA) fields for thin-body acoustics, however, has been reported.A continuum-based shape DSA method is presented for the radiated noise from the thin-body. The normal derivative integral equation is employed as an analysis formulation. And, for the acoustic shape design sensitivity formulation, the equation is differentiated directly by using material derivative concept. To solve the normal derivative integral equation, the normal velocities on the surface should be calculated. In the acoustic shape sensitivity formulation, not only the normal velocities on the surface are required but also derivative coefficients of the normal velocities (structural shape design sensitivity) are also required as the input. Hence, the shape design sensitivity of structural velocities on the surface, with respect to the shape change, should be calculated. In this research, the structural shape design sensitivities are also obtained by using a continuum approach. And both a modified interpolation function and the Cauchy principle value are used to regularize the singularities generated from the acoustic shape design sensitivity formulation.A simple annular disk is considered as a numerical example to validate the accuracy and efficiency of the shape design sensitivity equations derived in this research. The commercial BEM code, SYSNOISE, is utilized to confirm the results of the developed in-house code based on a normal derivative integral equation. To validate the calculated design sensitivity results, central finite difference method (FDM) is employed. The error between FDM and the analytical result are less than 3%. This comparison demonstrates that the proposed design sensitivities of the radiated pressure are very accurate.