The ultraspherical spectral element method

Abstract

We introduce a novel spectral element method based on the ultraspherical spectral method and the hierarchical Poincaré-Steklov scheme for solving second-order linear partial differential equations on polygonal domains with unstructured quadrilateral or triangular meshes. Properties of the ultraspherical spectral method lead to almost banded linear systems, allowing the element method to be competitive in the high-polynomial regime ($p>5$). The hierarchical Poincaré-Steklov scheme enables precomputed solution operators to be reused, allowing for fast elliptic solves in implicit and semi-implicit time-steppers. The resulting spectral element method achieves an overall computational complexity of $\mathcal{O}(p^4/h^3)$ for mesh size $h$ and polynomial order $p$, enabling $hp$-adaptivity to be efficiently performed. We develop an open-source software system, ultraSEM, for flexible, user-friendly spectral element computations in MATLAB.

References in corpus (5)

• Firedrake: automating the finite element method by composing abstractions
• MFEM: a modular finite element methods library
• Nektar++: enhancing the capability and application of high-fidelity spectral/$hp$ element methods
• The automatic solution of partial differential equations using a global spectral method
• HPS Accelerated Spectral Solvers for Time Dependent Problems

Cited by in corpus (9)

• Learning developmental mode dynamics from single-cell trajectories
• A contour method for time-fractional PDEs and an application to fractional viscoelastic beam equations
• Asymmetric transport computations in Dirac models of topological insulators
• Fast global spectral methods for three-dimensional partial differential equations
• SlabLU: A Two-Level Sparse Direct Solver for Elliptic PDEs
• A spectral element method for modelling streamer discharges in low-temperature atmospheric-pressure plasmas
• Towards a Machine-Learned Poisson Solver for Low-Temperature Plasma Simulations in Complex Geometries
• Quantitative bounds on Impedance-to-Impedance operators with applications to fast direct solvers for PDEs
• Change of Basis between Classical Orthogonal Polynomials