Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Center for Computing Research
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Center for Computing Research; Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Mathematics and Statistics
- Technical Univ. of Dortmund, Dortmund (Germany). Inst. of Applied Mathematics (LS III)
The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. In conclusion, for time discretization, Crank–Nicolson scheme and backward Euler scheme are utilized.
- Research Organization:
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC04-94AL85000; NA0003525
- OSTI ID:
- 1464187
- Alternate ID(s):
- OSTI ID: 1548780
- Report Number(s):
- SAND-2017-8659J; 663067; TRN: US1902374
- Journal Information:
- Journal of Computational Physics, Vol. 361, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Monotonicity-preserving finite element schemes with adaptive mesh refinement for hyperbolic problems
|
journal | September 2020 |
On differentiable local bounds preserving stabilization for Euler equations | text | January 2019 |
Similar Records
Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations
On differentiable local bounds preserving stabilization for Euler equations