小阻尼封闭空间声压级灵敏度分析与优化设计

本文研究了小阻尼界面封闭空间低频声学有限元分析、灵敏度分析和优化设计问题。分别用模态法和直接法计算了封闭空间内声压级响应，并推导了声压级响应对声空间边界形状控制参数的灵敏度分析公式，在此基础上建立了小阻尼空间声学问题的优化模型，同时给出了优化求解算法，并在JIFEX软件中进行了程序实现。本文提出的灵敏度分析和优化设计方法可以使声场的边界布局更为合理，从而达到改进小阻尼界面封闭空间声学性能的目的。数值算例验证了本文提出的灵敏度分析和优化算法的有效性。     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

基于声粒子分布积分的无网格声场计算方法

小尺度封闭空间内部声场的数值计算是声学设计、噪声控制等领域的关键技术。由于波动声学及几何声学方法计算频率上的限制,中频段声场计算问题一直是个难点。本文以声学无网格法为基础,提出了一种基于声粒子分布积分的无网格声场数值计算方法。文中利用声线跟踪理论计算声场中的声粒子分布,并以某个时间点上的声粒子作为蒙特卡罗法中的积分点,将其应用于无网格法中,从而获得声场中的节点声压。利用该方法对一个矩形封闭空间的中低频声场进行了计算,并与模态叠加法、商用声场计算软件、经典无网格法的结果进行了对比,证明基于声粒子分布积分的无网格声场数值计算方法在中低频段相较于传统基于网格的方法具有更高的精度。     

小尺度封闭空间声场的无网格数值算法

无网格法是一种新兴的数值计算方法,具有不需要网格支持的特点。本文将该方法引入室内声学建模,推导了无网格声场数值计算模型,并将其应用于典型小尺度封闭空间内部声场的数值分析。针对声传递函数,将本方法与理论解和SYSNOISE计算结果进行了比较,并将计算的混响时间与实验测量结果作了对比,表明本方法具有良好的精度。     

适用于任意形状封闭空间的声场预测相干模型

针对现有几何声学的方法对封闭空间内声场进行预测时在中低频段出现较大误差的问题,该文提出一种近似圆锥声束追踪法和相干反射场理论相结合的声场预测新模型。在近似圆锥声束追踪法基础上,考虑声束轴线在边界多次反射时声压和相位的改变,最后计算不同声波之间的干涉效应,建立一种适用于任意形状封闭空间的声场预测相干模型。利用该模型对某一矩形封闭空间进行声场预测,通过对边界元法、Raynoise软件相干和非相干算法的预测结果和本模型的数值模拟结果对比。结果表明,文中提出的方法和边界元法的计算结果在中低频段非常吻合,两者的计算结果平均绝对误差为1.1 d B。本模型在中低频率下与同样考虑了相位的Raynoise相干算法相比有更好的准确性,在较高频率上,本模型计算结果与Raynoise相干算法计算结果非常吻合。     

有限长导管声场预报的一种边界元方法

为了研究介质流动情况下有限长导管中点声源产生的声场,提出了一种以源势密度作为未知数的间接边界元方法.首先将整个求解空间划分为管内和管外两个相互独立的封闭域,然后利用压力和速度在开口处的连续条件和壁面边界条件将内外域的声场方程联立为一个矩阵方程,求解出源势密度后由源势密度计算出任意位置的声场.与Myers的基准数据比较表...     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法，该方程求解声辐射和散射问题具有相同的影响矩阵，能够同时求解薄壳体气动和振动噪声的辐射和散射现象，以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上，可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

边界元法计算已知振速封闭面的声辐射

采用边界元法求解封闭面在无限域声媒质中的辐射声场具有内存小、计算精度高、速度快等优点。但需处理被积函数在边界面上的奇异积分及表面Helmholtz方程在特征频率下无唯一解的问题。本文提出把内部Helmholtz方程与它关于内点坐标取导后的式子构成补充方程式,经与表面Helmholtz方程相结合,可求解任意频率下的声辐射问题;而在奇点附近区域,则提出用极坐标变换消除积分的奇异性。文中以轴对称形封闭面为例,计算了具有已知表面振速分布下的辐射声场。     

阻尼受击板的时域建模及声音合成

为满足声源辨识中对合成冲击声的迫切需求,建立了球-板撞击的时域模型,提出一种时域快速求解方法,并进行了实验验证.首先给出一种将时域有限差分法(FDTD)和模态展开法(MEM)相结合的时域混合方法,求解板的振动方程,并解决了混合方法中MEM的模态截断和初值问题,及两种方法中阻尼的一致性问题;随后,给出了简支矩形板的冲击声计算结果,通过与FDTD方法的运算量进行对比,验证了混合方法的高效性;最后,针对自由边界下的L形板进行了实验验证.结果表明,与传统FDTD方法相比,时域混合方法在保证合成冲击声精度的前提下可将计算效率提高100至1200倍。     

多约束纳米结构的声子热导率模型研究

纳米技术的快速发展使得对微纳尺度导热机理的深入研究变得至关重要. 理论和实验都表明, 在纳米尺度下声子热导率将表现出尺寸效应. 基于声子玻尔兹曼方程和修正声子平均自由程的方法得到了多约束纳米结构的声子热导率模型, 可以描述多个几何约束共同作用下热导率的尺寸效应. 不同几何约束对声子输运的限制作用可以分开计算, 总体影响则通过马西森定则进行耦合. 对于热流方向的约束, 采用扩散近似的方法求解声子玻尔兹曼方程; 对于侧面边界约束, 采用修正平均自由程的方法计算边界散射对热导率的影响. 得到的模型能够预测纳米薄膜(法向和面向)及有限长度方形纳米线的热导率随相应特征尺寸的变化. 与蒙特卡罗模拟及硅纳米结构热导率实验值的对比验证了模型的正确性.     

基于联合波叠加法的浅海信道下圆柱壳声辐射研究

针对浅海信道下弹性结构声辐射预报尚无高效可靠的研究方法,提出了一种浅海信道下弹性结构声辐射快速预报的联合波叠加法.该方法结合了浅海信道传输函数、多物理场耦合数值计算法和波叠加法理论,运用该方法可对浅海信道下弹性结构辐射声场进行快速预报.经数值法和解析解法验证后,从信道下辐射源、环境影响和辐射声场测量的角度研究分析了浅海信道下弹性圆柱壳的声辐射特性,阐释了进行浅海信道下结构声辐射研究的必要性.研究结果表明,仅在低频浅海信道下弹性结构可近似等效为点源,信道上下边界对声场产生显著的耦合影响,高频段的空间声场指向性分布尤为明显,垂直线列阵进行信道下结构辐射声功率测量时,测量结果受到信道环境边界和潜深的影响较大.     

求解水中加肋板声辐射特性的虚拟声源法*

发展一种利用虚拟声源离散声场的方法求解加肋板在水中的声振耦合问题。由波叠加原理和单元体积速度匹配的原则,根据离散的结构单元满足的动力方程和结构与介质的交界相容性条件,确定虚拟声源强度,计算结构的声辐射功率。本文以简支矩形加肋板为例,在不获得结构表面振速和声压的情况下,计算了结构在水中的声辐射功率,并与解析方法计算的结果进行了比较,表明了该方法具有较好的计算精度。     

Sound propagation above an inhomogeneous plane: Boundary integral equation methods

A boundary integral equation method is used to compute the sound pressure emitted by a harmonic source above an inhomogeneous plane. First, the theoretical aspects of the problem (behaviour of the pressure around the discontinuities,…) are studied. Then, a comparison between theoretical levels and experimental levels obtained in an anechoic room is presented. It shows that the boundary integral equation (BIE) method is quite convenient for solving this kind of problem. Two interesting results are pointed out: (i) if only a prediction of maximum sound levels is needed, the attenuation is the same for a cylindrical source, a spherical source and N spherical sources, and so it is possible to transform some three-dimensional problems into two-dimensional ones; (ii) a numerical method of computation of the sound field above an inhomogeneous plane does not provide a correct prediction if each part of the plane is not accurately described by the boundary condition chosen.     

小阻尼界面房间声传输函数和声脉冲响应函数的有限元素法

有限元法可用于以声波动方程为基础通过数值计算求解室内声场,适用于分析界面阻抗非均匀分布和复杂形状房间内声场的低频特性。本文首先介绍了小阻尼界面条件下室内声场简正方式、声衰变系数、混响时间的ＦＥＭ计算方法。在此基础上导出了房间内两点之间声传输函数和声脉冲响应函数的ＦＥＭ计算模型,并以矩形房间为例详细讨论了有关细节。本文所讨论的计算模型可以映房间内不同的声源点,接收点位置上的声压频谱特性和脉冲响应的时     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

A computational method of evaluating noncompact sound based on vortex sound theory

A numerical investigation is made of the production of sound by turbulence interacting with a noncompact body. The problem is formulated in the frequency domain by extending the theory of vortex sound proposed by Howe. The anomalous "numerical" generation of sound by the sudden termination of Lighthill's stress tensor at the outer boundary of a finite computational domain is avoided by identification of "scattered" sound sources that generate sound principally by interaction with the solid surface. It is argued that the boundary element method is the most efficient means of computing the aeroacoustic Green's function for the problem, because it requires a minimum of CPU time, is not prone to numerical errors such as dispersion and dissipation during propagation, and the radiation condition is easily applied at the outer boundary. The method is applied to the problem of sound generation by high Reynolds number flow past a circular cylinder. The "scattered" sources are shown to be confined to the vicinity of the cylinder surface. At low frequencies the radiation has a dipole-like directivity in agreement with the compact approximation. However, the directivity is quite different at high frequencies, where our noncompact method predicts a more complicated "leaf-like" radiation pattern.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

快速多极边界元法用于扩散光学断层成像研究

在求解扩散光学断层成像中的正向问题时, 目前普遍采用有限元法, 但是随着实际模型规模的增大, 有限元法的计算量问题日益显著, 而边界元法则由于可以降低计算维度使计算量减少而备受关注. 本文以均匀的高散射介质为模型, 研究了将快速多极边界元法用于扩散光学断层成像的正向问题. 快速多极边界元法利用核函数的多极展开, 将常规边界元法中系数矩阵和迭代矢量的乘积项等价为相应四叉树结构的一次递归, 再结合广义最小残量法进行迭代求解. 将计算结果和蒙特卡罗法的模拟结果进行了比较, 表明利用快速多极边界元法的模拟结果和蒙特卡罗法的结果有很好的一致性. 研究结果验证了快速多极边界元法可以用于扩散光学断层成像, 为其大规模和实时成像带来可观的前景. 关键词： 扩散光学断层成像 边界元法 快速多极边界元法     

On boundary conditions for the diffusion equation in room-acoustic prediction: Theory, simulations, and experiments

This paper proposes a modified boundary condition to improve the room-acoustic prediction accuracy of a diffusion equation model. Previous boundary conditions for the diffusion equation model have certain limitations which restrict its application to a certain number of room types. The boundary condition employing the Sabine absorption coefficient [V. Valeau et al., J. Acoust. Soc. Am. 119, 1504-1513 (2006)] cannot predict the sound field well when the absorption coefficient is high, while the boundary condition employing the Eyring absorption coefficient [Y. Jing and N. Xiang, J. Acoust. Soc. Am. 121, 3284-3287 (2007); A. Billon et al., Appl. Acoust. 69, (2008)] has a singularity whenever any surface material has an absorption coefficient of 1.0. The modified boundary condition is derived based on an analogy between sound propagation and light propagation. Simulated and experimental data are compared to verify the modified boundary condition in terms of room-acoustic parameter prediction. The results of this comparison suggest that the modified boundary condition is valid for a range of absorption coefficient values and successfully eliminates the singularity problem.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题