×

Augmented AMG-shifted Laplacian preconditioners for indefinite Helmholtz problems. (English) Zbl 1374.65039

Summary: Discrete representations of the Helmholtz operator generally give rise to extremely difficult linear systems from an iterative solver perspective. This is due in part to the large oscillatory near null space of the linear system. Typical iterative methods do not effectively reduce error components in the subspace associated with this near null space. Traditional coarse grids used within multilevel solvers also cannot capture these components because of their oscillatory nature. While the shifted Laplacian is a popular method allowing multigrid preconditioners to achieve improved convergence, it is still unsatisfactory for higher-frequency problems. This paper analyzes the shifted Laplacian through polynomial smoothers to show its effect on the spectrum of the discrete Helmholtz operator. This analysis reveals desirable features of the shifted Laplacian preconditioners as well as limitations. Shifted Laplacian techniques are first extended to smoothed aggregation algebraic multigrid. Motivated by analysis, we propose augmenting the algebraic multigrid (AMG) shifted Laplacian with a two-grid error correction step. This additional step consists of the generalized minimal residual iterations applied to a projected version of the unshifted equations on a single auxiliary coarse grid. Results indicate that augmentation can improve the convergence significantly and can reduce the overall solve time on two-dimensional and three-dimensional problems.

MSC:

65F08 Preconditioners for iterative methods
65F10 Iterative numerical methods for linear systems
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs
65N12 Stability and convergence of numerical methods for boundary value problems involving PDEs

Software:

Matlab; Salinas; PSP; clique

References:

[1] VirieuxJ, OpertoS. An overview of full waveform inversion in exploration geophysics. Geophysics2009; 74(6):WCC127-WCC152.
[2] ErnstOG, GanderMJ. Why it is difficult to solve Helmholtz problems with classical iterative methods. In Numerical Analysis of Multiscale Problems, Vol. 83, GrahamI (ed.), HouT (ed.), LakkisO (ed.), ScheichlR (ed.) (eds). Springer‐Verlag: Heidelberg, Germany, 2011; 325-361.
[3] GanderMJ, NatafF. An incomplete LU preconditioner for problems in acoustics. Journal of Computational Acoustics2005; 13:455-476. · Zbl 1189.76362
[4] BenamouJ, DespresB. A domain decomposition method for the Helmholtz equation and related optimal control problems. Journal of Computational Physics1997; 136(1):68-82. · Zbl 0884.65118
[5] BrandtA, LivshitsI. Wave‐ray multigrid methods for standing wave equations. Electronic Transactions on Numerical Analysis1997; 6:162-181. · Zbl 0891.65127
[6] EngquistB, YingL. Sweeping preconditioner for the Helmholtz equation: hierarchical matrix representation. Communications on Pure and Applied Mathematics2011; 64:697-735. · Zbl 1229.35037
[7] EngquistB, YingL. Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers. SIAM Multiscale Modeling and Simulation2011; 9:686-710. · Zbl 1228.65234
[8] PoulsonJ, EngquistB, LiS, YingL. A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations. SIAM Journal on Scientific Computing2013; 35(3):194-212. · Zbl 1275.65073
[9] BaylissA, GoldsteinCI, TurkelE. An iterative method for the Helmholtz equation. Journal of Computational Physics1983; 49(3):443-457. · Zbl 0524.65068
[10] ErlanggaYA, VuikC, OosterleeCW. On a class of preconditioners for solving the Helmholtz equation. Applied Numerical Mathematics2004; 50(3‐4):409-425. · Zbl 1051.65101
[11] vanGijzenM, ErlanggaY, VuikC. Spectral analysis of the discrete Helmholtz operator preconditioned with a shifted Laplacian. SIAM Journal on Scientific Computing2007; 29(5):1942-1958. · Zbl 1155.65088
[12] CoolsS, VanrooseW. Local Fourier analysis of the complex shifted Laplacian preconditioner for Helmholtz problems. Numerical Linear Algebra with Applications2013; 20(4):575-597. · Zbl 1313.65320
[13] GanderM, GrahamI, SpenceE. How should one choose the shift for the shifted Laplacian to be a good preconditioner for the Helmholtz equation?, 2014.
[14] McCormickS, RugeJ. Unigrid for multigrid simulation. Mathematics of Computation1983; 41(163):43-62. · Zbl 0563.65065
[15] ManteuffelTA, McCormickSF, RöhrleO, RugeJ. Projection multilevel methods for quasilinear elliptic partial differential equations: numerical results. SIAM Journal on Numerical Analysis2006; 44(1):120-138. · Zbl 1145.65108
[16] FredericksonP, McBryanO. Parallel superconvergent multigrid. In Multigrid Methods: Theory, Applications, and Supercomputing, Vol. 110, McCormickS (ed.) (ed.). CRC Press: Boca Raton, FL, 1988; 195-210. · Zbl 0651.65028
[17] ChanT, TuminaroR. Analysis of a parallel multigrid algorithm. In Proceedings of the Fourth Copper Mountain Conference on Multigrid Methods, MandelJ (ed.) (ed.). SIAM: Philadelphia, 1989; 66-86.
[18] ElmanHC, ErnstOG, O’LearyDP. A multigrid method enhanced by Krylov subspace iteration for discrete Helmholtz equations. SIAM Journal on Scientific Computing2001; 23(4):1291-1315. · Zbl 1004.65134
[19] RepsB, VanrooseW, ZubairH. GMRES‐based multigrid for the complex scaled preconditioner for the indefinite Helmholtz equation. arXiv preprint arXiv:1012.5379 2010.
[20] OlsonL, SchroderJ. Smoothed aggregation for Helmholtz problems. Numerical Linear Algebra with Applications2010; 17(2‐3):361-386. · Zbl 1240.65359
[21] AdamsM. Algebraic multigrid methods for direct frequency response analyses in solid mechanics. Computational Mechanics2007; 39:497-507. · Zbl 1163.74043
[22] UmetaniN, MaclachlanS, OosterleeC. A multigrid‐based shifted Laplacian preconditioner for a fourth‐order Helmholtz discretization. Numerical Linear Algebra with Applications2009; 16(8):603-626. · Zbl 1224.65078
[23] CoolsS, RepsB, VanrooseW. A new level‐dependent coarse grid correction scheme for indefinite Helmholtz problems. Numerical Linear Algebra with Applications2014; 21(4):513-533. · Zbl 1340.65294
[24] ErlanggaY, NabbenR. On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian. Electronic Transactions on Numerical Analysis2008; 31:403-424. · Zbl 1190.65050
[25] SheikhA, LahayeD, VuikC. On the convergence of shifted Laplace preconditioner combined with multilevel deflation. Numerical Linear Algebra with Applications2013; 20(4):645-662. · Zbl 1313.65293
[26] ChewW, LiuQ. Perfectly matched layers for elastodynamics: a new absorbing boundary condition. Journal of Computational Acoustics1996; 4(04):341-359.
[27] BourgeoisA, BourgetM, LaillyP, PouletM, RicarteP, VersteegR. Marmousi, model and data. The Marmousi Experience: Proceedings of the 1990 EAEG Workshop on Practical Aspects of Seismic Data Inversion. In 52nd EAEG Meeting, Copenhagen, 1991; 5-16.
[28] BhardwajM, PiersonK, ReeseG, WalshT, DayD, AlvinK, LesoinneM. Salinas: a scalable software for high‐performance structural and solid mechanics simulations. Proceedings of the 2002 ACM/IEEE Conference on Supercomputing, IEEE Computer Society Press: Baltimore, Maryland, 2002; 1-19.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.