Commun. Comput. Phys., 11 (2012), pp. 594-609.


A Preconditioned 3-D Multi-Region Fast Multipole Solver for Seismic Wave Propagation in Complex Geometries

S. Chaillat 1*, J. F. Semblat 2, M. Bonnet 3

1 Georgia Tech, College of Computing, Atlanta, USA; and POems (CNRS-ENSTA-INRIA), Appl. Math. Dept., ENSTA, Paris, France.
2 Universite Paris Est, IFSTTAR, Paris, France.
3 POems (CNRS-ENSTA-INRIA), Appl. Math. Dept., ENSTA, Paris, France.

Received 23 December 2009; Accepted (in revised version) 3 January 2011
Available online 24 October 2011
doi:10.4208/cicp.231209.030111s

Abstract

The analysis of seismic wave propagation and amplification in complex geological structures requires efficient numerical methods. In this article, following up on recent studies devoted to the formulation, implementation and evaluation of 3-D single- and multi-region elastodynamic fast multipole boundary element methods (FM-BEMs), a simple preconditioning strategy is proposed. Its efficiency is demonstrated on both the single- and multi-region versions using benchmark examples (scattering of plane waves by canyons and basins). Finally, the preconditioned FM-BEM is applied to the scattering of plane seismic waves in an actual configuration (alpine basin of Grenoble, France), for which the high velocity contrast is seen to significantly affect the overall efficiency of the multi-region FM-BEM.

AMS subject classifications: 74J20, 74S15, 86-08

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Fast multipole method, preconditioning strategy, 3-D elastodynamics, seismic wave propagation.

*Corresponding author.
Email: stephanie.chaillat@ensta-paristech.fr (S. Chaillat), semblat@lcpc.fr (J. F. Semblat), bonnet@lms.polytechnique.fr (M. Bonnet)
 

The Global Science Journal