pith. sign in

arxiv: 1907.04115 · v1 · pith:Z7EMBJQTnew · submitted 2019-07-09 · 🧮 math.NA · cs.NA

Shock Capturing by Bernstein Polynomials for Scalar Conservation Laws

Jan Glaubitz This is my paper

Pith reviewed 2026-05-25 00:13 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords shock capturingBernstein polynomialsspectral element methodsscalar conservation lawstotal variation diminishingdiscontinuity sensorpositivity preservation
0
0 comments X

The pith

A convex combination with Bernstein polynomial reconstructions stabilizes high-order spectral approximations near shocks while preserving total variation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a shock-capturing procedure that replaces polluted high-order approximations with a convex blend of the original value and its Bernstein reconstruction inside elements flagged by a discontinuity sensor. Because the blend rests on classical Bernstein operators, the resulting scheme is total variation diminishing and exactly preserves monotone shock profiles for scalar conservation laws. The same construction can be adjusted to enforce solution bounds such as positivity. Unlike artificial-viscosity approaches, the procedure leaves the original CFL condition unchanged and requires only local modifications inside existing codes.

Core claim

The Bernstein procedure, formed by a sensor-controlled convex combination of a high-order spectral element approximation and its Bernstein reconstruction, is total variation diminishing and preserves monotone profiles; the same mechanism can be tuned to enforce bounds while leaving the time-step restriction untouched.

What carries the argument

The sensor-activated convex combination of the original approximation and its Bernstein polynomial reconstruction, applied only inside troubled elements.

If this is right

  • Spurious oscillations near discontinuities are suppressed without introducing new time-step restrictions.
  • Monotone shock profiles remain exactly monotone after each update.
  • Solution bounds such as positivity can be enforced by a simple adjustment of the same convex combination.
  • The method integrates into any existing spectral-element code by modifying only the reconstruction step inside flagged elements.
  • Global stability follows once the local TVD property holds in every troubled element.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same blending idea could be tested on systems of equations if an appropriate discontinuity sensor and component-wise Bernstein reconstruction are supplied.
  • Accuracy in smooth regions might be recovered by driving the blending coefficient to zero faster than the sensor threshold suggests.
  • The procedure's independence from the underlying time integrator suggests it could be paired with implicit or adaptive time-stepping schemes without further analysis of stability limits.
  • Because Bernstein operators converge to the identity as the polynomial degree rises, the method may recover full high-order accuracy once the sensor deactivates.

Load-bearing premise

The discontinuity sensor must correctly locate every element that contains a shock or discontinuity.

What would settle it

A scalar conservation-law test with a known discontinuity in which the sensor misses at least one troubled element and the computed solution then exhibits either persistent oscillations or an increase in total variation.

Figures

Figures reproduced from arXiv: 1907.04115 by Jan Glaubitz.

Figure 1
Figure 1. Figure 1: Illustration of the Bernstein reconstruction [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Parameter function α(S) as defined in (43). For the later numerical tests we investigated different other parameter functions as well, of which some have been discussed in [27]. Yet, we obtained the best results with (43). It should be noted that the above revisited PA sensor is only recommended for high orders N ≥ 4. In our numerical test, we have observed some miss-identifications for N = 3. For N ≤ 3 ot… view at source ↗
Figure 3
Figure 3. Figure 3: Parameter study for κ for the linear advection equation (65) at time t = 1. [22] J.-L. Guermond, R. Pasquetti, and B. Popov. Entropy viscosity method for nonlinear conservation laws. Journal of Computational Physics, 230(11):4248–4267, 2011. [23] H. Gzyl and J. L. Palacios. On the approximation properties of Bernstein polynomials via probabilistic tools. Boletın de la Asociaci´on Matem´atica Venezolana, 10… view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p019_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p022_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p023_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of the Bernstein procedure in a DG method with a usual filtering technique and no filtering [PITH_FULL_IMAGE:figures/full_fig_p024_9.png] view at source ↗
read the original abstract

A main disadvantage of many high-order methods for hyperbolic conservation laws lies in the famous Gibbs-Wilbraham phenomenon, once discontinuities appear in the solution. Due to the Gibbs-Wilbraham phenomenon, the numerical approximation will be polluted by spurious oscillations, which produce unphysical numerical solutions and might finally blow up the computation. In this work, we propose a new shock capturing procedure to stabilise high-order spectral element approximations. The procedure consists of going over from the original (polluted) approximation to a convex combination of the original approximation and its Bernstein reconstruction, yielding a stabilised approximation. The coefficient in the convex combination, and therefore the procedure, is steered by a discontinuity sensor and is only activated in troubled elements. Building up on classical Bernstein operators, we are thus able to prove that the resulting Bernstein procedure is total variation diminishing and preserves monotone (shock) profiles. Further, the procedure can be modified to not just preserve but also to enforce certain bounds for the solution, such as positivity. In contrast to other shock capturing methods, e.g. artificial viscosity methods, the new procedure does not reduce the time step or CFL condition and can be easily and efficiently implemented into any existing code. Numerical tests demonstrate that the proposed shock-capturing procedure is able to stabilise and enhance spectral element approximations in the presence of shocks.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes a shock-capturing procedure for high-order spectral element approximations of scalar conservation laws. In elements flagged by a discontinuity sensor, the polluted polynomial is replaced by a convex combination of the original approximation and its Bernstein reconstruction. Building on classical Bernstein operator properties, the authors claim to prove that the resulting procedure is total variation diminishing (TVD) and preserves monotone shock profiles; they also discuss a modification to enforce bounds such as positivity. The method is asserted to require no CFL reduction and to be straightforward to implement. Numerical tests are presented to illustrate stabilization of spectral element solutions near discontinuities.

Significance. If the global TVD claim holds, the approach would supply a non-intrusive, provably monotone shock-capturing technique that preserves high-order accuracy away from discontinuities and avoids the time-step penalties of artificial viscosity. The explicit reduction to classical Bernstein-operator theory for the local monotonicity proof is a clear methodological strength that distinguishes the work from purely heuristic limiters.

major comments (2)
  1. [Abstract and TVD proof section] Abstract and the section presenting the TVD proof: the global TVD and monotonicity-preservation claims are derived from local properties of the Bernstein operator and the convex combination, yet total variation is a global quantity. The argument does not address how inter-element numerical fluxes couple the locally modified elements to the rest of the domain when the sensor is imperfect.
  2. [Abstract and sensor definition] Abstract and discontinuity-sensor description: the entire global TVD guarantee is transferred to the auxiliary discontinuity sensor, whose precise definition, threshold, and error analysis (false-positive rate in smooth regions or false-negative rate at shocks) are not supplied. Without such analysis the local TVD property does not automatically extend to the assembled solution.
minor comments (1)
  1. [Numerical tests] Numerical experiments section: quantitative tables comparing L1 or L2 errors against standard TVD limiters or artificial-viscosity methods, together with explicit sensor-threshold values used in each test, would strengthen the validation.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments. Below we respond point-by-point to the major concerns, clarifying the local nature of the TVD proof and the dependence on the sensor. Partial revisions will be made to qualify the claims and add discussion of the sensor's role.

read point-by-point responses

  1. Referee: [Abstract and TVD proof section] Abstract and the section presenting the TVD proof: the global TVD and monotonicity-preservation claims are derived from local properties of the Bernstein operator and the convex combination, yet total variation is a global quantity. The argument does not address how inter-element numerical fluxes couple the locally modified elements to the rest of the domain when the sensor is imperfect.

    Authors: We agree that total variation is a global quantity and that our proof relies on local properties of the Bernstein operator applied element-wise. The manuscript establishes that the convex combination in a flagged element is TVD and monotone-preserving. However, when the sensor is imperfect, inter-element numerical fluxes can indeed couple the modified elements to the global solution and potentially affect the total variation. We will revise the abstract and TVD section to explicitly state that the TVD property holds locally within modified elements and to discuss the dependence on sensor accuracy for any global guarantee. This qualification addresses the coupling issue under the assumption of reliable sensor detection. revision: partial

  2. Referee: [Abstract and sensor definition] Abstract and discontinuity-sensor description: the entire global TVD guarantee is transferred to the auxiliary discontinuity sensor, whose precise definition, threshold, and error analysis (false-positive rate in smooth regions or false-negative rate at shocks) are not supplied. Without such analysis the local TVD property does not automatically extend to the assembled solution.

    Authors: The sensor employed in the numerical tests is defined in Section 3.2 as a modal-decay indicator with a fixed threshold. We acknowledge that the manuscript provides no rigorous error analysis of false-positive rates in smooth regions or false-negative rates at shocks. The global TVD claim therefore rests on the empirical performance of this sensor rather than a proven guarantee. In the revision we will add an explicit remark noting this limitation and stating that the local TVD property extends to the assembled solution only when the sensor correctly identifies troubled elements. A full sensor error analysis lies outside the present scope. revision: partial

standing simulated objections not resolved
  • Rigorous error analysis (false-positive and false-negative rates) of the discontinuity sensor

Circularity Check

0 steps flagged

No circularity; TVD proof derives from classical Bernstein operators

full rationale

The paper's central derivation states that the Bernstein procedure is TVD and preserves monotone profiles by building directly on properties of classical Bernstein operators. The discontinuity sensor is an auxiliary activation mechanism whose correctness is assumed but not part of the operator-based proof chain. No equations reduce a claimed result to a fitted parameter by construction, no self-citations are load-bearing for the uniqueness or ansatz of the result, and the derivation remains independent of the present paper's inputs. This is the normal case of a self-contained argument resting on external classical theory.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the procedure relies on classical properties of Bernstein operators and on the existence of a reliable local discontinuity sensor whose precise definition is not supplied.

free parameters (1)
  • discontinuity sensor threshold or scaling
    The coefficient of the convex combination is steered by a sensor whose internal parameters are not stated in the abstract.
axioms (1)
  • standard math Classical Bernstein operators are variation-diminishing and map monotone data to monotone polynomials
    The TVD and monotonicity claims are explicitly built on these known properties of Bernstein operators.

pith-pipeline@v0.9.0 · 5760 in / 1192 out tokens · 26346 ms · 2026-05-25T00:13:23.197795+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

60 extracted references · 60 canonical work pages

  1. [1]

    Ainsworth, G

    M. Ainsworth, G. Andriamaro, and O. Davydov. Bernstein–B´ ezier finite elements of arbitrary order and optimal assembly procedures. SIAM Journal on Scientific Computing , 33(6):3087–3109, 2011

    work page 2011

  2. [2]

    Anderson, V

    R. Anderson, V. Dobrev, T. Kolev, D. Kuzmin, M. Q. de Luna, R. Rieben, and V. Tomov. High-order local maxi- mum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation. Journal of Computational Physics, 334:102–124, 2017

    work page 2017

  3. [3]

    Archibald, A

    R. Archibald, A. Gelb, and J. Yoon. Polynomial fitting for edge detection in irregularly sampled signals and images. SIAM Journal on Numerical Analysis , 43(1):259–279, 2005

    work page 2005

  4. [4]

    G. E. Barter and D. L. Darmofal. Shock capturing with PDE-based artificial viscosity for DGFEM: Part i. Formulation. Journal of Computational Physics , 229(5):1810–1827, 2010

    work page 2010

  5. [5]

    C. E. Baumann and J. T. Oden. A discontinuous hp finite element method for convection—diffusion problems. Computer Methods in Applied Mechanics and Engineering , 175(3-4):311–341, 1999

    work page 1999

  6. [6]

    Beisiegel and J

    N. Beisiegel and J. Behrens. Quasi-nodal third-order Bernstein polynomials in a discontinuous Galerkin model for flooding and drying. Environmental Earth Sciences, 74(11):7275–7284, 2015

    work page 2015

  7. [7]

    D. Berman. On the impossibility of constructing a linear polynomial operator furnishing an approximation within the order of the best approximation. DOKLADY AKADEMII NAUK SSSR , 120(6):1175–1177, 1958

    work page 1958

  8. [8]

    Bernstein

    S. Bernstein. D´ emonstration du th´ eoreme de Weierstrass fond´ ee sur le calcul des probabilities.Comm. Soc. Math. Kharkov, 13(1):1–2, 1912

    work page 1912

  9. [9]

    Bernstein

    S. Bernstein. On the best approximation of continuous functions by polynomials of a given degree. Comm. Soc. Math. Kharkow, 2(13):49–194, 1912

    work page 1912

  10. [10]

    Burbeau, P

    A. Burbeau, P. Sagaut, and C.-H. Bruneau. A problem-independent limiter for high-order Runge–Kutta discontinuous Galerkin methods. Journal of Computational Physics , 169(1):111–150, 2001

    work page 2001

  11. [11]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu. Tvb Runge-Kutta local projection discontinuous Galerkin finite element method for conser- vation laws. ii. general framework. Mathematics of Computation , 52(186):411–435, 1989

    work page 1989

  12. [12]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu. The Runge-Kutta local projection P 1-discontinuous-Galerkin finite element method for scalar conservation laws. ESAIM: Mathematical Modelling and Numerical Analysis , 25(3):337–361, 1991

    work page 1991

  13. [13]

    C. M. Dafermos. Generalized characteristics and the structure of solutions of hyperbolic conservation laws. Indiana University Mathematics Journal , 26(6):1097–1119, 1977

    work page 1977

  14. [14]

    Dervieux, D

    A. Dervieux, D. Leservoisier, P.-L. George, and Y. Coudi` ere. About theoretical and practical impact of mesh adaptation on approximation of functions and PDE solutions. International Journal for Numerical Methods in Fluids , 43(5):507–516, 2003

    work page 2003

  15. [15]

    L. Evans. Measure theory and fine properties of functions . Routledge, 2018

    work page 2018

  16. [16]

    G. Faber. ¨Uber die interpolatorische Darstellung stetiger Funktionen. Jber. Deutsch. Math. Verein , 23:192–210, 1914

    work page 1914

  17. [17]

    Glaubitz and A

    J. Glaubitz and A. Gelb. High order edge sensors with ℓ1 regularization for enhanced discontinuous Galerkin methods. SIAM Journal on Scientific Computing , 41(2):A1304–A1330, 2019

    work page 2019

  18. [18]

    Glaubitz, A

    J. Glaubitz, A. Nogueira, J. Almeida, R. Cant˜ ao, and C. Silva. Smooth and compactly supported viscous sub-cell shock capturing for discontinuous Galerkin methods. Journal of Scientific Computing , 79:249–272, 2019

    work page 2019

  19. [19]

    Glaubitz, P

    J. Glaubitz, P. ¨Offner, and T. Sonar. Application of modal filtering to a spectral difference method. Mathematics of Computation, 87(309):175–207, 2018

    work page 2018

  20. [20]

    Gottlieb and C.-W

    D. Gottlieb and C.-W. Shu. On the Gibbs phenomenon and its resolution. SIAM Review, 39(4):644–668, 1997

    work page 1997

  21. [21]

    Gottlieb and C.-W

    S. Gottlieb and C.-W. Shu. Total variation diminishing Runge-Kutta schemes. Mathematics of computation of the American Mathematical Society, 67(221):73–85, 1998. 17 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref = 0.6 = 0.7 = 0.8 = 0.9 (a) N = 3 & I = 10. 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref = 0.6 = 0.7 = 0.8 = 0.9 (b) N = 4 & I ...

    work page 1998

  22. [22]

    Guermond, R

    J.-L. Guermond, R. Pasquetti, and B. Popov. Entropy viscosity method for nonlinear conservation laws. Journal of Computational Physics, 230(11):4248–4267, 2011

    work page 2011

  23. [23]

    Gzyl and J

    H. Gzyl and J. L. Palacios. On the approximation properties of Bernstein polynomials via probabilistic tools. Boletın de la Asociaci´ on Matem´ atica Venezolana, 10(1):5–13, 2003. 18 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref no filtering mean-value Bernstein (a) N = 3 & I = 20. 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref no filteri...

    work page 2003

  24. [24]

    A. Harten. On a class of high resolution total-variation-stable finite-difference schemes. SIAM Journal on Numerical Analysis, 21(1):1–23, 1984

    work page 1984

  25. [25]

    J. S. Hesthaven and T. Warburton. Nodal discontinuous Galerkin methods: Algorithms, analysis, and applications . Springer Science & Business Media, 2007

    work page 2007

  26. [26]

    Hewitt and R

    E. Hewitt and R. E. Hewitt. The Gibbs-Wilbraham phenomenon: An episode in Fourier analysis. Archive for History of Exact Sciences, 21(2):129–160, 1979

    work page 1979

  27. [27]

    Huerta, E

    A. Huerta, E. Casoni, and J. Peraire. A simple shock-capturing technique for high-order discontinuous Galerkin methods. International Journal for Numerical Methods in Fluids , 69(10):1614–1632, 2012

    work page 2012

  28. [28]

    H. T. Huynh. A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods. AIAA Paper, 4079:2007, 2007

    work page 2007

  29. [29]

    Karniadakis and S

    G. Karniadakis and S. Sherwin. Spectral/hp element methods for computational fluid dynamics . Oxford University Press, 2013

    work page 2013

  30. [30]

    R. C. Kirby. Fast simplicial finite element algorithms using Bernstein polynomials. Numerische Mathematik, 117(4):631– 652, 2011

    work page 2011

  31. [31]

    Kl¨ ockner, T

    A. Kl¨ ockner, T. Warburton, and J. S. Hesthaven. Viscous shock capturing in a time-explicit discontinuous Galerkin method. Mathematical Modelling of Natural Phenomena , 6(3):57–83, 2011

    work page 2011

  32. [32]

    Krivodonova, J

    L. Krivodonova, J. Xin, J.-F. Remacle, N. Chevaugeon, and J. E. Flaherty. Shock detection and limiting with discontinuous 19 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref no filtering mean-value Bernstein (a) N = 3 & I = 20. 0.0 0.2 0.4 0.6 0.8 1.0 x 0.0 0.2 0.4 0.6 0.8 1.0u ref no filtering mean-value Bernstein (b) N = 4 & I = 20. 0.0 0.2 0.4 0....

    work page 2004

  33. [33]

    Kuzmin, R

    D. Kuzmin, R. L¨ ohner, and S. Turek. Flux-Corrected Transport: Principles, Algorithms, and Applications . Scientific Computation. Springer Berlin Heidelberg, 2006

    work page 2006

  34. [34]

    P. D. Lax. Hyperbolic systems of conservation laws ii. Communications on Pure and Applied Mathematics, 10(4):537–566, 1957

    work page 1957

  35. [35]

    P. D. Lax. Hyperbolic systems of conservation laws and the mathematical theory of shock waves . SIAM, 1973

    work page 1973

  36. [36]

    R. J. LeVeque. Finite volume methods for hyperbolic problems , volume 31. Cambridge University Press, 2002

    work page 2002

  37. [37]

    T.-P. Liu. The Riemann problem for general systems of conservation laws. Journal of Differential Equations , 18(1):218– 234, 1975

    work page 1975

  38. [38]

    Lohmann, D

    C. Lohmann, D. Kuzmin, J. N. Shadid, and S. Mabuza. Flux-corrected transport algorithms for continuous Galerkin methods based on high order Bernstein finite elements. Journal of Computational Physics , 344:151–186, 2017

    work page 2017

  39. [39]

    G. G. Lorentz. Bernstein polynomials. American Mathematical Soc., 2012

    work page 2012

  40. [40]

    O. Oleinik. The Cauchy problem for nonlinear equations in a class of discontinuous functions. Amer. Math. Soc. Transl, 42:7–12, 1964

    work page 1964

  41. [41]

    O. A. Oleinik. Discontinuous solutions of non-linear differential equations. Uspekhi Matematicheskikh Nauk , 12(3):3–73, 1957

    work page 1957

  42. [42]

    A. T. Patera. A spectral element method for fluid dynamics: Laminar flow in a channel expansion. Journal of Computa- 20 0.0 0.2 0.4 0.6 0.8 1.0 x 0.92 0.94 0.96 0.98 1.00 1.02 1.04 1.06 1.08u ref no filtering mean-value Bernstein (a) N = 3 & I = 20. 0.0 0.2 0.4 0.6 0.8 1.0 x 0.92 0.94 0.96 0.98 1.00 1.02 1.04 1.06 1.08u ref no filtering mean-value Bernstein...

    work page 1984

  43. [43]

    Persson and J

    P.-O. Persson and J. Peraire. Sub-cell shock capturing for discontinuous Galerkin methods. In 44th AIAA Aerospace Sciences Meeting and Exhibit , page 112, 2006

    work page 2006

  44. [44]

    K. Petras. On the minimal norms of polynomial projections. Journal of Approximation Theory , 62(2):206–212, 1990

    work page 1990

  45. [45]

    G. M. Phillips. Interpolation and approximation by polynomials , volume 14. Springer Science & Business Media, 2003

    work page 2003

  46. [46]

    Qiu and C.-W

    J. Qiu and C.-W. Shu. A comparison of troubled-cell indicators for Runge–Kutta discontinuous Galerkin methods using weighted essentially nonoscillatory limiters. SIAM Journal on Scientific Computing , 27(3):995–1013, 2005

    work page 2005

  47. [47]

    J. L. Randall. Numerical methods for conservation laws. Lectures in Mathematics ETH Z¨ urich, 1992

    work page 1992

  48. [48]

    Ranocha, J

    H. Ranocha, J. Glaubitz, P. ¨Offner, and T. Sonar. Stability of artificial dissipation and modal filtering for flux reconstruc- tion schemes using summation-by-parts operators. Applied Numerical Mathematics, 128:1–23, 2018

    work page 2018

  49. [49]

    Richards

    F. Richards. A Gibbs phenomenon for spline functions. Journal of Approximation Theory , 66(3):334–351, 1991

    work page 1991

  50. [50]

    Riemann et al

    B. Riemann et al. ¨Uber die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite . Verlag der Dieterichschen Buchhandlung, 1860

    work page

  51. [51]

    D. G. Schaeffer. A regularity theorem for conservation laws. Advances in Mathematics, 11(3):368–386, 1973

    work page 1973

  52. [52]

    Shu and S

    C.-W. Shu and S. Osher. Efficient implementation of essentially non-oscillatory shock-capturing schemes. Journal of Computational Physics, 77(2):439–471, 1988

    work page 1988

  53. [53]

    Shu and S

    C.-W. Shu and S. Osher. Efficient implementation of essentially non-oscillatory shock-capturing schemes, ii. Journal of 21 0.0 0.2 0.4 0.6 0.8 1.0 x 0.900 0.925 0.950 0.975 1.000 1.025 1.050 1.075 1.100u ref no filtering mean-value Bernstein (a) N = 3 & I = 20. 0.0 0.2 0.4 0.6 0.8 1.0 x 0.925 0.950 0.975 1.000 1.025 1.050 1.075 1.100u ref no filtering mean-...

    work page 1989

  54. [54]

    Tadmor and T

    E. Tadmor and T. Tassa. On the piecewise smoothness of entropy solutions to scalar conservation laws. Communications in Partial Differential Equations , 18(9-10):1631–1652, 1993

    work page 1993

  55. [55]

    E. F. Toro. Riemann solvers and numerical methods for fluid dynamics: A practical introduction . Springer Science & Business Media, 2013

    work page 2013

  56. [56]

    L. N. Trefethen and D. Bau III. Numerical linear algebra, volume 50. Siam, 1997

    work page 1997

  57. [57]

    VonNeumann and R

    J. VonNeumann and R. D. Richtmyer. A method for the numerical calculation of hydrodynamic shocks. Journal of Applied Physics, 21(3):232–237, 1950

    work page 1950

  58. [58]

    Weierstrass

    K. Weierstrass. ¨Uber die analytische Darstellbarkeit sogenannter willk¨ urlicher Functionen einer reellen Ver¨ anderlichen. Sitzungsberichte der K¨ oniglich Preußischen Akademie der Wissenschaften zu Berlin, 2:633–639, 1885

    work page

  59. [59]

    B. Wendroff. The Riemann problem for materials with nonconvex equations of state i: Isentropic flow. Journal of Mathematical Analysis and Applications , 38(2):454–466, 1972

    work page 1972

  60. [60]

    B. Wendroff. The Riemann problem for materials with nonconvex equations of state: Ii: General flow. Journal of Mathematical Analysis and Applications , 38(3):640–658, 1972. 22 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 x 0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 u ref no filtering mean-value Bernstein (a) N = 3 & I = 20. 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00...

    work page 1972