Abstract:An efficient finite-difference method for solving the isotropic Helmholtz equation relies on a discretization scheme and an appropriate solver. Accordingly, we have adopted an average-derivative optimal scheme that has two advantages: (1) it can be applied to unequal directional sampling intervals and (2) it requires less than four grid points of sampling per wavelength. Direct methods are not of interest for industry-sized problems due to the high memory requirements; Krylov subspace methods such as the biconjugate gradient stabilized method and the flexible generalized minimal residual method that combine a multigrid-based preconditioner are better alternatives. However, standard geometric multigrid algorithms fail to converge when there exist unequal directional sampling intervals; this is called anisotropic grids in terms of the multigrid. We first review our previous research on 2D anisotropic grids: the semicoarsening strategy, line-wise relaxation operator, and matrix-dependent interpolation were used to modify the standard V-cycle multigrid algorithms, resulting in convergence. Although directly extending to the 3D case by substituting line relaxation for plane relaxation deteriorates the convergence rate considerably, we then find that a multilevel generalized minimal residual preconditioner-combined semicoarsening strategy is more suitable for anisotropic grids and the convergence rate is faster in the 2D and 3D cases. The results of the numerical experiments indicate that the standard geometric multigrid does not work for anisotropic grids, whereas our method demonstrates a faster convergence rate than the previous method.

Multigrid in H(div) on Axisymmetric Domains

The Hodge Laplacian on Axisymmetric Domains and its Discretization

Convergence Analysis of Hss-Multigrid Methods for Second-Order Nonselfadjoint Elliptic Problems

Notes on Multigrid Methods

Uniform Convergence of Adaptive Multigrid Methods for Elliptic Problems and Maxwell's Equations

Homogeneous multigrid for hybrid discretizations: application to HHO methods

Stand-alone Multigrid for Helmholtz Revisited: Towards Convergence Using Standard Components

A multigrid reduction framework for domains with symmetries

Uniform Convergence of a Multigrid Method for Elliptic Equations with Anisotropic Coefficients

Multigrid Methods for Hellan–Herrmann–Johnson Mixed Method of Kirchhoff Plate Bending Problems

An anisotropic multilevel preconditioner for solving the Helmholtz equation with unequal directional sampling intervals

An efficient multigrid method with preconditioned smoother for two-dimensional anisotropic space-fractional diffusion equations

Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems

Multigrid Solutions Based on the Implicit Schemes of the Three-Dimensional Wave Equation

Multigrid as an exact solver

Convergence Analysis of a Multigrid Method for a Nonlocal Model.

A robust solver for H(curl) convection-diffusion and its local Fourier analysis

Multigrid for two-sided fractional differential equations discretized by finite volume elements on graded meshes

Generalizing Reduction-Based Algebraic Multigrid

$hp$-Multigrid preconditioner for a divergence-conforming HDG scheme for the incompressible flow problems

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation