A Weak Form of the Conjugate Gradient FFT Method for Two-Dimensional Elastodynamics

Author: Pelekanos G.   Kleinman R.E.   van den Berg P.M.  

Publisher: Academic Press

ISSN: 0021-9991

Source: Journal of Computational Physics, Vol.160, Iss.2, 2000-05, pp. : 597-611

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

The problem of two-dimensional scattering of elastic waves by an elastic inclusion can be formulated in terms of a domain integral equation, in which the grad-div operator acts on a vector potential. The vector potential is the spatial convolution of a Green's function with the product of the density and the displacement over the domain of interest. A weak form of the integral equation for the unknown displacement is obtained by testing it with rooftop functions. This method shows excellent numerical performance.