Towards Accuracy and Scalability: Combining Isogeometric Analysis with Deflation to Obtain Scalable Convergence for the Helmholtz Equation
arXiv:2010.10232 · doi:10.1016/j.cma.2021.113694
Abstract
Finding fast yet accurate numerical solutions to the Helmholtz equation remains a challenging task. The pollution error (i.e. the discrepancy between the numerical and analytical wave number k) requires the mesh resolution to be kept fine enough to obtain accurate solutions. A recent study showed that the use of Isogeometric Analysis (IgA) for the spatial discretization significantly reduces the pollution error. However, solving the resulting linear systems by means of a direct solver remains computationally expensive when large wave numbers or multiple dimensions are considered. An alternative lies in the use of (preconditioned) Krylov subspace methods. Recently, the use of the exact Complex Shifted Laplacian Preconditioner (CSLP) with a small complex shift has shown to lead to wave number independent convergence while obtaining more accurate numerical solutions using IgA. In this paper, we propose the use of deflation techniques combined with an approximated inverse of the CSLP using a geometric multigrid method. Numerical results obtained for both one- and two-dimensional model problems, including constant and non-constant wave numbers, show scalable convergence with respect to the wave number and approximation order p of the spatial discretization. Furthermore, when kh is kept constant, the proposed approach leads to a significant reduction of the computational time compared to the use of the exact inverse of the CSLP with a small shift.
References in corpus (10)
- Robust Multigrid for Isogeometric Analysis Based on Stable Splittings of Spline Spaces
- Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption
- A Robust Multigrid Method for Isogeometric Analysis using Boundary Correction
- L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation
- An improved sweeping domain decomposition preconditioner for the Helmholtz equation
- Robust multigrid solvers for the biharmonic problem in isogeometric analysis
- A p-multigrid method enhanced with an ILUT smoother and its comparison to h-multigrid methods within Isogeometric Analysis
- Two-level preconditioners for the Helmholtz equation
- The surrogate matrix methodology: Accelerating isogeometric analysis of waves
- Isogeometric solution of Helmholtz equation with Dirichlet boundary condition: numerical experiences