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 - Electrical Engineering, Mathematics and Computer Science)

R. Astudillo (TU Delft - Electrical Engineering, Mathematics and Computer Science)

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

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

M.B. van Gijzen (TU Delft - Electrical Engineering, Mathematics and Computer Science)

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

Research Group
Numerical Analysis
DOI related publication
https://doi.org/10.1007/s10596-017-9667-7 Final published version
More Info
expand_more
Publication Year
2018
Language
English
Research Group
Numerical Analysis
Journal title
Computational Geosciences: modeling, simulation and data analysis
Issue number
1
Volume number
22
Pages (from-to)
43-61
Downloads counter
353
Collections
Institutional Repository
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.