An MSSS-preconditioned matrix equation approach for the time-harmonic elastic wave equation at multiple frequencies

Journal Article (2018)
Author(s)

M. Baumann (TU Delft - Numerical Analysis)

R.A. Astudillo Rengifo (TU Delft - Numerical Analysis)

Y Qiu (Max Planck Institute for Dynamics of Complex Technical Systems)

Y.M.E. Ang (Nanyang Technological University)

Martin B. Van Gijzen (TU Delft - Numerical Analysis)

R.E. Plessix (Shell Global Solutions International B.V.)

Research Group
Numerical Analysis
Copyright
© 2018 M.M. Baumann, R.A. Astudillo Rengifo, Y. Qiu, Y.M.E. Ang, M.B. van Gijzen, R.E. Plessix
DOI related publication
https://doi.org/10.1007/s10596-017-9667-7
More Info
expand_more
Publication Year
2018
Language
English
Copyright
© 2018 M.M. Baumann, R.A. Astudillo Rengifo, Y. Qiu, Y.M.E. Ang, M.B. van Gijzen, R.E. Plessix
Research Group
Numerical Analysis
Issue number
1
Volume number
22
Pages (from-to)
43-61
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

In this work, we present a new numerical framework for the efficient solution of the time-harmonic elastic wave equation at multiple frequencies. We show that multiple frequencies (and multiple right-hand sides) can be incorporated when the discretized problem is written as a matrix equation. This matrix equation can be solved efficiently using the preconditioned IDR(s) method. We present an efficient and robust way to apply a single preconditioner using MSSS matrix computations. For 3D problems, we present a memory-efficient implementation that exploits the solution of a sequence of 2D problems. Realistic examples in two and three spatial dimensions demonstrate the performance of the new algorithm.