Insights on Aliasing Driven Instabilities for Advection Equations with Application to Gauss-Lobatto Discontinuous Galerkin Methods

摘要

We analyse instabilities due to aliasing errors when solving one dimensional non-constant advection speed equations and discuss means to alleviate these types of errors when using high order discontinuous Galerkin (DG) schemes. First, we compare analytical bounds for the continuous and discrete version of the PDEs. Whilst traditional L-2 norm energy bounds applied to the discrete PDE do not always predict the physical behaviour of the continuous version of the equation, more strict elliptic norm bounds correctly bound the behaviour of the continuous PDE. Having derived consistent bounds, we analyse the effectiveness of two stabilising techniques: over-integration and split form variations (conservative, nonconservative and skew-symmetric). Whilst the former is shown to not alleviate aliasing in general, the latter ensures an aliasing-free solution if the splitting form of the discrete PDE is consistent with the continuous equation. The success of the split form de-aliasing is restricted to DG schemes with the summation-by-parts simultaneous-approximation-term properties (e.g. DG with Gauss-Lobatto points). Numerical experiments are included to illustrate the theoretical findings.