Vanka-smoothed shifted Laplacian multigrid preconditioners for the Helmholtz equations

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Yovel, Rachel, He, Yunhui, Treister, Eran
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866917096480309248
author Yovel, Rachel
He, Yunhui
Treister, Eran
author_facet Yovel, Rachel
He, Yunhui
Treister, Eran
contents We present an improved multigrid preconditioner for the acoustic Helmholtz equation with enhanced scalability. Standard multigrid fails to converge for the Helmholtz equation, and the well-known complex shifted Laplacian method overcomes it by adding a complex shift and using the shifted system as a preconditioner. However, the added complex shift grows with the frequency and interferes with the preconditioner's scalability. In this work, we present an additive Vanka smoother that requires a much lower shift than point-wise smoothers, and thereby enhances the scalability. By carefully designing different ingredients of the multigrid cycle, the presented method enables deep V-cycles with a small and bounded shift, even when many levels are used. We validate our method theoretically by local Fourier analysis, and hold numerical experiments for homogeneous and heterogeneous media. We show that our method outperforms plain shifted Laplacian in terms of runtimes and performs well on challenging geophysical media in 2D and 3D.
format Preprint
id arxiv_https___arxiv_org_abs_2511_16808
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Vanka-smoothed shifted Laplacian multigrid preconditioners for the Helmholtz equations
Yovel, Rachel
He, Yunhui
Treister, Eran
Numerical Analysis
65F10, 65N55, 35J05 65F10, 65N55, 35J05 65F10, 65N55, 35J05
We present an improved multigrid preconditioner for the acoustic Helmholtz equation with enhanced scalability. Standard multigrid fails to converge for the Helmholtz equation, and the well-known complex shifted Laplacian method overcomes it by adding a complex shift and using the shifted system as a preconditioner. However, the added complex shift grows with the frequency and interferes with the preconditioner's scalability. In this work, we present an additive Vanka smoother that requires a much lower shift than point-wise smoothers, and thereby enhances the scalability. By carefully designing different ingredients of the multigrid cycle, the presented method enables deep V-cycles with a small and bounded shift, even when many levels are used. We validate our method theoretically by local Fourier analysis, and hold numerical experiments for homogeneous and heterogeneous media. We show that our method outperforms plain shifted Laplacian in terms of runtimes and performs well on challenging geophysical media in 2D and 3D.
title Vanka-smoothed shifted Laplacian multigrid preconditioners for the Helmholtz equations
topic Numerical Analysis
65F10, 65N55, 35J05 65F10, 65N55, 35J05 65F10, 65N55, 35J05
url https://arxiv.org/abs/2511.16808