Open AccessNumerical analysis of Krylov multigrid methods for stationary advection-diffusion equation
Author(s) -
Н. М. Евстигнеев
Publication year - 2019
Publication title -
journal of physics. conference series
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.21
H-Index - 85
eISSN - 1742-6596
pISSN - 1742-6588
DOI - 10.1088/1742-6596/1391/1/012080
Subject(s) - multigrid method , krylov subspace , discretization , mathematics , preconditioner , advection , generalized minimal residual method , dirichlet boundary condition , convection–diffusion equation , mathematical analysis , boundary value problem , linear system , partial differential equation , physics , thermodynamics
This paper outlines the problem of applying multigrid methods for the stationary scalar advection-diffusion equation for the given advection filed. We assume that all functions are smooth enough to be represented by discretization methods. Two main discretization are considered: finite difference method for Dirichlet boundary conditions and pseudo-spectral method for the periodic boundary. We test the following smoothers for multigrid methods: Jacobi smoother, Gauss-Seidel smoother, Krylov subspace smoother (GMRES method with and without preconditioners). The analysis is performed in the space formed by the cross product of discretization parameters, diffusion coefficient values, multigrid levels and smoothers. We demonstrate that the most efficient strategy depends on parameter value and given velocity field. Best variants include Gauss-Seidel smoothers which is optimal for advection-dominated problem while multigrid method is used as a preconditioner for a Krylov method. Such methods can be used for spectral or pseudo-spectral methods where explicit dense matrix storage is impossible.