Picture of smart phone in human hand

World leading smartphone and mobile technology research at Strathclyde...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by University of Strathclyde researchers, including by Strathclyde researchers from the Department of Computer & Information Sciences involved in researching exciting new applications for mobile and smartphone technology. But the transformative application of mobile technologies is also the focus of research within disciplines as diverse as Electronic & Electrical Engineering, Marketing, Human Resource Management and Biomedical Enginering, among others.

Explore Strathclyde's Open Access research on smartphone technology now...

Dispersive and dissipative behavior of the spectral element method

Ainsworth, Mark and Wajid, Hafiz Abdul (2009) Dispersive and dissipative behavior of the spectral element method. SIAM Journal on Numerical Analysis, 47 (5). pp. 3910-3937. ISSN 0036-1429

Full text not available in this repository. (Request a copy from the Strathclyde author)

Abstract

If the nodes for the spectral element method are chosen to be the Gauss-Legendre-Lobatto points and a Lagrange basis is used, then the resulting mass matrix is diagonal and the method is sometimes then described as the Gauss-point mass lumped finite element scheme. We study the dispersive behavior of the scheme in detail and provide both a qualitative description of the nature of the dispersive and dissipative behavior of the scheme along with precise quantitative statements of the accuracy in terms of the mesh-size and the order of the scheme. We prove that (a) the Gauss-point mass lumped scheme (i.e., spectral element method) tends to exhibit phase lag whereas the (consistent) finite element scheme tends to exhibit phase lead; (b) the absolute accuracy of the spectral element scheme is $1/p$ times better than that of the finite element scheme despite the use of numerical integration; (c) when the order $p$, the mesh-size $h$, and the frequency of the wave $\omega$ satisfy $2p+1 \approx \omega h$ the true wave is essentially fully resolved. As a consequence, one obtains a proof of the general rule of thumb sometimes quoted in the context of spectral element methods: $\pi$ modes per wavelength are needed to resolve a wave.