摘要
针对边界元法在声辐射多频分析中的困难,提出了一种基于子空间投影的快速插值算法。该方法将所求解的声学量视为一个线性系统的传递函数,然后采用Krylov投影降阶方法的无矩阵形式对原系统进行模型降阶,并通过降阶系统的传递函数来实现对所求声学量的频率插值计算。文中分别采用等距频率插值和自适应插值方案,计算了三种模型算例的辐射声压和声功率。数值结果表明,本算法可以在保证一定计算精度的前提下,显著地提高声辐射多频分析的计算速度。
An interpolation algorithm based on the Krylov subspace projection is proposed to resolve the computational difficulty in multi-frequency acoustic radiation analysis by the boundary element method. It recasts acoustic quantities to be solved into transfer functions of a linear system. Then a matrix-free model reduction formulation is employed to reduce the original system in the frequency band of interest. Interpolation of the acoustic quantities is realized by the transfer functions of the reduced order system. Equally-spacing and adaptive interpolation schemes are presented to calculate the acoustic radiation pressure and power for three model examples. Numerical results demonstrate that the proposed algorithm can significantly speed up the multi-frequency analysis of acoustic radiation problems with a guaranteed computational accuracy.
出处
《声学学报》
EI
CSCD
北大核心
2010年第3期321-328,共8页
Acta Acustica
基金
国家自然科学基金资助项目(10602012)
关键词
子空间投影
插值算法
声辐射
频分
线性系统
模型降阶
传递函数
插值计算
Acoustic emissions
Acoustic radiators
Acoustic wave propagation
Acoustic wave transmission
Algorithms
Atmospheric pressure
Frequency bands
Interpolation
Linear systems
Transfer functions
