阶梯圆柱形耦合声场的特征正交-里兹法建模及声学特性分析

针对阶梯圆柱形耦合声场建模问题,提出基于特征正交-里兹能量原理的声学建模方法.该方法利用二维特征正交多项式和周向傅里叶级数表征阶梯圆柱形耦合声场子分段的声压函数,从能量角度考虑邻近子声场间声学连续性条件,并结合里兹法获得耦合声场的声学特性.基于本建模方法对不同分段的耦合声场开展声学特性分析,结果表明,本建模方法在保证计算准确性的基础上有效提高了计算效率,且对任意阶梯分段的圆柱形耦合声场普遍适用;圆柱形耦合声场固有频率会随着腔体外径增大而普遍增大,而腔深的影响规律相反;降低声学边界阻抗可抑制声学响应幅值,为此类声场的噪声控制提供了设计依据.     

可渗透壁对翼型声涡相互作用下流场及声场的影响

由仿生学原理构建的可渗透翼型对湍流气动噪声抑制作用已展现良好的应用前景。对NACA 0012可渗透翼型和实体翼型进行了数值计算,得到了声涡相互作用下气动噪声声场和流场,分析了可渗透壁对翼型流场和声场的影响。研究表明,相对实体翼型,可渗透壁通过减小声源强度降低了主纯音噪声声压级幅值和远场总声压级,消除了高阶离散纯音,但对噪声的指向性没有较大改变。进一步的流场分析表明,可渗透壁对翼型气动性能影响不大的情况下能够降低边界层扰动和翼型后缘大尺度涡旋强度,并推迟分离泡转捩和再附位置。     

Predicting the wind noise from the pantograph cover of a train

Finite element and boundary element calculations are combined to predict the flow noise radiated from a 1/10th-scale model of an aerodynamic cover used around the pantograph on a train at 250 km h⁻¹. The solutions of the unsteady air flow over the cover and the resulting sound propagation are divided into two parts in order to keep the problem tractable. First the unsteady fluid flow is solved using large-eddy simulation (LES). The pressure histories on the cover are then used to predict the radiated sound, using a boundary element method to solve the Helmholtz equation. The result thus leans heavily on assumptions about the coupling of the two solutions, the propagation of sound in a disturbed medium and the efficacy of LES. The predicted sound pressure levels are compared with experimental measurements made in an anechoic wind tunnel. © 1997 John Wiley & Sons, Ltd.     

Dirichlet and Neumann boundary conditions for the pressure poisson equation of incompressible flow

In a recent paper Gresho and Sani showed that Dirichlet and Neumann boundary conditions for the pressure Poisson equation give the same solution. The purpose of this paper is to confirm this (for one case at least) by numerically solving the pressure equation with Dirichlet and Neumann boundary conditions for the inviscid stagnation point flow problem. The Dirichlet boundary condition is obtained by integrating the tangential component of the momentum equation along the boundary. The Neumann boundary condition is obtained by applying the normal component of the momentum equation at the boundary. In this work solutions for the Neumann problem exist only if a compatibility condition is satisfied. A consistent finite difference procedure which satisfies this condition on non-staggered grids is used for the solution of the pressure equation with Neumann conditions. Two test cases are computed. In the first case the velocity field is given from the analytical solution and the pressure is recovered from the solution of the associated Poisson equation. The computed results are identical for both Dirichlet and Neumann boundary conditions. However, the Dirichlet problem converges faster than the Neumann case. In the second test case the velocity field is computed from the momentum equations, which are solved iteratively with the pressure Poisson equation. In this case the Neumann problem converges faster than the Dirichlet problem.     

平板湍流边界层的声辐射

基于平板湍流边界层的壁压起伏波数—频率谱 ,给出了一种湍流边界层声辐射的估算方法 ,并对光滑平板湍流边界层和平板表面粗糙度引起的湍流边界层声辐射进行了分析。结果表明 :湍流边界层声辐射是一种四极子声辐射 ,且其辐射声能集中于平板表面粗糙度引起的湍流边界层声辐射 ;光滑平板湍流边界层的声辐射也不可忽略。     

Influence of acoustic perturbations on the flow structure in a boundary layer with adverse pressure gradient

An experimental Investigation was made of the flow in a laminar boundary layer in a region of adverse pressure gradient in the presence of an acoustic field. There is a significant rearrangement of the laminar flow under the influence of the sound. The strong influence of the sound on the average flow is due to the excitation of vorticity perturbations, but this does not completely explain the phenomenon.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 48–52, March–April, 1983.     

Diffraction of a sound wave by a nonmoving plate

The problem of determining the velocity field excited by a sound wave impinging on a plate at rest is analyzed as an initial- and boundary-value problem with a movable boundary for the two-dimensional wave equation. The latter problem is solved by the formulation and inversion of integral equations of the Volterra type. The solution is obtained in closed form for any angle of inclination of the incident wave relative to the plate surface and is represented by recursion relations allowing for the influence of any number of diffracted waves generated in succession at the plate boundary.Moscow. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 123–130, March–April, 1972.     

Impact on the boundary of a compressible two-layer fluid

Within the framework of the acoustic approximation a solution of the plane nonstationary problem of impact on a fluid boundary is found. The fluid occupies the lower half-plane and consists of two layers with given speeds of sound and densities. The upper layer has a constant depth and is bounded above by a plate with a given normal velocity. The solution is constructed using the Fourier and Laplace integral transforms. Numerical calculations are performed for piston impact across a rigid screen and the impact of a jet with an aerated head on a rigid wall. It is shown that the presence of an interlayer with reduced speed of sound and/or density considerably changes the evolution of the hydrodynamic pressure distribution over the impacting surface: the absolute pressure maximum decreases but pressures of significant amplitude are maintained for a longer time than for a homogeneous fluid.     

二维含裂纹结构与声耦合问题研究

应用分形有限元方法结合边界元方法研究了二维含裂纹结构和声耦合问题.采用二级分形有限元方法对含裂纹的弹性结构体进行离散处理,这样可以使得自由度数大大地减少;无限大外域声场的计算使用边界元方法,可以自动满足无穷远辐射条件.数值仿真算例结果表明:结构声耦合系统的共振频率随着裂纹深度的增加而下降;裂纹附近的声场所受的影响较为明显.     

Variational principles for hydrodynamic impact problems

We first establish the rigorous field equations of the two continuous stages before and after entering water. Then correspondently, we obtain the specific variational principles, bounded theorems, and boundary integral equations of the second stage problems. The existence of solutions are proved and the scheme of solving the solutions are provided. Finally, as a numerical example, the ship's wave resistence problem is used to demonstrate the specific application of the second stage problems and its accuracy. Then we provide a rigorous and sound theoretical basis of variational finite element method and boundary element method for calculating the accurately fundamental equations.     

An efficient discontinuous Galerkin method for aeroacoustic propagation

An efficient discontinuous Galerkin formulation is applied to the solution of the linearized Euler equations and the acoustic perturbation equations for the simulation of aeroacoustic propagation in two‐dimensional and axisymmetric problems, with triangular and quadrilateral elements. To improve computational efficiency, a new strategy of variable interpolation order is proposed in addition to a quadrature‐free approach and parallel implementation. Moreover, an accurate wall boundary condition is formulated on the basis of the solution of the Riemann problem for a reflective wall. Time discretization is based on a low dissipation formulation of a fourth‐order, low storage Runge–Kutta scheme. Along the far‐field boundaries a perfectly matched layer boundary condition is used. For the far‐field computations, the integral formulation of Ffowcs Williams and Hawkings is coupled with the near‐field solver. The efficiency and accuracy of the proposed variable order formulation is assessed for realistic geometries, namely sound propagation around a high‐lift airfoil and the Munt problem. Copyright © 2011 John Wiley & Sons, Ltd.     

Scattered noise prediction using acoustic velocity formulations V1A and KV1A

Aeroacoustic scattering prediction generally relies on boundary integral methods which require evaluation of the impermeability condition on the scattering surface. The boundary condition implies zero normal velocity relative to the scattering surface. This condition has been expressed by relating the acoustic velocity to the acoustic pressure gradient, allowing indirect evaluation of the boundary condition by existent acoustic pressure gradient formulations. In the present paper, a direct evaluation of the hardwall boundary condition in scattering problems is demonstrated by time-domain analytic acoustic velocity formulae. Acoustic velocity formulations V1A and KV1A are implemented for acoustic scattering prediction, by hybrid approaches based on the FW–H equation and the Kirchhoff method These formulations can be coupled to any scattering solver, allowing time-domain prediction of the incident acoustic field when broadband noise generation is concerned. Formulation V1A offers mathematical simplicity and computational efficiency, which can be advantageous for realistic scattering applications. Implementation of formula KV1A enables acoustic scattering prediction by existing solvers based on the Kirchhoff method. The validity of the suggested methodology is assessed through the analytical test case of harmonic sound scattered by a rigid sphere. Sound propagation and scattering effects are analyzed by examination of the acoustic velocity field characteristics.     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.     

密闭腔体声-结构耦合系统的动力灵敏度分析

以密闭空腔为对象,开展了声-结构耦合系统的动力分析和灵敏度计算,为系统性态优化设计提供理论和算法基础。分别把结构和声场进行离散化,推导了声-结构耦合系统的有限元方程,求解了耦合系统的频率和声压级响应。在此基础上,以结构尺寸为设计变量,计算了耦合系统的固有频率和声压级响应的灵敏度,解决了声-结构耦合系统动力灵敏度的数值算法问题。     

Acoustic radiation induced by bubble motion in compressible fluid

Based on the theory of compressible fluid, a three-dimension boundary element method is utilized to research the motion of bubble. The far-field noise radiation during the growth and contraction is calculated by the Kirchhoff formula and the Ffowcs Williams-Hawkings （FW-H） formula with a fixed radiation surface being arranged at the near-field of bubble as a new acoustic source. The results show that the amplitude of the sound pressure induced by non-spherical bubble is lower than that of spherical bubble in the contraction phase. The retardance effect is more obvious when the observer is farther away from the bubble. In the anaphase of contraction, the observer with the maximum amplitude of sound pressure moves up with the obvious jet. Larger buoyance parameters will generate lower sound pressure amplitudes in the anaphase, while larger intensive parameters will cause higher sound pressure amplitudes in the whole procedure of bubble motion.     

一种求解三维声腔声压分布问题的无网格法

利用传统有限元法求解声压分布问题常常受到污染误差和色散误差的困扰。加权最小二乘无网格法（MWLS）是一种基于移动最小二乘（MLS）近似的无网格方法,求解声腔声压分布问题具有低色散、高精度的特点。然而传统的MLS近似有时容易产生病态矩阵,利用加权正交基函数构建改进的移动最小二乘（IMLS）近似,得到的系统方程为非病态的。本文基于改进的加权最小二乘无网格法（IMWLS）求解三维声腔内部声压分布。计算得到的声压分布和声压频响曲线都与参考值十分吻合,峰值误差和污染误差都比FEM的小,计算成本相比无单元伽辽金法显著降低。计算结果表明IMWLS相比传统的FEM,能在更高的频段内达到高精度,并且相比EFGM能大幅提高计算效率。     

一种模拟大规模高频声场的双层奇异边界法

大规模高频声场的数值模拟是一项非常有计算挑战性的课题. 为了解决传统边界型离散方法由于全局支撑的满阵限制, 不易应用于大规模高频声场模拟的计算瓶颈, 本文提出了一种用于模拟大规模高频声场的双层奇异边界法. 在该方法中, 通过引入双层结构, 细网格上的全局支撑的满阵被转化为局部支撑的大规模稀疏矩阵, 传统奇异边界法模拟大规模问题时所面临的高计算量以及过度存储需求遂得以解决. 其次, 双层奇异边界法仅通过粗网格评估远场作用, 且独立于特定的插值核函数. 相较于快速多级方法, 该方法具有更强的适应性和灵活性, 且多层结构使该方法具有一定的预调节作用, 非常适合求解具有大规模、高秩、高条件数特点的高频波矩阵. 在其后的散射球模型算例中, 双层奇异边界法配置10万个节点, 成功模拟了无量纲波数高达160的声散射问题. 在对于人头模型的声散射特性分析中, 双层奇异边界法比COMSOL软件计算速度快了约78.13\mathrm{\%}. 当配置8万个节点时, 双层奇异边界法成功模拟了频率高达25 kHz 的工况, 该频率已远远超出了人耳的听力极限.      

基于奇异分解和自适应BEM积分算法的水下航行体机械噪声预报

提出一种改进的声学边界元法（M-BEM）用于准确计算水下航行体发动机振动引起的近场辐射噪声。分别采用奇异分解技术和自适应边界元积分算法解决了Helmholtz积分方程在求解近场声压时出现的超奇异积分和奇异积分问题。采用一脉动球源的声辐射算例对方法进行验证,数值解与精确解误差小于1．5dB。结合有限元方法并考虑流固耦合作...     

管道振动引起辐射噪声场的有限元分析

对于管束振动引起的辐射声场的分析,利用解析法较为困难．运用有限元分析法对管道耦合振动噪声进行了比较系统的分析,得到一系列有益结果,为振动噪声控制提供了理论依据．     

Review on underwater sound absorption materials and mechanisms

随着我国加速实施海洋强国战略,对先进水下吸声材料的需求日益迫切.与空气吸声不同,水下的高静水压力和复杂的海洋环境对水下吸声材料提出了更为苛刻的要求.吸声问题的本质是如何将弹性能高效地转化为热能或其他形式能量.本文综述了主要以聚合物分子内摩擦机制及界面耗能机制为基础的传统水下吸声材料.传统水下吸声材料面临的主要是其在低频及高静水压力下吸声性能差的问题.这是因为：一方面受质量密度定律的限制,有限厚度的水下吸声材料无法有效吸收水中传来的低频声波;另一方面,在高静水压力下,弹性材料如高分子聚合物会变“硬”,从而大大降低了声波弹性能的转换效率.随着局域共振理论及超材料概念的提出,发展出了一系列新型水下吸声材料,为解决水下吸声材料遇到的难题提供了新思路.局域共振理论的特点是可以用小尺度结构控制长波声波的传播,从而可以解决低频吸声问题.本文重点综述了局域共振理论,以及由此发展出的声子木堆、声子玻璃等新型水下吸声材料.声子玻璃材料在局域共振理论基础上,通过引入多孔金属骨架结构提高了材料的抗压性能,从而解决了高静水压力下材料吸声性能变差的问题.本文最后对水下吸声材料未来发展方向进行了展望.