pith. sign in

arxiv: 2605.00189 · v1 · submitted 2026-04-30 · 🧮 math.AP

Local Asymptotic Patterns for Viscous Approximations of Conservation Laws

Alberto Bressan , Laura Caravenna , Wen Shen This is my paper

Pith reviewed 2026-05-09 19:57 UTC · model grok-4.3

classification 🧮 math.AP
keywords viscous approximationsconservation lawsvanishing viscosityeternal solutionslocal rescalingshock singularitieshyperbolic equationsasymptotic patterns
0
0 comments X

The pith

Viscous approximations to conservation laws converge under local rescaling to eternal solutions near shocks.

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

The paper studies the local behavior of viscous approximations to solutions of hyperbolic conservation laws near points where the limiting solution develops singularities. For the vanishing-viscosity case it proves that, given a sequence of ε-approximate solutions, a suitable local rescaling of space and time produces, as ε tends to zero, a unique limit that is an eternal solution to the viscous equation, defined for all times and all positions. This holds when the underlying hyperbolic solution has a shock, an interaction of two shocks, or the formation of a new shock. A reader cares because the result supplies an explicit description of how the smoothing effect of viscosity unfolds at discontinuities, replacing vague convergence statements with concrete local patterns.

Core claim

Given a sequence of ε-approximate solutions, in the presence of a hyperbolic solution possessing one of the specified singularities (a point along a shock, the interaction of two shocks, or the formation of a new shock), a local rescaling of coordinates yields, as ε approaches zero, a well-defined limit that is an eternal solution to the viscous conservation law, globally defined in both space and time.

What carries the argument

Local rescaling of space-time coordinates around the singularity that extracts a globally defined eternal solution as the viscosity parameter vanishes.

Load-bearing premise

The hyperbolic solution has one of the three listed singularity types and the given sequence of approximations is compatible with the local rescaling construction.

What would settle it

Construct a concrete viscous approximation near a known shock and compute its rescaled profile; the claim fails if the limit either does not exist or fails to satisfy the viscous equation for all positive and negative times.

Figures

Figures reproduced from arXiv: 2605.00189 by Alberto Bressan, Laura Caravenna, Wen Shen.

Figure 1
Figure 1. Figure 1: The three types of generic singularities for a solution to a conservation law. A: a point along a shock. B: a point where a new shock forms. C: a point where two shocks merge. Here the dashed lines denote characteristics, while the solid curves denote shocks. For such generic solutions, singularities can only be of the above types (I)–(III). For each type, there is a natural asymptotic rescaling that chara… view at source ↗
Figure 2
Figure 2. Figure 2: Three examples of backward domains for a point P. Left: a point where the solution u = u(t, x) is continuous. Center: a point where u has a shock. Right: a point where two shocks merge. Definition 1.1 Let u = u(t, x) be a solution to the inviscid equation (1.3). We say that Ω ⊂ R 2 is a backward domain for the point P = (τ, ξ) if Ω contains all backward characteristics starting from a neighborhood of P. De… view at source ↗
Figure 3
Figure 3. Figure 3: An initial data satisfying the conditions (4.1). We seek an estimate of how long it takes for the solution to become close to a traveling profile. Notice that (4.1) does not require the existence of the limits limx→±∞ u¯(x). There￾fore, as t → +∞, the solution may not approach asymptotically any viscous shock profile. Differently from most literature on the topic, rather than on the asymptotic limit here w… view at source ↗
Figure 4
Figure 4. Figure 4: The initial phase in which a viscous shock profile is approached. Left: a curve γ corresponding to general initial data, as in (4.1). Right: after a first phase, the curve γ is entirely contained in a vertical strip, where u ∈ [u + − δ1, u− + δ1]. _ _ _ _ γ γ + P P + P P 2 3 u + u u + u t = T t = T [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The second and third phase in which a viscous shock profile is approached. Left: for t ≥ T2, the curve γ is entirely contained in a neighborhood of the convex closure of the graph of f, between u + and u −. Right: at the end of a third phase, the curve γ is contained in a thin strip around the segment with endpoints P +, P −. This implies that for t ≥ T3 the function u(t, ·) is well approximated by a visco… view at source ↗
Figure 6
Figure 6. Figure 6: Left: proof of Lemma 3.1, part (i). Right: the construction to prove part (iii). 1 2 ~ γ u _ γ γ ~ γ w u + u f _u w q(u) w (t,u) + u − δ1 u + u _ _ δ 1 u + u + f [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Left: the upper solution w +(t, u) constructed at (4.12). Right: the construction used in the proof of Lemma 4.3. The curve γ corresponding to the graph of u(T2, ·), as defined at (3.9), and the segments γe1, γe2, γb, corresponding to the traveling profiles ψe1, ψe2, ψb. In (4.30), the u-coordinates of the endpoints of the segments γe1, γe2, γb are denoted by ψe± 1 , ψe± 2 , ψb± respectively. In the settin… view at source ↗
Figure 8
Figure 8. Figure 8: Left: in the setting shown in [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Left: at time t = T2, the crossings between the graphs of u and ψe1 remain to the left of the crossings between the graphs of u and ψe2. Right: for t − T2 large enough, the graphs of u(t, ·) and ψe1 cross at a single point. By the second identity in (4.27), we have the implication ( u(T2, x) ≥ ψe 1(x − eλT2), u(T2, y) ≤ ψe 2(y − eλT2), =⇒ x < y. (4.35) Indeed, since ψe+ 1 − ψe+ 2 = 3δ1, by (4.27) this foll… view at source ↗
Figure 10
Figure 10. Figure 10: An initial data at τ << 0, obtained by interpolating two viscous traveling profiles. Lemma 5.1 There exists a unique solution W = W(t, x) of (5.1), defined for all (t, x) ∈ R 2 , such that lim τ→−∞ (Z λ ∗τ −∞ [PITH_FULL_IMAGE:figures/full_fig_p025_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Left: the function ϕ +(s) defined as (6.14). Right: the assumption (6.7) requires that characteristic curves cross the shaded domain transversally. 2. We shall estimate W on the region where x > λt by establishing a bound on the integrated function Φ(t, x) .= Z +∞ x W(t, y) dy. By (6.9), this function satisfies Φt − A(U, V )Φx = Φxx , (6.12) with initial and boundary conditions Φ(0, x) = 0 for x ≥ 0, Φx(t… view at source ↗
Figure 12
Figure 12. Figure 12: The comparison argument used in the proof of (6.16). 3. By a comparison argument, we claim that −Φ +(t, x) ≤ Z +∞ x [PITH_FULL_IMAGE:figures/full_fig_p029_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Left: the construction used in the proof of Theorem 1.1. For |x − γ(t)| > ε1/3 , by (7.3) the viscous approximation uε is close to the inviscid solution u. Right: the domain Dε in (7.19). 7.1 A point along a single shock. 1. Assume that the inviscid solution has a shock along the curve x = γ(t), passing through the point P = (τ, ξ). Call u ± .= u(τ, ξ±) the left and right states at time t = τ . W.l.o.g. w… view at source ↗
Figure 14
Figure 14. Figure 14: Left: the shock curves Γ1ε, Γ2ε, in the rescaled coordinates. Based on the bound (7.3), transversality of characteristics is ensured only at a distance > ε−2/3 . Right: the argument at (7.31) yields the transversality of characteristics up to time ε −1/15 . and consider the straight lines (see [PITH_FULL_IMAGE:figures/full_fig_p036_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: The solution to Burgers’ equation implicitly defined by (8.3), for t ≤ 0. Let z = z(t, x) be the backward solution to Burgers’ equation zt + zzx = 0, t ≤ 0, (8.1) with terminal data z(0, x) = −x 1/3 . (8.2) For t < 0 this solution is smooth (see [PITH_FULL_IMAGE:figures/full_fig_p038_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: The construction used in the proof of Theorem 1.2. Here the trapezoidal domain has the form Dε =  (t, x) ; t ∈ [−ε −α, ε−γ ], |x| ≤ rε(t) [PITH_FULL_IMAGE:figures/full_fig_p043_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Left: the function xe(τ, ·) with the properties (P1) is obtained by interpolating the graphs of x(τ, ·) and Y (u) = −u 3 (dashed line). Right: the flux function fewith the properties (P2) is obtained by interpolating the graphs of f and B(u) = u 2/2 (dashed line). (P1) The profile x 7→ ue(τ, x) is monotone decreasing. Moreover, its inverse function u 7→ xe(τ, u) satisfies (see [PITH_FULL_IMAGE:figures/fu… view at source ↗
Figure 18
Figure 18. Figure 18: The construction showing that the modifications of the data in [PITH_FULL_IMAGE:figures/full_fig_p045_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: The trapezoidal domain Dε at (10.30). The darker shaded parallelograms are the regions D− ε , D+ ε defined at (10.35). for all ε > 0 sufficiently small. This proves that the integral e2(t) in (10.33) approaches zero as ε → 0, for all t ∈ [ε −α, ε−γ ]. A similar argument shows that the same holds for the first integral e1(t). 7. We now argue as in the proof of Lemma 6.1. Since the domains Dε invade the ent… view at source ↗
read the original abstract

Solutions to hyperbolic conservation laws can be approximated in many different ways: by vanishing viscosity, relaxations, discrete or semi-discrete numerical schemes, approximation with a nonlocal flux, etc$\ldots$ For some of these methods, general ${\bf L}^1$ convergence results are available. Aim of this paper is to understand the local behavior of these approximations, in a neighborhood of point where the hyperbolic solution has a singularity. Specifically: a point along a shock, or where two shocks interact, or where a new shock is formed. Given a sequence of $\epsilon$-approximate solutions, a general expectation is that, by a suitable local rescaling of coordinates, as $\epsilon\to 0$ a well defined limit is obtained. This corresponds to a specific ``eternal solution" (globally defined both in space and in time) to the approximating equation. Precise results this direction are here given, in the case of vanishing viscosity.

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

1 major / 1 minor

Summary. The manuscript studies local asymptotic behavior of vanishing-viscosity approximations to hyperbolic conservation laws near singularities (points along shocks, shock interactions, or shock formation). It claims that, given a sequence of ε-approximate solutions, a suitable local rescaling yields convergence as ε→0 to a specific eternal solution (globally defined in space and time) of the viscous equation, with precise results provided for the vanishing-viscosity case.

Significance. If the convergence statements hold with the claimed precision, the results supply detailed local patterns for how viscous regularizations resolve singularities, extending general L¹ convergence theory. This could inform stability analysis and the design of numerical schemes that capture shock structures accurately. The focus on three distinct singularity classes and the identification of corresponding eternal solutions would be a concrete advance if the proofs are complete.

major comments (1)
  1. Abstract and main convergence statements: the headline claim that 'a well defined limit is obtained' (i.e., the entire sequence converges) is stronger than what standard compactness (BV bounds or compensated compactness on the rescaled parabolic problem) typically delivers, which is only that every subsequential limit satisfies the eternal viscous equation. Establishing that all limit points coincide requires a uniqueness or stability result for the rescaled eternal problem in each of the three singularity classes. If the manuscript only classifies possible limits without proving uniqueness (or without a direct, subsequence-free convergence argument), the central assertion does not follow from the compactness step alone and needs explicit treatment in the main theorems.
minor comments (1)
  1. Abstract: 'Precise results this direction' is missing 'in'.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and for identifying this key point about the strength of the convergence claim. We clarify below that the manuscript does establish uniqueness of the eternal solutions, which upgrades subsequential convergence to convergence of the full sequence. We will revise the manuscript to make this explicit.

read point-by-point responses

  1. Referee: Abstract and main convergence statements: the headline claim that 'a well defined limit is obtained' (i.e., the entire sequence converges) is stronger than what standard compactness (BV bounds or compensated compactness on the rescaled parabolic problem) typically delivers, which is only that every subsequential limit satisfies the eternal viscous equation. Establishing that all limit points coincide requires a uniqueness or stability result for the rescaled eternal problem in each of the three singularity classes. If the manuscript only classifies possible limits without proving uniqueness (or without a direct, subsequence-free convergence argument), the central assertion does not follow from the compactness step alone and needs explicit treatment in the main theorems.

    Authors: We agree that compactness alone yields only subsequential limits and that uniqueness (or a direct argument) is required to conclude convergence of the entire sequence. In the manuscript this is achieved: for each of the three singularity classes we prove that any subsequential limit must coincide with a specific eternal solution of the viscous equation (the viscous shock profile for a single shock, the interaction profile for shock collisions, and the formation profile for shock formation). These uniqueness statements are obtained by combining the parabolic maximum principle with precise asymptotic matching to the far-field hyperbolic states in the rescaled variables; see Theorems 2.3, 3.2 and 4.1 together with the uniqueness lemmas in Sections 2.4, 3.3 and 4.2. We will revise the abstract, the statements of the main theorems, and add a short paragraph after the compactness argument to separate the compactness step from the uniqueness step, thereby addressing the referee's request for explicit treatment. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation is a standard convergence proof

full rationale

The paper establishes local convergence of rescaled viscous approximations to specific eternal solutions near singularities of the hyperbolic limit. No quoted step reduces a claimed prediction or limit to a fitted parameter, self-definition, or unverified self-citation chain by construction. Compactness plus identification of limit points is the expected structure for such results; any prior citations by the authors are not shown to be the sole load-bearing justification for uniqueness or existence in the rescaled problems. The derivation is therefore self-contained as an independent mathematical argument.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard domain assumptions from conservation-law theory; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (2)
  • domain assumption Hyperbolic conservation laws admit solutions with isolated shocks, shock interactions, and shock formation points.
    Implicit in the statement that the hyperbolic solution has a singularity of one of the listed types.
  • domain assumption ε-approximate solutions exist and satisfy suitable L1 bounds compatible with local rescaling.
    Stated directly as the given sequence of ε-approximate solutions.

pith-pipeline@v0.9.0 · 5460 in / 1253 out tokens · 76462 ms · 2026-05-09T19:57:08.743105+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

34 extracted references · 34 canonical work pages

  1. [1]

    Preprint 2025

    J.Anderson, S.Chaturvedi and C.Graham, Shock formation in 1D conservation laws II: vanishing viscosity. Preprint 2025. arXiv:2506.17156

    work page arXiv 2025

  2. [2]

    Pure Appl

    S.Bianchini, Hyperbolic limit of the Jin-Xin relaxation model.Comm. Pure Appl. Math.59(2006), 688–753

    work page 2006

  3. [3]

    S.Bianchini and A.Bressan, On a Lyapunov functional relating shortening curves and viscous conservation laws,Nonlinear Analysis T.M.A.51(2002), 649–662

    work page 2002

  4. [4]

    S.Bianchini and A.Bressan, Vanishing viscosity solutions of nonlinear hyperbolic systems,Annals of Mathematics161(2005), 223–342

    work page 2005

  5. [5]

    The one-dimensional Cauchy problem

    A.Bressan,Hyperbolic systems of conservation laws. The one-dimensional Cauchy problem. Oxford University Press, Oxford, 2000

    work page 2000

  6. [6]

    A.Bressan and G.Chen, Generic regularity of conservative solutions to a nonlinear wave equation.Ann. Inst. H. Poincaré34(2017), 335–354

    work page 2017

  7. [7]

    A.Bressan and C.Donadello, On the formation of scalar viscous shocks,Int. J. Dynam. Syst. Differ. Equat.1(2007), 1–11. 55

    work page 2007

  8. [8]

    A.Bressan and C.Donadello, On the convergence of viscous approximations after shock interactions.Discr. Cont. Dynam. Syst.23(2009), 29–48

    work page 2009

  9. [9]

    Dynamical Systems3(1997), 35–58

    A.Bressan and G.Guerra, Shift differentiability of the flow generated by a conser- vation law,Discrete Cont. Dynamical Systems3(1997), 35–58

    work page 1997

  10. [10]

    A.Bressan, T.Huang and F.Yu, Structurally stable singularities for a nonlinear wave equation.Bull. Inst. Math. Acad. Sinica,10(2015), 449–478

    work page 2015

  11. [11]

    A.Bressan and A.Marson, A variational calculus for discontinuous solutions of con- servative systems,Comm. Part. Diff. Equat.20(1995), 1491–1552

    work page 1995

  12. [12]

    Rational Mech

    A.Bressan and W.Shen, On traffic flow with nonlocal flux: a relaxation represen- tation,Arch. Rational Mech. Anal.237(2020), 1213–1236

    work page 2020

  13. [13]

    S.ChaturvediandC.Graham, TheinviscidlimitofviscousBurgersatnondegenerate shock formation.Annals of PDE(2023) 9:1

    work page 2023

  14. [14]

    M.Chiani, D.Dardari and M.K.Simon, New exponential bounds and approxima- tions for the computation of error probability in fading channels.IEEE Transactions on Wireless Communications2(2003) 840–845

    work page 2003

  15. [15]

    J.Choi, C-Y.Jung and H.Lee, On boundary layers for the Burgers equations in a bounded domainComm. Nonlin. Sci. Num. Simul.67(2019), 637–657

    work page 2019

  16. [16]

    Rational Mech

    M.Colombo, G.Crippa, andLV.Spinolo, Onthesingularlocallimitforconservation laws with nonlocal fluxes.Arch. Rational Mech. Anal.233(2019), 1131–1167

    work page 2019

  17. [17]

    C.Dafermos,Hyperbolic Conservation Laws in Continuum Physics, 4-th ed., Springer-Verlag, Berlin, 2016

    work page 2016

  18. [18]

    P. C. Fife,Mathematical Aspects of Reacting and Diffusing Systems.Lecture Notes in Biomathematics, Springer, 1979

    work page 1979

  19. [19]

    Pure Appl

    H.Freistühler and D.Serre,L1 stability of shock waves in scalar viscous conservation laws.Comm. Pure Appl. Math.51(1998), 291–301

    work page 1998

  20. [20]

    Rational Mech

    J.Goodman and Z.Xin, Viscous limits for piecewise smooth solutions to systems of conservation laws.Arch. Rational Mech. Anal.121(1992), 235–265

    work page 1992

  21. [21]

    Math.15(1975), 65-71

    M.Golubitsky and D.Schaeffer, Stability of shock waves for a single conservation law.Adv. Math.15(1975), 65-71

    work page 1975

  22. [22]

    Dynamical systems - War- wick 1974

    J.Guckenheimer, Solving a single conservation law. In “Dynamical systems - War- wick 1974”, pp. 108–134. Springer Lecture Notes in Math,468

    work page 1974

  23. [23]

    Springer-Verlag, Berlin, 2015

    H.Holden and N.Risebro,Front Tracking for Hyperbolic Conservation Laws.Second Edition. Springer-Verlag, Berlin, 2015

    work page 2015

  24. [24]

    A.Keimer and L.Pflug, On approximation of local conservation laws by nonlocal conservation laws.J. Math. Anal. Appl.475(2019), 1927–1955

    work page 2019

  25. [25]

    USSR Sbornik10(1970), 217–243

    S.Kruzhkov, First order quasilinear equations with several space variables,Math. USSR Sbornik10(1970), 217–243. 56

    work page 1970

  26. [26]

    A.Kurganov and E.Tadmor, Stiff systems of hyperbolic conservation laws: conver- gence and error estimates.SIAM J. Math. Anal.28(1997), 1446–1456

    work page 1997

  27. [27]

    R.J.LeVeque,Numerical Methods for Conservation Laws.Birkhäuser, Basel, 1992

    work page 1992

  28. [28]

    M.Li and Q.Zhang, Generic regularity of conservative solutions to Camassa-Holm type equations.SIAM J. Math. Anal.49(2017), 2920–2949

    work page 2017

  29. [29]

    T.P.Liu, Hyperbolic conservation laws with relaxation.Commun. Math. Phys.108 (1987),153–175

    work page 1987

  30. [30]

    Soc., Providence, RI, 2021

    T.P.Liu,Shock Waves.American Math. Soc., Providence, RI, 2021

    work page 2021

  31. [31]

    in Math.11(1973), 368–386

    D.Schaeffer, A regularity theorem for conservation laws,Adv. in Math.11(1973), 368–386

    work page 1973

  32. [32]

    Differential Equations244(2008), 1692–1711

    W.Shen and Z.Xu, Vanishing viscosity approximations to hyperbolic conservation laws,J. Differential Equations244(2008), 1692–1711

    work page 2008

  33. [33]

    A.Van Harten and R.R.Van Hassel, A quasilinear, singular perturbation problem of hyperbolic type,SIAM J. Math. Anal.16(1985), 1258–1267

    work page 1985

  34. [34]

    Differential Equations155(1999), 89–132

    Wen-An Yong, Singular perturbations of first-order hyperbolic systems with stiff source terms,J. Differential Equations155(1999), 89–132. 57

    work page 1999