Sweeping preconditioners for elastic wave propagation with spectral element methods
1 Sandia National Laboratories, Org.
1442: Numerical Analysis and Applications, Livermore, CA
2 Georgia Institute of Technology, School of Computational Science and Engineering, Atlanta, GA 30332, USA.
3 University of Texas at Austin, Department of Mathematics, Austin, TX 78712, USA.
4 Stanford University, Department of Mathematics, Stanford, CA 94305, USA.
We present a parallel preconditioning method for the iterative solution of the time-harmonic elastic wave equation which makes use of higher-order spectral elements to reduce pollution error. In particular, the method leverages perfectly matched layer boundary conditions to efficiently approximate the Schur complement matrices of a block LDLT factorization. Both sequential and parallel versions of the algorithm are discussed and results for large-scale problems from exploration geophysics are presented.
Mathematics Subject Classification: 65F08 / 65N22 / 65N80
Key words: Elastic wave / seismic wave / time-harmonic / frequency domain / spectral elements / parallel preconditioner / iterative solver / sparse-direct / perfectly matched layers / full waveform inversion
© EDP Sciences, SMAI, 2014