Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems
- Univ. of Zurich (Switzerland). Inst. of Mathematics
- Linköping Univ. (Sweden); Univ. of Johannesburg (South Africa)
- Johannes Gutenberg Univ., Mainz (Germany). Inst. of Mathematics
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE; SNF; UZH
- Grant/Contract Number:
- 89233218CNA000001; 200021_175784; FK-19-104; 200021_153604
- OSTI ID:
- 1823748
- Report Number(s):
- LA-UR-19-32410
- Journal Information:
- Journal of Scientific Computing, Vol. 85, Issue 2; ISSN 0885-7474
- Publisher:
- SpringerCopyright Statement
- Country of Publication:
- United States
- Language:
- English
A Broad Class of Conservative Numerical Methods for Dispersive Wave Equations | text | January 2020 |
Stabilizing Radial Basis Function Methods for Conservation Laws Using Weakly Enforced Boundary Conditions
|
journal | March 2021 |
Similar Records
Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)
Higher-order compositional modeling of three-phase flow in 3D fractured porous media based on cross-flow equilibrium