Krylov-Subspace Preconditioners for Discontinuous Galerkin Finite Element Methods
نویسندگان
چکیده
Standard (conforming) finite element approximations of convection-dominated convectiondiffusion problems often exhibit poor stability properties that manifest themselves as nonphysical oscillations polluting the numerical solution. Various techniques have been proposed for the stabilisation of finite element methods (FEMs) for convection-diffusion problems, such as the popular streamline upwind Petrov-Galerkin (SUPG) method, and its variants. During the last decade, families of discontinuous Galerkin finite element methods (DGFEMs) have been proposed for the numerical solution of convection-diffusion problems, due to the many attractive properties they exhibit. In particular, DGFEMs admit good stability properties, they offer flexibility in the mesh design (irregular meshes are admissible) and in the imposition of boundary conditions (Dirichlet boundary conditions are weakly imposed), and they are increasingly popular in the context of hp-adaptive algorithms. The increase in popularity for DGFEMs has created a corresponding demand for developing corresponding linear solvers. This work aims to provide an overview of the current state of affairs in the solution of DGFEM-linear problems and present some recent results on the preconditioning of stiffness matrices arising from DGFEM discretisations of steady-state convection-diffusion boundaryvalue problems. More specifically, preconditioners are derived for which theoretical results indicate GMRES performance independent of discretisation parameters.
منابع مشابه
Preconditioned MHSS iteration methods for a class of block two-by-two linear systems with applications to distributed control problems
We construct a preconditioned modified Hermitian and skew-Hermitian splitting (PMHSS) iteration scheme for solving and preconditioning a class of block two-by-two linear systems arising from the Galerkin finite element discretizations of a class of distributed control problems. The convergence theory of this class of PMHSS iteration methods is established and the spectral properties of the PMHS...
متن کاملTruncated Hierarchical Preconditioning for the Stochastic Galerkin Fem
Stochastic Galerkin finite element discretizations of partial differential equations with coefficients characterized by arbitrary distributions lead, in general, to fully block dense linear systems. We propose two novel strategies for constructing preconditioners for these systems to be used with Krylov subspace iterative solvers. In particular, we present a variation of the hierarchical Schur ...
متن کاملSolving log-transformed random diffusion problems by stochastic Galerkin mixed finite element methods
Stochastic Galerkin finite element discretisations of PDEs with stochastically nonlinear coefficients lead to linear systems of equations with block dense matrices. In contrast, stochastic Galerkin finite element discretisations of PDEs with stochastically linear coefficients lead to linear systems of equations with block sparse matrices which are cheaper to manipulate and precondition in the f...
متن کاملPreconditionning Techniques for the Solution of the Helmholtz Equation by the Finite Element Method
This paper discusses 2D and 3D solutions of the harmonic Helmholtz equation by finite elements. It begins with a short survey of the absorbing and transparent boundary conditions associated with the DtN technique. The solution of the discretized system by means of a standard Galerkin or Galerkin least-squares (GLS) scheme is obtained by a preconditioned Krylov subspace technique, specifically a...
متن کاملOn Non-overlapping Domain Decomposition Preconditioners for Discontinuous Galerkin Finite Element Methods in H-type Norms
Abstract. We analyse the spectral bounds of non-overlapping domain decomposition additive Schwarz preconditioners for hp-version discontinuous Galerkin finite element methods in H-type norms. Using original approximation results for discontinuous finite element spaces, it is found that these preconditioners yield a condition number bound of order 1 + Hp/hq, where H and h are respectively the co...
متن کامل