INVARIANT DOMAINS AND FIRST-ORDER CONTINUOUS FINITE ELEMENT APPROXIMATION FOR HYPERBOLIC SYSTEMS
Academic Article
Overview
Research
Identity
Additional Document Info
Other
View All
Overview
abstract
2016 Societ y for Industrial and Applied Mathematics. We propose a numerical method for solving general hyperbolic systems in any space dimension using forward Euler time stepping and continuous finite elements on nonuniform grids. The properties of the method are based on the introduction of an artificial dissipation that is defined so that any convex invariant set containing the initial data is an invariant domain for the method. The invariant domain property is proved for any hyperbolic system provided a CFL condition holds. The solution is also shown to satisfy a discrete entropy inequality for every admissible entropy of the system. The method is formally first-order accurate in space and can be made high-order in time by using strong stability preserving algorithms. This technique extends to continuous finite elements the work of [D. Hoff, Math. Comp., 33(1979), pp. 1171-1193], [D. Hoff, Trans. Amer. Math. Soc., 289(1985), pp. 591-610], and [H. Frid, Arch. Ration. Mech. Anal., 160(2001), pp. 245-269].