Print Email Facebook Twitter A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation Title A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation Author Chen, J. (TU Delft Numerical Analysis) Dwarka, V.N.S.R. (TU Delft Numerical Analysis) Vuik, Cornelis (TU Delft Delft Institute of Applied Mathematics) Date 2024 Abstract The Helmholtz equation is related to seismic exploration, sonar, antennas, and medical imaging applications. It is one of the most challenging problems to solve in terms of accuracy and convergence due to the scalability issues of the numerical solvers. For 3D large-scale applications, high-performance parallel solvers are also needed. In this paper, a matrix-free parallel iterative solver is presented for the three-dimensional (3D) heterogeneous Helmholtz equation. We consider the preconditioned Krylov subspace methods for solving the linear system obtained from finite-difference discretization. The Complex Shifted Laplace Preconditioner (CSLP) is employed since it results in a linear increase in the number of iterations as a function of the wavenumber. The preconditioner is approximately inverted using one parallel 3D multigrid cycle. For parallel computing, the global domain is partitioned blockwise. The matrix-vector multiplication and preconditioning operator are implemented in a matrix-free way instead of constructing large, memory-consuming coefficient matrices. Numerical experiments of 3D model problems demonstrate the robustness and outstanding strong scaling of our matrix-free parallel solution method. Moreover, the weak parallel scalability indicates our approach is suitable for realistic 3D heterogeneous Helmholtz problems with minimized pollution error. Subject Helmholtz equationparallel computationmatrix-freegeometric multigridpreconditionerscalability To reference this document use: http://resolver.tudelft.nl/uuid:01b73100-531d-4923-8c36-bcbc4b4b1855 DOI https://doi.org/10.1553/etna_vol59s270 Embargo date 2024-07-15 ISSN 1068-9613 Source Electronic Transactions on Numerical Analysis, 59, 270 - 294 Bibliographical note Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public. Part of collection Institutional Repository Document type journal article Rights © 2024 J. Chen, V.N.S.R. Dwarka, Cornelis Vuik Files file embargo until 2024-07-15