TY - JOUR
T1 - Algebraic entropy fixes and convex limiting for continuous finite element discretizations of scalar hyperbolic conservation laws
AU - Kuzmin, Dmitri
AU - Quezada de Luna, Manuel
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The work of Dmitri Kuzmin was supported by the German Research Association (DFG) under grant KU 1530/23-1. The authors would like to thank Hennes Hajduk (TU Dortmund University) for careful proofreading of the manuscript and helpful feedback.
PY - 2020/8/25
Y1 - 2020/8/25
N2 - In this work, we modify a continuous Galerkin discretization of a scalar hyperbolic conservation law using new algebraic correction procedures. Discrete entropy conditions are used to determine the minimal amount of entropy stabilization and constrain antidiffusive corrections of a property-preserving low-order scheme. The addition of a second-order entropy dissipative component to the antidiffusive part of a nearly entropy conservative numerical flux is generally insufficient to prevent violations of local bounds in shock regions. Our monolithic convex limiting technique adjusts a given target flux in a manner which guarantees preservation of invariant domains, validity of local maximum principles, and entropy stability. The new methodology combines the advantages of modern entropy stable/entropy conservative schemes and their local extremum diminishing counterparts. The process of algebraic flux correction is based on inequality constraints which provably provide the desired properties. No free parameters are involved. The proposed algebraic fixes are readily applicable to unstructured meshes, finite element methods, general time discretizations, and steady-state residuals. Numerical studies of explicit entropy-constrained schemes are performed for linear and nonlinear test problems.
AB - In this work, we modify a continuous Galerkin discretization of a scalar hyperbolic conservation law using new algebraic correction procedures. Discrete entropy conditions are used to determine the minimal amount of entropy stabilization and constrain antidiffusive corrections of a property-preserving low-order scheme. The addition of a second-order entropy dissipative component to the antidiffusive part of a nearly entropy conservative numerical flux is generally insufficient to prevent violations of local bounds in shock regions. Our monolithic convex limiting technique adjusts a given target flux in a manner which guarantees preservation of invariant domains, validity of local maximum principles, and entropy stability. The new methodology combines the advantages of modern entropy stable/entropy conservative schemes and their local extremum diminishing counterparts. The process of algebraic flux correction is based on inequality constraints which provably provide the desired properties. No free parameters are involved. The proposed algebraic fixes are readily applicable to unstructured meshes, finite element methods, general time discretizations, and steady-state residuals. Numerical studies of explicit entropy-constrained schemes are performed for linear and nonlinear test problems.
UR - http://hdl.handle.net/10754/662417
UR - https://linkinghub.elsevier.com/retrieve/pii/S0045782520305557
UR - http://www.scopus.com/inward/record.url?scp=85089801499&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2020.113370
DO - 10.1016/j.cma.2020.113370
M3 - Article
SN - 0045-7825
VL - 372
SP - 113370
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
ER -