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

DOI10.1007/S10915-021-01632-7zbMATH Open1480.65268arXiv2103.16158OpenAlexW3200699243WikidataQ114225587 ScholiaQ114225587MaRDI QIDQ2233969FDOQ2233969

Authors: Sixtine Michel, Davide Torlo, M. Ricchiuto, R. Abgrall 

Publication date: 12 October 2021

Published in: Journal of Scientific Computing (Search for Journal in Brave)

Abstract: 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.

Full work available at URL: https://arxiv.org/abs/2103.16158

Recommendations

• Spectral analysis of high order continuous FEM for hyperbolic PDEs on triangular meshes: influence of approximation, stabilization, and time-stepping
• From \(h\) to \(p\) efficiently: optimal implementation strategies for explicit time-dependent problems using the spectral/\(hp\) element method
• The CFL Condition for Spectral Approximations to Hyperbolic Initial-Boundary Value Problems
• Analysis of the SBP-SAT stabilization for finite element methods. I: Linear problems
• Some Observations on the Interaction Between Linear and Nonlinear Stabilization for Continuous Finite Element Methods Applied to Hyperbolic Conservation Laws

zbMATH Keywords

continuous Galerkin methodmass lumpinghigh order accuracydispersion analysisstabilization techniquesnonstandard elements

Mathematics Subject Classification ID

Numerical quadrature and cubature formulas (65D32) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Numerical methods for discrete and fast Fourier transforms (65T50)

Cites Work

• Efficient implementation of essentially nonoscillatory shock-capturing schemes
• Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations
• Error analysis of the SUPG finite element discretization of evolutionary convection-diffusion-reaction equations
• Consistent SUPG-method for transient transport problems: stability and convergence
• Explicit Runge-Kutta schemes and finite elements with symmetric stabilization for first-order linear PDE systems
• Finite element methods for linear hyperbolic problems
• Unified Stabilized Finite Element Formulations for the Stokes and the Darcy Problems
• Global optimization of explicit strong-stability-preserving Runge-Kutta methods
• Title not available (Why is that?)
• An algorithm for computing Fekete points in the triangle
• Stabilized methods for compressible flows
• Entropy viscosity method for nonlinear conservation laws
• Semi-implicit spectral deferred correction methods for ordinary differential equations
• Spectral deferred correction methods for ordinary differential equations
• Stabilization of incompressibility and convection through orthogonal sub-scales in finite element methods
• Higher order triangular finite elements with mass lumping for the wave equation
• Title not available (Why is that?)
• A finite element formulation for the Stokes problem allowing equal velocity-pressure interpolation
• On the accuracy of the finite element method plus time relaxation
• Edge stabilization for Galerkin approximations of convection-diffusion-reaction problems
• Stabilization strategies for high order methods for transport dominated problems
• Title not available (Why is that?)
• A diagonal-mass-matrix triangular-spectral-element method based on cubature points
• Interior penalty continuous and discontinuous finite element approximations of hyperbolic equations
• Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling
• High order schemes for hyperbolic problems using globally continuous approximation and avoiding mass matrices
• Spatial eigenanalysis of spectral/\textit{hp} continuous Galerkin schemes and their stabilisation via DG-mimicking spectral vanishing viscosity for high Reynolds number flows
• Skeleton-stabilized immersogeometric analysis for incompressible viscous flow problems
• High-order residual distribution scheme for the time-dependent Euler equations of fluid dynamics
• A cut finite element method for a model of pressure in fractured media
• Arbitrary high-order, conservative and positivity preserving Patankar-type deferred correction schemes
• Stabilization of high order cut finite element methods on surfaces
• IDENT: identifying differential equations with numerical time evolution
• High order asymptotic preserving deferred correction implicit-explicit schemes for kinetic models
• A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes
• Analysis of the SBP-SAT stabilization for finite element methods. I: Linear problems
• A finite element time relaxation method

Cited In (6)

• Relaxation deferred correction methods and their applications to residual distribution schemes
• Novel well-balanced continuous interior penalty stabilizations
• A new efficient explicit deferred correction framework: analysis and applications to hyperbolic PDEs and adaptivity
• On the convergence of residual distribution schemes for the compressible Euler equations via dissipative weak solutions
• ``\textit{A posteriori} limited high order and robust schemes for transient simulations of fluid flows in gas dynamics
• Spectral analysis of high order continuous FEM for hyperbolic PDEs on triangular meshes: influence of approximation, stabilization, and time-stepping

Uses Software

• SciPy
• Python

This page was built for publication: Spectral analysis of continuous FEM for hyperbolic PDEs: influence of approximation, stabilization, and time-stepping

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2233969)