An application of graph-based algebraic multigrid method

Anssi Pennanen

Department of Mathematical Information Technology
University of Jyväskylä Finland


Abstract

In this study, we describe an Algebraic Multigrid Method (AMG) which is a modification of F. Kickinger's graph-based AMG [1]. In Kickinger's AMG, coarsening is based on the graph of the system matrix only, instead of strongly connected components in the matrix, as it is in the "classical" AMG by Ruge and Stüben [2]. This leads to fast computation of both coarse grid matrices and interpolation operators between coarse and fine levels. By some minor changes in the selection of coarse grid variables and processing of Dirichlet boundary conditions, we have obtained more robustness to the method. Additionally, we have extended the method in a way presented e.g. in [3] to be more suitable for vector-valued problems.

In numerical tests, we have used our method as well as a solver as a preconditioner for linear systems of equations arising from discretization of various mathematical models of physical phenomena, such as fluid flow, acoustics, linear elasticity, etc. Here we will show numerical results for two different application areas:

1) incompressible, viscous flow problems in 3-D, and,

2) scattering of sound waves in 2-D.

In the first case, AMG is used as a solver for the Oseen problem arising from Picard-type linearization of steady state Navier-Stokes. Linear Lagrangian finite elements are used as a discretization, stabilized by the well-known scheme by Franca, Frey and Hughes. Incomplete LU factorization with relaxation in smoothing steps is employed as a smoother. In the second case, our method is used to approximate inverse of the shifted-Laplacian operator which is used as a preconditioner for the Helmholz equation. In [4], this was done by using a geometric multigrid. The equations are discretized by linear, quadratic, and cubic Lagrangian finite elements. For both cases we have obtained good convergence results, but the stabilization scheme in the first case restricts the number of coarse levels.

References

[1] F. Kickinger, Algebraic Multi-grid Solver for Discrete Elliptic Second-Order Problems, in: Multigrid Methods V (Stuttgart, 1996), Springer, Berlin, 1998, 157-172.
[2] J.W. Ruge and K. Stüben, Algebraic Multigrid (AMG), in: S.F. McCormick (Ed.), Multigrid Methods, Frontiers in Applied Mathematics, SIAM, Philadelphia, Pennsylvania, 1987, pp. 73-130.
[3] M. Wabro, Algebraic Multigrid Methods for the Numerical Solution of the Incompressible Navier-Stokes Equations, Ph.D. thesis, Johannes Kepler Universität, Linz, 2003.
[4] Y. A. Erlangga, C. W. Oosterlee, C. Vuik, A novel multigrid based preconditioner for heterogeneous Helmholtz problems, SIAM J. Sci. Comput. 27 (4) (2006) 1471 1492.