World Library  
Flag as Inappropriate
Email this Article

Spectral element method

Article Id: WHEBN0010504376
Reproduction Date:

Title: Spectral element method  
Author: World Heritage Encyclopedia
Language: English
Subject: Finite element method, Neumann–Neumann methods, Abstract additive Schwarz method, Mortar methods, Lax–Wendroff method
Publisher: World Heritage Encyclopedia

Spectral element method

In the numerical solution of partial differential equations, a topic in mathematics, the spectral element method (SEM) is a formulation of the finite element method (FEM) that uses high degree piecewise polynomials as basis functions. The spectral element method was introduced in a 1984 paper[1] by A. T. Patera.


The spectral method expands the solution in trigonometric series, a chief advantage being that the resulting method is of very high order. This approach relies on the fact that trigonometric polynomials are an orthonormal basis for L^2(\Omega) . The spectral element method chooses instead a high degree piecewise polynomial basis functions, also achieving a very high order of accuracy. Such polynomials are usually orthogonal Chebyshev polynomials or very high order Legendre polynomials over non-uniformly spaced nodes. In SEM computational error decreases exponentially as the order of approximating polynomial, therefore a fast convergence of solution to the exact solution is realized with fewer degrees of freedom of the structure in comparison with FEM. In structural health monitoring FEM can be used for detecting large flaws in a structure, but as the size of the flaw is reduced there is a need to use a high frequency wave with a small wavelength. Therefore, the FEM mesh must be much finer, resulting in increased computational time and an inexact solution. SEM, with fewer degrees of freedom per node, can be useful for detecting small flaws. Non-uniformity of nodes helps to make the mass matrix diagonal, which saves time and memory and is also useful for adopting central difference method (CDM). The disadvantages of SEM include modeling of complex geometry, compared to FEM.

A-priori error estimate

The classic analysis of Galerkin methods and Céa's lemma holds here and it can be shown that, if u is the solution of the weak equation, uN is the approximate solution and u \in H^{s+1}(\Omega):

\|u-u_N\|_{H^1(\Omega)} \leqq C_s N^{-s} \| u \|_{H^{s+1}(\Omega)}

where C is independent from N and s is no larger than the degree of the piecewise polynomial basis. As we increase N, we can also increase the degree of the basis functions. In this case, if u is an analytic function:

\|u-u_N\|_{H^1(\Omega)} \leqq C \exp( - \gamma N )

where \gamma depends only on u.

Related methods

  • G-NI or SEM-NI are the most used spectral methods. The Galerkin formulation of spectral methods or spectral element methods, for G-NI or SEM-NI respectively, is modified and Gaussian numerical integration is used instead of integrals in the definition of the bilinear form a(\cdot,\cdot) and in the functional F. These methods are a family of Petrov–Galerkin methods. Their convergence is a consequence of Strang's lemma.
  • The spectral element method uses a tensor product space spanned by nodal basis functions associated with Gauss–Lobatto points. In contrast, the p-version finite element method spans a space of high order polynomials by nodeless basis functions, chosen approximately orthogonal for numerical stability. Since not all interior basis functions need to be present, the p-version finite element method can create a space that contains all polynomials up to a given degree with fewer degrees of freedom.[2] However, some speedup techniques possible in spectral methods due to their tensor-product character are no longer available. The name p-version means that accuracy is increased by increasing the order of the approximating polynomials (thus, p) rather than decreasing the mesh size, h.
  • The hp finite element method (hp-FEM) combines the advantages of the h and p refinements to obtain exponential convergence rates.[3]


  1. ^ A. T. Patera. A spectral element method for fluid dynamics - Laminar flow in a channel expansion. Journal of Computational Physics, 54:468--488, 1984.
  2. ^ Barna Szabó and Ivo Babuška, Finite element analysis, John Wiley & Sons, Inc., New York, 1991. ISBN 0-471-50273-1
  3. ^ P. Šolín, K. Segeth, I. Doležel: Higher-order finite element methods, Chapman & Hall/CRC Press, 2003. ISBN 1-58488-438-X
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.

Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.