Skip to main content Accessibility help
×
Home
Hostname: page-component-7ccbd9845f-s2vjv Total loading time: 0.465 Render date: 2023-01-27T09:32:17.588Z Has data issue: true Feature Flags: { "useRatesEcommerce": false } hasContentIssue true

Sweeping preconditioners for elastic wave propagation with spectral element methods

Published online by Cambridge University Press:  20 February 2014

Paul Tsuji
Affiliation:
Sandia National Laboratories, Org. 1442: Numerical Analysis and Applications, Livermore, CA 94550, USA.
Jack Poulson
Affiliation:
Georgia Institute of Technology, School of Computational Science and Engineering, Atlanta, GA 30332, USA.
Björn Engquist
Affiliation:
University of Texas at Austin, Department of Mathematics, Austin, TX 78712, USA.
Lexing Ying
Affiliation:
Stanford University, Department of Mathematics, Stanford, CA 94305, USA.. lexing@math.stanford.edu
Get access

Abstract

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.

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

F. Aminzadeh, J. Brac and T. Kunz, 3-D Salt and Overthrust models. In SEG/EAGE 3-D Modeling Series 1. Tulsa, OK (1997).
Berenger, J.P., A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. 114 (1994) 185200. Google Scholar
Bramble, J. and Pasciak, J., A note on the existence and uniqueness of solutions of frequency domain elastic wave problems: a priori estimates in H 1. J. Math Anal. Appl. 345 (2008) 396404. Google Scholar
Chew, W.C. and Weedon, W.H., A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates. Microwave Optical Tech. Lett. 7 (1994) 599604. Google Scholar
J. Choi, J.J. Dongarra, R. Pozo and D.W. Walker, ScaLAPACK: A scalable linear algebra library for distributed memory concurrent computers, in Proc. of the Fourth Symposium on the Frontiers of Massively Parallel Computation, IEEE Comput. Soc. Press(1992) 120–127.
Engquist, B. and Ying, L., Sweeping preconditioner for the Helmholtz equation: hierarchical matrix representation. Commun. Pure Appl. Math. 64 (2011) 697735. Google Scholar
Engquist, B. and Ying, L., Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers. Multiscale Model. Simul. 9 (2011) 686710. Google Scholar
Erlangga, Y.A., Vuik, C. and Oosterlee, C.W., On a class of preconditioners for solving the Helmholtz equation. Appl. Numer. Math. 50 (2004) 409425. Google Scholar
O.G. Ernst and M.J. Gander, Why it is difficult to solve Helmholtz problems with classical iterative methods, in vol. 83 of Numerical Analysis of Multiscale Problems. Edited by I. Graham, T. Hou, O. Lakkis and R. Scheichl. Springer-Verlag (2011) 325–361.
Grasedyck, L. and Hackbusch, W., Construction and arithmetics of ℋ-matrices. Computing 70 (2003) 295334. Google Scholar
Gupta, A., Karypis, G. and Kumar, V., A highly scalable parallel algorithm for sparse matrix factorization. IEEE Transactions on Parallel and Distributed Systems 8 (1997) 502520. Google Scholar
A. Gupta, S. Koric and T. George, Sparse matrix factorization on massively parallel computers, in Proc. of the Conference on High Performance Computing, Networking, Storage and Analysis. Portland, OR (2009).
Harari, I. and Albocher, U., Studies of FE/PML for exterior problems of time-harmonic elastic waves. Comput. Methods Appl. Mech. Eng. 195 (2006) 38543879. Google Scholar
T. Hughes, The finite element method: linear static and dynamic finite element analysis. Prentice-Hall, Inc. (1987).
G. Karniadakis, Spectral/hp element methods for CFD. Oxford University Press (1999).
Komatitsch, D. and Tromp, J., Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophys. J. Int. 139 (1999) 806822. Google Scholar
Liu, J.. The multifrontal method for sparse matrix solution: theory and practice. SIAM Review 34 (1992) 82109. Google Scholar
Patera, A., A spectral element method for fluid dynamics: laminar flow in a channel expansion. J. Comput. Phys. 54 (1984) 468488. Google Scholar
J. Poulson, B. Engquist, S. Li and L. Ying, A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations e-prints ArXiv (2012).
Poulson, J., Marker, B., van de Geijn, R.A., Hammond, J.R. and Romero, N.A., Elemental: A new framework for distributed memory dense matrix computations. ACM Trans. Math. Software 39.
Raghavan, P., Efficient parallel sparse triangular solution with selective inversion. Parallel Proc. Lett. 8 (1998) 2940. Google Scholar
Saad, Y. and Schultz, M.H., A generalized minimum residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 7 (1986) 856869. Google Scholar
Schreiber, R., A new implementation of sparse Gaussian elimination. ACM Trans. Math. Software 8 (1982) 256276. Google Scholar
Tsuji, P., Engquist, B. and Ying, L., A sweeping preconditioner for time-harmonic Maxwell’s equations with finite elements. J. Comput. Phys. 231 (2012) 37703783. Google Scholar
Tsuji, P. and Ying, L., A sweeping preconditioner for Yee’s finite difference approximation of time-harmonic Maxwell’s equations. Frontiers of Mathematics in China 7 (2012) 347363. Google Scholar

Save article to Kindle

To save this article to your Kindle, first ensure coreplatform@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about saving to your Kindle.

Note you can select to save to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

Sweeping preconditioners for elastic wave propagation with spectral element methods
Available formats
×

Save article to Dropbox

To save this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Dropbox account. Find out more about saving content to Dropbox.

Sweeping preconditioners for elastic wave propagation with spectral element methods
Available formats
×

Save article to Google Drive

To save this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Google Drive account. Find out more about saving content to Google Drive.

Sweeping preconditioners for elastic wave propagation with spectral element methods
Available formats
×
×

Reply to: Submit a response

Please enter your response.

Your details

Please enter a valid email address.

Conflicting interests

Do you have any conflicting interests? *