### Abstract

Language | English |
---|---|

Pages | 3910-3937 |

Number of pages | 28 |

Journal | SIAM Journal on Numerical Analysis |

Volume | 47 |

Issue number | 5 |

Early online date | 16 Dec 2009 |

DOIs | |

Publication status | Published - 2009 |

### Fingerprint

### Keywords

- mass lumped scheme
- numerical dispersion
- spectral element method

### Cite this

*SIAM Journal on Numerical Analysis*,

*47*(5), 3910-3937. https://doi.org/10.1137/080724976

}

*SIAM Journal on Numerical Analysis*, vol. 47, no. 5, pp. 3910-3937. https://doi.org/10.1137/080724976

**Dispersive and dissipative behavior of the spectral element method.** / Ainsworth, Mark; Wajid, Hafiz Abdul.

Research output: Contribution to journal › Article

TY - JOUR

T1 - Dispersive and dissipative behavior of the spectral element method

AU - Ainsworth, Mark

AU - Wajid, Hafiz Abdul

PY - 2009

Y1 - 2009

N2 - 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.

AB - 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.

KW - mass lumped scheme

KW - numerical dispersion

KW - spectral element method

U2 - 10.1137/080724976

DO - 10.1137/080724976

M3 - Article

VL - 47

SP - 3910

EP - 3937

JO - SIAM Journal on Numerical Analysis

T2 - SIAM Journal on Numerical Analysis

JF - SIAM Journal on Numerical Analysis

SN - 0036-1429

IS - 5

ER -