Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
8256798 | Wave Motion | 2018 | 40 Pages |
Abstract
We investigate efficient methods to simulate the multiple scattering of obstacles in homogeneous media. With a large number of small obstacles on a large domain, optimized pieces of software based on spatial discretization such as Finite Element Method (FEM) or Finite Difference lose their robustness. As an alternative, we work with an integral equation method, which uses single-layer potentials and truncation of Fourier series to describe the approximate scattered field. In the theoretical part of the paper, we describe in detail the linear systems generated by the method for impenetrable obstacles, accompanied by a well-posedness study. For the numerical performance study, we limit ourselves to the case of circular obstacles. We first compare and validate our codes with the highly optimized FEM-based software Montjoie. Secondly, we investigate the efficiency of different solver types (direct and iterative of type GMRES) in solving the dense linear system generated by the method. We observe the robustness of direct solvers over iterative ones for closely-spaced obstacles, and that of GMRES with Lower-Upper Symmetric Gauss-Seidel and Symmetric Gauss-Seidel preconditioners for far-apart obstacles.
Keywords
Related Topics
Physical Sciences and Engineering
Earth and Planetary Sciences
Geology
Authors
Hélène Barucq, Juliette Chabassier, Ha Pham, Sébastien Tordeux,