We propose a distributed parallel algorithm for the solution of block circulant linear systems arising from acoustic radiation problems with rotationally symmetric boundary surfaces. When large structural acoustics problems are solved using a coupling finite element/boundary element formulation, the most time consuming part of the analysis is the solution of the linear system of equations for the boundary element computation. In general, the problem is solved frequency by frequency, and the coefficient matrix for the boundary element analysis is fully populated and exhibits no exploitable structure. This typically limits the number of acoustic degrees of freedom to 10–20 thousand. Because acoustic boundary element calculations require approximately six elements per wavelength to produce accurate solutions, the formation is limited to relatively low frequencies. However, when the outer surface of the structure is rotationally symmetric, the system of linear equations becomes block circulant. Building upon a known inversion formula for block circulant matrices, a parallel algorithm for the efficient solution of linear systems arising from acoustic radiation problems with rotationally symmetric boundary surfaces is developed. We show through a runtime, speedup, and efficiency analysis that the reductions in computation time are significant for an increasing number of processors.
On performance of SOR method for solving nonsymmetric linear systems