A new Fast Multipole formulation for the elastodynamic half-space Green's tensor - Archive ouverte HAL Accéder directement au contenu
Article Dans Une Revue Journal of Computational Physics Année : 2014

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

Résumé

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.
Fichier principal
Vignette du fichier
jcp13_final.pdf (434.26 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-00904660 , version 1 (14-11-2013)

Identifiants

Citer

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

Altmetric

Partager

Gmail Facebook X LinkedIn More