Ulrik S. Fjordholm

Ulrik S. Fjordholm, Siddhartha Mishra, Eitan Tadmor

We prove stability estimates for the ENO reconstruction and ENO interpolation procedures. In particular, we show that the jump of the reconstructed ENO pointvalues at each cell interface has the same sign as the jump of the underlying cell averages across that interface. We also prove that the jump of the reconstructed values can be upper-bounded in terms of the jump of the underlying cell averages. Similar sign properties hold for the ENO interpolation procedure. These estimates, which are shown to hold for ENO reconstruction and interpolation of arbitrary order of accuracy and on non-uniform meshes, indicate a remarkable rigidity of the piecewise-polynomial ENO procedure.

Ulrik S. Fjordholm, Susanne Solem

We prove that a class of monotone, \emph{$W_1$-contractive} schemes for scalar conservation laws converge at a rate of $Δx^2$ in the Wasserstein distance ($W_1$-distance), whenever the initial data is decreasing and consists of a finite number of piecewise constants. It is shown that the Lax--Friedrichs, Enquist--Osher and Godunov schemes are $W_1$-contractive. Numerical experiments are presented to illustrate the main result. To the best of our knowledge, this is the first proof of second-order convergence of any numerical method for discontinuous solutions of nonlinear conservation laws.

Ulrik S. Fjordholm, Deep Ray

We propose a third-order WENO reconstruction which satisfies the sign property, required for constructing high resolution entropy stable finite difference scheme for conservation laws. The reconstruction technique, which is termed as SP-WENO, is endowed with additional properties making it a more robust option compared to ENO schemes of the same order. The performance of the proposed reconstruction is demonstrated via a series of numerical experiments for linear and nonlinear scalar conservation laws. The scheme is easily extended to multi-dimensional conservation laws.

Ulrik S. Fjordholm, Roger Käppeli, Siddhartha Mishra et al.

Entropy solutions have been widely accepted as the suitable solution framework for systems of conservation laws in several space dimensions. However, recent results in \cite{CDL1,CDL2} have demonstrated that entropy solutions may not be unique. In this paper, we present numerical evidence that demonstrates that state of the art numerical schemes \emph{may not} necessarily converge to an entropy solution of systems of conservation laws as the mesh is refined. Combining these two facts, we argue that entropy solutions may not be suitable as a solution framework for systems of conservation laws, particularly in several space dimensions. Furthermore, we propose a more general notion, that of \emph{entropy measure valued solutions}, as an appropriate solution paradigm for systems of conservation laws. To this end, we present a detailed numerical procedure, which constructs stable approximations to entropy measure valued solutions and provide sufficient conditions that guarantee that these approximations converge to an entropy measure valued solution as the mesh is refined, thus providing a viable numerical framework for systems of conservation laws in several space dimensions. A large number of numerical experiments that illustrate the proposed schemes are presented and are utilized to examine several interesting properties of the computed entropy measure valued solutions.