High-order accurate entropy stable numercial schemes for hyperbolic conservation laws

The role of entropy for the stability of hyperbolic conservation laws is well-understood, and to some extent, also for firstand second-order accurate numerical schemes. The stability of higher-order accurate schemes, however, is to a large degree an open problem. In this thesis we adopt the framework of entropy stability as a design principle. We combine entropy conservative fluxes with appropriate diffusion operators, and to obtain high order accuracy, we consider diffusion operators using the ENO reconstruction procedure. We show that the ENO procedure satisfies the so-called sign property, which ensures entropy stability of our (arbitrarily) high-order accurate finite difference scheme. Moreover, we pose and argue for a conjecture on the total variation of the ENO reconstruction, which would imply a weak total variation bound for our scheme. For hyperbolic systems, the reconstruction is performed using a novel, computationally efficient characteristic-wise decomposition. The scheme is easily generalized to multi-dimensional systems on Cartesian meshes. To study the convergence properties of our scheme for scalar equations, we consider the framework of compensated compactness. Under the assumption that our conjecture on the total variation bound of ENO holds, we prove convergence to the entropy solution. For multidimensional hyperbolic systems we are able to prove convergence to an entropy measurevalued solution. The robustness, accuracy and computational efficiency of the scheme are demonstrated in a series of numerical experiments.

[1]  Eitan Tadmor,et al.  Arbitrarily High-order Accurate Entropy Stable Essentially Nonoscillatory Schemes for Systems of Conservation Laws , 2012, SIAM J. Numer. Anal..

[2]  E. Hopf On the Right Weak Solution of the Cauchy Problem for a Quasilinear Equation of First Order Final Report , 1969 .

[3]  Steven Schochet,et al.  Examples of measure-valued solutions , 1989 .

[4]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[5]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[6]  Nira Dyn,et al.  Approximation of piecewise smooth functions and images by edge-adapted (ENO-EA) nonlinear multiresolution techniques , 2008 .

[7]  J. Ball A version of the fundamental theorem for young measures , 1989 .

[8]  A. Harten ENO schemes with subcell resolution , 1989 .

[9]  E. Panov Erratum to: Existence and Strong Pre-compactness Properties for Entropy Solutions of a First-Order Quasilinear Equation with Discontinuous Flux , 2010 .

[10]  P. Lax Hyperbolic systems of conservation laws II , 1957 .

[11]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[12]  Philippe G. LeFloch,et al.  Fully Discrete, Entropy Conservative Schemes of ArbitraryOrder , 2002, SIAM J. Numer. Anal..

[13]  Yann Brenier,et al.  Weak-Strong Uniqueness for Measure-Valued Solutions , 2009, 0912.1028.

[14]  P. Lax Shock Waves and Entropy , 1971 .

[15]  Ulrik Skre Fjordholm,et al.  Structure preserving finite volume methods for the shallow water equations , 2009 .

[16]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier—Stokes equations and the second law of thermodynamics , 1986 .

[17]  Timothy J. Barth,et al.  Numerical Methods for Gasdynamic Systems on Unstructured Meshes , 1997, Theory and Numerics for Conservation Laws.

[18]  M. Mock,et al.  Systems of conservation laws of mixed type , 1980 .

[19]  Eitan Tadmor,et al.  ENO Reconstruction and ENO Interpolation Are Stable , 2011, Found. Comput. Math..

[20]  M. Crandall,et al.  Monotone difference approximations for scalar conservation laws , 1979 .

[21]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[22]  Eitan Tadmor,et al.  Numerical Viscosity and the Entropy Condition for Conservative Difference Schemes , 1984 .

[23]  J. Málek Weak and Measure-valued Solutions to Evolutionary PDEs , 1996 .

[24]  Eitan Tadmor,et al.  Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems , 2003, Acta Numerica.

[25]  Chi-Wang Shu Numerical experiments on the accuracy of ENO and modified ENO schemes , 1990 .

[26]  P. Lax Hyperbolic systems of conservation laws , 2006 .

[27]  B. Piccoli,et al.  Well-posedness of the Cauchy problem for × systems of conservation laws , 2000 .

[28]  Michael Westdickenberg,et al.  Convergence of Approximate Solutions of Conservation Laws , 2015, 1502.00798.

[29]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[30]  A. Harten On the symmetric form of systems of conservation laws with entropy , 1983 .

[31]  Chi-Wang Shu,et al.  High Order ENO and WENO Schemes for Computational Fluid Dynamics , 1999 .

[32]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[33]  Nira Dyn,et al.  Interpolation and Approximation of Piecewise Smooth Functions , 2005, SIAM J. Numer. Anal..

[34]  Luc Tartar,et al.  Compensated compactness and applications to partial differential equations , 1979 .

[35]  Eitan Tadmor,et al.  Energy Preserving and Energy Stable Schemes for the Shallow Water Equations , 2009 .

[36]  Tadmor Entropy functions for symmetric systems of conservation laws. Final Report , 1987 .

[37]  O. Oleinik Discontinuous solutions of non-linear differential equations , 1963 .

[38]  Pierre Degond,et al.  Polynomial upwind schemes for hyperbolic systems , 1999 .

[39]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[40]  J. Glimm Solutions in the large for nonlinear hyperbolic systems of equations , 1965 .

[41]  James M. Hyman,et al.  On Finite-Difference Approximations and Entropy Conditions for Shocks , 2015 .

[42]  R. Bartle The elements of integration and Lebesgue measure , 1995 .

[43]  R. J. Diperna,et al.  Measure-valued solutions to conservation laws , 1985 .

[44]  Eitan Tadmor,et al.  The numerical viscosity of entropy stable schemes for systems of conservation laws. I , 1987 .

[45]  Applications of the Theory of Compensated Compactness , 1986 .

[46]  E. Panov Existence and Strong Pre-compactness Properties for Entropy Solutions of a First-Order Quasilinear Equation with Discontinuous Flux , 2010 .

[47]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[48]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[49]  S. Osher Riemann Solvers, the Entropy Condition, and Difference , 1984 .

[50]  A. Bressan,et al.  L1 Stability Estimates for n×n Conservation Laws , 1999 .

[51]  S. Kružkov FIRST ORDER QUASILINEAR EQUATIONS IN SEVERAL INDEPENDENT VARIABLES , 1970 .

[52]  Philip L. Roe,et al.  Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks , 2009, J. Comput. Phys..

[53]  E. Yu. Panov,et al.  Uniqueness of the solution of the Cauchy problem for a first order quasilinear equation with one admissible strictly convex entropy , 1994 .

[54]  Tony F. Chan,et al.  ENO-Wavelet Transforms for Piecewise Smooth Functions , 2002, SIAM J. Numer. Anal..