A fast Fourier transform based direct solver for the Helmholtz problem

This paper is devoted to the efficient numerical solution of the Helmholtz equation in a two- or three-dimensional rectangular domain with an absorbing boundary condition (ABC). The Helmholtz problem is discretized by standard bilinear and trilinear finite elements on an orthogonal mesh yielding a separable system of linear equations. The main key to high performance is to employ the Fast Fourier transform (FFT) within a fast direct solver to solve the large separable systems. The computational complexity of the proposed FFT based direct solver is O(N log N) operations. Numerical results for both two- and three-dimensional problems are presented confirming the efficiency of the method discussed.

[1]  Tuomo Rossi,et al.  Fast direct solution of the Helmholtz equation with a perfectly matched layer or an absorbing boundary condition , 2003 .

[2]  Tuomo Rossi,et al.  A Nonstandard Cyclic Reduction Method, Its Variants and Stability , 1999, SIAM J. Matrix Anal. Appl..

[3]  Patrick Joly,et al.  Second-order absorbing boundary conditions for the wave equation: a solution for the corner problem , 1990 .

[4]  Kazufumi Ito,et al.  A Parallel Domain Decomposition Method for the Helmholtz Equation in Layered Media , 2019, SIAM J. Sci. Comput..

[5]  Kazufumi Ito,et al.  A domain decomposition solver for acoustic scattering by elastic objects in layered media , 2008, J. Comput. Phys..

[6]  Yuri A. Kuznetsov,et al.  Fictitious Domain Methods For The Numerical Solution Of Three-Dimensional Acoustic Scattering Proble , 1999 .

[7]  Charbel Farhat,et al.  Iterative solution of large‐scale acoustic scattering problems with multiple right hand‐sides by a domain decomposition method with Lagrange multipliers , 2001 .

[8]  P. Swarztrauber.  THE METHODS OF CYCLIC REDUCTION, FOURIER ANALYSIS AND THE FACR ALGORITHM FOR THE DISCRETE SOLUTION OF POISSON'S EQUATION ON A RECTANGLE* , 1977 .

[9]  B. L. Buzbee,et al.  The direct solution of the discrete Poisson equation on irregular regions , 1970 .

[10]  Alexandra Banegas,et al.  Fast Poisson solvers for problems with sparsity , 1978 .

[11]  John R. Rice,et al.  Direct solution of partial difference equations by tensor product methods , 1964 .

[12]  Yuri A. Kuznetsov,et al.  Fictitious Domain Methods for the Numerical Solution of Two-Dimensional Scattering Problems , 1998 .

[13]  Oliver G. Ernst,et al.  A Domain Decomposition Approach to Solving the Helmholtz Equation with a Radiation Boundary Conditio , 1992 .

[14]  Timothy A. Davis,et al.  Direct methods for sparse linear systems , 2006, Fundamentals of algorithms.

[15]  Yu. A. KUZNETSOV,et al.  On partial solution of systems of linear algebraic equations , 1989 .

[16]  Tuomo Rossi,et al.  A Parallel Fictitious Domain Method for the Three-Dimensional Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[17]  Oliver G. Ernst,et al.  A finite-element capacitance matrix method for exterior Helmholtz problems , 1996 .

[18]  Tuomo Rossi,et al.  A Parallel Fast Direct Solver for Block Tridiagonal Systems with Separable Matrices of Arbitrary Dimension , 1999, SIAM J. Sci. Comput..

[19]  M. Guddati,et al.  Modified integration rules for reducing dispersion error in finite element methods , 2004 .