Title
A new Fast Multipole formulation for the elastodynamic half-space Green's tensor
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=6x10^5 degrees of freedom. In particular, a complexity significantly lower than that of the non-multipole version is shown to be achieved.
Year
DOI
Venue
2014
10.1016/j.jcp.2013.11.010
Journal of Computational Physics
Keywords
Field
DocType
elastodynamic half-space Green,half-space Green,free surface,frequency-domain elastodynamic,new Fast Multipole formulation,efficient quadrature,full-space Green,single-layer elastodynamic,non-multipole version,Gaussian quadrature rule,classical Gaussian quadratures
Discretization,Mathematical optimization,Tensor,Mathematical analysis,Fourier transform,Gaussian,Fast multipole method,Fundamental solution,Tensor contraction,Gaussian quadrature,Mathematics
Journal
Volume
ISSN
Citations 
258,
0021-9991
2
PageRank 
References 
Authors
0.54
4
2
Name
Order
Citations
PageRank
Stéphanie Chaillat193.05
Marc Bonnet232.65