Parametric Interpolation Framework for 1-D Scalar Conservation Laws with Non-Convex Flux Functions

arXiv1911.13174MaRDI QIDQ6330108FDOQ6330108

Authors: Geoffrey McGregor, Jean-Christophe Nave 

Publication date: 29 November 2019

Abstract: In this paper we present a novel framework for obtaining high order numerical methods for 1-D scalar conservation laws with non-convex flux functions. When solving Riemann problems, the Oleinik entropy condition, [16], is satisfied when the resulting shocks and rarefactions correspond to correct portions of the appropriate (upper or lower) convex envelope of the flux function. We show that the standard equal-area principle fails to select these solutions in general, and therefore we introduce a generalized equal-area principle which always selects the weak solution corresponding to the correct convex envelope. The resulting numerical scheme presented here relies on the area-preserving parametric interpolation framework introduced in [14] and locates shock position to fifth order in space, conserves area exactly and admits weak solutions which satisfy the Oleinik entropy condition numerically regardless of the initial states.

This page was built for publication: Parametric Interpolation Framework for 1-D Scalar Conservation Laws with Non-Convex Flux Functions

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