2017/05/03 by Juan Manzanero, Manzanero, Juan, Gonzalo Rubio +7 · 2 citations
Engineering · Mathematics · #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #FOS: Mathematics #Numerical Analysis (math.NA) #Numerical methods for differential equations
paper · pdf · doi:10.48550/arxiv.1705.01503
openalex publication_date 2017/05/03 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/01
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 L2 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, non-conservative 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 split form de-aliasing is restricted to DG schemes with the summation-by-parts simultaneous-approximation-term (SBP-SAT) properties (e.g. DG with Gauss-Lobatto points). Numerical experiments are included to illustrate the theoretical findings.