基于子空间投影的声辐射多频分析算法

针对边界元法在声辐射多频分析中的困难,提出了一种基于子空间投影的快速插值算法。该方法将所求解的声学量视为一个线性系统的传递函数,然后采用Krylov投影降阶方法的无矩阵形式对原系统进行模型降阶,并通过降阶系统的传递函数来实现对所求声学量的频率插值计算。文中分别采用等距频率插值和自适应插值方案,计算了三种模型算例的辐射声压和声功率。数值结果表明,本算法可以在保证一定计算精度的前提下,显著地提高声辐射多频分析的计算速度。

A projection-based algorithm for multi-frequency acoustic radiation analysis

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.