Spectral analysis of continuous FEM for hyperbolic PDEs: influence of approximation, stabilization, and time-stepping

03/30/2021
by   Sixtine Michel, et al.
0

We study continuous finite element dicretizations for one dimensional hyperbolic partial differential equations. The main contribution of the paper is to provide a fully discrete spectral analysis, which is used to suggest optimal values of the CFL number and of the stabilization parameters involved in different types of stabilization operators. In particular, we analyze the streamline-upwind Petrov-Galerkin (SUPG) stabilization technique, the continuous interior penalty (CIP) stabilization method and the local projection stabilization (LPS). Three different choices for the continuous finite element space are compared: Bernstein polynomials, Lagrangian polynomials on equispaced nodes, and Lagrangian polynomials on Gauss-Lobatto cubature nodes. For the last choice, we only consider inexact quadrature based on the formulas corresponding to the degrees of freedom of the element, which allows to obtain a fully diagonal mass matrix. We also compare different time stepping strategies, namely Runge-Kutta (RK), strong stability preserving RK (SSPRK) and deferred correction time integration methods. The latter allows to alleviate the computational cost as the mass matrix inversion is replaced by the high order correction iterations. To understand the effects of these choices, both time-continuous and fully discrete Fourier analysis are performed. These allow to compare all the different combinations in terms of accuracy and stability, as well as to provide suggestions for optimal values discretization parameters involved. The results are thoroughly verified numerically both on linear and non-linear problems, and error-CPU time curves are provided. Our final conclusions suggest that cubature elements combined with SSPRK and CIP or LPS stabilization are the most promising combinations.

READ FULL TEXT

page 13

page 14

research
06/25/2019

A High-Order Lower-Triangular Pseudo-Mass Matrix for Explicit Time Advancement of hp Triangular Finite Element Methods

Explicit time advancement for continuous finite elements requires the in...
research
05/29/2022

Continuous finite element subgrid basis functions for Discontinuous Galerkin schemes on unstructured polygonal Voronoi meshes

We propose a new high order accurate nodal discontinuous Galerkin (DG) m...
research
09/07/2019

Subcell flux limiting for high-order Bernstein finite element discretizations of scalar hyperbolic conservation laws

This work extends the concepts of algebraic flux correction and convex l...
research
08/14/2019

Stability of explicit Runge-Kutta methods for high order finite element approximation of linear parabolic equations

We study the stability of explicit Runge-Kutta methods for high order La...
research
01/25/2022

Variationally consistent mass scaling for explicit time-integration schemes of lower- and higher-order finite element methods

In this paper, we propose a variationally consistent technique for decre...

Please sign up or login with your details

Forgot password? Click here to reset