A CG-FFT approach to the solution of a stress-velocity formulation of three-dimensional elastic scattering problems

Authors:
Jiaqi Yang;Aria Abubakar;Peter M. van den Berg;Tarek M. Habashy;Fernando Reitich
Affiliations:
University of Minnesota, School of Mathematics, 206 Church Street SE, 127 Vincent Hall, Minneapolis, MN 55455, USA;Schlumberger-Doll Research, Cambridge, USA;Delft University of Technology, Delft, The Netherlands;Schlumberger-Doll Research, Cambridge, USA;University of Minnesota, School of Mathematics, 206 Church Street SE, 127 Vincent Hall, Minneapolis, MN 55455, USA
Venue:
Journal of Computational Physics
Year:
2008

Citing 6
Cited 0

BI-CGSTAB: a fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems

SIAM Journal on Scientific and Statistical Computing
Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing

Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing
Iterative methods for solving linear systems

Iterative methods for solving linear systems
A weak form of the conjugate gradient FFT method for two-dimensional elastodynamics

Journal of Computational Physics
Spatial Parallelism of a 3D Finite Difference Velocity-Stress Elastic Wave Propagation Code

SIAM Journal on Scientific Computing
Iterative forward and inverse algorithms based on domain integral equations for three-dimensional electric and magnetic objects

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

In this paper we introduce a Conjugate Gradient Fast Fourier Transform (CG-FFT) scheme for the numerical solution of the integral equation formulating three-dimensional elastic scattering problems. The formulation is in terms of the stress tensor and particle velocities as the unknown field variables. In contrast with the formulation based on particle displacements, this approach leads to integral representations that do not involve derivatives of the unknown fields, thus resulting in simplified and more stable numerics. The numerical procedure is based on suitable quadrature formulas that provide (second order) accurate approximations while retaining the convolution nature of the relevant integrals that make them amenable to efficient evaluation via FFTs. The scheme is further improved through the introduction of (approximation-based) pre-conditioners that are shown to accelerate the convergence of the CG iterations. Numerical results are presented that demonstrate the accuracy and efficiency of the proposed methodology.