Sweeping preconditioners for stratified media in the presence of reflections

DOI10.1007/S42985-020-00019-XMaRDI QIDQ2218207zbMATH OpenOpenAlexFDO

Authors Janosch Preuß, Thorsten Hohage, Christoph Lehrenfeld

Publication date 15 January 2021

Published in SN Partial Differential Equations and Applications (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/2002.12016

zbMATH Keywords

computational seismologyDirichlet-to-Neumann operatordomain decomposition methodsHelmholtz equationperfectly matched layerspreconditioning

Mathematics Subject Classification ID

Wave equation (35L05) Preconditioners for iterative methods (65F08) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Computational methods for problems pertaining to geophysics (86-08) Seismology (including tsunami modeling), earthquakes (86A15)

Abstract: In this paper we consider sweeping preconditioners for time harmonic wave propagation in stratified media, especially in the presence of reflections. In the most famous class of sweeping preconditioners Dirichlet-to-Neumann operators for half-space problems are approximated through absorbing boundary conditions. In the presence of reflections absorbing boundary conditions are not accurate resulting in an unsatisfactory performance of these sweeping preconditioners. We explore the potential of using more accurate Dirichlet-to-Neumann operators within the sweep. To this end, we make use of the separability of the equation for the background model. While this improves the accuracy of the Dirichlet-to-Neumann operator, we find both from numerical tests and analytical arguments that it is very sensitive to perturbations in the presence of reflections. This implies that even if accurate approximations to Dirichlet-to-Neumann operators can be devised for a stratified medium, sweeping preconditioners are limited to very small perturbations.

Recommendations

• Sweeping preconditioners for the iterative solution of quasiperiodic Helmholtz transmission problems in layered media
• Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers
• Sweeping preconditioner for the Helmholtz equation: hierarchical matrix representation
• Preconditioning Methods for Thin Scattering Structures Based on Asymptotic Results
• Sweeping preconditioners for elastic wave propagation with spectral element methods
• Preconditioning techniques for the iterative solution of scattering problems
• Preconditioned Krylov subspace methods for solving radiative transfer problems with scattering and reflection
• An efficient, preconditioned, high-order solver for scattering by two-dimensional inhomogeneous media
• A Physically-Based Preconditioner for Quasi-Planar Scattering Problems
• A double-sweeping preconditioner for the Helmholtz equation

Cites work

• A class of iterative solvers for the Helmholtz equation: factorizations, sweeping preconditioners, source transfer, single layer potentials, polarized traces, and optimized Schwarz methods
• Effective algorithms with circulant-block matrices
• Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers?
• NETGEN: An advancing front 2D/3D-mesh generator based on abstract rules
• Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers
• Sweeping preconditioners for the iterative solution of quasiperiodic Helmholtz transmission problems in layered media
• The method of polarized traces for the 2D Helmholtz equation
• Why it is difficult to solve Helmholtz problems with classical iterative methods

Cited in

(2)

• Restrictions on the use of sweeping type preconditioners for Helmholtz problems
• Schwarz methods by domain truncation

Describes a project that uses

Uses Software

• Netgen
• NGSolve

This page was built for publication: Sweeping preconditioners for stratified media in the presence of reflections

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2218207)