NUMERICAL ANALYSIS OF LARGE-SCALE SOUND FIELDS USING ITERATIVE METHODS PART I: APPLICATION OF KRYLOV SUBSPACE METHODS TO BOUNDARY ELEMENT ANALYSIS∗

The convergence behavior of the Krylov subspace iterative solvers towards the systems with the 3D acoustical BEM is investigated through numerical experiments. The fast multipole BEM, which is an efficient BEM based on the fast multipole method, is used for solving problems with up to about 100,000 DOF. It is verified that the convergence behavior of solvers is much affected by the formulation of the BEM (singular, hypersingular, and Burton-Miller formulation), the complexity of the shape of the problem, and the sound absorption property of the boundaries. In BiCG-like solvers, GPBiCG and BiCGStab2 have more stable convergence than others, and these solvers are useful when solving interior problems in basic singular formulation. When solving exterior problems with greatly complex shape in Burton-Miller formulation, all solvers hardly converge without preconditioning, whereas the convergence behavior is much improved with ILU-type preconditioning.

[1]  S. Marburg,et al.  Performance of iterative solvers for acoustic problems. Part I. Solvers and effect of diagonal preconditioning , 2003 .

[2]  R. Coifman,et al.  The fast multipole method for the wave equation: a pedestrian prescription , 1993, IEEE Antennas and Propagation Magazine.

[3]  Tetsuya Sakuma,et al.  A TECHNIQUE FOR PLANE-SYMMETRIC SOUND FIELD ANALYSIS IN THE FAST MULTIPOLE BOUNDARY ELEMENT METHOD , 2005 .

[4]  Martin H. Gutknecht,et al.  Variants of BICGSTAB for Matrices with Complex Spectrum , 1993, SIAM J. Sci. Comput..

[5]  V. Rokhlin.  Rapid Solution of Integral Equations of Scattering Theory , 1990 .

[6]  Yijun Liu,et al.  A unified boundary element method for the analysis of sound and shell-like structure interactions. II. Efficient solution techniques , 2000 .

[7]  Weng Cho Chew,et al.  Calculation of acoustical scattering from a cluster of scatterers , 1998 .

[8]  N. Nishimura.  Fast multipole accelerated boundary integral equation methods , 2002 .

[9]  K. H. Chen,et al.  Applications of the dual integral formulation in conjunction with fast multipole method in large-scale problems for 2D exterior acoustics , 2004 .

[10]  Toru Otsuru,et al.  Basic concept , accuracy and application of large-scale finite element sound field analysis of rooms , 2022 .

[11]  Tetsuya Sakuma,et al.  Fast multipole boundary element method for large-scale steady-state sound field analysis, Part I setup and validation , 2002 .

[12]  H. Hong,et al.  Review of Dual Boundary Element Methods With Emphasis on Hypersingular Integrals and Divergent Series , 1999 .

[13]  David E. Keyes,et al.  Preconditioned Krylov solvers for BEA , 1994 .

[14]  Michael A. Epton,et al.  Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations , 1995, SIAM J. Sci. Comput..

[15]  G. Golub,et al.  Gmres: a Generalized Minimum Residual Algorithm for Solving , 2022 .

[16]  Stefan Schneider.  Application of Fast Methods for Acoustic Scattering and Radiation Problems , 2003 .

[17]  Tetsuya Sakuma,et al.  AN EFFECTIVE SETTING OF HIERARCHICAL CELL STRUCTURE FOR THE FAST MULTIPOLE BOUNDARY ELEMENT METHOD , 2005 .

[18]  V. Rokhlin.  Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions , 1993 .

[19]  P. Sonneveld.  CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear systems , 1989 .

[20]  Steffen Marburg,et al.  Performance of iterative solvers for acoustic problems. Part II. Acceleration by ILU-type preconditioner , 2003 .

[21]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[22]  Shao-Liang Zhang,et al.  GPBi-CG: Generalized Product-type Methods Based on Bi-CG for Solving Nonsymmetric Linear Systems , 1997, SIAM J. Sci. Comput..

[23]  Y. Saad,et al.  Iterative solution of linear systems in the 20th century , 2000 .

[24]  S. Amini,et al.  PRECONDITIONED KRYLOV SUBSPACE METHODS FOR BOUNDARY ELEMENT SOLUTION OF THE HELMHOLTZ EQUATION , 1998 .

[25]  Jeng-Tzong Chen,et al.  Dual integral formulation for determining the acoustic modes of a two-dimensional cavity with a degenerate boundary , 1998 .

[26]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[27]  L. Greengard,et al.  Accelerating fast multipole methods for the Helmholtz equation at low frequencies , 1998 .

[28]  G. F. Miller,et al.  The application of integral equation methods to the numerical solution of some exterior boundary-value problems , 1971, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[29]  S. Amini,et al.  Multi-level fast multipole solution of the scattering problem , 2003 .

[30]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[31]  T. Terai.  On calculation of sound fields around three dimensional objects by integral equation methods , 1980 .

[32]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

[33]  Zydrunas Gimbutas,et al.  A wideband fast multipole method for the Helmholtz equation in three dimensions , 2006, J. Comput. Phys..

[34]  V. Rokhlin.  Rapid solution of integral equations of classical potential theory , 1985 .