Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 41, Issue 5, Pages S2-S18Publisher
SIAM PUBLICATIONS
DOI: 10.1137/18M1194997
Keywords
spectral element method; finite element method; high-order; low-order; preconditioner; algebraic multigrid
Categories
Ask authors/readers for more resources
Low-order finite element (FE) systems are considered as preconditioners for spectral element (SE) discretizations of the Poisson problem in canonical and complex domains. The FE matrices are based on the same mapped set of Gauss-Lobatto-Legendre (GLL) points as the SE discretization. Three different versions of the preconditioner based on combinations of the low-order stiffness and mass matrices are tested for 2D and 3D geometries. When building the preconditioning operators, a new meshing approach that allows elements to overlap without needing to fill out the volume of the mesh is explored and shown to be better than traditional schemes. These preconditioners are robust with respect to cell aspect ratio, demonstrate bounded iteration counts with hand p-refinement, and have lower iteration counts than scalable hybrid-Schwarz multigrid schemes currently used in production-level SE simulations. Overall costs for large-scale parallel applications are dependent on fast and robust solvers for the sparse FE systems. Algebraic multigrid is shown to offer a pathway to realizing a robust and fast preconditioning strategy in this context.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available