# A new Fast Multipole formulation for the elastodynamic half-space Green's tensor

1 POEMS - Propagation des Ondes : Étude Mathématique et Simulation
Inria Saclay - Ile de France, UMA - Unité de Mathématiques Appliquées, CNRS - Centre National de la Recherche Scientifique : UMR7231
Abstract : In this article, a version of the frequency-domain elastodynamic Fast Multipole-Boundary Element Method (FM-BEM) for semi-infinite media, based on the half-space Green's tensor (and hence avoiding any discretization of the planar traction-free surface), is presented. The half-space Green's tensor is often used (in non-multipole form until now) for computing elastic wave propagation in the context of soil-structure interaction, with applications to seismology or civil engineering. However, unlike the full-space Green's tensor, the elastodynamic half-space Green's tensor cannot be expressed using derivatives of the Helmholtz fundamental solution. As a result, multipole expansions of that tensor cannot be obtained directly from known expansions, and are instead derived here by means of a partial Fourier transform with respect to the spatial coordinates parallel to the free surface. The obtained formulation critically requires an efficient quadrature for the Fourier integral, whose integrand is both singular and oscillatory. Under these conditions, classical Gaussian quadratures would perform poorly, fail or require a large number of points. Instead, a version custom-tailored for the present needs of a methodology proposed by Rokhlin and coauthors, which generates generalized Gaussian quadrature rules for specific types of integrals, has been implemented. The accuracy and efficiency of the proposed formulation is demonstrated through numerical experiments on single-layer elastodynamic potentials involving up to about $N=6 10^5$ degrees of freedom. In particular, a complexity significantly lower than that of the non-multipole version is shown to be achieved.
Document type :
Journal articles

Cited literature [36 references]

https://hal.archives-ouvertes.fr/hal-00904660
Contributor : Marc Bonnet Connect in order to contact the contributor
Submitted on : Thursday, November 14, 2013 - 6:32:16 PM
Last modification on : Wednesday, May 11, 2022 - 12:06:04 PM
Long-term archiving on: : Saturday, February 15, 2014 - 4:37:37 AM

### File

jcp13_final.pdf
Files produced by the author(s)

### Citation

Stéphanie Chaillat, Marc Bonnet. A new Fast Multipole formulation for the elastodynamic half-space Green's tensor. Journal of Computational Physics, Elsevier, 2014, 258, pp.787-808. ⟨10.1016/j.jcp.2013.11.010⟩. ⟨hal-00904660⟩

Record views