Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

✓ 31 results · all verified · 25 also independently AI-judged
Every result on this page is machine-checked by a proof checker and read in full and owner-audited; the judge is an additional, independent cross-model AI review of the proofs. The 6 not AI-judged were verified by owner audit (typically over a confirmed judge false positive), not failures.

Scalar Conservation Laws and Entropy Solutions

1 · Prerequisites

2 · Summary

This page develops the theory of scalar conservation laws ut+div⁡xf(u)=0 from the definitions through existence, uniqueness and wave structure, with the vanishing-viscosity method as its backbone. It fixes the flux and the distributional weak formulation, verifies that classical solutions are weak solutions, derives the constancy of u along characteristics together with the gradient-catastrophe formula for one-dimensional solutions, and sets up piecewise C1 shocks with one-sided traces and the space--time Rankine--Hugoniot condition. Non-uniqueness of the weak formulation is demonstrated before any selection principle is introduced. Convex entropy--entropy-flux pairs and the viscous entropy-dissipation identity lead to the Kruzhkov notion: bounded solutions satisfying the entropy inequalities for ηk(s)=∣s−k∣, with the strong local L1 initial trace.

The constructive half of the page proves the global classical solvability of the viscous Cauchy problem with smooth data, its uniform L∞, mass and energy bounds, and the uniform contraction of spatial translates; the doubling-variables (Kato) inequality for two entropy solutions yields the local L1 contraction, from which uniqueness, comparison, order preservation and finite propagation follow. Uniform bounds and the translate and time moduli give local precompactness of vanishing-viscosity families, the global L1 contraction, and finally the existence of a bounded Kruzhkov entropy solution for L1∩L∞ data by passing the weak equation and the viscous entropy balance to the limit. The structure theory then characterises admissible jumps by the flux-chord inequality, solves the Riemann problem for strictly convex fluxes with its shock and centred-rarefaction profiles, proves Oleinik's one-sided estimate ∂xu≤(κt)−1 as an equivalent entropy condition under uniform convexity on the state range, derives the Lax shock inequalities, and establishes the one-dimensional Hamilton--Jacobi correspondence between entropy solutions and primitives of viscosity solutions.

The closing items record quantitative consequences: mass conservation for compactly supported integrable data, invariance of the entropy inequality under additive constants in the entropy flux, the L∞ maximum bound, the entropy solution semigroup on L1∩L∞ with its extension to L1 for globally Lipschitz flux with f(0)=0, and strong L1 continuity of the orbits. Countable Choice and Dependent Choice are declared where the analytic interfaces require them, including inheritance by consumers. The Hamilton--Jacobi correspondence is proved through viscous primitives and localized Hopf--Lax formulas; no viscosity/entropy equivalence is silently imported.

3 · Logical flowchart

4 · Definitions, theorems and proofs

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

Scalar conservation laws, fluxes and Cauchy data

Definition

Let n≥1 and let Π=O×(a,b)⊆Rn×R be an open space--time cylinder. The equation ut+div⁡xf(u)=0in Π is the scalar conservation law in conservation form, where the unknown u ⁣:Π→R is the conserved quantity and the flux is a map f ⁣:R→Rn. The classical expression is used when u and f∘u are C1; the distributional expression is used whenever u,f(u)∈Lloc1(Π). In particular, if u is bounded measurable and f is continuous, then f(u) is locally integrable and its distributional divergence is defined (A locally integrable function on Rn, The space Lp(μ) as the quotient by null functions).

A Cauchy problem is posed separately on ΠT=Rn×(0,T) and has initial datum u0∈L∞(Rn)∩Lloc1(Rn), understood as an equivalence class; the weak formulation records it in the initial boundary term and the entropy formulation uses a strong local L1 trace. If f∈C1(R;Rn) and u∈C1, the chain rule gives the equivalent quasilinear equation ut+f′(u)⋅∇u=0 (The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c), Ck maps and multi-index derivative notation in Euclidean space, The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder, Directional derivatives and partial derivatives of a map U⊆Rm→Rn, Divergence and curl of a C1 vector field). This equivalence is not asserted for discontinuous u: the product of f′(u) with a distributional gradient is not generally defined, while div⁡xf(u) is defined distributionally whenever f(u)∈Lloc1.

The one-dimensional case is ut+f(u)x=0 with scalar flux f ⁣:R→R. Whenever the Riemann theory, the Rankine--Hugoniot condition or the characteristic formula below is invoked, f is assumed at least C1 (and C2 where f′ or f′′ is differentiated); adding a constant vector to f does not change the equation because its divergence is zero.

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

Distributional weak solutions of the Cauchy problem

Definition

Let n≥1, T>0, ΠT=Rn×(0,T), let f ⁣:R→Rn be continuous and let u0∈L∞(Rn)∩Lloc1(Rn). A bounded measurable u ⁣:ΠT→R is a distributional weak solution of the Cauchy problem ut+div⁡xf(u)=0, u(⋅,0)=u0, if for every φ∈Cc∞(Rn×(−∞,T)) --- test functions whose support may meet the initial plane t=0 --- one has ∫ΠT(u φt+f(u)⋅∇xφ) dx dt+∫Rnu0(x)φ(x,0) dx=0. Equivalently ut+div⁡xf(u)=0 in D′(ΠT), together with the displayed initial boundary term (Distribution, Distributional derivative, Test function space d of an open set). The integrals are absolutely convergent: u is bounded, f(u) is bounded on the bounded range of u, and φ has compact support, so the pairings are Lloc1 pairings and are representative-independent (A locally integrable function on Rn, The space Lp(μ) as the quotient by null functions; the product-space identities are those of Fubini's theorem for L^1 functions on a sigma-finite product and Tonelli's theorem for nonnegative measurable functions on a sigma-finite product). The condition is an equality in the Lloc1 classes of u and f(u); it involves no pointwise assignment of u on {t=0}, and the initial datum enters only through the boundary term (Scalar conservation laws, fluxes and Cauchy data). For the entropy formulation of this page the initial condition is instead imposed as a strong local L1 trace (Kruzhkov entropy solutions).

PropositionStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Classical solutions are distributional weak solutions, and conversely

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the injectivity interface below. Let n≥1, T>0, f∈C1(R;Rn) and u0∈L∞(Rn)∩Lloc1(Rn).

(i) If u∈C1(ΠT)∩C0([0,T];Lloc1(Rn))∩L∞(ΠT) satisfies ut+div⁡xf(u)=0 pointwise in ΠT and has initial trace u(⋅,t)→u0 in Lloc1 as t↓0, then u is a distributional weak solution in the sense of Distributional weak solutions of the Cauchy problem.

(ii) Conversely, if u is a bounded distributional weak solution, u∈C1(ΠT), and u(⋅,t) has an Lloc1-continuous trace uˉ0 as t↓0, then ut+div⁡xf(u)=0 pointwise in ΠT and uˉ0=u0 as Lloc1 equivalence classes.

The initial-trace conclusion is an almost-everywhere class equality; no pointwise representative on t=0 is asserted (Scalar conservation laws, fluxes and Cauchy data).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, f∈C1(R;Rn), u0∈L∞∩Lloc1, a classical solution u∈C1(ΠT) satisfying ut+div⁡xf(u)=0 pointwise with an Lloc1 initial trace (i), and, in (ii), a bounded distributional weak solution u∈C1(ΠT) with an Lloc1-continuous initial trace uˉ0.

[F2]

Integration by parts on a box follows coordinatewise from The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a) and Fubini's theorem for L^1 functions on a sigma-finite product applied to the C1 products uφ and fi(u)φ. Compact support kills spatial and terminal faces. No surface divergence theorem is used.

[F3]

A continuous function on an open set whose integral against every compactly supported smooth test function vanishes is identically zero there; testing against a translate of a fixed nonzero smooth compactly supported bump gives a nonzero pairing where the function does not vanish (Explicit compactly supported smooth cutoffs).

[F4]

Under Countable Choice, two locally integrable functions with equal pairings against every smooth compactly supported test represent the same almost-everywhere class (Locally integrable functions embed in distributions).

Proof

technique · direct
1.1F1F2

Fix φ∈Cc∞(Rn×(−∞,T)) and choose QR=(−R,R)n and τ with supp⁡φ⊂QR×(−∞,τ), τ<T. Multiplying the pointwise equation by φ and integrating over QR×(δ,τ), 0<δ<τ, [F2] gives 0=∫QRuφ∣δτ dx−∫δτ ⁣ ⁣∫QRu φt−∫δτ ⁣ ⁣∫QRf(u)⋅∇xφ, since φ vanishes on the lateral and terminal faces.

1.2F1F3given

For (ii), test the weak identity with φ∈Cc∞(ΠT) (so φ=0 near t=0): integration by parts over the support of φ gives 0=∫ΠT(uφt+f(u)⋅∇xφ)=−⟨ut+div⁡xf(u),φ⟩, where the residual R:=ut+div⁡xf(u) is continuous by [F1]. By [F3] applied on the open set ΠT, R≡0, so the equation holds pointwise.

2.1F2givenstep 1.1

In step 1.1 the terminal term vanishes and φ(⋅,δ)→φ(⋅,0) uniformly on the compact spatial support while u(⋅,δ)→u0 in Lloc1; letting δ↓0 gives ∫ΠT(uφt+f(u)⋅∇xφ) dx dt+∫Rnu0φ(⋅,0) dx=0, which is the weak formulation for this test function. As φ was arbitrary, (i) holds.

3.1F2F4givenstep 1.2step 2.1∎

Identification of the trace. Fix ψ∈Cc∞(Rn) and β∈Cc∞((−∞,T)) with β(0)=1. The pointwise equation from step 1.2, integrated on a box times (δ,τ) containing the positive-time support of ψβ, gives the calculation of steps 1.1 and 2.1 with trace uˉ0. Hence ∫ΠT(uφt+f(u)⋅∇φ)+∫uˉ0ψ=0 for φ=ψβ. Subtract the given weak identity, whose bottom term is ∫u0ψ, to get ∫(uˉ0−u0)ψ=0. By [F4], uˉ0=u0 almost everywhere.

PropositionStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Characteristics and the Riccati equation for the spatial derivative

Statement

Let T>0, f∈C2(R), and let u∈C1(R×(0,T)) be a classical solution of ut+f(u)x=0 (Scalar conservation laws, fluxes and Cauchy data, Ck maps and multi-index derivative notation in Euclidean space).

(i) Along every characteristic x(t) solving x˙=f′(u(x(t),t)), the value u(x(t),t) is constant.

(ii) If in addition u∈C2(R×(0,T)), then p(t)=ux(x(t),t) satisfies along each characteristic p˙(t)=−f′′(u(x(t),t))p(t)2. Consequently, if f′′≥0 on the range of u, then p is nonincreasing along characteristics. More precisely, fix t0∈(0,T) and a characteristic through (x0,t0), and set p0=ux(x0,t0). If p0<0 and f′′(u(x0,t0))≥κ>0, then, as long as the classical solution exists along that characteristic, p(t)=p01+f′′(u(x0,t0))p0(t−t0). If the solution exists along this characteristic up to that time, its derivative tends to −∞ at t∗=t0+1/(f′′(u(x0,t0))∣p0∣)≤t0+1/(κ∣p0∣); hence a C2 solution cannot persist through t∗ along this characteristic (The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c), Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, Linear transport equations and their characteristic flow).

Facts & Assumptions

Given: T>0, f∈C2(R), a classical solution u∈C1(R×(0,T)) of ut+f(u)x=0, together with a characteristic t↦x(t) solving x˙=f′(u(x(t),t)); in part (ii) additionally u∈C2(R×(0,T)) and p=ux.

[F1]

A classical solution satisfies ut+(f(u))x=0 pointwise; since f∈C2⊂C1 and u∈C1, the composition f(u) is C1 with (f(u))x=f′(u)ux (Scalar conservation laws, fluxes and Cauchy data, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c)).

[F2]

For the transport field a(x,t):=f′(u(x,t)), which is C1 because f′ is C1 and u is C1, the solution u satisfies ut+a⋅ux=0, and the curve x(⋅) is a characteristic of that transport equation, so ddtu(x(t),t)=ut(x(t),t)+ux(x(t),t)x˙(t) (Linear transport equations and their characteristic flow, A transport equation restricts to a linear ODE along each characteristic).

[F3]

Sums, scalar multiples and products of differentiable functions are differentiable, with the usual sum and product rules; for u∈C2 and p=ux, the functions p and f′(u)p are C1, while pt, px, and f′′(u) are continuous, so pt+f′(u)px+f′′(u)p2 is continuous (Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, Ck maps and multi-index derivative notation in Euclidean space, The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder). Equality utx=uxt is supplied by Clairaut--Schwarz theorem for continuous second partial derivatives.

[F4]

A continuous field locally Lipschitz in its state has a unique local and maximal ODE solution (Picard-Lindelöf local existence and uniqueness for first-order systems, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval). The field f′(u(x,t)) is locally state-Lipschitz because its x derivative f′′(u)ux is continuous and bounded on compact boxes.

Proof

technique · direct
1.1F1F2

Part (i). With a(x,t)=f′(u(x,t)), [F2] gives ddtu(x(t),t)=ut+uxx˙ along the characteristic. Substituting the characteristic ODE x˙=f′(u(x(t),t)) and then the pointwise equation of [F1] gives ddtu(x(t),t)=ut+f′(u)ux=ut+(f(u))x=0. Hence t↦u(x(t),t) is constant along every characteristic.

1.2F1F3

Part (ii): differentiation in x. Assume now u∈C2, so that p=ux is C1 and f′(u)p is C1 by [F1] and [F3]. Differentiating the pointwise equation ut+(f(u))x=0 in x gives pt+(f′(u)p)x=0, and the product and chain rules give (f′(u)p)x=f′′(u)uxp+f′(u)px=f′′(u)p2+f′(u)px, so that pt+f′(u)px+f′′(u)p2=0.

2.1F2step 1.2

The Riccati equation along characteristics. Along the characteristic of step 1.1, [F2] applied to the C1 function p gives p˙=pt+pxx˙=pt+f′(u)px. By step 1.2 this equals −f′′(u(x(t),t))p2, which is the asserted Riccati equation.

3.1step 1.1step 2.1

Monotonicity. By step 1.1, u(x(t),t)≡u∗ is constant along the characteristic, so along that curve f′′(u(x(t),t))=f′′(u∗) is a constant c and step 2.1 reads p˙=−cp2. If f′′≥0 on the range of u, then c≥0, so p˙≤0 and p is nonincreasing along the characteristic.

4.1step 3.1F4algebra

Exact Riccati solution. Suppose p0=p(t0)<0 and c:=f′′(u(x0,t0))≥κ>0; by step 3.1 the constant is c=f′′(u∗) along the whole characteristic. The scalar ODE p˙=−cp2 has a locally Lipschitz right-hand side, so uniqueness in [F4] and the zero solution imply that p cannot reach zero on its interval of existence. Since p(t0)<0, it remains negative there. Hence one has ddt(1p)=−p˙/p2=c by step 3.1, hence 1p(t)=1p0+c(t−t0) and therefore p(t)=p01+cp0(t−t0)=p01−c∣p0∣(t−t0).

5.1step 1.1step 4.1F3F4∎

Blow-up and the persistence bound. Since u(x(t),t)≡u∗ by step 1.1, the speed x˙=f′(u∗) is constant, so the characteristic is the straight line x(t)=x0+f′(u∗)(t−t0), on its maximal interval. If that interval ended at an interior time t1∈(0,T), the straight line would have a finite endpoint x1; continuity would give u(x1,t1)=u∗, and [F4] would extend the characteristic with initial condition x(t1)=x1, contradicting maximality. Thus its interval is (0,T). The denominator in step 4.1 is positive exactly for t<t∗=t0+1/(c∣p0∣)≤t0+1/(κ∣p0∣) and tends to 0 as t↑t∗, while p0<0, so p(t)→−∞. Were a C2 solution defined on R×(0,T) with T>t∗, then p=ux would be continuous, hence finite, on the compact rectangle [−A,A]×[t0,12(t∗+T)], with A>1+∣x0∣+∣f′(u∗)∣(T−t0), and along the straight characteristic it would equal the explicit solution of step 4.1, which is unbounded on that interval near t∗: a contradiction. Hence the C2 solution cannot persist through t∗ along this characteristic.

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

Piecewise smooth shocks and one-sided traces

Definition

Let n≥1, T>0, and f∈C1(R;Rn) (Scalar conservation laws, fluxes and Cauchy data). Let Γ be a C1 hypersurface in ΠT=Rn×(0,T) with a two-sided open neighbourhood U⊂ΠT, so that U∖Γ=U−∪˙U+ for disjoint open sides U± (Ck Euclidean maps and diffeomorphisms, Interior, closure, boundary, limit point, isolated point and dense subset of a metric space, The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder). Fix the unit normal ν=(νt,νx) on Γ oriented from U− toward U+, with νt2+∣νx∣2=1.

A piecewise C1 weak solution near Γ is a distributional weak solution u of ut+div⁡xf(u)=0 on ΠT (Distributional weak solutions of the Cauchy problem) such that at each ζ∈Γ there is a neighbourhood Vζ⊂U and functions u~ζ±∈C1(Vζ) representing u on Vζ∩U±, respectively. The one-sided traces at ζ are u±(ζ):=u~ζ±(ζ). They are independent of the chosen local extensions, since continuous extensions agreeing almost everywhere on an open side agree throughout that side and at its interface points. Thus genuine traces are specified at every point of Γ⊂U; arbitrary representative values on Γ do not affect the weak-solution class or these traces. Write [u](ζ)=u+(ζ)−u−(ζ),[f](ζ)=f(u+(ζ))−f(u−(ζ)). A point is a shock point when [u](ζ)≠0.

In one space dimension, this includes a C1 graph x=s(t), with sides x<s(t) and x>s(t); the definition is not restricted to that case, and the existence of the strong one-sided traces is part of the piecewise-smooth hypothesis, not a conclusion for arbitrary L∞ weak solutions.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The Rankine--Hugoniot jump condition in space--time normal form

Statement

Let n≥1, T>0, f∈C1(R;Rn) (Scalar conservation laws, fluxes and Cauchy data), and let u be a piecewise C1 distributional weak solution (Distributional weak solutions of the Cauchy problem) with two-sided C1 interface Γ and traces as in Piecewise smooth shocks and one-sided traces. Orient the unit space--time normal ν=(νt,νx) from the minus side to the plus side. Then at every ζ∈Γ⊂U, [u](ζ)νt(ζ)+[f](ζ)⋅νx(ζ)=0, where [u]=u+−u− and [f]=f(u+)−f(u−).

In one space dimension, for a graph x=s(t) with minus side x<s(t) and plus side x>s(t), ν=(−s′(t),1)/1+s′(t)2, so the condition is s′(t)[u](t)=[f](t) at every graph point. If [u]≠0 in one dimension, then s′=[f]/[u]; if [u]=0, then [f]=0 and the relation is 0=0, with no speed constraint.

Facts & Assumptions

Given: n≥1, a piecewise C1 weak solution u with interface Γ⊂U and traces u±, a point ζ0∈Γ, and one-sided C1 local extensions u~± on a neighbourhood V⊂U of ζ0.

[F1]

The weak identity reads ∫ΠT(uφt+f(u)⋅∇xφ)=0 for every φ∈Cc∞(ΠT), i.e. div⁡t,xF(u)=0 in distributions, where F(u)=(u,f(u)) (Distributional weak solutions of the Cauchy problem, Scalar conservation laws, fluxes and Cauchy data).

[F2]

Locally about a point of a C1 hypersurface, after permuting coordinates, a patch is a graph zk=γ(z′) over the remaining n coordinates z′, with γ∈C1; the unnormalised normal N=ek−∑j≠k(∂jγ)ej points from the region below the graph to the region above it, and the unit normal of [F1]'s orientation is ν=σN/∣N∣ with σ=1 if the minus side lies below the graph and σ=−1 otherwise; also [F]=([u],[f]), and a continuous function vanishing against all nonnegative smooth bumps on an open set vanishes there (Piecewise smooth shocks and one-sided traces, Explicit compactly supported smooth cutoffs).

Proof

technique · direct
1.1F1F2

Side extensions solve the equation classically. On each side of Γ the function u agrees with a C1 extension u~±; testing away from t=0 and against bumps supported in a single side, the weak identity [F1] shows that div⁡t,xF(u~±) vanishes as a distribution on that side. Since F(u~±) is C1, its divergence is continuous, and by [F2] it vanishes pointwise; consequently, for smooth compactly supported φ supported in the side, div⁡t,x(F(u±)φ)=F(u±)⋅Dφ.

2.1F3step 1.1

Graph computation on one side. After a permutation of coordinates, write the graph locally as zk=γ(z′) and take φ supported in a box B′×(a,b) in which the graph stays in (a,b). With Gj=Fj(uL)φ on the lower side, [F3] gives ∫lowerF(uL)⋅Dφ=∫lowerdiv⁡(F(uL)φ)=∫B′GkL(z′,γ(z′)) dz′−∑j≠k∫B′GjL(z′,γ(z′))∂jγ(z′) dz′=∫B′GL(z′,γ(z′))⋅N(z′) dz′, because the k-derivative integrates to the trace at the graph and each tangential derivative of the moving-endpoint integral Aj(z′)=∫aγ(z′)Gj(z′,r) dr contributes −∂jγ Gj at the graph, the integral of ∂jAj vanishing by compact support.

3.1F2step 1.1step 2.1

Upper side and the interface term. The same computation on the upper side, with the graph as its lower boundary, gives ∫upperF(uR)⋅Dφ=−∫B′GR(z′,γ(z′))⋅N(z′) dz′; adding with step 2.1, and noting that F(u+)−F(u−)=([u],[f]) with u± the traces from the plus and minus sides, the weak identity becomes 0=∫ΠTF(u)⋅Dφ=−σ∫B′φ(z′,γ(z′)) [F](z′,γ(z′))⋅N(z′) dz′ for every φ supported in the box.

4.1F2step 3.1

Continuity and vanishing of the bracket. The function z′↦[F](z′,γ(z′))⋅N(z′) is continuous, being a composition of continuous data; if it were nonzero at the point corresponding to ζ0, it would keep one sign on a smaller patch, and a nonnegative smooth bump supported there and positive at ζ0 would make the integral of step 3.1 nonzero, a contradiction. Hence [F]⋅N=0 at ζ0.

5.1step 4.1F2∎

Normal form and the one-dimensional case. Since ν=σN/∣N∣ with σ=±1 by [F2], 0=[F]⋅N=[F]⋅ν σ∣N∣, and ∣N∣=∣ek−∑j≠k∂jγ ej∣>0, so [F]⋅ν=[u]νt+[f]⋅νx=0 at every point of Γ. For a one-dimensional graph x=s(t) with minus side x<s(t), the graph function is γ(t)=s(t), so N=(−s′(t),1) and ν=(−s′,1)/1+s′2; the condition becomes (−s′[u]+[f])/1+s′2=0, that is, s′[u]=[f]. If [u]≠0 this determines s′=[f]/[u], and if [u]=0 the relation reads 0=[f], so [f]=0 and no speed is constrained.

PropositionStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Distributional weak solutions of the Cauchy problem are not unique

Statement

Take n=1, f(u)=u2 and u0≡0. For every δ>0 define, for t>0, uδ(t,x)={0,x<−δt,−δ,−δt<x<0,+δ,0<x<δt,0,x>δt. This is a bounded distributional weak solution of ut+(u2)x=0 on ΠT with strong local L1 initial trace u0=0, and it is not identically zero. At its three jumps the speeds s=[u2]/[u] are, respectively, −δ, 0 and +δ, so each jump coefficient [u2]−s[u] in the weak equation vanishes. The middle stationary jump violates the Kruzhkov entropy inequality for k=0: with η0(u)=∣u∣ and q0(u)=sgn⁡(u)u2 its entropy production is [q0]−0⋅[η0]=δ2−(−δ2)=2δ2>0, whereas the entropy inequality requires this coefficient to be nonpositive (Kruzhkov entropy solutions). Hence uδ is a weak solution but not an entropy solution, and the zero solution is a distinct weak solution with the same initial data: distributional weak solutions are not unique (Distributional weak solutions of the Cauchy problem, Piecewise smooth shocks and one-sided traces).

Facts & Assumptions

Given: n=1, f(u)=u2, u0≡0, T>0, δ>0, and the piecewise constant function uδ above, whose jump rays are Γ−={x=−δt}, Γ0={x=0} and Γ+={x=δt} in ΠT.

[F1]

On each of the four regions the function is constant and f is smooth, so uδ solves the equation classically there; across a jump ray x=s(t) of a piecewise C1 weak solution of ut+(u2)x=0, the distributional identity holds iff the jump coefficient [u2]−s [u] vanishes, where [⋅] denotes the right minus left trace across the ray (Piecewise smooth shocks and one-sided traces, Distributional weak solutions of the Cauchy problem, Scalar conservation laws, fluxes and Cauchy data).

[F2]

The Kruzhkov entropy pair for k=0 is η0(u)=∣u∣, q0(u)=sgn⁡(u)u2 with sgn⁡(0)=0; an entropy solution must satisfy ∂tη0(u)+∂xq0(u)≤0 in D′(ΠT), so across a jump ray the entropy production coefficient [q0]−s[η0] must be nonpositive (Kruzhkov entropy solutions).

[F3]

Basic computation with the explicit states and speeds: for the ray x=−δt the left state is 0, the right state is −δ, and the rightward speed is s=−δ; for x=0 the states are −δ (left) and +δ (right) with s=0; for x=δt the states are +δ (left) and 0 (right) with s=δ; directly [u2]−s[u]=0 in all three cases.

Proof

technique · direct
1.1F1F3

Weak equation. The profile is constant on its four regions. At x=−δt the right-minus-left jumps are [u]=−δ, [u2]=δ2, and s=−δ, so [u2]−s[u]=δ2−(−δ)(−δ)=0. At x=0, [u]=2δ, [u2]=0, and s=0. At x=δt, [u]=−δ, [u2]=−δ2, and s=δ, again giving zero. To verify the distributional equation, integrate uφt+u2φx in each region using the moving-endpoint FTC formula: each interface contributes (s[u]−[u2])∫φ(t,st) dt, which vanishes.

1.2F1given

Initial trace. For every compact K⊆R and 0<t<δ0, the set where uδ(t,⋅) differs from 0 is contained in [−δt,δt], so ∫K∣uδ(t,x)∣ dx≤2δ2t→0 as t↓0; hence uδ has the strong local L1 initial trace 0. The function is bounded, hence a distributional weak solution of the Cauchy problem with datum u0≡0 in the sense of [F1].

2.1F2step 1.2

Failure of the entropy condition. At the middle ray x=0 the left and right states are −δ and +δ. The entropy production coefficient is [q0]−s[η0] with [η0]=∣+δ∣−∣−δ∣=0, s=0, and [q0]=q0(+δ)−q0(−δ)=δ2−(−δ2)=2δ2>0. By [F2] the required entropy inequality fails: the distribution ∂tη0(uδ)+∂xq0(uδ) carries the positive coefficient 2δ2 on the ray x=0.

3.1step 1.1step 1.2step 2.1∎

Non-uniqueness. By steps 1.1 and 1.2, both uδ and the zero function are bounded distributional weak solutions of the same Cauchy problem with initial datum u0≡0; they differ on a set of positive measure for every δ>0. By step 2.1, uδ is not a Kruzhkov entropy solution, so the non-uniqueness occurs strictly within the class of distributional weak solutions.

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

Convex entropy--entropy flux pairs

Definition

Let n≥1, let f ⁣:R→Rn be C1, and let η ⁣:R→R be a finite convex locally Lipschitz function (Convex and strictly convex functions on Euclidean convex sets). An entropy flux is specified by the coordinatewise formula qi(s)=Ci+∫0sη′(r)fi′(r) dr. One may use the left derivative of η in this formula. Convex secant inequalities show that it is bounded and nondecreasing on compact intervals, and that it equals the ordinary derivative except at at most countably many points: assign a distinct rational to each nonempty gap between left and right slopes. A bounded monotone function is Riemann integrable, since the difference of upper and lower sums on an equal mesh is at most the mesh size times its total increase. The same holds after multiplication by the continuous fi′: uniform continuity controls the additional oscillation in the product sums. Thus the displayed integrals exist and give locally Lipschitz q. The convex secant bounds also squeeze the telescoping sum η(b)−η(a) between the left and right derivative sums, proving η(b)−η(a)=∫abη′ without a choice principle. At every continuity point of η′, averaging the integrand over a shrinking interval gives q′(s)=η′(s)f′(s); the exceptional points are countable. This includes nonsmooth entropies such as ∣s−k∣. The constants Ci are the only normalization freedom in this construction. For η∈C1, ordinary FTC makes q∈C1 (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G, The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)); second-derivative identities require η∈C2 with η′′≥0. Conversely, under Countable Choice and Dependent Choice, any locally Lipschitz q satisfying q′=η′f′ almost everywhere has this formula, by the fundamental theorem for absolutely continuous functions (Fundamental theorem of calculus for absolutely continuous functions, The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain). This converse analytic interface is separate from defining the explicit integral pairs used here.

Such a pair (η,q) is an entropy--entropy flux pair. For bounded measurable u its entropy inequality is ∂tη(u)+div⁡xq(u)≤0in D′(ΠT), equivalently ∫ΠT(η(u)φt+q(u)⋅∇φ)≥0 for every nonnegative φ∈Cc∞(ΠT) (Distribution, Distributional derivative). All compositions are locally integrable on this bounded range. Adding a constant vector to q leaves the inequality unchanged. The weak conservation law implies equality for the affine pairs (η,q)=(s,f(s)) and (−s,−f(s)); conversely their two entropy inequalities together imply that equality. The inequality for (s,f(s)) alone does not imply the weak equation (Scalar conservation laws, fluxes and Cauchy data).

PropositionStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The viscous entropy dissipation identity

Statement

Let n≥1, ε>0, f∈C1(R;Rn) and let uε∈C1,2(ΠT) be a classical solution of the viscous conservation law utε+div⁡xf(uε)=εΔuεin ΠT (Scalar conservation laws, fluxes and Cauchy data, The Laplacian of a C2 function and of a C2 vector field). For every convex η∈C2(R) and entropy flux q with q′(s)=η′(s)f′(s), the pointwise viscous entropy balance is η(uε)t+div⁡xq(uε)=εΔxη(uε)−εη′′(uε)∣∇xuε∣2≤εΔxη(uε). The nonpositive term is the entropy dissipation; for fixed ε>0 this is a balance with diffusion, not the first-order entropy inequality (Convex entropy--entropy flux pairs, Divergence and curl of a C1 vector field).

Facts & Assumptions

Given: n≥1, ε>0, f∈C1(R;Rn), a classical solution uε∈C1,2(ΠT) of the viscous conservation law, and a convex η∈C2(R) with entropy flux q, q′=η′f′.

Proof

technique · direct
1.1F1F2

By [F1] the two left-hand terms are η(uε)t=η′(uε)utε and div⁡xq(uε)=η′(uε)f′(uε)⋅∇xuε, and by [F2] the diffusion term is Δxη(uε)=η′(uε)Δuε+η′′(uε)∣∇xuε∣2.

2.1step 1.1F1algebra

Multiplying the viscous equation pointwise by η′(uε) and adding the second identity of step 1.1 gives η(uε)t+div⁡xq(uε)=η′(uε)(utε+div⁡xf(uε)−εΔuε)+εΔxη(uε)−εη′′(uε)∣∇xuε∣2=εΔxη(uε)−εη′′(uε)∣∇xuε∣2, which is the asserted balance.

3.1step 2.1F2∎

The convexity assumption gives η′′≥0, so the dissipation term −εη′′(uε)∣∇xuε∣2 is nonpositive pointwise and the balance implies the stated inequality; the term εΔxη(uε) may change sign and cannot be dropped pointwise for fixed ε>0.

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

Kruzhkov entropy solutions

Definition

Let n≥1, T>0, f∈C1(R;Rn), and ΠT=Rn×(0,T). The unknown and initial datum are equivalence classes [u]∈L∞(ΠT) and [u0]∈L∞(Rn)∩Lloc1(Rn), where equality is Lebesgue-a.e. (The space Lp(μ) as the quotient by null functions, A locally integrable function on Rn).

For k∈R define the Kruzhkov entropy pair ηk(s)=∣s−k∣,qk(s)=sgn⁡(s−k)(f(s)−f(k)), with sgn⁡(0)=0 (Absolute value in an ordered field). Then ηk is convex and qk′=ηk′f′ almost everywhere, so (ηk,qk) is a locally Lipschitz convex entropy--entropy flux pair (Convex entropy--entropy flux pairs).

The class [u] is a Kruzhkov entropy solution of ut+div⁡xf(u)=0 with initial trace [u0] if:

(i) for every k∈R and every nonnegative φ∈Cc∞(ΠT), ∫ΠT(ηk(u) φt+qk(u)⋅∇xφ) dx dt≥0, that is, ∂tηk(u)+div⁡xqk(u)≤0 in D′(ΠT); and

(ii) for every compact K⊆Rn, lim⁡δ↓0ess sup⁡0<t<min⁡{δ,T}∫K∣u(t,x)−u0(x)∣ dx=0.

The integral in (i) is independent of the representative because its integrand is unchanged almost everywhere. For (ii), Fubini's theorem gives locally integrable spatial sections for almost every t (Fubini's theorem for L^1 functions on a sigma-finite product); the displayed slice integral is defined for those times and its essential supremum ignores the exceptional null set. If u or u0 is changed on a null set in its respective space, Fubini's theorem shows that the slice-integral function changes only for a null set of times, so the trace condition is well defined on the equivalence classes: this is the strong local L1 initial trace.

Taking k above and below the essential range of u makes ηk equal k−u and u−k, whose t- and x-derivatives cancel the constant terms against compactly supported test functions, so the entropy inequalities imply the weak conservation law of Distributional weak solutions of the Cauchy problem tested against nonnegative test functions, hence by linearity against all test functions.

To recover the Cauchy boundary term, apply the interior weak identity to φ(t,x)ζ(t/δ), where ζ is smooth and nondecreasing, ζ=0 on (−∞,1/2] and ζ=1 on [1,∞). For δ>0 this product is supported away from t=0, so it is an admissible interior test. The term containing ζ′/δ converges to ∫u0(x)φ(0,x) dx by (ii), while the other terms converge on the compact support. This proves the full weak formulation with initial datum u0, rather than only its interior equation.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The viscous scalar Cauchy problem with smooth data has a global classical solution

Statement

Assume Countable Choice (CC) for the heat-kernel and L1 interfaces below. Let n≥1, ε>0, let f ⁣:R→Rn be C2 with f(0)=0, and let u0∈Cc∞(Rn). For every T>0 there is a mild solution u of ut+div⁡xf(u)=εΔuon Rn×(0,T),u(x,0)=u0(x), with u∈C([0,T];Cb(Rn))∩C([0,T];L1(Rn))∩C1,2(Rn×(0,T)). The solution is global in the sense that these solutions are compatible on finite time intervals; for every t∈[0,T], sup⁡x∈Rnu(t,x)≤sup⁡xu0(x),inf⁡x∈Rnu(t,x)≥inf⁡xu0(x). In particular it is bounded and remains in the initial range. No uniqueness beyond the constructed mild solution is asserted.

Facts & Assumptions

Given: Countable Choice, n≥1, ε>0, f∈C2(R;Rn) with f(0)=0, u0∈Cc∞(Rn), and T>0.

[F1]

The heat evolution Ht is the convolution with the heat kernel Γ: Htg=Γ(⋅,t)∗g is defined for g∈Lp with ∥Htg∥p≤∥g∥p, and Γ(⋅,t) has unit mass, is strictly positive, is C∞ with ∂tΓ=ΔΓ, and satisfies the Gaussian derivative bounds ∣DαΓ(x,t)∣≤Cn,αt−(∣α∣+n)/2e−∣x∣2/(8t) (The heat evolution Ht of initial data, Normalisation, parabolic scaling, heat equation and derivative bounds for the heat kernel).

[F2]

For every multi-index α and 1≤p≤q≤∞ there is Cα,p,q with ∥DαHtg∥q≤Cα,p,qt−∣α∣/2−n2(1/p−1/q)∥g∥p for all g∈Lp and t>0; moreover the spatial and time derivatives pass through heat convolution for positive time, and Young's inequality bounds convolutions (Spatial derivative estimates for the heat flow, Lp to Lq smoothing estimate for the heat flow, Spatial and time derivatives pass through heat convolution for positive time, Young's convolution inequality under Countable Choice).

[F3]

Htu0 is bounded and uniformly continuous with ∥Htu0∥∞≤∥u0∥∞, and C([0,τ];Cb(Rn)) with the supremum norm is complete: a uniformly Cauchy sequence of bounded functions converges pointwise in R, its Cauchy bound then gives uniform convergence and boundedness of the limit; uniform limits preserve spatial continuity. Applying the same argument to continuous paths with values in this complete Cb space gives a uniform limit continuous in time (the triangle inequality with one approximating path proves continuity). Thus this path space is complete, the uniform limit of continuous functions is continuous, and the Banach fixed point theorem applies to a contraction of a nonempty complete metric space (The heat Cauchy problem for bounded uniformly continuous data, The uniform limit of continuous real-valued functions on a metric space is continuous, Closed subspaces of complete metric spaces are complete; the converse under countable choice, A contraction of a nonempty complete metric space into itself has exactly one fixed point, the limit of the iterates from any starting point).

[F6]

Write BUC(Rn) for the bounded uniformly continuous functions. The Gaussian kernels form an L1 approximate identity, so Hεtg→g in L1 for g∈L1 by the approximate-identity theorem. For g∈BUC, define ωg(ρ):=sup⁡∣z∣≤ρ∥g(⋅−z)−g∥∞, which tends to zero with ρ. Unit mass and the Gaussian tail give ∥Hεtg−g∥∞≤ωg(ρ)+2∥g∥∞∫∣z∣>ρΓ(z,εt) dz, so Hεtg→g uniformly as t↓0. The Gaussian convolution identity Γt∗Γs=Γt+s follows by completing the square in its integrand and using Gaussian unit mass. The contraction and this identity give strong continuity of t↦Hεtg in both L1 and BUC on [0,τ] (Normalisation, parabolic scaling, heat equation and derivative bounds for the heat kernel, Every L1 approximate identity converges to the identity in Lp for 1≤p<∞). The product-space rearrangements use Fubini's theorem for L^1 functions on a sigma-finite product, kernel time integrations use The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a), and the flux Lipschitz bound uses The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a).

Proof

technique · direct
1.1F1F2F6

Setup and kernel estimates. Put M0:=∥u0∥∞, R:=2M0+1 and L:=sup⁡∣s∣≤R∣f′(s)∣<∞. Since f(0)=0 and f is C1, the mean value theorem gives ∣f(s)∣≤L∣s∣ for ∣s∣≤R. By [F2] there are constants C1,C2 with ∥∇Hσg∥p≤C1σ−1/2∥g∥p for p∈{1,∞}. The Gaussian first moment gives ∥Hσg−g∥∞≤C2σ1/2∥∇g∥∞ for g∈C1∩L∞ with ∇g bounded. Also, ∥∇Γσ(⋅+h)−∇Γσ∥1≤min⁡{2C1σ−1/2,C2∣h∣σ−1} and ∥Γσ(⋅+h)−Γσ∥1≤min⁡{2,C2∣h∣σ−1/2} for every h∈Rn and σ>0. For the scaled kernel, ∥∂σ∇Γεσ∥1≤Cεσ−3/2, so for Δ>0 and σ>0, ∥∇Γε(σ+Δ)−∇Γεσ∥1≤Cεmin⁡{σ−1/2,Δσ−3/2}.

2.1F1F3F6step 1.1

A local mild solution for bounded uniformly continuous L1 data. Fix a starting datum g∈BUC(Rn)∩L1(Rn) with ∥g∥∞≤M0. Choose τ0∈(0,T] with q:=2C1Lτ0/ε≤12 and let Xτ0:=C([0,τ0];Cb(Rn)) with the supremum norm (Closed subspaces of complete metric spaces are complete; the converse under countable choice). For v∈Xτ0 with ∥v∥Xτ0≤R define (Tgv)(t):=Hεtg−∫0t∇Hε(t−s)∗f(v(s)) ds. By [F6], t↦Hεtg is continuous in the supremum norm at t=0; by [F1] the Duhamel term is continuous in Xτ0, since its integrand is norm continuous for s<t and has the integrable majorant C1(ε(t−s))−1/2∥f(v(s))∥∞. The bounds ∥Hεtg∥∞≤M0 and ∥f(v(s))∥∞≤LR give ∥Tgv∥Xτ0≤M0+qR≤R,∥Tgu−Tgv∥Xτ0≤q∥u−v∥Xτ0. Thus Tg is a contraction of the closed radius-R ball; the Banach fixed point theorem gives a unique fixed point u∈Xτ0 satisfying u(t)=Hεtg−∫0t∇Hε(t−s)∗f(u(s))ds. The original datum u0∈Cc∞ satisfies these assumptions.

2.2F1F2step 1.1

Hölder regularity on strips. Fix 0<δ<τ0 and let t,t′∈[δ,τ0]. Since f(u) is bounded with ∥f(u)∥∞≤LR, the mild identity and the kernel bounds of step 1.1 give, for every α∈(0,1), ∥u(t,⋅+h)−u(t,⋅)∥∞≤Cδ∣h∣α,∥u(t,⋅)−u(t′,⋅)∥∞≤Cδ∣t−t′∣1/2. For the spatial estimate, ∥∇Hεtg∥∞≤Cε,δ∥g∥∞ on the strip. The Duhamel term is bounded by LR∫0tmin⁡{Cεσ−1/2,Cε∣h∣σ−1} dσ, with t≤τ0; splitting at σ=∣h∣2 gives at most Cε(∣h∣+∣h∣log⁡+(τ0/∣h∣2))≤Cε,τ0,α∣h∣α for 0<∣h∣≤1, where log⁡+(r):=max⁡{log⁡r,0}. Larger ∣h∣ are covered by boundedness of u. For the temporal estimate assume t=t′+Δ with Δ>0. The initial heat term is Lipschitz in t on [δ,τ0] by the positive-time estimate for ∂tHεtg=εΔHεtg. On the common Duhamel interval, with σ=t′−s, the kernel difference is a time difference and satisfies ∥∇Γε(σ+Δ)−∇Γεσ∥1≤Cεmin⁡{σ−1/2,Δσ−3/2}. Splitting its integral at σ=Δ gives ∫0t′∥∇Γε(σ+Δ)−∇Γεσ∥1 dσ≤CεΔ. The remaining time slab has norm at most LR∫t′tC1(ε(t−s))−1/2ds≤2C1LRΔ/ε. Thus both estimates hold with Cδ depending only on δ,τ0,ε,n,∥g∥∞,LR, and u is continuous in (t,x) for t>0.

3.1F2F5F6step 2.1

L1 continuity and identification of the limits. For the arbitrary starting datum g∈BUC∩L1 of step 2.1, take Picard iterates u(0):=g, u(k+1):=Tgu(k). By [F6], u(0)∈C([0,τ0];L1). If u(k)∈C([0,τ0];L1), the integral defining u(k+1) is defined by norm limits of Riemann sums on truncated intervals s≤t−η. Completeness supplies these limits, and the bound on the omitted interval is at most 2C1L∥u(k)∥CtL1η/ε, which tends to zero. Its integrand has integrable majorant C1(εs)−1/2L∥u(k)∥CtL1, so u(k+1) is L1-continuous and ∥u(k+1)∥CtL1≤∥g∥1+q∥u(k)∥CtL1≤2∥g∥1. Successive differences obey ∥u(k+1)−u(k)∥CtL1≤q∥u(k)−u(k−1)∥CtL1, so the iterates converge in the complete space C([0,τ0];L1) (Riesz-Fischer completeness of Lp for 1≤p≤∞) to some w. They also converge uniformly on [0,τ0]×Rn to the fixed point u of step 2.1. On every bounded ball Bm, this uniform convergence implies convergence to u in C([0,τ0];L1(Bm)), while the global L1 convergence gives convergence to w in the same local space. Uniqueness of the local L1 limit yields u=w a.e. on every Bm, hence a.e. on Rn; thus u∈C([0,τ0];L1) and ∥u(t)∥1≤2∥g∥1.

3.2F1F2step 2.1

The equation in distributions. Testing the mild identity and using Fubini (the bound (t−s)−1/2 is integrable on each finite time triangle) gives ∫Π(uφt+f(u)⋅∇φ)+∫gφ(⋅,0)=−ε∫ΠuΔφ. Indeed the initial heat term pairs as −∫gφ(⋅,0) against φt+εΔφ, and the divergence Duhamel term pairs as −∫f(u)⋅∇φ against the same expression, by integrating ∂tHε(t−s)=εΔHε(t−s) in t. Thus ut+div⁡f(u)=εΔu distributionally, with datum g.

3.3F1F2F4step 1.1step 2.2

Hölder continuity and boundedness of the spatial gradient. Fix 0<δ0<δ<τ0 and use the mild identity restarted at δ0: u(t)=Hε(t−δ0)u(δ0)−div⁡W(t),W(t):=∫δ0tHε(t−s)f(u(s)) ds. By step 2.2, f(u) is spatially Cα and temporally Cα/2 on each positive strip (the C1/2 bound implies the weaker Cα/2 bound). Componentwise, differentiation of the divergence heat potential gives ∂ku(t,x)=∂kHε(t−δ0)u(δ0,x)−∑i∫0t−δ0∫RnDkiΓεσ(z)(fi(u(t−σ,x−z))−fi(u(t,x))) dz dσ, where the frozen value is subtracted using ∫DkiΓεσ=0. Gaussian scaling gives ∥D2Γεσ∥1≤Cεσ−1, ∥∇D2Γεσ∥1≤Cεσ−3/2, and ∥∂σD2Γεσ∥1≤Cεσ−2. With the parabolic Cα,α/2 modulus of f(u), a spatial increment h is bounded after splitting at σ=∣h∣2 by C∫0∣h∣2σ−1+α/2 dσ+C∣h∣∫∣h∣2τ0−δ0σ−3/2+α/2 dσ≤C∣h∣α, and a time increment Δ=∣t−t′∣ is bounded after splitting at σ=Δ by C∫0Δσ−1+α/2 dσ+CΔ∫Δτ0−δ0σ−2+α/2 dσ≤CΔα/2. If either split point exceeds the finite integration horizon, the same bounds follow by increasing C. The heat initial term is smooth with bounded derivatives on t≥δ. Thus for each 0<α<1, [∇u]Cxα(Rn×[δ,τ0))+[∇u]Ctα/2(Rn×[δ,τ0))≤Cδ0,δ. These are seminorm bounds; in particular they establish continuity of ∇u before any absolute bound is used. The absolute gradient bound follows from boundedness of u. If ∇u(t,x)≠0, put e=∇u(t,x)/∣∇u(t,x)∣. Along the unit segment x+se, 0≤s≤1, the fundamental theorem of calculus and the spatial seminorm Mδ:=[∇u]Cxα give ∣∇u(t,x)∣≤∣u(t,x+e)−u(t,x)∣+∫01∣∇u(t,x+se)−∇u(t,x)∣ ds≤2R+Mδ1+α. Thus ∇u is bounded on the positive-time strip. Since f′′ is bounded on [−R,R], f′(u) has the same parabolic Hölder regularity as u there; hence h:=−f′(u)⋅∇u is bounded and belongs to Cα,α/2(Rn×[δ,τ0)).

4.1F1F4F5step 2.1step 3.3

Heat-potential cancellation and C1,2. Fix 0<δ<τ0, use the bounds of step 3.3 on a slightly larger positive strip, and restart the heat equation at δ. Then u(t)=Hε(t−δ)u(δ)+w(t),w(t,x):=∫δtHε(t−s)h(s)(x) ds. For σ>0, ∫RnDijΓεσ(z) dz=0, so cancellation gives Dijw(t,x)=∫0t−δ∫RnDijΓεσ(z)(h(t−σ,x−z)−h(t,x)) dz dσ. The parabolic Hölder bound from step 3.3 yields ∣h(t−σ,x−z)−h(t,x)∣≤C(∣z∣α+σα/2). Gaussian scaling therefore gives ∫Rn∣D2Γεσ(z)∣(∣z∣α+σα/2) dz≤Cε,n,ασ−1+α/2, which is integrable at σ=0. Truncating the heat-potential integral at σ=η>0, differentiating, and letting η↓0 with this integrable dominator shows that Dijw exists and is continuous; the heat-kernel identity also gives wt−εΔw=h in distributions. Since h and Δw are continuous, this identity gives a continuous classical time derivative. The first term Hε(t−δ)u(δ) is smooth for t>δ, so u∈C1,2(Rn×(δ,τ0)) and satisfies ut=εΔu−f′(u)⋅∇u, hence ut+div⁡xf(u)=εΔu pointwise. Since δ>0 is arbitrary, u∈C1,2(Rn×(0,τ0)).

5.1F3F4step 2.1step 4.1

The range bound by a barrier. Let a<b in [0,τ0) and put Ma:=sup⁡xu(a,x). Write B:=f′(u(t,x)), which is bounded on [a,b]×Rn by L. For δ1>0 and λ>0 set Φ(t,x):=δ1eλ(t−a)(1+∣x∣2); then, using steps 4.1 and [F4], (∂t−εΔ+B⋅∇)Φ=δ1eλ(t−a)(λ(1+∣x∣2)−2εn+2B⋅x)>0 for λ large enough, because ∣B∣≤L and −εΔ(1+∣x∣2)=−2εn. The function z:=u−Ma−Φ is negative at t=a and, by step 2.1, ∣u∣≤R while Φ≥δ1(1+∣x∣2), so z is negative on the lateral boundary of every sufficiently large ball intersected with [a,b]; if z had a positive maximum in the closed cylinder, then at that point zt≥0, Dz=0, Δz≤0 (A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value, Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line), contradicting (∂t−εΔ+B⋅∇)z<0 there. Letting δ1↓0 gives u≤Ma on [a,b]; applying the same argument to −u, whose flux is f~(s)=−f(−s) and whose range bound is the same, gives u≥ma:=inf⁡xu(a,x). In particular, with a=0, sup⁡xu(t,x)≤sup⁡xu0 and inf⁡xu(t,x)≥inf⁡xu0 on [0,τ0), the range of u is contained in I0:=[inf⁡u0,sup⁡u0].

6.1F3F6step 2.1step 2.2step 3.1step 4.1step 5.1∎

Global continuation and conclusion. For every t1>0, the positive-strip spatial modulus of step 2.2 makes u(t1,⋅) bounded and uniformly continuous, step 3.1 gives u(t1)∈L1, and step 5.1 places its values in I0. Thus g:=u(t1) is an admissible BUC∩L1 starting datum for the local fixed-point and Picard iteration of steps 2.1 and 3.1. The strong heat continuity required for this restart is [F6]: Gaussian approximate-identity convergence gives Hεtg→g in the supremum norm for BUC data and in L1 for L1 data. With the same global constants R,L, the local solution on [t1,t1+τ0) agrees with u at t1; local uniqueness in step 2.1 makes the pieces agree on overlaps, and ∥u(t)∥∞≤max⁡{∣inf⁡u0∣,∣sup⁡u0∣}≤R throughout. Iterating finitely many times covers [0,T], and the same local uniqueness shows that two solutions obtained with different terminal times T,T′ agree on their common interval. The range bounds of step 5.1 hold on all of [0,T]; L1 continuity holds by step 3.1 and C1,2 regularity on Rn×(0,T) by step 4.1. This is the asserted global mild solution.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Uniform L-infinity, mass and energy bounds for the viscous approximations

Statement

Assume Countable Choice (CC). Let n≥1, 0<ε≤1, f∈C2(R;Rn) with f(0)=0, and u0∈Cc∞(Rn). Let uε be the viscous solution constructed in The viscous scalar Cauchy problem with smooth data has a global classical solution, so that uε∈C([0,T];Cb(Rn))∩C([0,T];L1(Rn))∩C1,2(Rn×(0,T)) solves utε+div⁡xf(uε)=εΔuε pointwise. Then for every t∈[0,T]:

(i) ∥uε(t)∥∞≤∥u0∥∞;

(ii) signed mass is conserved, ∫Rnuε(t,x) dx=∫Rnu0(x) dx, while the L1 norm is nonincreasing, ∥uε(t)∥1≤∥u0∥1 (in general it is not constant);

(iii) the integrated energy identity 12∥uε(t)∥22+ε∫0t ⁣ ⁣∫Rn∣∇uε∣2 dx ds=12∥u0∥22, so that ε∫0T ⁣ ⁣∫∣∇uε∣2≤12∥u0∥22, uniformly for 0<ε≤1.

Facts & Assumptions

Given: Countable Choice, n≥1, 0<ε≤1, f∈C2 with f(0)=0, u0∈Cc∞, the viscous solution uε, a time t∈[0,T], and the constant LM:=sup⁡∣s∣≤∥u0∥∞∣f′(s)∣.

[F1]

The constructed solution obeys the range bound sup⁡xuε(t,x)≤sup⁡xu0 and inf⁡xuε(t,x)≥inf⁡xu0, has AT:=sup⁡0≤s≤T∥uε(s)∥1<∞ by its L1-continuous orbit, and solves utε+div⁡f(uε)=εΔuε pointwise (The viscous scalar Cauchy problem with smooth data has a global classical solution).

[F2]

The viscous entropy balance holds pointwise for every convex C2 entropy: ∂tη(uε)+div⁡q(uε)=εΔη(uε)−εη′′(uε)∣∇uε∣2 with q′=η′f′ (The viscous entropy dissipation identity, Convex entropy--entropy flux pairs).

[F3]

There are smooth radial cutoffs 0≤χR≤1 with χR=1 on BR, χR=0 outside B2R, χR↑1 as R→∞, ∣DχR∣≤C/R and ∣ΔχR∣≤C/R2: use the explicit profile χR(x)=σ((4−∣x∣2/R2)/3) from the cited construction. On 0<t<1, differentiating its defining quotient gives σ′(t)=σ(t)(1−σ(t))(t−2+(1−t)−2)>0, so the profile is nonincreasing in radius and χR increases with R. The derivative scaling gives the stated gradient and Laplacian bounds (Explicit compactly supported smooth cutoffs, The Laplacian of a C2 function and of a C2 vector field).

[F4]

Spatial integration by parts against a compactly supported smooth χ follows from the one-dimensional theorem, not from a theorem on balls: enclose supp⁡χ in the interior of a box, fix all coordinates except xi, and apply Integration by parts for continuous factors with Riemann-integrable extensions of their interior derivatives on that coordinate interval. The integrands are continuous, their Riemann and Lebesgue integrals agree by A bounded Riemann integrable function on a closed bounded interval is Lebesgue measurable and has the same integral, and Fubini (Fubini's theorem for L^1 functions on a sigma-finite product) integrates the identity over the remaining coordinates. Boundary terms vanish since χ vanishes near the box boundary. Summing gives ∫χdiv⁡Q=−∫Q⋅∇χ; applying the same argument twice gives ∫χΔv=∫vΔχ for Q∈C1 and v∈C2. Dominated and monotone convergence justify the indicated cutoff and nonnegative limits (Dominated convergence, Monotone convergence for the integral, The space Lp(μ) as the quotient by null functions).

Proof

technique · direct
1.1F1

Supremum bound. Statement (i) is exactly the range bound of the construction theorem [F1]: ∥uε(t)∥∞≤max⁡{∣sup⁡xu0∣,∣inf⁡xu0∣}=∥u0∥∞.

1.2F1F3F4

Mass identity with cutoffs. Multiply the pointwise equation by a cutoff χR of [F3], integrate first from a positive lower time, and then pass that time to zero by C([0,T];L1) continuity: since χR vanishes outside a compact set, integration by parts is legitimate and gives ∫uε(t)χR−∫u0χR=∫0t ⁣ ⁣∫(f(uε)⋅∇χR+εuεΔχR); by [F1] the right side is bounded in absolute value by (CLMR−1+εCR−2)∫0t∥uε(s)∥1ds≤CT(R−1+εR−2).

1.3F1F2F3F4

Positive-time energy identity. Put M=∥u0∥∞. Fix 0<s<t≤T, take the convex entropy η(r)=r2/2 and its flux q(r)=∫0raf′(a) da, and integrate the balance [F2] against χR over [s,t]×Rn. This is legitimate on each compact support because uε∈C1,2 for positive times. The resulting identity is 12∫∣uε(t)∣2χR+ε∫st ⁣ ⁣∫∣∇uε∣2χR=12∫∣uε(s)∣2χR+∫st ⁣ ⁣∫q(uε)⋅∇χR+ε∫st ⁣ ⁣∫η(uε)ΔχR. On the solution range, ∣q(uε)∣≤LMM∣uε∣/2 and η(uε)≤M∣uε∣/2, so both cutoff errors tend to zero by [F1, F3]; the endpoint energies converge by dominated convergence. Since χR↑1, monotone convergence applies to the nonnegative dissipation and gives 12∥uε(t)∥22+ε∫st ⁣ ⁣∫∣∇uε∣2=12∥uε(s)∥22. In particular the dissipation is finite on every positive-time interval.

2.1F1F2F4step 1.2

Signed mass and L1 bound. In step 1.2 the right side tends to 0 when R→∞, while uε(t)χR→uε(t) and u0χR→u0 pointwise with ∣uε(t)∣χR≤∣uε(t)∣ and ∣u0∣χR≤∣u0∣, and both majorants are integrable by [F1]; dominated convergence gives ∫uε(t)=∫u0, which is signed-mass conservation. For the L1 bound let ηδ(s)=s2+δ2−δ, a convex C2 function with 0≤ηδ≤∣⋅∣, ηδ↑∣⋅∣ as δ↓0, and let qδ′=ηδ′f′ with qδ(0)=0. Testing the balance [F2] with χR, integrating first on a positive-time interval and passing its lower endpoint to zero by the Lipschitz entropy and L1 continuity and using ∣qδ(s)∣≤LM∣s∣ and ηδ(s)≤∣s∣ on the range of [F1], the cutoff terms vanish in the limit R→∞ exactly as in step 1.2, and dropping the nonpositive dissipation gives ∫ηδ(uε(t))≤∫ηδ(u0); monotone convergence in δ↓0 yields ∥uε(t)∥1≤∥u0∥1.

3.1F1step 1.3F4∎

Passage to the initial time. The construction gives uε∈C([0,T];L1) and ∥uε(s)∥∞≤M=∥u0∥∞. Hence ∥uε(s)−u0∥22≤2M∥uε(s)−u0∥1⟶0(s↓0), so the endpoint energy in step 1.3 converges to 12∥u0∥22. Letting s↓0, monotone convergence for the nonnegative space-time dissipation gives the identity in (iii) for every t>0; at t=0 it is immediate. Dropping the nonnegative final energy yields the stated uniform bound for 0<ε≤1.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Viscous solutions contract spatial translates in L-one

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the analytic prerequisites used below.

Let n≥1, 0<ε≤1, let f∈C2(R;Rn) satisfy f(0)=0, and let u0∈Cc∞(Rn). Let uε be the mild viscous solution with initial datum u0 from The viscous scalar Cauchy problem with smooth data has a global classical solution. Then for every h∈Rn and t∈[0,T], ∥uε(⋅+h,t)−uε(⋅,t)∥1≤∥u0(⋅+h)−u0(⋅)∥1. The estimate depends on f only through a bound for ∣f′∣ on the solution range. Thus it is uniform for any family of such fluxes with a common derivative bound on that range and the same initial datum (Absolute value in an ordered field, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable Choice, n≥1, 0<ε≤1, f∈C2 with f(0)=0, u0∈Cc∞, the viscous solution uε, a shift h∈Rn, and a spatial cutoff family χR(x)=χ(x/R) with χ∈Cc∞, 0≤χ≤1, χ=1 on B1 and χR→1 pointwise.

[F1]

The viscous solution is a classical solution with uε∈C([0,T];Cb)∩C([0,T];L1)∩C1,2(Rn×(0,T)), its range is the initial range, and it is bounded in L1 uniformly on [0,T] (The viscous scalar Cauchy problem with smooth data has a global classical solution, Uniform L-infinity, mass and energy bounds for the viscous approximations).

[F3]

Dominated convergence on compact sets and in the cutoff limit, and continuity in C([0,T];L1) permitting the limits R→∞ and s↓0 (Dominated convergence, The space Lp(μ) as the quotient by null functions). The smooth cutoffs with ∣DχR∣≤C/R and ∣ΔχR∣≤C/R2 are supplied by Explicit compactly supported smooth cutoffs.

Proof

technique · direct
1.1F1F2

The translate difference solves a linear equation. Fix h and put w(t,x)=uε(t,x+h)−uε(t,x), so that w∈C([0,T];L1)∩C1,2 by [F1]. Define a(t,x)=∫01f′(suε(t,x+h)+(1−s)uε(t,x)) ds; then a is C1 and bounded by LM=sup⁡∣s∣≤∥u0∥∞∣f′(s)∣, and f(uε(t,x+h))−f(uε(t,x))=a(t,x)w(t,x). Subtracting the two pointwise viscous equations and using the chain rule gives wt+div⁡(aw)=εΔw.

2.1F2step 1.1

The modulus balance. For δ>0 let ηδ(r)=r2+δ2, so that ηδ∈C∞, ηδ≥∣r∣, ∣ηδ′∣≤1, and ηδ(r)−rηδ′(r)=δ2/ηδ(r), ηδ′′(r)=δ2/ηδ(r)3. On compact subsets of Rn×(0,T), using step 1.1 and [F2], ∂tηδ(w)+div⁡(aηδ(w))−εΔηδ(w)=ηδ′(w)(wt+a⋅∇w−εΔw)+(div⁡a)ηδ(w)−εηδ′′(w)∣∇w∣2=δ2ηδ(w)div⁡a−εδ2∣∇w∣2ηδ(w)3.

3.1F1F3step 2.1

The limit δ↓0. The second term of step 2.1 is nonpositive, and the first is bounded in absolute value by δ ∣div⁡a∣, which tends to 0 in Lloc1 as δ↓0 because div⁡a is bounded on compact sets; since ηδ(w)→∣w∣ pointwise and ∣ηδ(w)∣≤∣w∣+δ, dominated convergence gives the distributional inequality ∂t∣w∣+div⁡(a∣w∣)≤εΔ∣w∣ on ΠT.

4.1F1F2F3step 3.1

Cutoff estimate. Test step 3.1 with χR(x) times a nonnegative smooth time test. In distributions in time this gives ddt∫∣w(t)∣χR≤C(LMR−1+εR−2)∥w(t)∥1. Approximating the indicator of (s,t) by smooth time cutoffs and using the L1 continuity of w gives, for all 0<s<t<T, ∫∣w(t)∣χR≤∫∣w(s)∣χR+C(LMR−1+εR−2)∫st∥w(r)∥1dr. The time integral is finite by [F1]. Dominated convergence as R→∞ yields ∥w(t)∥1≤∥w(s)∥1; continuity includes t=T.

5.1step 4.1F1F3∎

Conclusion. Letting s↓0 in step 4.1 and using the C([0,T];L1) continuity of w and w(0,⋅)=u0(⋅+h)−u0(⋅) gives ∥uε(⋅+h,t)−uε(⋅,t)∥1=∥w(t)∥1≤∥w(0)∥1=∥u0(⋅+h)−u0(⋅)∥1 for every t∈[0,T]. Every constant used depends on f only through LM, the bound for ∣f′∣ on the solution range, so the estimate is uniform over such flux families.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The Kruzhkov doubling inequality for two entropy solutions

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)). Let n≥1, T>0, f∈C1(R;Rn), and let u,v be bounded Kruzhkov entropy solutions on ΠT in the sense of Kruzhkov entropy solutions. Then, in the sense of distributions on ΠT, ∂t∣u−v∣+div⁡x(sgn⁡(u−v)(f(u)−f(v)))≤0. Equivalently, for every nonnegative φ∈Cc∞(ΠT), ∫ΠT(∣u−v∣ φt+sgn⁡(u−v)(f(u)−f(v))⋅∇xφ) dx dt≥0. Here sgn⁡(0)=0. Together with the weak equation this is the doubling-variables inequality from which uniqueness and the local L1 contraction are read off (Distribution, Distributional derivative).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, f∈C1, bounded entropy solutions u,v on ΠT, a nonnegative test function φ∈Cc∞(ΠT), and nonnegative unit-mass even mollifiers θh on R and ρh on Rn with ∫ρh=1, supp⁡θh⊂(−h,h) and supp⁡ρh⊂Bh (A radial mollifier family in Rn, Convolution of a distribution with a test function).

[F1]

For every k∈R the pair (ηk,qk) with ηk(s)=∣s−k∣, qk(s)=sgn⁡(s−k)(f(s)−f(k)) is a convex entropy pair with qk′=ηk′f′, and each of u,v satisfies the corresponding distributional inequality against every nonnegative test function (Kruzhkov entropy solutions, Convex entropy--entropy flux pairs, Absolute value in an ordered field).

[F2]

Fubini and dominated convergence apply on compact supports (Fubini's theorem for L^1 functions on a sigma-finite product, Dominated convergence). After multiplication by a fixed cutoff, u,v lie in L1(Rn+1) and their translations are norm continuous under Countable Choice (∥τhf−f∥p→0 in Lp(Rn) as h→0, for 1≤p<∞, The Axiom of Countable Choice (ACω), The space Lp(μ) as the quotient by null functions). This is the L1 diagonal interface; distributional mollifier convergence alone would not supply it.

Proof

technique · direct
1.1F2given

The doubled test function. For h>0 smaller than half the distance of the temporal support of φ from {0,T}, set gh(t,x,τ,y)=φ(t+τ2,x+y2)θh(t−τ)ρh(x−y); it is nonnegative and smooth with compact support in each pair of variables, and ∂tgh+∂τgh=(∂tφ)(t+τ2,x+y2)θh(t−τ)ρh(x−y), ∇xgh+∇ygh=(∇φ)(t+τ2,x+y2)θh(t−τ)ρh(x−y).

1.2F2given

The inequality for u with state k=v(τ,y). For almost every (τ,y) the constant k=v(τ,y) is admissible in [F1], and testing the entropy inequality for u by the nonnegative function gh(⋅,⋅,τ,y) gives ∫ΠT(∣u(t,x)−v(τ,y)∣ ∂tgh+sgn⁡(u−v(τ,y))(f(u)−f(v(τ,y)))⋅∇xgh) dx dt≥0; the integrand is bounded by a constant times the compactly supported smooth gh and its derivatives, so the left side is a bounded measurable function of (τ,y).

2.1F2step 1.2

Adding the symmetric inequality. Integrating the inequality of step 1.2 over (τ,y)∈ΠT, and likewise testing the entropy inequality for v with k=u(t,x) and integrating over (t,x), then adding, Fubini's theorem gives 0≤∫ ⁣ ⁣∫ΠT×ΠT[∣u−v∣ (∂tgh+∂τgh)+sgn⁡(u−v)(f(u)−f(v))⋅(∇xgh+∇ygh)]; the two integrals have the same bounded integrand because of the symmetry of gh.

3.1F2step 1.1step 2.1∎

Passing to the diagonal. Put z=(t,x) and z′=(τ,y)=z−h′. On a fixed compact set containing the doubled supports, translation continuity gives ∥v(⋅−h′)−v∥1→0 uniformly for ∣h′∣≤2h. The map Q(a,b)=sgn⁡(a−b)(f(a)−f(b)) is Lipschitz in each variable on the common bounded range: when a variable crosses the other one, split the interval at that point and use Q(b,b)=0 and the flux Lipschitz bound. Hence replacing v(z′) by v(z) changes the doubled integral by at most Csup⁡∣h′∣≤2h∥v(⋅−h′)−v∥1. Replacing Dφ((z+z′)/2) by Dφ(z) has error O(h) by smoothness, boundedness and unit kernel mass. Step 2.1 therefore converges to ∫ΠT(∣u−v∣φt+Q(u,v)⋅∇φ)≥0.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Local L1 contraction for two entropy solutions

Statement

Assume Countable Choice. Let n≥1, T>0, and f ⁣:R→Rn be C1. Let M,L≥0 and assume ∣f(a)−f(b)∣≤L∣a−b∣ for all a,b∈[−M,M]. Let u,v be bounded Kruzhkov entropy solutions on ΠT in the sense of Kruzhkov entropy solutions, with ∣u∣,∣v∣≤M almost everywhere and initial data u0,v0∈L∞∩Lloc1. For each fixed x0∈Rn and R>0, for almost every t∈(0,T) satisfying Lt<R, ∫B(x0,R−Lt)∣u(t,x)−v(t,x)∣ dx≤∫B(x0,R)∣u0(x)−v0(x)∣ dx. In particular, for each fixed x0,R, if u0=v0 almost everywhere on B(x0,R), then u(t,⋅)=v(t,⋅) almost everywhere on B(x0,R−Lt) for almost every t with Lt<R. The exceptional null set may depend on x0 and R.

Facts & Assumptions

Given: n≥1, T>0, f∈C1(R;Rn), constants M,L≥0 with ∣f(a)−f(b)∣≤L∣a−b∣ on [−M,M], bounded Kruzhkov entropy solutions u,v on ΠT with ∣u∣,∣v∣≤M almost everywhere and initial data u0,v0∈L∞∩Lloc1(Rn), a centre x0∈Rn and radius R>0, and the abbreviations w=∣u−v∣, q=sgn⁡(u−v)(f(u)−f(v)).

[F1]

Kato's inequality: for every nonnegative φ∈Cc∞(ΠT), ∫ΠT(w φt+q⋅∇xφ) dx dt≥0. Since ∣u∣,∣v∣≤M almost everywhere and f is L-Lipschitz on [−M,M], also ∣q∣≤Lw almost everywhere (The Kruzhkov doubling inequality for two entropy solutions, Lipschitz map, α-Hölder map for rational 0<α≤1, and contraction, Kruzhkov entropy solutions).

[F2]

Strong local L1 initial traces: for every compact K⊆Rn, lim⁡δ↓0ess sup⁡0<t<δ∫K∣u(t,x)−u0(x)∣ dx=0, and the same holds for v and v0 (Kruzhkov entropy solutions).

[F3]

Cutoff profiles: for every σ>0 there is a smooth nonincreasing βσ ⁣:R→[0,1] with βσ=1 on (−∞,R−σ] and βσ=0 on [R,∞), obtained by integrating a nonnegative smooth bump supported in (R−σ,R); then βσ′≤0 and Φσ(t,x)=βσ(∣x−x0∣+Lt) satisfies ∂tΦσ+L∣∇xΦσ∣=0 on the region where ∣x−x0∣+Lt>R−σ, and Φσ=1 on the region where ∣x−x0∣+Lt≤R−σ; hence Φσ is smooth on the slab 0<t<t3 whenever Lt3<R−σ, has compact spatial support contained in B(x0,R), and vanishes identically for t≥R/L if L>0 (A Euclidean bump for a compact set inside an open set, Open ball, closed ball and sphere in a metric space).

[F4]

Slice functions: Fσ(t)=∫Rnw(t,x)Φσ(t,x) dx is well defined for almost every t and locally integrable on its interval of definition, because w is bounded and Φσ is bounded with compact spatial support; hence almost every point is a Lebesgue point of Fσ, and the intersection of countably many full-measure sets is again full measure. Dominated and monotone convergence justify limits of integrals with uniformly bounded integrands against fixed integrable functions (Lebesgue differentiation theorem on Rn, Dominated convergence, The space Lp(μ) as the quotient by null functions, The Axiom of Countable Choice (ACω)).

Proof

technique · direct
1.1F1F3F4

Cutoff inequalities on a time slab. Fix t3∈(0,T) with Lt3<R — for L>0 such t3 exist by taking t3<min⁡{T,R/L}, and for L=0 every t3∈(0,T) works — and fix σ∈(0,R−Lt3). Let Φσ be as in [F3] and set Fσ(t)=∫Rnw(t,x)Φσ(t,x) dx and Gσ(t)=∫Rn(w ∂tΦσ+q⋅∇xΦσ)(t,x) dx for t∈(0,t3). Because βσ′≤0 and [F3] holds, Gσ≤∫Rnw(∂tΦσ+L∣∇xΦσ∣) dx=0 by [F1]. For nonnegative η∈Cc∞((0,t3)) the function φ=ηΦσ is an admissible nonnegative test function in [F1], since Φσ is smooth on the slab and compactly supported in x; hence 0≤∫0t3(Fση′+Gση)≤∫0t3Fση′, that is, ∫0t3Fση′≥0.

2.1F4step 1.1

Monotonicity in time. Fix a nonnegative smooth bump ζ supported in (0,1) with ∫01ζ=1 and put H(r)=∫−∞rζ. For 0<s<t<t3 and small ε>0, the function ηε(τ)=H(τ−sε)−H(τ−tε) is admissible in step 1.1 and ηε′→δs−δt as ε↓0. At Lebesgue points s<t of Fσ, step 1.1 gives Fσ(s)−Fσ(t)=lim⁡ε↓0∫0t3Fσηε′≥0, so Fσ(s)≥Fσ(t) for all Lebesgue points 0<s<t<t3 of Fσ, a full-measure set of pairs by [F4].

3.1F2F4step 2.1

The limit as s↓0. We claim ess lim⁡s↓0Fσ(s)=∫Rn∣u0(x)−v0(x)∣ βσ(∣x−x0∣) dx. Indeed, the difference is bounded by ∫B(x0,R)∣u(s)−u0∣ βσ+∫B(x0,R)∣v(s)−v0∣ βσ+∫Rn∣u0−v0∣ ∣βσ(∣x−x0∣+Ls)−βσ(∣x−x0∣)∣; the first two terms tend to 0 by [F2], since βσ is supported in B(x0,R), and the third tends to 0 because βσ has bounded derivative and ∣u0−v0∣∈Lloc1. Combining with step 2.1 and letting s↓0 through Lebesgue points of Fσ, for almost every t∈(0,t3), Fσ(t)≤∫Rn∣u0−v0∣ βσ(∣x−x0∣) dx.

4.1F4step 3.1∎

Removing the cutoff. Let σj↓0 with σj<R−Lt3 and intersect the full-measure sets of step 3.1 over all j using [F4]: for almost every t∈(0,t3) the inequality of step 3.1 with σ=σj holds for every j. For such t, βσj(∣x−x0∣+Lt)→1{∣x−x0∣<R−Lt}(x) pointwise away from the sphere ∣x−x0∣=R−Lt, and the corresponding integrands are dominated by w(t,⋅)1B(x0,R) respectively ∣u0−v0∣1B(x0,R), which are integrable; hence dominated convergence gives ∫B(x0,R−Lt)∣u−v∣(t,x) dx≤∫B(x0,R)∣u0−v0∣(x) dx. Every t∈(0,T) with Lt<R lies in (0,t3) for some admissible t3 — put H=min⁡{T,R/L} for L>0, and H=T for L=0, and use the explicit sequence t3(j)=H(1−1/(j+1))↑H — so the estimate holds for almost every such t, with exceptional set depending on x0,R. If u0=v0 almost everywhere on B(x0,R) the right-hand side vanishes, so u(t,⋅)=v(t,⋅) almost everywhere on B(x0,R−Lt) for almost every such t.

CorollaryStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Uniqueness, comparison and order preservation of entropy solutions

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the analytic prerequisites used below.

Let n≥1, T>0 and let f ⁣:R→Rn be C1 (hence Lipschitz on bounded intervals).

(i) If u,v are bounded Kruzhkov entropy solutions on ΠT in the sense of Kruzhkov entropy solutions with ∣u∣,∣v∣≤M and u0≤v0 almost everywhere, then u≤v almost everywhere on ΠT.

(ii) There is at most one bounded Kruzhkov entropy solution with a given initial datum u0∈L∞∩Lloc1; if u0=v0 almost everywhere on the whole of Rn, then u=v almost everywhere on ΠT.

(iii) The positive part contracts: ∫Rn(u(t,⋅)−v(t,⋅))+ dx≤∫Rn(u0−v0)+ dx for almost every t∈(0,T) whenever both sides are finite (Absolute value in an ordered field, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, a C1 flux f, bounded Kruzhkov entropy solutions u,v on ΠT with common essential bound M, and a constant L≥0 with ∣f(a)−f(b)∣≤L∣a−b∣ for a,b∈[−M,M].

[F1]

Entropy solutions are distributional weak solutions: ∂t(u−v)+div⁡x(f(u)−f(v))=0 in D′(ΠT); moreover each of u,v has the strong local L1 initial trace: for every compact K⊆Rn, ess sup⁡0<t<δ∫K∣u(t,x)−u0(x)∣ dx→0 as δ↓0, and likewise for v (Kruzhkov entropy solutions).

[F2]

Kato's inequality: ∂t∣u−v∣+div⁡x(sgn⁡(u−v)(f(u)−f(v)))≤0 in D′(ΠT) (The Kruzhkov doubling inequality for two entropy solutions).

[F3]

Positive-part identities: for every r∈R, (r)+=12(∣r∣+r) and 1{r>0}=12(sgn⁡(r)+1) for r≠0; at u=v the flux difference is zero, so the flux identity remains valid with sgn⁡(0)=0, so adding [F1] and [F2] gives ∂t(u−v)++div⁡x(1{u>v}(f(u)−f(v)))≤0 in D′(ΠT); writing w=(u−v)+ and q=1{u>v}(f(u)−f(v)), the Lipschitz hypothesis gives ∣q∣≤Lw almost everywhere, since ∣u∣,∣v∣≤M and f is L-Lipschitz on [−M,M] (Absolute value in an ordered field).

[F4]

Cutoff machinery of Local L1 contraction for two entropy solutions: for t3∈(0,T) with Lt3<R and 0<σ<R−Lt3 there is a smooth nonincreasing βσ with βσ=1 on (−∞,R−σ] and βσ=0 on [R,∞), and Φσ(t,x)=βσ(∣x−x0∣+Lt) satisfies ∂tΦσ+L∣∇xΦσ∣=0, is compactly supported in x, and is admissible as a test factor on (0,t3); for w≥0 with strong local L1 initial trace w0 and ∣q∣≤Lw, testing ∂tw+div⁡xq≤0 against ηΦσ and letting σ↓0 along a decreasing sequence yields ∫B(x0,R−Lt)w(t,x) dx≤∫B(x0,R)w0(x) dx for almost every t∈(0,T) with Lt<R; the argument uses Lebesgue points of t↦∫wΦσ, monotone and dominated convergence (Dominated convergence, Monotone convergence for the integral, The space Lp(μ) as the quotient by null functions).

Proof

technique · direct
1.1F1F2F3

The positive-part inequality. By [F1] and [F2], the sum of the weak equation for (u−v) and Kato's inequality is the distributional inequality ∂tw+div⁡xq≤0 with w=(u−v)+ and q=1{u>v}(f(u)−f(v)) by the identities of [F3], and ∣q∣≤Lw almost everywhere.

2.1F4step 1.1

Local positive-part estimate. Apply the cutoff computation [F4] to the pair (w,q) of step 1.1 with any centre x0 and radius R: ∫B(x0,R−Lt)(u−v)+(t,x) dx≤∫B(x0,R)(u0−v0)+(x) dx for almost every t∈(0,T) with Lt<R. Indeed the structural hypotheses of [F4] are met: the strong local L1 trace of w at 0 is w0=(u0−v0)+ because ∣w(t,x)−w0(x)∣≤∣u(t,x)−u0(x)∣+∣v(t,x)−v0(x)∣ almost everywhere, the positive part being 1-Lipschitz, and ∣q∣≤Lw holds by step 1.1; the cutoff, Lebesgue-point, initial-trace and σ↓0 steps are those of the proof of Local L1 contraction for two entropy solutions with ∣u−v∣ replaced by (u−v)+.

3.1F4step 2.1

Order preservation. Assume u0≤v0 almost everywhere, so (u0−v0)+=0 almost everywhere and the right-hand side of step 2.1 vanishes for every centre and radius. Take centres x0=0 and radii Rm=LT+m, m≥1, so that B(0,Rm−Lt)⊇B(0,m−LT) for every t∈(0,T); intersecting the countably many full-measure sets of times supplied by step 2.1, for almost every t∈(0,T) one has ∫B(0,m−LT)(u−v)+(t,x) dx=0 for every m with m>LT, hence (u−v)+(t,⋅)=0 almost everywhere on the union ⋃mB(0,m−LT)=Rn. By Fubini, (u−v)+=0 almost everywhere on ΠT, that is, u≤v almost everywhere.

4.1step 3.1

Uniqueness. If u,v are bounded entropy solutions with the same datum u0=v0, then both u0≤v0 and v0≤u0 hold almost everywhere, so step 3.1 gives u≤v and v≤u almost everywhere on ΠT, whence u=v almost everywhere. This proves both assertions of (ii).

5.1F4step 2.1∎

Positive-part contraction. Let u,v be any two bounded entropy solutions with both integrals finite; step 2.1 with centre 0 and radii Rm=LT+m gives ∫B(0,Rm−Lt)(u−v)+≤∫B(0,Rm)(u0−v0)+≤∫Rn(u0−v0)+<∞ for almost every t outside a null set Nm. On the complement of the null set ⋃mNm, all these inequalities hold, and monotone convergence over the increasing balls B(0,Rm−Lt)↑Rn gives ∫Rn(u(t,⋅)−v(t,⋅))+ dx≤∫Rn(u0−v0)+ dx, which is (iii).

CorollaryStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Finite propagation for scalar conservation laws

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the analytic prerequisites used below.

Let n≥1, T>0, and let f ⁣:R→Rn be C1 (hence Lipschitz on bounded intervals), with Lipschitz constant L≥0 on the common essential range of the solutions below. Let u,v be bounded Kruzhkov entropy solutions on ΠT in the sense of Kruzhkov entropy solutions. Fix x0∈Rn and R>0. If u0=v0 almost everywhere on B(x0,R), then for almost every t∈(0,T) with Lt<R, u(t,⋅)=v(t,⋅)almost everywhere on B(x0,R−Lt). For L=0 this holds for almost every t∈(0,T) on the stationary ball B(x0,R). In particular, if u0=0 almost everywhere outside B(x0,R), then u(t,x)=0 for almost every (t,x)∈ΠT with ∣x−x0∣>R+Lt. The conclusions are almost-everywhere statements for each fixed cone; an every-time claim requires a chosen time-continuous representative (Open ball, closed ball and sphere in a metric space).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, f∈C1, bounded Kruzhkov entropy solutions u,v on ΠT with common essential bound M and ∣u∣,∣v∣≤M almost everywhere, a constant L≥0 with ∣f(a)−f(b)∣≤L∣a−b∣ for a,b∈[−M,M], a centre x0∈Rn, and R>0.

[F1]

Local L1 contraction: for almost every t∈(0,T) with Lt<R, ∫B(x0,R−Lt)∣u(t,x)−v(t,x)∣ dx≤∫B(x0,R)∣u0(x)−v0(x)∣ dx, with an exceptional null set depending on x0,R; when L=0 the time condition is vacuous and the ball is stationary (Local L1 contraction for two entropy solutions, Lipschitz map, α-Hölder map for rational 0<α≤1, and contraction, Kruzhkov entropy solutions).

[F2]

The function identically 0 on ΠT is a Kruzhkov entropy solution with initial datum 0 for every flux f: for each k∈R the functions ηk(0)=∣k∣ and qk(0)=sgn⁡(−k)(f(0)−f(k)) are constant, so ∂tηk(0)+div⁡xqk(0)=0 in distributions, and the strong local L1 initial condition holds because ∫K∣0−0∣ dx=0 for every compact K (Kruzhkov entropy solutions).

[F3]

A nonnegative function in Lloc1 has vanishing integral over an open ball if and only if it vanishes almost everywhere there; almost-everywhere statements are statements about equivalence classes (The space Lp(μ) as the quotient by null functions).

[F4]

Countable measure bookkeeping: Q is countable and dense in R (Q is countably infinite, Both Q and R∖Q are dense in R, and every nonempty open subset of R is uncountable), hence finite products are countable (A product of two at most countable sets is at most countable). To approximate a point of Rn within δ, approximate each coordinate by a rational within δ/n; this proves density of Qn. a countable intersection of full-measure sets of times is full measure, and a countable union of null subsets of ΠT is null; Fubini gives the section-to-product nullity implication (Fubini's theorem for L^1 functions on a sigma-finite product). A countable union of null sets is null: finite-union indicators are bounded by the finite sums of the null-set indicators, and monotone convergence passes to their increasing union. Limits of integrals over expanding balls are covered by monotone and dominated convergence (Monotone convergence for the integral, Dominated convergence, Open ball, closed ball and sphere in a metric space).

Proof

technique · direct
1.1F1F3

The local estimate with vanishing right-hand side. Fix x0,R and suppose first that u0=v0 almost everywhere on B(x0,R). By [F1], for almost every t∈(0,T) with Lt<R, ∫B(x0,R−Lt)∣u(t,x)−v(t,x)∣ dx≤∫B(x0,R)∣u0(x)−v0(x)∣ dx=0, so u(t,⋅)=v(t,⋅) almost everywhere on B(x0,R−Lt) by [F3]. If L=0 the condition Lt<R reads 0<R and is automatic, the ball B(x0,R−Lt)=B(x0,R) is stationary, and the statement holds for almost every t∈(0,T).

2.1F2step 1.1

The support claim for the atomic ball. Suppose now that u0=0 almost everywhere outside B(x0,R) and let q∈Rn, r0>0 with B(q,r0)⊆Rn∖B‾(x0,R). Then u0=0 almost everywhere on B(q,r0), so step 1.1 applied to the pair (u,0) with centre q and radius r0 — legitimate because 0 is an entropy solution with datum 0 by [F2] — gives u(t,⋅)=0 almost everywhere on B(q,r0−Lt) for almost every t with Lt<r0.

3.1F4step 2.1

Covering the exterior cone. Let Q consist of rational pairs (q,r0)∈Qn×Q>0 satisfying r0<∣q−x0∣−R, and let C(q,r0)={(t,x):Lt<r0, ∣x−q∣<r0−Lt}. If ∣x−x0∣>R+Lt, choose rational q sufficiently close to x that Lt+∣x−q∣<∣q−x0∣−R, then a rational r0 between these bounds. Thus B(q,r0) lies in the strict exterior of the initial ball and (t,x)∈C(q,r0). The countable family of these cones covers the strict exterior cone.

4.1F4step 2.1step 3.1∎

Conclusion of the support claim. For each (q,r0)∈Q, step 2.1 exhibits a null set of times t with Lt<r0 such that u(t,⋅)≠0 on a positive-measure subset of B(q,r0−Lt); hence the set N(q,r0)={(t,x)∈C(q,r0) ⁣:u(t,x)≠0} is a null subset of ΠT, by Fubini: its sections are null for almost every time, and bounded spatial sections at the exceptional null set of times contribute zero. A countable union of null sets is null by [F4], so N=⋃(q,r0)∈QN(q,r0) is null, and by step 3.1 the set {(t,x)∈ΠT ⁣:∣x−x0∣>R+Lt, u(t,x)≠0} is contained in N. Therefore u=0 almost everywhere in the exterior cone.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Vanishing-viscosity families are locally precompact in L1

Statement

Assume Dependent Choice. Let n≥1, T>0, u0∈Cc∞(Rn), and put M=∥u0∥∞. Let (fj)j≥1 be C2 fluxes fj ⁣:R→Rn with fj(0)=0 and sup⁡jsup⁡∣s∣≤M+1∣fj′(s)∣<∞. For 0<εj≤1 with εj↓0, let uj be the global mild classical viscous solution with flux fj and datum u0, as supplied by The viscous scalar Cauchy problem with smooth data has a global classical solution. Then every subsequence has a further subsequence converging in L1(K) for every compact K⋐ΠT to u∈L∞∩Lloc1(ΠT) with ∥u∥∞≤M; a further subsequence converges almost everywhere on ΠT. The limit has a representative in C([0,T];Lloc1(Rn)) with u(0)=u0 in Lloc1. In particular this applies to a C1-convergent smooth approximation of one flux. The extraction uses Dependent Choice; energy dissipation alone does not give this strong compactness (Open ball, closed ball and sphere in a metric space, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Dependent Choice, n≥1, T>0, u0∈Cc∞(Rn), M=∥u0∥∞, C2 fluxes fj with fj(0)=0 and L:=sup⁡jsup⁡∣s∣≤M+1∣fj′(s)∣<∞, 0<εj≤1 with εj↓0, and the viscous solutions uj of utj+div⁡xfj(uj)=εjΔuj on ΠT with uj(0,⋅)=u0.

[F1]

Each uj is a classical global solution with uj∈C([0,T];Cb(Rn))∩C([0,T];L1(Rn))∩C1,2(Rn×(0,T)), and uj(t,⋅)→u0 in L1(Rn) as t↓0 (The viscous scalar Cauchy problem with smooth data has a global classical solution).

[F2]

Uniform bounds: ∣uj(t,x)∣≤M and ∥uj(t,⋅)∥1≤∥u0∥1 for all j≥1 and t∈[0,T] (Uniform L-infinity, mass and energy bounds for the viscous approximations).

[F3]

Uniform spatial modulus: for all j, t∈[0,T] and z∈Rn, ∥uj(t,⋅+z)−uj(t,⋅)∥1≤∥u0(⋅+z)−u0(⋅)∥1≤ω0(∣z∣), where ω0(r)=sup⁡∣z∣≤r∥u0(⋅+z)−u0(⋅)∥1↓0 as r↓0 by uniform continuity of the compactly supported u0; the estimate depends on fj only through the derivative bound L on the range, so it is uniform in j (Viscous solutions contract spatial translates in L-one).

[F4]

Mollification and cutoffs: for an even mollifier ϱh with support in Bh and a bounded compactly supported F∈L∞, the convolution g=ϱh∗F is smooth with ∇g=(∇ϱh)∗F, Δg=(Δϱh)∗F (differentiation under the integral sign via difference quotients and dominated convergence) and ∥Dαg∥∞≤Ch−∣α∣ for ∣α∣≤2; Fubini's theorem gives ∫gw=∫F(ϱh∗w) dx whenever w∈Lloc1 and F is bounded with compact support, and the mollification error obeys ∥ϱh∗w−w∥L1(BR+ρ/2)≤sup⁡∣y∣≤h∥w(⋅+y)−w(⋅)∥L1(BR+ρ) (A radial mollifier family in Rn, Fubini's theorem for L^1 functions on a sigma-finite product, Dominated convergence).

[F5]

For every R>0, ρ∈(0,1) and every compact cylinder there exist smooth cutoffs 0≤χ∈Cc∞(BR+ρ/2) with χ=1 on BR, and ψ∈Cc∞(Rn), θ∈Cc∞((0,T)) equal to 1 on the spatial and temporal projections of that cylinder (Open ball, closed ball and sphere in a metric space, A Euclidean bump for a compact set inside an open set).

[F6]

Fréchet–Kolmogorov criterion: a family in L1(Rd) that is uniformly bounded, has uniformly small tails, and is uniformly translation-continuous is relatively compact, and every sequence in it has a subsequence converging in L1(Rd); the criterion is used with the Axioms of Countable and Dependent Choice, and L1 convergence yields an almost-everywhere convergent subsequence, while each L1 space is complete (The Fr'echet--Kolmogorov compactness criterion in Lp(Rn), Convergence in L^1(mu) has an almost-everywhere convergent subsequence, Riesz-Fischer completeness of Lp for 1≤p≤∞, The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain, The Axiom of Countable Choice (ACω)).

Proof

technique · direct
1.1F2F3

Uniform bounds and spatial modulus. By [F2], ∣uj(t,x)∣≤M and ∥uj(t,⋅)∥1≤∥u0∥1 for all j,t; by [F3], ∥uj(t,⋅+z)−uj(t,⋅)∥1≤ω0(∣z∣) with ω0(r)↓0, uniformly in j and t.

1.2F1F2F4F5

The local time modulus. Fix R>0, 0<ρ<1, j, and t,τ with t,t+τ∈[0,T]; put w=uj(t+τ,⋅)−uj(t,⋅). Choose χ as in [F5] and an even mollifier ϱh with 0<h<ρ/4; set z=ϱh∗w and g=ϱh∗(χ sgn⁡z)∈Cc∞(BR+ρ), so that ∥Dαg∥∞≤Ch−∣α∣ for ∣α∣≤2 and ∫Rngw dx=∫Rnχ∣z∣ dx by [F4]. Assume τ>0; negative increments follow by reversing the two times. Since uj solves the viscous equation pointwise at positive times and is C1 there, dds∫g uj(t+s,⋅) dx=∫∇g⋅fj(uj(t+s,⋅)) dx+εj∫Δg uj(t+s,⋅) dx for s∈[0,τ]; integrating over s and using ∥fj(v)∥1≤L∥v∥1 for ∣v∣≤M and [F2], ∣∫gw∣≤∫0τ(∥∇g∥∞L∥u0∥1+εj∥Δg∥∞∥u0∥1)ds≤CRτh−2.

2.1F3F4step 1.2

The uniform modulus in Lloc1. With χ,w,z as in step 1.2 and [F4], ∫BR∣w∣≤∫χ∣w∣≤∫χ∣z∣+∫χ∣w−z∣≤CRτh−2+2ω0(h), because ∫χ∣z∣=∣∫gw∣ and ∥w−z∥L1(BR+ρ/2)≤2ω0(h) by [F3] applied at the two times t and t+τ. Taking h=τ1/3 (for 0<τ<(ρ/4)3) gives ∫BR∣uj(t+τ)−uj(t)∣≤CR(ω0(τ1/3)+τ1/3)=:ηR(τ) with ηR(τ)→0 as τ↓0, uniformly in j and in t∈[0,T−τ]; the endpoint t=0 follows by first integrating from a positive time and then using the L1 continuity uj(t)→u0.

3.1F5F6F7step 2.1

Relative compactness on cylinders. Fix l≥1 and cutoffs ψl,θl as in [F5] equal to 1 on the projections of the cylinder Kl=B‾(0,l)×[1/l,T−1/l], and set Flj(t,x)=θl(t)ψl(x) uj(t,x), extended by zero. The family (Flj)j is uniformly bounded in L1(Rn+1) by T∥θlψl∥∞∥u0∥1, has common compact support (uniform tails), and is uniformly translation-continuous: for a shift (y,s) one has ∥Flj(⋅+(y,s))−Flj∥1≤∥Flj(t+s,x+y)−Flj(t,x+y)∥1+∥Flj(t,x+y)−Flj(t,x)∥1, where the spatial part is bounded by T∥(θlψl)(⋅+y)−(θlψl)∥∞∥u0∥1+T∥θlψl∥∞ω0(∣y∣) and the temporal part by T∥θl(⋅+s)−θl∥∞∥ψl∥∞∥u0∥1+T∥θlψl∥∞ηl′(∣s∣) for a radius l′ containing supp⁡ψl, both tending to 0 as ∣(y,s)∣→0 uniformly in j by step 1.1 and step 2.1. By [F6] every subsequence of (Flj)j has a further subsequence converging in L1(Rn+1), hence, after diagonal extraction over the countably many l≥1 (Dependent Choice), some subsequence of (uj) converges in L1(Kl) for every l, and therefore in L1(K) for every compact K⋐ΠT, since each such K lies in some Kl.

3.2F6F7step 2.1

Continuous representative and the initial trace. Pass to a further subsequence converging almost everywhere on ΠT by [F6]. Since ∣ujk∣≤M, this gives ∣u∣≤M almost everywhere; by Fubini and dominated convergence, there is a common full-measure set S⊂(0,T) on which ujk(t)→u(t) in L1(Bm) for every integer m≥1. Thus S is dense. Fix m and choose an integer m′>m. For s,t∈S sufficiently close that step 2.1 applies, lower semicontinuity on the local ball, followed by its time-modulus estimate on the larger ball, gives ∥u(t)−u(s)∥L1(Bm)≤lim inf⁡k∥ujk(t)−ujk(s)∥L1(Bm)≤ηm′(∣t−s∣). For t∈S sufficiently close to 0, the same local comparison with ujk(0,⋅)=u0 gives ∥u(t)−u0∥L1(Bm)≤ηm′(t). This modulus tending to zero at zero makes the map on S uniformly continuous; completeness of L1(Bm) extends it uniquely to a continuous map on [0,T], with value u0 at t=0. These extensions agree on nested balls because they agree on the dense set S. Consequently u has a representative in C([0,T];Lloc1(Rn)) with u(0)=u0 in Lloc1.

4.1F2F6step 3.1step 3.2

Limit properties. The diagonal limit of step 3.1 lies in Lloc1(ΠT) and the selected subsequence converges in L1(K) for every compact K⋐ΠT; the further almost-everywhere subsequence in step 3.2 preserves these convergences and passes the uniform bound ∣ujk∣≤M to u. This establishes the asserted limit and almost-everywhere convergence.

5.1F2F6step 1.1step 2.1step 3.1∎

Applicability and hypotheses. If fj→f in C1 on compact sets with sup⁡jsup⁡∣s∣≤M+1∣fj′(s)∣<∞ — a C1-convergent smooth approximation of one flux — then the hypotheses above hold, so the conclusions apply. The extraction used Dependent Choice in the diagonal step 3.1 (and Countable Choice inside [F6]); the uniform energy dissipation εj∫0T ⁣ ⁣∫∣∇uj∣2≤12∥u0∥22 supplies only a uniform gradient bound and would not by itself give the compactness in L1 obtained from the uniform bounds and translation moduli of steps 1.1 and 2.1.

CorollaryStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Global L1 contraction from the local estimate

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the analytic prerequisites used below.

Let n≥1, T>0, let f ⁣:R→Rn be C1 and Lipschitz on the common essential range of the solutions below with constant L≥0, and let u,v be bounded Kruzhkov entropy solutions on ΠT in the sense of Kruzhkov entropy solutions, with initial data u0,v0∈L1(Rn)∩L∞(Rn). Then for almost every t∈(0,T), ∥u(t,⋅)−v(t,⋅)∥1≤∥u0−v0∥1. If u,v have representatives continuous in Lloc1 on [0,T], the same inequality holds for every t∈[0,T] (Open ball, closed ball and sphere in a metric space, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, f Lipschitz with constant L≥0 on the common essential range of u and v, bounded Kruzhkov entropy solutions u,v on ΠT with ∣u∣,∣v∣≤M almost everywhere, and initial data u0,v0∈L1(Rn)∩L∞(Rn). Set L∗:=sup⁡∣z∣≤M∣f′(z)∣<∞; in the proof use L∗, irrespective of the given constant on the common essential range.

[F1]

Since f∈C1, L∗<∞, and coordinatewise FTC gives f(b)−f(a)=∫abf′(z) dz and ∣f(b)−f(a)∣≤L∗∣b−a∣ for a,b∈[−M,M] (The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)). Local L1 contraction with this interval constant: for every centre x0 and radius R>0, for almost every t∈(0,T) with L∗t<R, ∫B(x0,R−L∗t)∣u(t,x)−v(t,x)∣ dx≤∫B(x0,R)∣u0(x)−v0(x)∣ dx, with an exceptional null set depending on x0,R (Local L1 contraction for two entropy solutions, Kruzhkov entropy solutions).

[F2]

Monotone convergence for integrals of nonnegative functions over increasing sets: if 0≤gm↑g pointwise then ∫gm↑∫g; in particular the integrals of a fixed nonnegative function over the balls B(0,ρ)↑Rn increase to its integral over Rn, finite or infinite (Monotone convergence for the integral).

[F3]

Almost-everywhere assertions concern equivalence classes: a countable union of null sets in (0,T) is null, and members of L1 are defined up to modification on null sets (The space Lp(μ) as the quotient by null functions, Open ball, closed ball and sphere in a metric space).

Proof

technique · direct
1.1F1

Ball estimates along an exhausting sequence. For m≥1 put Rm=L∗T+m, so that Rm>L∗t for every t∈(0,T) and Rm−L∗t=m+L∗(T−t)>0. Applying [F1] with centre 0 and radius Rm gives, for every m, an exceptional null set Nm⊆(0,T) such that for all t∈(0,T)∖Nm ∫B(0,Rm−L∗t)∣u(t,x)−v(t,x)∣ dx≤∫B(0,Rm)∣u0(x)−v0(x)∣ dx≤∥u0−v0∥1.

2.1F2F3step 1.1

Intersection and monotone limit. The set N=⋃m≥1Nm is null by [F3]. Fix t∈(0,T)∖N, so that the estimates of step 1.1 hold for every m. The balls B(0,Rm−L∗t) increase to Rn as m↑∞, hence the integrals of the fixed nonnegative function ∣u(t,⋅)−v(t,⋅)∣ over them increase to ∫Rn∣u(t,x)−v(t,x)∣ dx, while ∫B(0,Rm)∣u0−v0∣↑∥u0−v0∥1 by monotone convergence. Passing to the limit in step 1.1 gives ∥u(t,⋅)−v(t,⋅)∥1≤∥u0−v0∥1<∞ for every t∈(0,T)∖N, which is the almost-everywhere assertion and shows that the slice integrals are finite for almost every t.

3.1F2step 2.1∎

The every-time assertion under Lloc1 continuity. Assume now that u and v have representatives on [0,T] such that tj→t implies u(tj)→u(t) and v(tj)→v(t) in L1(K) for every compact K⊆Rn (with one-sided sequences at t=0,T). Fix t∈[0,T] and choose tj∈(0,T)∖N with tj→t, possible because N is null. For every fixed ball B(0,ρ), step 2.1 gives ∫B(0,ρ)∣u(tj)−v(tj)∣≤∥u(tj)−v(tj)∥1≤∥u0−v0∥1, and ∣u(tj)−v(tj)∣→∣u(t)−v(t)∣ in L1(B(0,ρ)) because u(tj)→u(t) and v(tj)→v(t) there; the inequality ∣∥a∥1−∥b∥1∣≤∥a−b∥1 gives convergence of these integrals directly. Hence ∫B(0,ρ)∣u(t)−v(t)∣≤lim inf⁡j∫B(0,ρ)∣u(tj)−v(tj)∣≤∥u0−v0∥1. Letting ρ↑∞ and using monotone convergence once more gives ∥u(t)−v(t)∥1≤∥u0−v0∥1.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Existence of bounded Kruzhkov entropy solutions

Statement

Assume Countable Choice and Dependent Choice (The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) for the heat-kernel, L1 completeness and vanishing-viscosity extraction interfaces used below. Let n≥1, T>0, let f ⁣:R→Rn be locally Lipschitz and C1, and let u0∈L∞(Rn)∩L1(Rn). Then there exists a Kruzhkov entropy solution u∈L∞(ΠT)∩C0([0,T];Lloc1(Rn)) of ut+div⁡xf(u)=0 with u(⋅,0)=u0 in the strong Lloc1 sense, satisfying ∥u(t,⋅)∥∞≤∥u0∥∞ for every t. By Uniqueness, comparison and order preservation of entropy solutions this solution is unique. The route is vanishing viscosity, and the compactness comes from the translation lemmas, not from the energy dissipation alone (Scalar conservation laws, fluxes and Cauchy data, Kruzhkov entropy solutions).

Facts & Assumptions

Given: Countable and Dependent Choice, n≥1, T>0, a locally Lipschitz C1 flux f, and u0∈L∞(Rn)∩L1(Rn) with M=∥u0∥∞.

[F1]

The weak formulation: u is a distributional weak solution of ut+div⁡xf(u)=0 on ΠT iff ∫ΠT(uφt+f(u)⋅∇φ)=0 for every φ∈Cc∞(ΠT); subtracting the constant f(0) from the flux changes neither the divergence term nor the Kruzhkov fluxes sgn⁡(u−k)(f(u)−f(k)), so all existence and entropy statements may be proved for the normalized flux f~=f−f(0) and transferred back (Scalar conservation laws, fluxes and Cauchy data, Kruzhkov entropy solutions).

[F2]

Viscous solutions: for every C2 flux g with g(0)=0, every ε∈(0,1] and every datum in Cc∞, there is a global classical solution of ut+div⁡xg(u)=εΔu with u∈C([0,T];Cb)∩C([0,T];L1)∩C1,2(Rn×(0,T)), range contained in the initial range, and ∥u(t,⋅)∥1≤∥u0∥1 (The viscous scalar Cauchy problem with smooth data has a global classical solution, Uniform L-infinity, mass and energy bounds for the viscous approximations).

[F3]

Viscous entropy balance: for every convex C2 entropy η with flux q′=η′g′, ∂tη(u)+div⁡xq(u)=εΔη(u)−εη′′(u)∣∇u∣2 pointwise, and the Laplacian term integrates by parts against compactly supported tests (The viscous entropy dissipation identity).

[F4]

Mollification: convolving a locally integrable function with a radial mollifier gives a smooth function; the mollified derivatives are the convolutions of the derivatives, and on compact sets the mollified flux and its derivative converge uniformly to the original for C1 data; approximate identities converge in L1, and the classes of Lp are equivalence classes (A unit-mass smooth bump generates an L1 approximate identity, Convolution with a mollifier is smooth, and derivatives pass under the integral sign, The space Lp(μ) as the quotient by null functions). The norm convergence assertion is Every L1 approximate identity converges to the identity in Lp for 1≤p<∞; the compact cutoffs are A Euclidean bump for a compact set inside an open set.

[F5]

Vanishing-viscosity compactness: for smooth compactly supported data, C2 fluxes fm with fm(0)=0 and uniformly bounded derivatives on the range [−M−1,M+1], and 0<εm≤1, εm↓0, the viscous solutions have a subsequence converging in L1(K) for every compact K⋐ΠT and almost everywhere on ΠT to u∈L∞∩Lloc1(ΠT) with ∥u∥∞≤M, having a representative in C([0,T];Lloc1) with u(0)=u0 in Lloc1; the extraction uses Dependent Choice (Vanishing-viscosity families are locally precompact in L1).

[F6]

Contraction and uniqueness: for data in L1∩L∞ the global L1 difference of two entropy solutions is bounded by the L1 difference of their data, at every time if the representatives are Lloc1-continuous; two bounded Kruzhkov entropy solutions with the same data coincide almost everywhere (Global L1 contraction from the local estimate, Uniqueness, comparison and order preservation of entropy solutions).

[F7]

Limits and completeness: dominated and monotone convergence for integrals over fixed compact sets with uniformly bounded integrands, and completeness of L1(K) for compact K (Dominated convergence, Monotone convergence for the integral, Riesz-Fischer completeness of Lp for 1≤p≤∞).

Proof

technique · direct
1.1F1F4

Normalization and smooth flux approximation. By [F1] it suffices to treat f~=f−f(0). Choose a radial mollifier and, for m≥1, set fm(s)=χ(s) (f~∗ρ1/m)(s)−(f~∗ρ1/m)(0), where χ∈Cc∞ equals 1 on [−M−12,M+12] and is supported in [−M−1,M+1]. Then each fm is C∞ (hence C2), fm(0)=0, and on [−M−12,M+12] one has fm→f~ and fm′=f~′∗ρ1/m→f~′ uniformly by [F4]; in particular sup⁡msup⁡∣s∣≤M+1∣fm′(s)∣<∞.

1.2F2F5

Smooth-datum case: extraction. Let u0∈Cc∞(Rn) with ∥u0∥∞=M, choose εm↓0 with 0<εm≤1, and let um be the global classical solution with flux fm and datum u0 given by [F2]. Then ∣um(t,x)∣≤∥u0∥∞=M by the range bound of [F2], so the hypotheses of [F5] are met; passing to a subsequence (relabelled) there is u∈L∞∩Lloc1(ΠT) with ∥u∥∞≤M, um→u in L1(K) for every compact K⋐ΠT and almost everywhere, and u has a representative in C([0,T];Lloc1) with u(0)=u0 in Lloc1. At every time t this representative obeys ∣u(t,⋅)∣≤M almost everywhere: for t in the full-measure set where the construction of [F5] passes the bound, slicewise almost-everywhere convergence preserves it, and general t follows by Lloc1-continuity.

2.1F1F7step 1.1step 1.2

The weak equation passes to the limit. For φ∈Cc∞(ΠT), testing the pointwise viscous equation of um gives ∫ΠT(umφt+fm(um)⋅∇φ)=−εm∫ΠTumΔφ, whose right side is bounded by εmM∥Δφ∥1→0. The left side converges to ∫ΠT(uφt+f~(u)⋅∇φ) by [F7], because um→u almost everywhere with uniform bounds and fm→f~ uniformly on [−M,M]; hence u is a weak solution for f~, and therefore for f by [F1].

2.2F3F4F7step 1.1step 1.2

Entropy inequalities for ∣k∣≤M. Fix k∈[−M,M] and δ>0, and put ηδ(r)=r2+δ2−δ, a convex C2 function with 0≤ηδ≤∣⋅∣, ∣ηδ(r)−∣r∣∣≤δ, and ηδ′′≥0; let qδ,k,m(s)=∫ksηδ′(r−k)fm′(r) dr, so qδ,k,m′=ηδ′(⋅−k)fm′. Let φ∈Cc∞(Rn×[0,T)) be nonnegative. Multiplying the exact viscous balance [F3] for the pair (ηδ(⋅−k),qδ,k,m) by φ, integrating over ΠT and integrating the Laplacian by parts gives Im,δ:=∫ΠT(ηδ(um−k)φt+qδ,k,m(um)⋅∇φ)+∫Rnηδ(u0−k)φ(x,0) dx=−εm∫ΠTηδ(um−k)Δφ+εm∫ΠTηδ′′(um−k)∣∇um∣2φ — the only boundary term is the one at t=0, displayed with the initial datum; the last term is nonnegative and the first is bounded by εm⋅2M∥Δφ∥1→0, so lim inf⁡mIm,δ≥0. On the other hand Im,δ converges by [F7]: um→u almost everywhere, ∣ηδ(um−k)∣≤2M, and qδ,k,m→qδ,k uniformly on [−M,M], where qδ,k(s)=∫ksηδ′(r−k)f~′(r) dr, so the limit obeys Iδ=∫ΠT(ηδ(u−k)φt+qδ,k(u)⋅∇φ)+∫Rnηδ(u0−k)φ(x,0) dx≥0.

3.1F7step 2.1step 2.2

The Kruzhkov inequalities. Letting δ↓0 in step 2.2, ηδ(r−k)→∣r−k∣ uniformly on [−M,M] and qδ,k(s)→∫kssgn⁡(r−k)f~′(r) dr=sgn⁡(s−k)(f~(s)−f~(k))=qk(s) uniformly on [−M,M] by dominated convergence, since ∣ηδ′∣≤1 and f~′ is continuous there; hence dominated convergence gives ∫ΠT(∣u−k∣φt+qk(u)⋅∇φ)+∫∣u0−k∣φ(x,0) dx≥0 for every nonnegative φ∈Cc∞(Rn×[0,T)) and every k∈[−M,M]. If k>M, then u−k<0 almost everywhere, ηk is affine on the range [−M,M], and qk(s)=f~(k)−f~(s), so the identity ∂tηk(u)+div⁡xqk(u)=−(ut+div⁡xf~(u))=0 holds by step 2.1, and similarly for k<−M; thus all Kruzhkov inequalities hold and, with the trace of step 1.2, u is a bounded Kruzhkov entropy solution with datum u0.

4.1F4F6step 2.1step 2.2step 3.1

General datum: approximation and Cauchy property. Now let u0∈L∞∩L1 be arbitrary. By [F4] choose u0r∈Cc∞ with ∥u0r∥∞≤∥u0∥∞ and u0r→u0 in L1 (truncate u0 to a large ball and mollify). Steps 1.1–3.1 applied to each smooth datum u0r give bounded Kruzhkov entropy solutions ur for the flux f, with Lloc1-continuous representatives and ∥ur(t,⋅)∥∞≤∥u0r∥∞≤M. By [F6], for every t∈[0,T], ∥ur(t)−uℓ(t)∥1≤∥u0r−u0ℓ∥1, so sup⁡t∈[0,T]∥ur(t)−uℓ(t)∥L1(K)≤∥u0r−u0ℓ∥1 for every compact K, and the right side tends to 0 as r,ℓ→∞.

5.1F7step 3.1step 4.1

The limit for general datum. By completeness of L1(K) [F7] and the uniform-in-time contraction in step 4.1, ur converges in C([0,T];L1(K)) for each compact ball K, consistently on nested balls, to u∈C([0,T];Lloc1(Rn)). Since ur(0)=u0r→u0 in L1, this representative has initial trace u0. For every t and ball K, the inequality (∣u(t)∣−M)+≤∣u(t)−ur(t)∣+(∣ur(t)∣−M)+ and ∥ur(t)∥∞≤M show, after integration on K and passage to the L1(K) limit, that (∣u(t)∣−M)+=0 almost everywhere there. Thus ∥u(t)∥∞≤M for every t. The weak equation passes to the limit because f is Lipschitz on [−M,M] and ur→u in local L1. For each fixed k∈R, s↦∣s−k∣ is 1-Lipschitz and qk(s)=sgn⁡(s−k)(f(s)−f(k)) is Lipschitz on [−M,M] with constant at most sup⁡∣s∣≤M∣f′(s)∣. Therefore the entropy and flux terms in the inequality of step 3.1 converge in L1 on every test support; the initial entropy term converges by u0r→u0 in L1 and the same Lipschitz bound for ∣⋅−k∣. Passing to the limit proves every Kruzhkov inequality without requiring an almost-everywhere subsequence for the general-data approximation. Hence u is a bounded Kruzhkov entropy solution with datum u0.

6.1F6step 1.2step 5.1∎

Uniqueness and conclusion. If v is another bounded Kruzhkov entropy solution with the same datum u0, then u=v almost everywhere on ΠT by [F6], so the constructed solution is the unique bounded Kruzhkov entropy solution with this datum; this completes the proof.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The convex entropy condition for a single shock is the chord condition

Statement

Let n=1, f∈C1(R), and let u be a piecewise C1 weak solution with a single jump from the left state u− to the right state u+ across a C1 curve x=s(t) whose speed satisfies the Rankine--Hugoniot condition s′=(f(u+)−f(u−))/(u+−u−). Put F(z)=f(z)−f(u−)−s′(z−u−), so that F(u−)=F(u+)=0. Then the entropy inequality η(u)t+q(u)x≤0 of Convex entropy--entropy flux pairs holds for every convex entropy pair (η,q) if and only if F(z) (u+−u−) ≥ 0for every z between u− and u+; equivalently, in the case u−<u+ the graph of f on [u−,u+] lies above the chord joining (u−,f(u−)) and (u+,f(u+)), while in the case u−>u+ it lies below that chord, both in the non-strict sense (Piecewise smooth shocks and one-sided traces, Convex and strictly convex functions on Euclidean convex sets).

Facts & Assumptions

Given: n=1, f∈C1, a single-jump piecewise C1 weak solution with states u−≠u+ and speed s′ satisfying Rankine--Hugoniot, and an arbitrary convex entropy pair (η,q) with η∈C2 and q′=η′f′.

[F1]

The jump configuration and Rankine--Hugoniot condition are as in Piecewise smooth shocks and one-sided traces and The Rankine--Hugoniot jump condition in space--time normal form: in one dimension the interface is a graph x=s(t) with minus side x<s(t), plus side x>s(t), unit normal ν=(−s′,1)/1+s′2, and s′(u+−u−)=f(u+)−f(u−).

[F2]

The graph integration of The Rankine--Hugoniot jump condition in space--time normal form, applied to (η(u),q(u)), gives the interface production ([q]−s′[η])δ(x−s(t)), where the latter distribution pairs by ∫φ(t,s(t))dt. For smooth pairs the production vanishes in the classical side regions by the chain rule. Smooth nonnegative bumps can be placed on any interface patch (Explicit compactly supported smooth cutoffs). Thus the entropy inequality is equivalent to nonpositive jump production at every point (Convex entropy--entropy flux pairs, Kruzhkov entropy solutions).

[F3]

Primitives: since η′ and f′ are continuous, q(z)=∫0zη′(r)f′(r)dr up to a constant, and increments of C1 functions are integrals of their derivatives; the fundamental theorem of calculus, its use under limits, and the primitive construction are as in Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G and The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a).

[F4]

Approximation tools: monotone bounded convergence for limits of test functions and dominated convergence for the passing of inequalities (Monotone convergence for the integral, Dominated convergence, Absolute value in an ordered field).

Proof

technique · direct
1.1F2F3

Jump entropy production. By [F2] the entropy inequality for (η,q) is equivalent to [q]−s′[η]≤0. Using [F3] in the orientation of the jump, [q]=∫u−u+q′(z) dz=∫u−u+η′(z)f′(z) dz and [η]=∫u−u+η′(z) dz, so the condition is ∫u−u+η′(z)(f′(z)−s′) dz≤0.

2.1F1F3step 1.1

Smooth-pair sufficiency. At a fixed interface point put a=u−, b=u+ and s′=[f]/[u]. Integration by parts, with F(a)=F(b)=0, gives [q]−s′[η]=∫abη′(z)F′(z) dz=−∫abη′′(z)F(z) dz. If a<b and F≥0, this is nonpositive since η′′≥0. If a>b and F≤0, reversal of the integral gives the same conclusion.

3.1F1F2F3step 2.1

Necessity. If a<b and F(z0)<0 for some interior z0, continuity supplies an interval on which F<0. Choose a smooth nonnegative bump β supported there and not identically zero, and define η′(z)=∫0zβ(r)dr, η(z)=∫0zη′(r)dr. Then η′′=β≥0, so this is a smooth convex entropy; step 2.1 gives strictly positive production, a contradiction. If a>b and F(z0)>0, the same bump and reversed integral again give positive production. Thus all smooth convex inequalities force F(z)(b−a)≥0. For a Kruzhkov pair with k between the states, a direct subtraction gives [qk]−s′[ηk]=−2sgn⁡(b−a)F(k); outside the interval it is zero. This also proves exact equivalence with the Kruzhkov jump criterion.

4.1F3F4step 2.1step 3.1∎

Nonsmooth pairs and chord interpretation. A finite convex entropy is uniformly approximated on compact intervals by its convolution with a nonnegative smooth unit-mass bump at scale δ. These convolutions are smooth and convex (average the convexity inequality), and their derivatives converge at each differentiability point of η, while remaining bounded by a common local Lipschitz constant. The integral fluxes therefore converge uniformly by dominated convergence, so the smooth entropy inequalities of step 2.1 pass to every locally Lipschitz convex pair, both in the side regions and at the jump. Together with step 3.1 this proves the equivalence. Finally F(z)≥0 means f(z) lies above f(a)+s′(z−a) for a<b; F(z)≤0 means it lies below for a>b. This line is the chord through the two states.

DefinitionDefinition: Literature-sourcedProof: Not applicableprecheck passjudge pass (gpt-6.1-sol)Open item page →

The self-similar Riemann problem

Definition

Let n=1 and f∈C1(R) (Scalar conservation laws, fluxes and Cauchy data). The Riemann problem for ut+f(u)x=0 prescribes the two-state initial datum u0(x)={uL,x<0,uR,x>0,uL,uR∈R. One looks for self-similar solutions u(t,x)=v(x/t), t>0, that is, solutions invariant under the scaling (t,x)↦(λt,λx), λ>0; such a u is determined by the single function v ⁣:R→R and is constant on each ray x=ξt.

The initial condition is read as the strong local L1 trace u(t,⋅)→u0 as t↓0 (Kruzhkov entropy solutions), and any jump or corner of v occurs on a ray; admissibility is the entropy condition of Kruzhkov entropy solutions, not a further restriction on the self-similar ansatz. Self-similarity is an ansatz to be justified by the uniqueness theorem rather than an additional hypothesis. No choice principle occurs.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The Riemann solver for a strictly convex flux

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the analytic prerequisites used below.

Let f∈C2(R) be strictly convex and let uL,uR∈R. The unique Kruzhkov entropy solution of the Riemann problem (The self-similar Riemann problem) is:

(i) if uL>uR, the shock u(t,x)={uL,x<st,uR,x>st,s=f(uR)−f(uL)uR−uL;

(ii) if uL<uR, the centred rarefaction u(t,x)={uL,x/t≤f′(uL),(f′)−1(x/t),f′(uL)<x/t<f′(uR),uR,x/t≥f′(uR).

Here f′ is continuous and strictly increasing, so its inverse on [f′(uL),f′(uR)] is continuous; it need not be differentiable, and the rarefaction may have a cusp when f′′ vanishes. Both profiles satisfy the weak conservation law, all Kruzhkov entropy inequalities, and the strong local L1 initial trace. Uniqueness in the bounded Kruzhkov class follows from Uniqueness, comparison and order preservation of entropy solutions. If uL=uR, the constant solution is the unique one (Scalar conservation laws, fluxes and Cauchy data, Kruzhkov entropy solutions).

Facts & Assumptions

Given: Countable Choice, a strictly convex flux f∈C2(R), states uL,uR∈R, the Riemann datum u0, and the self-similar profiles of the statement.

[F1]

The interior weak equation and the self-similar ansatz: the interior distributional equation is equivalent to ∫ΠT(uφt+f(u)φx)=0 for every φ∈Cc∞(ΠT), and a self-similar solution has the form u(t,x)=U(x/t), constant along rays (Scalar conservation laws, fluxes and Cauchy data, The self-similar Riemann problem).

[F2]

Interface computation at a single jump: for a piecewise C1 function with one interface and speed s, the weak residual against a test function supported near the interface equals −∫Γφ ([u]νt+[f(u)]νx)dS; in one dimension with the graph x=st, [u]νt+[f]νx=([f]−s[u])/1+s2. Thus Rankine--Hugoniot s[u]=[f] makes the residual vanish there, and across a continuous interface (equal traces of u, hence of f(u)) the contribution vanishes identically (The Rankine--Hugoniot jump condition in space--time normal form).

[F3]

Chord criterion at a jump: a piecewise C1 weak solution with a single nontrivial jump of speed s satisfies the entropy inequality for every convex C2 pair if and only if F(z)(u+−u−)≥0 for all z between the states, where F(z)=f(z)−f(u−)−s(z−u−) (The convex entropy condition for a single shock is the chord condition).

[F4]

Strict convexity: f′ is strictly increasing by the secant argument in The Lax shock inequalities for convex scalar laws, so f′(uL)<f′(uR) when uL<uR and the inverse (f′)−1 ⁣:[f′(uL),f′(uR)]→[uL,uR] is continuous, strictly increasing; the graph of f lies strictly below every chord on the interior of its interval (Convex and strictly convex functions on Euclidean convex sets, A differentiable function on an open interval is convex if and only if its derivative is nondecreasing).

[F5]

Calculus and regularization: the chain rule and algebra of derivatives compute the classical residual of a C1 self-similar profile; for fϵ(r)=f(r)+ϵ2r2 one has fϵ′′≥ϵ, so fϵ′ is a C1 diffeomorphism of [uL,uR] onto its image and ψϵ=(fϵ′)−1 is C1 on [fϵ′(uL),fϵ′(uR)] by the inverse function theorem; on the fixed interval [uL,uR] the derivatives fϵ′ converge uniformly to f′ as ϵ↓0 (The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c), Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, The Euclidean inverse function theorem).

[F6]

Limits: dominated convergence and uniform convergence on compact sets justify passing to the limit in the weak and entropy test integrals, and Lp membership is a property of equivalence classes (Dominated convergence, The space Lp(μ) as the quotient by null functions).

Proof

technique · direct
1.1F1F2given

The shock profile solves the equation and has the right trace. Assume uL>uR and let u be the profile (i) with speed s=(f(uR)−f(uL))/(uR−uL), so s[u]=[f] with [h]=h(uR)−h(uL): this is the Rankine--Hugoniot condition. By [F2] the weak residual of a piecewise constant profile with a single jump reduces to the interface integral −∫φ(t,st)([f]−s[u]) dt, which vanishes; hence u is a distributional weak solution. Moreover u(t,⋅) differs from u0 only on the interval between 0 and st, whose length is ∣s∣t, and ∣u−u0∣≤∣uL−uR∣ there, so ∫K∣u(t,x)−u0(x)∣ dx≤∣uL−uR∣ ∣s∣t→0 as t↓0: the strong local L1 trace holds.

1.2F2F5

Regularized fans are weak solutions with vanishing entropy production. For ϵ>0 put fϵ(r)=f(r)+ϵ2r2, let ψϵ=(fϵ′)−1 on [fϵ′(uL),fϵ′(uR)], and define Uϵ by the same three-branch formula with ξLϵ=fϵ′(uL), ξRϵ=fϵ′(uR), and uϵ(t,x)=Uϵ(x/t). Each branch is C1 by [F5], and on the open middle region the chain rule gives utϵ+fϵ(uϵ)x=1t ψϵ′(ξ)(fϵ′(ψϵ(ξ))−ξ)=0 since fϵ′(ψϵ(ξ))=ξ; on the outer regions the profile is constant, so the residual vanishes pointwise there as well. At the two interfaces the traces of uϵ, hence of fϵ(uϵ), agree from both sides, so by [F2] no interface term arises and uϵ is a distributional weak solution of ut+fϵ(u)x=0. The same computation applied to a convex C2 pair (η,qϵ) with qϵ′=η′fϵ′ gives ∂tη(uϵ)+∂xqϵ(uϵ)=η′(uϵ)(utϵ+fϵ′(uϵ)uxϵ)=0 on each open branch and continuous traces η(uϵ),qϵ(uϵ) at the interfaces, so the entropy residual is identically 0 for every convex C2 pair.

2.1F3F4step 1.1

The shock is entropic. With u−=uL, u+=uR and speed s, the function F(z)=f(z)−f(uL)−s(z−uL) satisfies F(uL)=F(uR)=0, and by strict convexity [F4] the graph of f lies strictly below the chord through (uL,f(uL)), (uR,f(uR)) on (uR,uL); that chord has slope s, so F(z)<0 for z∈(uR,uL). Since uR−uL<0, F(z)(uR−uL)≥0 on [uR,uL], and the chord criterion [F3] gives the entropy inequality for every convex C2 pair.

2.2F1F4F5

The rarefaction profile and its inverse. Assume uL<uR and set ξL=f′(uL)<ξR=f′(uR); by [F4] the inverse ψ=(f′)−1 ⁣:[ξL,ξR]→[uL,uR] is continuous and strictly increasing. Define U(ξ)=uL for ξ≤ξL, U(ξ)=ψ(ξ) for ξL≤ξ≤ξR, and U(ξ)=uR for ξ≥ξR, and set u(t,x)=U(x/t) for t>0. Then u takes values in [uL,uR], differs from u0 only between min⁡{0,ξL}t and max⁡{0,ξR}t, an interval of length (max⁡{0,ξR}−min⁡{0,ξL})t, and satisfies the strong local L1 initial trace by the same estimate as in step 1.1.

2.3F6step 1.2

Passage to the limiting fan. As ϵ↓0, fϵ′→f′ uniformly on [uL,uR] by [F5], so the clamped inverse profiles converge uniformly on R. Indeed, if D=max⁡{∣uL∣,∣uR∣}, then ∣fϵ′−f′∣≤ϵD on the state interval, so U(ξ−ϵD)≤Uϵ(ξ)≤U(ξ+ϵD). The extended continuous profile U is uniformly continuous (it is constant outside a compact interval), yielding sup⁡ξ∣Uϵ−U∣→0; also fϵ→f and qϵ→q uniformly on the compact range [uL,uR], where q′=η′f′ is the flux of the same convex pair for f. Passing to the limit in the weak residual: ∣∫ΠT(Uϵ(x/t)−U(x/t))φt∣≤sup⁡∣Uϵ−U∣ ∥φt∥1→0 and ∥fϵ(Uϵ)−f(U)∥∞≤∥fϵ−f∥∞,[uL,uR]+Lip(f)sup⁡∣Uϵ−U∣→0, so ∫ΠT(Uφt+f(U)φx)=0 and u is a weak solution. The entropy residual passes similarly, so for every convex C2 pair (η,q) with q′=η′f′ one has ∫ΠT(η(u)φt+q(u)φx)≥0 for every nonnegative φ∈Cc∞(ΠT).

3.1F6step 2.1step 2.3

Kruzhkov pairs by smoothing. Fix k∈R and let ηδ(r)=(r−k)2+δ2−δ, a smooth convex function with 0≤ηδ(r)≤∣r−k∣ and ∣ηδ(r)−∣r−k∣∣≤δ, and let qδ(s)=∫ks(ηδ)′(r)f′(r) dr, so that (qδ)′=(ηδ)′f′. Since z↦sgn⁡(z−k) is bounded and f′ is continuous on the compact range of the profiles, dominated convergence gives qδ(s)→∫kssgn⁡(r−k)f′(r) dr=sgn⁡(s−k)(f(s)−f(k))=qk(s) uniformly for s in the profile range; and ηδ→ηk uniformly there. Applying the entropy inequality of step 2.1 (shock case, via [F3]) or of step 2.3 (rarefaction case) to (ηδ,qδ) and passing to the limit using uniform convergence and φ∈Cc∞ gives ∫ΠT(ηk(u)φt+qk(u)φx)≥0 for every nonnegative φ; hence both profiles satisfy all Kruzhkov entropy inequalities.

4.1step 1.1step 2.1step 2.2step 1.2step 2.3step 3.1∎

The constant case and uniqueness. If uL=uR, the constant u≡uL is a distributional weak solution with the exact trace, and its entropy production vanishes, so it is a Kruzhkov entropy solution; any bounded Kruzhkov entropy solution with the same constant datum equals it by order preservation applied in both directions. In the cases (i) and (ii), the profiles are bounded Kruzhkov entropy solutions with datum u0 by steps 1.1, 2.1, 2.2–2.3 and 3.1, and any bounded Kruzhkov entropy solution with datum u0 coincides with the profile almost everywhere on ΠT by Uniqueness, comparison and order preservation of entropy solutions applied in both directions. This proves existence, uniqueness, and the asserted profile in each case.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Oleinik's one-sided estimate characterizes bounded entropy solutions

Statement

Assume Countable Choice and Dependent Choice, as used by the mollification and vanishing-viscosity compactness interfaces. Let T>0, f∈C2(R), and let I⊂R be a bounded closed interval with f′′≥κ>0 on I. Let u0∈L∞(R) and let u∈L∞((0,T)×R) be a distributional weak solution of ut+∂xf(u)=0, with u0 and the essential range of u in I, and with the strong local L1 initial trace u0 in Kruzhkov entropy solutions. The following are equivalent: (i) u is a Kruzhkov entropy solution with that trace; (ii) for almost every t∈(0,T) and almost every pair x<y, u(t,y)−u(t,x)≤y−xκt; (iii) for almost every t∈(0,T), Dxu(t,⋅)≤(κt)−1L1 in distributional order. The initial trace and weak equation are hypotheses of the equivalence; the slope bound alone is not a definition of an entropy solution. For piecewise C1 solutions, the bound in particular excludes upward jumps, and the convex chord criterion makes the remaining shocks entropy-admissible. The state-slope constant 1/κ requires uniform convexity on the solution range.

Facts & Assumptions

Given: Countable and Dependent Choice, T>0, f∈C2(R), a bounded closed interval I with f′′≥κ>0 on I, a bounded weak solution u with datum u0∈L∞, both taking values in I, and a nonnegative test function φ in the arguments below.

[F1]

The weak equation and the trace: ∫ΠT(uφt+f(u)φx)=0 for every φ∈Cc∞(ΠT), and u has the strong local L1 trace u0; a Kruzhkov entropy solution is a bounded weak solution satisfying ∂t∣u−k∣+∂xqk(u)≤0 for all k, qk(s)=sgn⁡(s−k)(f(s)−f(k)) (Kruzhkov entropy solutions).

[F2]

The viscous construction supplies bounded C1,2 solutions for C2 fluxes and smooth compactly supported data (The viscous scalar Cauchy problem with smooth data has a global classical solution). Its heat-potential cancellation estimates apply on every positive-time strip (The heat evolution Ht of initial data, Normalisation, parabolic scaling, heat equation and derivative bounds for the heat kernel). For a smooth flux, these also make p=vx classical: p satisfies the differentiated divergence equation with source g′(v)p, which is parabolically Hölder by the gradient estimate in the construction. Applying its second-kernel cancellation first gives Hölder px; then the source −g′(v)px−g′′(v)p2 is Hölder, and the nondifferentiated heat-potential estimate gives p∈C1,2 locally. Chain and product rules are The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c) and Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0. Smooth flux regularization is used before this differentiation.

[F3]

Vanishing viscosity: for each smooth compactly supported datum the viscous solutions have a subsequence converging in Lloc1(ΠT) and almost everywhere to the unique bounded Kruzhkov entropy solution of that datum, and the entropy inequalities pass to the limit (Existence of bounded Kruzhkov entropy solutions, Vanishing-viscosity families are locally precompact in L1).

[F4]

Local contraction: two bounded Kruzhkov entropy solutions with data in Lloc1∩L∞ satisfy ∫B(x0,R−L′t)∣v−w∣≤∫B(x0,R)∣v0−w0∣ for almost every t with L′t<R, where L′ is a Lipschitz constant of the (shifted) flux on the common range (Local L1 contraction for two entropy solutions).

[F5]

Mollification and distributional calculus: convolutions with radial mollifiers are smooth and converge in L1 (or Lloc1) to the original function; derivatives may be taken inside the convolution; approximate identities converge in Lp; distributional derivatives commute with convolution against test functions; almost-everywhere convergence is available along subsequences of L1-convergent sequences (A radial mollifier family in Rn, A unit-mass smooth bump generates an L1 approximate identity, Convolution with a mollifier is smooth, and derivatives pass under the integral sign, Distributional derivatives commute with test-function convolution, Every L1 approximate identity converges to the identity in Lp for 1≤p<∞, Convergence in L^1(mu) has an almost-everywhere convergent subsequence, Dominated convergence, Fubini's theorem for L^1 functions on a sigma-finite product, The space Lp(μ) as the quotient by null functions).

[F6]

Piecewise C1 interpretation: at a single shock satisfying the Rankine--Hugoniot condition, the entropy inequalities for all convex C2 pairs hold iff the chord residual F(z)=f(z)−f(u−)−s(z−u−) satisfies F(z)(u+−u−)≥0 between the states; in particular only compressive jumps u−>u+ are admissible (The convex entropy condition for a single shock is the chord condition).

Proof

technique · direct
1.1F2

Shifting and positive-strip gradient bounds. Fix c∈I and put J:=I−c, which contains 0. Set f~(s):=f(s)−f(0) and g(v):=f~(v+c)−f~(c)=f(v+c)−f(c). Then g(0)=0 and g′′≥κ on J. For a smooth compactly supported datum v0s with values in J, apply [F2] to flux g and datum v0s, obtaining vε with values in J. The unshifted profile Uε:=vε+c has datum v0s+c, equals c outside a compact set, and solves the equation with normalized flux f~; thus the compact datum to which [F2] is applied is v0s, not the generally noncompact function v0s+c or a nonzero-tail shift. The Oleinik slope is unchanged by adding c. Set Mv:=∥v0s∥∞ and G0:=sup⁡v∈J∣g′(v)∣<∞. By variation of constants, for 0<d<t, vε(t)=Hε(t−d)vε(d)−∫dt∂xHε(t−s)g(vε(s)) ds. Indeed, differentiating Hε(t−s)vε(s) for d<s<t and using the viscous equation gives −∂xHε(t−s)g(vε(s)); the integral converges at s=t because ∥∂xΓεr∥1≤Cεr−1/2. The Gaussian derivative bounds [F2] also give ∥∂xΓεr(⋅+h)−∂xΓεr∥1≤min⁡(Cεr−1/2,Cε∣h∣r−1),∥∂r∂xΓεr∥1≤Cεr−3/2. Splitting the spatial integral at r=∣h∣2 and the time integral at r=∣t−t′∣ in the restarted identity shows that, for every 0<α<1 and every τ>0, ∣vε(t,x+h)−vε(t,x)∣≤Cτ,T∣h∣α,∣vε(t,x)−vε(t′,x)∣≤Cτ,T∣t−t′∣1/2 on [τ/2,T]×R. Fix τ>0, put d=τ/2, and differentiate the restarted identity for t≥τ, first omitting its final η-length of integration and then letting η↓0. Since ∫∂xxΓεr=0, cancellation gives p(t,x)=∂xHε(t−d)vε(d,x)−∫0t−d∫R∂xxΓεr(z)(g(vε(t−r,x−z))−g(vε(t,x))) dz dr. The range bound and the positive-strip estimates imply ∣g(vε(t−r,x−z))−g(vε(t,x))∣≤G0Cτ,T(∣z∣α+r1/2). By scaling, ∫∣∂xxΓεr(z)∣(∣z∣α+r1/2) dz≤Cε,α(r−1+α/2+r−1/2), which is integrable at r=0; the first term is bounded by Cε,τMv because t−d≥τ/2. Thus sup⁡(t,x)∈[τ,T]×R∣p(t,x)∣<∞. This bound is global in space and holds on every positive-time strip.

1.2F5

Equivalence of the pointwise and distributional forms. Fix t∈(0,T), assume first that (ii) holds at this t, and put w(x)=u(t,x)−x/(κt), so that w(y)≤w(x) for almost every x<y. For n≥1 set a:=1/n and wn(x)=n∫xx+aw(z) dz. If x<x′, put d:=x′−x>0. When d≤a, cancellation of the overlap gives wn(x)−wn(x′)=n(∫xx+dw(z) dz−∫x+ax+a+dw(ζ) dζ)≥0, because the two intervals have equal length d and every interior pair z∈[x,x+d], ζ∈[x+a,x+a+d] satisfies z<ζ; explicitly, d(∫xx+dw−∫x+ax+a+dw)=∫[x,x+d]×[x+a,x+a+d](w(z)−w(ζ)) dz dζ≥0 by the assumed a.e. pair inequality. When d≥a, the intervals [x,x+a] and [x′,x′+a] are ordered and have equal length, so wn(x)−wn(x′)=n(∫xx+aw(z) dz−∫x′x′+aw(ζ) dζ)≥0, since a(∫xx+aw−∫x′x′+aw)=∫[x,x+a]×[x′,x′+a](w(z)−w(ζ)) dz dζ≥0. Thus wn is nonincreasing and Dxwn≤0 distributionally. Since wn→w in Lloc1 as n→∞, distributional differentiation passes to the limit, so Dxw≤0, which is (iii) at this t. Conversely, assume (iii) at some t and let ρδ be a nonnegative spatial mollifier; then wδ=w∗ρδ is smooth with Dxwδ=(Dxw)∗ρδ≤0, so wδ is nonincreasing and wδ(y)≤wδ(x) for all x<y; since wδ→w in Lloc1, passing to an almost-everywhere convergent subsequence gives the two-point inequality of (ii) at this t for almost every pair. Hence (ii) and (iii) hold for the same full-measure set of times, proving the equivalence.

1.3F5

Mollification commutator. Assume (iii). Fix a nonnegative test function φ supported in [a,b]×[−R,R] with 0<a<b<T, and let ρδ be a nonnegative unit-mass mollifier on R2 supported in the ball of radius δ<12min⁡{a,T−b}. Extend u boundedly to all of R2 by a fixed value in I outside (0,T)×R, and set uδ=u∗ρδ, Fδ=f(u)∗ρδ, rδ=f(uδ)−Fδ. By [F5], distributional derivatives commute with convolution, so on a neighbourhood of the support of φ the identity ∂tuδ+∂xFδ=0 holds; since the mollification averages only over times ≥a−δ>a/2, the distributional bound (iii) gives ∂xuδ≤2/(κa) there; and uδ still takes values in I. The tangent inequality f(s)≥f(uδ)+f′(uδ)(s−uδ) for convex f, averaged against ρδ, gives rδ≤0. Finally, localizing u and f(u) by a cutoff equal to 1 on a slightly larger compact set and applying approximate-identity convergence in L1 [F5] gives uδ→u, Fδ→f(u) and rδ→0 in L1 on the support of φ, the last two also using that f is Lipschitz on the bounded interval I.

2.1F2F5step 1.1

A classical barrier after flux regularization. Take smooth normalized gl converging to g in C2 on the compact state interval J, with gl(0)=0 and gl′′≥κl:=κ−1/l>0 there (mollify g at sufficiently small scales). Choose 0<εl≤1 tending to zero, and let vl have the fixed smooth datum of step 1.1 and flux gl. Its range lies in J. By [F2], p=vxl is classical at positive times, is bounded globally on positive strips by step 1.1, and satisfies pt+gl′(vl)px+gl′′(vl)p2=εlpxx. Fix τ>0, put q(t)=1/(κl(t−τ)), b=gl′(vl) and W=p−q. Where W>0, (∂t+b∂x−εl∂xx)W<0 since p>q>0. If Gl=sup⁡J∣gl′∣, choose A>Gl+2εl and Φ=eA(t−τ)(1+x2), so the same operator applied to Φ is strictly positive. For ρ>0, Z=W−ρΦ is negative at some t0>τ sufficiently close to τ, by the positive-strip bound on p, and negative on the sides of a sufficiently large rectangle. A positive maximum on that rectangle would have Zt≥0, Zx=0, Zxx≤0, contradicting the strict operator inequality there. Thus Z≤0. Let ρ↓0 and then τ↓0 to get vxl≤1/(κlt); integrating in x gives vl(t,y)−vl(t,x)≤(y−x)/(κlt). This uses a classical maximum argument, with no Sobolev positive-part test.

2.2F5step 1.3

The entropy production tends to a nonpositive limit. For any convex η∈C2 and q′=η′f′, the chain and product rules applied to ∂tuδ+∂x(f(uδ)−rδ)=0 give ∂tη(uδ)+∂xq(uδ)=∂x(η′(uδ)rδ)−η′′(uδ)(∂xuδ)rδ≤∂x(η′(uδ)rδ)+2∥η′′∥L∞(I)κa(−rδ), using rδ≤0 and ∂xuδ≤2/(κa) from step 1.3. Tested against the nonnegative φ, the first term is bounded by ∥η′∥L∞(I)∥φx∥∞∥rδ∥L1(supp⁡φ) after an integration by parts, and the second by 2∥η′′∥∞κa∥φ∥∞∥rδ∥L1(supp⁡φ); both tend to 0 as δ↓0. Since ∂tη(uδ)+∂xq(uδ)→∂tη(u)+∂xq(u) distributionally by [F5] (local L1 convergence of uδ and continuity of η,q), the limit satisfies ⟨∂tη(u)+∂xq(u),φ⟩≤0 for every nonnegative test function supported in (0,T)×R.

3.1F3F4F5step 2.1

Passage to the smooth-datum entropy solution. The varying-flux family vl of step 2.1 satisfies the common range and derivative bounds of the compactness lemma [F3]. Extract a locally L1 and almost-everywhere convergent subsequence. The existence proof passes its weak and entropy identities to the limit because gl→g in C1 on the range; the uniform local time modulus supplies the initial trace. Uniqueness identifies the limit as the entropy solution v for g with this smooth datum. Fubini gives slicewise convergence at almost every time, and passage to the bound in step 2.1, with κl→κ, gives v(t,y)−v(t,x)≤(y−x)/(κt) for almost every time and almost every pair x<y.

4.1F4F5step 3.1

Approximation for general data: (i) implies (ii). Now let u be the given entropy solution with datum u0 and range in I, and fix c∈I. For j≥1 set wj=((u0−c)1[−j,j])∗ρδj, where ρδj is a nonnegative mollifier of radius δj↓0: then wj∈Cc∞, its values lie in I−c (a convex combination of values of u0−c), and wj→u0−c in Lloc1(R). Let vj be the entropy solution with datum wj for the flux g: by step 3.1 each vj satisfies the two-point estimate, and by [F4] applied to vj and u~=u−c, the differences converge to 0 in L1 on every compact cylinder; a diagonal subsequence converges almost everywhere on ΠT. Passing the two-point estimate to that almost-everywhere limit proves (ii) for u.

5.1F1F5step 4.1step 1.2step 2.2

Kruzhkov pairs and conclusion of (iii) implies (i). For k∈R and m≥1 put ηm(s)=(s−k)2+m−2 and qm(s)=∫ksηm′(z)f′(z) dz, so (ηm,qm) is a smooth convex entropy pair. Step 2.2 gives ∂tηm(u)+∂xqm(u)≤0 for every m; letting m→∞, ηm→∣s−k∣ and qm→sgn⁡(s−k)(f(s)−f(k)) uniformly on the bounded interval I by dominated convergence, so the distributional inequality passes to the limit and every Kruzhkov inequality holds. Together with the weak equation and the strong local L1 trace (hypotheses), u is a Kruzhkov entropy solution, proving (iii) implies (i); the reverse implication (i) implies (ii) is step 4.1, and the equivalence of (ii) and (iii) is step 1.2.

6.1F6step 4.1step 1.2step 5.1∎

The piecewise C1 remark. If u is piecewise C1 with a single jump at a curve x=s(t) and satisfies the hypotheses, then (ii) forces the right trace not to exceed the left trace across an upward jump: taking x↑s(t), y↓s(t) in the two-point inequality and letting x,y→s(t) gives u+−u−≤0, so an upward jump u+>u− is excluded; for the remaining jumps with u−>u+ the Rankine--Hugoniot condition and the chord criterion [F6] make them entropy-admissible. This shows how the slope bound encodes admissibility in the piecewise smooth class, while the equivalence itself was proved for all bounded weak solutions.

CorollaryStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The Lax shock inequalities for convex scalar laws

Statement

Let f∈C2(R) be strictly convex. Consider a nontrivial one-dimensional jump from the left trace u− to the right trace u+ across x=s(t), with traces as in Piecewise smooth shocks and one-sided traces, satisfying the Rankine--Hugoniot condition of The Rankine--Hugoniot jump condition in space--time normal form. The jump is Kruzhkov entropy-admissible (Kruzhkov entropy solutions) if and only if u−>u+. In that case its speed is s′=f(u+)−f(u−)u+−u−, and it satisfies the Lax shock inequalities f′(u+)≤s′≤f′(u−). In particular, every nontrivial entropy-admissible jump is compressive; no admissible jump increases the state across the shock (Convex and strictly convex functions on Euclidean convex sets).

Facts & Assumptions

Given: a strictly convex f∈C2(R), a nontrivial single-jump piecewise C1 weak solution with left trace u−, right trace u+ across x=s(t), and speed σ=s′ satisfying the Rankine--Hugoniot condition.

[F1]

Rankine--Hugoniot and jump setup: the interface is the graph x=s(t) with minus side x<s(t) and plus side x>s(t), and σ(u+−u−)=f(u+)−f(u−), i.e. σ=(f(u+)−f(u−))/(u+−u−) since the jump is nontrivial (Piecewise smooth shocks and one-sided traces, The Rankine--Hugoniot jump condition in space--time normal form).

[F2]

Chord criterion: with F(z)=f(z)−f(u−)−σ(z−u−), the jump satisfies the Kruzhkov entropy inequalities for all convex entropy pairs if and only if F(z)(u+−u−)≥0 for every z between u− and u+; this is the notion of entropy admissibility at a single jump (The convex entropy condition for a single shock is the chord condition, Kruzhkov entropy solutions).

[F3]

Strict convexity: for a differentiable strictly convex f, the graph lies strictly below every chord on the interior of its interval; the derivative f′ is strictly increasing: it is nondecreasing by the cited theorem, and equality at a<b would make it constant on [a,b], so FTC would make f affine there, contradicting strict convexity; and for a<b the difference quotients satisfy f′(a)≤f(b)−f(a)b−a≤f′(b) with strict inequalities throughout, while the mean value theorem gives f(b)−f(a)b−a=f′(c) for some c∈(a,b) (Convex and strictly convex functions on Euclidean convex sets, A differentiable function on an open interval is convex if and only if its derivative is nondecreasing, The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a), The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)).

Proof

technique · direct
1.1F1F2F3

The forward case is never admissible. Suppose u−<u+. The chord criterion of [F2] requires F(z)≥0 for all z∈[u−,u+]. By strict convexity [F3] the graph of f lies strictly below the chord through (u−,f(u−)) and (u+,f(u+)) on (u−,u+), and that chord is z↦f(u−)+σ(z−u−) because its slope is σ; hence F(z)<0 for every z∈(u−,u+), contradicting the criterion. So a nontrivial Rankine--Hugoniot jump with u−<u+ is not entropy-admissible.

2.1F1F2F3

The backward case is admissible. Suppose u−>u+. On the interval between the states, strict convexity gives F(z)<0 for u+<z<u− and F(u+)=F(u−)=0. Since u+−u−<0, the product F(z)(u+−u−) is positive for interior z and vanishes at the endpoints, so the chord criterion of [F2] holds and the jump is entropy-admissible. Together with step 1.1 this shows that a nontrivial Rankine--Hugoniot jump is entropy-admissible if and only if u−>u+; in particular no admissible jump increases the state.

3.1F1F3step 1.1step 2.1∎

The Lax inequalities. Assume u−>u+ and write a=u+<b=u−. By [F1], σ=f(b)−f(a)b−a=f(u+)−f(u−)u+−u−, which is the stated speed. By the mean value theorem [F3] there is c∈(a,b) with f′(c)=σ, and since f′ is strictly increasing, f′(u+)=f′(a)<f′(c)=σ<f′(b)=f′(u−); a fortiori f′(u+)≤σ≤f′(u−), the Lax shock inequalities.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The Hamilton--Jacobi correspondence in one dimension

Statement

Assume Countable Choice and Dependent Choice (The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) for the heat-kernel, L1 completeness and vanishing-viscosity extraction interfaces used below. Let f∈C2(R) be strictly convex and superlinear, with f(p)/∣p∣→∞ as ∣p∣→∞.

(i) Let U0∈C0,1(R)∩L∞(R) and let u0=U0′ be its a.e. derivative. The Hopf--Lax function V(t,x)=QtU0(x)=inf⁡y∈R{U0(y)+tL ⁣(x−yt)},t>0,V(0,x)=U0(x), where L=f∗, is the unique viscosity solution of Vt+f(Vx)=0 with initial datum U0 among functions bounded and uniformly continuous on [0,S]×R for every finite S>0 (The Hamilton--Jacobi Cauchy problem and its classical solutions, Discontinuous viscosity solutions through the two envelopes). Its a.e. spatial derivative v=Vx is the bounded Kruzhkov entropy solution of vt+∂xf(v)=0 with initial datum u0.

(ii) Conversely, let u0∈L1(R)∩L∞(R) have compact support and let u be its bounded Kruzhkov entropy solution, using the strong local L1 initial trace. With U0(x)=∫−∞xu0(y) dy,U(t,x)=∫−∞xu(t,y) dy−tf(0), the function U is the unique viscosity solution of Ut+f(Ux)=0 with datum U0 in the same finite-slab class as in (i), and Ux=u almost everywhere. In the compactly supported datum class of (ii), differentiation and the normalized primitive are inverse correspondences (Kruzhkov entropy solutions, Existence of bounded Kruzhkov entropy solutions).

Facts & Assumptions

Given: Countable and Dependent Choice, a strictly convex superlinear C2 flux f, its conjugate L, and the two datum classes in the statement.

[F2]

The normalized flux g=f−f(0) gives the same conservation law. For smooth compactly supported data its viscous solutions are mild classical solutions, obey the range and L1 bounds, and are locally precompact in space--time L1, with limits continuous into local L1. The existence proof passes their weak and entropy identities to the unique entropy solution (The viscous scalar Cauchy problem with smooth data has a global classical solution, Uniform L-infinity, mass and energy bounds for the viscous approximations, Vanishing-viscosity families are locally precompact in L1, Existence of bounded Kruzhkov entropy solutions, Kruzhkov entropy solutions, Distributional weak solutions of the Cauchy problem).

[F3]

The heat kernels have unit mass, solve the heat equation, have Gaussian derivative estimates, and give the heat evolution; smooth cutoffs have derivatives O(R−1) and O(R−2) (The heat evolution Ht of initial data, Normalisation, parabolic scaling, heat equation and derivative bounds for the heat kernel, Explicit compactly supported smooth cutoffs). Fubini, dominated convergence and FTC justify the kernel calculations (Fubini's theorem for L^1 functions on a sigma-finite product, Dominated convergence, The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)).

[F4]

Entropy solutions contract in global L1 at every time for their continuous representatives, and locally on shrinking balls; they are unique in the bounded class. Compactly supported data stay supported in a common bounded interval on every finite horizon (Global L1 contraction from the local estimate, Local L1 contraction for two entropy solutions, Uniqueness, comparison and order preservation of entropy solutions, Finite propagation for scalar conservation laws).

Proof

technique · direct
1.1F1

Conjugate calculus and localization. Strict convexity makes f′ strictly increasing: convex secant inequalities give monotonicity, and equality at two distinct points would make f affine between them. Superlinearity makes its limits ±∞ (a finite derivative bound at either end would bound f linearly there). Thus p(q):=(f′)−1(q) is continuous, and is the unique maximiser of pq−f(p). The inequalities p(q)h≤L(q+h)−L(q)≤p(q+h)h for h>0, and their reversed versions for h<0, show that L′(q)=p(q). If W is M-Lipschitz and y minimises W(y)+tL((x−y)/t), perturb y in each direction and use the Lipschitz bound to obtain ∣L′((x−y)/t)∣≤M. Hence ∣x−y∣≤Ct, where C=max⁡∣p∣≤M∣f′(p)∣. Translating competitors shows that QtW is M-Lipschitz in x. Moreover L(v)−M∣v∣≥−C0, where C0=max⁡{f(M),f(−M)} by the conjugate definition, and the competitor y=x gives QhW−W≤hL(0). These bounds and the semigroup law give a uniform time Lipschitz bound on finite horizons. Also L(q)≥−f(0) for all q, with equality at q=f′(0), so Qt0=−tf(0). Sup contraction therefore gives ∥QtU0+tf(0)∥∞≤∥U0∥∞, proving boundedness on each finite slab, with no time-uniform bound asserted.

1.2F2F3

Viscous primitives for a smooth compact datum. Let a∈Cc∞, P0(x)=∫−∞xa, and let uε be [F2]'s viscous solution for g=f−f(0). Define Wε(t)=HεtP0−∫0tHε(t−s)f(uε(s)) ds. Differentiating in x gives exactly the mild identity for uε, so Wxε=uε. Heat-potential cancellation as in the viscous construction makes Wε classical at positive times, and its equation is Wtε+f(Wxε)=εWxxε. At x→−∞ its initial heat term tends to zero, while the integral tends to tf(0): g(uε(s))∈L1, and convolution of an L1 function with a bounded Gaussian tends to zero at spatial infinity; boundedness dominates the finite time integral. Thus Wε(t,x)=∫−∞xuε(t,y)dy−tf(0). Testing the smoothed ∣uε∣ balance with an exterior cutoff, then removing a second outer cutoff, gives ∫∣uε(t)∣(1−χR)≤∫∣a∣(1−χR)+CT(R−1+R−2)∥a∥1(0≤t≤T, 0<ε≤1). This is the cutoff calculation of the L1 bound in [F2], with ∣g(u)∣≤C∣u∣ and [F3]'s derivative bounds. It supplies uniform tails.

2.1F2F3step 1.2

A time modulus for the viscous primitives. The Gaussian convolution identity HrHs=Hr+s follows by completing the square in the kernel product and using unit mass and Fubini. Applying it to the definition in step 1.2 gives Wε(t+h)=HεhWε(t)−∫tt+hHε(t+h−s)f(uε(s)) ds. Put M=∥a∥∞ and B=max⁡∣z∣≤M∣f(z)∣. Since Wxε=uε, these primitives are M-Lipschitz in x, including at t=0. Gaussian scaling gives ∫∣z∣Γ(z,r) dz=Cr, with C<∞ by the Gaussian bound. Hence unit mass gives ∥HrWε(t)−Wε(t)∥∞≤CMr, and heat contraction bounds the time integral by Bh. Thus ∥Wε(t+h)−Wε(t)∥∞≤CMh+Bh for 0≤t<t+h≤T, uniformly in 0<ε≤1.

3.1F2F3step 1.2step 2.1

Uniform convergence of primitives. By [F2], choose a subsequence uε→u locally in space--time L1, where u is the entropy solution of datum a, continuous into local L1. A further subsequence converges on almost every time slice locally in L1. The uniform tails of step 1.2 pass to these slices by monotone exhaustion, and to every time by local continuity on bounded annuli followed by exhaustion. They imply u(t)∈L1 with uniformly small tails; local continuity then gives global L1 continuity on [0,T]. Moreover ∫0T∥uε(t)−u(t)∥1dt→0: the tails are uniformly small outside large intervals, the compact space--time convergence handles times away from 0,T, and the bound ∣uε∣,∣u∣≤M controls the remaining small time intervals on the fixed spatial interval. For P(t,x)=∫−∞xu(t,y) dy−tf(0), the primitive formula of step 1.2 gives ∥Wε(t)−P(t)∥∞≤∥uε(t)−u(t)∥1, so the integral in time of the left side tends to zero. The function P is continuous in time in the supremum norm by global L1 continuity. Together with the common modulus of step 2.1, this implies uniform convergence on [0,T]×R: a discrepancy of size d>0 at any time would persist with size at least d/2 on a one-sided interval of length bounded below independently of ε, contradicting that vanishing time integral.

4.1F1F2F3step 1.2step 3.1

The viscosity limit. At a strict local maximum of P−ϕ, with smooth ϕ, step 3.1 gives nearby local maxima of Wε−ϕ. The classical equation in step 1.2 gives ϕt+f(ϕx)≤εϕxx there; passing to the limit proves the subsolution inequality. Local minima give the supersolution inequality. Adding a fourth-power distance term makes a contact strict without changing its first derivatives; approximation in C1 on a compact contact neighbourhood reduces C1 tests to smooth tests. Thus P is a viscosity solution with initial datum P0. It is bounded by ∥a∥1+T∣f(0)∣, spatially M-Lipschitz, and uniformly continuous in time on [0,T] by step 3.1. Uniqueness in [F1] gives P=QtP0.

5.1F1F4F5step 4.1

Compactly supported bounded data. For compactly supported a∈L1∩L∞, choose smooth compactly supported aj→a in L1 with ∥aj∥∞≤∥a∥∞, using [F5]. Their primitives converge uniformly since sup⁡x∣∫−∞x(aj−a)∣≤∥aj−a∥1. By [F4], their entropy solutions converge uniformly in time in L1 to the solution of datum a; their normalized primitives therefore converge uniformly as well. The Hopf--Lax sup contraction in [F1] passes the identity of step 4.1 to P(t)=QtP0. In particular Px=u a.e. by [F5]. This proves (ii), including the normalization −tf(0).

6.1F1F4F5step 1.1step 5.1

A bounded Lipschitz primitive with nonintegrable derivative. Let U0 be as in (i), and set U0R(x)=U0(max⁡{−R,min⁡{x,R}}). It has the same sup and Lipschitz bounds, and derivative u0R=u01(−R,R) a.e. The difference between U0R and the normalized primitive of u0R is its constant value U0(−R); adding this constant commutes with Hopf--Lax. Thus step 5.1 shows that (QtU0R)x is an entropy solution with datum u0R. Step 1.1 places every minimiser for both U0 and U0R within Ct of x. Consequently QtU0R(x)=QtU0(x) whenever ∣x∣+Ct<R, since all those competitors see identical data. Every compact positive-time cylinder is contained in such a region for large R, so the a.e. derivative v=(QtU0)x is bounded by M and obeys the weak equation and every entropy inequality locally, hence globally. For a compact spatial set, fix R large enough that this equality holds throughout 0≤t≤T on that set; the strong local trace of the compact-data solution supplies the trace of v equal to u0. This proves (i), with entropy uniqueness from [F4].

7.1step 5.1step 6.1∎

Conclusion. Step 6.1 proves the derivative correspondence for the entire bounded Lipschitz primitive class, including nonintegrable derivatives, while step 5.1 proves the normalized primitive correspondence for compactly supported integrable data. In that latter class, a.e. differentiation returns u, and integration from −∞ with the time shift −tf(0) returns the prescribed viscosity potential. All arguments hold on an arbitrary finite horizon; [F1] gives viscosity uniqueness on each such slab, while [F4] gives compatibility of the entropy solutions on overlapping horizons. These are the global solutions and inverse correspondences asserted.

CorollaryStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Mass conservation for compactly supported entropy solutions

Statement

Assume Countable Choice and Dependent Choice (The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) for the heat-kernel, L1 completeness and vanishing-viscosity extraction interfaces used below. Let n≥1, let f ⁣:R→Rn be C1 with f(0)=0, let u0∈L1(Rn)∩L∞(Rn) be compactly supported, and let u be the entropy solution of Existence of bounded Kruzhkov entropy solutions. Then ∫Rnu(t,x) dx=∫Rnu0(x) dx for almost every t≥0, and the function t↦∫u(t,x) dx is constant on [0,∞) after choosing the continuous representative of the L1 orbit (Open ball, closed ball and sphere in a metric space, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable and Dependent Choice, n≥1, a C1 flux f with f(0)=0, a compactly supported datum u0∈L1∩L∞ with M=∥u0∥∞, the entropy solution u of Existence of bounded Kruzhkov entropy solutions with its representative in C0([0,T];Lloc1), and a centre x0∈Rn, R>0 with u0=0 almost everywhere outside B(x0,R).

[F1]

Existence and regularity: u is a bounded Kruzhkov entropy solution with ∣u∣≤M almost everywhere, ∥u(t,⋅)∥∞≤M for every t, strong local L1 trace u0, and an Lloc1-continuous representative (Existence of bounded Kruzhkov entropy solutions, Kruzhkov entropy solutions).

[F2]

Finite propagation: with L=sup⁡∣s∣≤M∣f′(s)∣, the solution u and the zero solution (which is a Kruzhkov entropy solution with datum 0) agree outside the cone: u(t,⋅)=0 almost everywhere outside B(x0,R+Lt) for almost every t≥0; the conclusion is an almost-everywhere statement at the level of the L1 classes (Finite propagation for scalar conservation laws, Open ball, closed ball and sphere in a metric space).

[F3]

Weak formulation: ∫ΠT(uφt+f(u)⋅∇φ)=0 for every φ∈Cc∞(ΠT) (Distributional weak solutions of the Cauchy problem, Kruzhkov entropy solutions).

[F4]

There is a smooth compactly supported β with β=1 on a neighbourhood of B‾(x0,R+LT) (Explicit compactly supported smooth cutoffs), and L1 functions are equivalence classes, so pointwise statements on full-measure sets determine the class (The space Lp(μ) as the quotient by null functions).

Proof

technique · direct
1.1F1F2

The cone support holds at every time for the chosen representative. By [F2] there is a full-measure set E⊆(0,∞) with u(t)=0 almost everywhere outside B(x0,R+Lt) for t∈E. Fix T>0, t0∈[0,T] and a compact set K⊆{x:∣x−x0∣>R+Lt0}. The positive distance of K from B‾(x0,R+Lt0) lets us choose tj∈E∩(0,T) with tj→t0 and K⊆{x:∣x−x0∣>R+Ltj} for all j; then u(tj)=0 in L1(K), and the Lloc1-continuity of the representative [F1] gives u(t0)=0 in L1(K). A countable exhaustion of the strict exterior of B(x0,R+Lt0) by compact sets gives u(t0)=0 almost everywhere there; hence for every t0∈[0,T], u(t0) is supported in B‾(x0,R+Lt0)⊆B‾(x0,R+LT).

2.1F3F4step 1.1

Truncated mass balance. By step 1.1, for every t∈[0,T] the function u(t,⋅) vanishes almost everywhere outside the fixed ball B(x0,R+LT); in particular u∈L1 at every time with ∥u(t)∥1≤M vol(B(x0,R+LT)), and the Lloc1 continuity of [F1] is continuity in L1(Rn). Choose β as in [F4] and, for ψ∈Cc∞((0,T)), test [F3] with φ(t,x)=ψ(t)β(x): since u(t,⋅) is supported where β=1 and f(0)=0 implies f(u)=0 wherever u=0, the flux term vanishes and only flat boundary terms contribute. The divergence theorem in the form of the weak identity then gives ∫0Tψ′(t) m(t) dt=0 for the function m(t)=∫Rnu(t,x) dx, that is, m′=0 in the sense of distributions on (0,T).

3.1F1F4step 1.1step 2.1∎

Conclusion. Since u is continuous in L1 on [0,T] and the support lies in the fixed ball, m is continuous on [0,T]; its distributional derivative vanishes on (0,T) by step 2.1 and m(0)=∫u0 by the strong L1 trace. To see constancy directly, convolve m locally in time with a smooth unit-mass bump: its derivative is zero by m′=0 tested against translated kernels, so FTC makes each convolution constant on every interior compact interval. Uniform continuity of m on such intervals makes the convolutions converge uniformly to m, hence m is constant on (0,T) and by continuity at its endpoints. Therefore m(t)=∫u0 for every t∈[0,T]; in particular the equality holds for almost every t and the chosen representative makes t↦∫u(t,x) dx constant on every [0,T]. As T>0 was arbitrary, the claims follow on [0,∞).

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

An additive constant in an entropy flux does not change the entropy inequality

Statement

Let n≥1, f∈C1(R;Rn), let (η,q) be an entropy pair as in Convex entropy--entropy flux pairs, and let qC=q+C with a constant vector C∈Rn. Then for every bounded measurable u ⁣:ΠT→R the distributional inequalities η(u)t+div⁡xq(u)≤0andη(u)t+div⁡xqC(u)≤0 are equivalent in D′(ΠT); the two divergences differ by the zero distribution, because the divergence of a constant vector field vanishes.

Facts & Assumptions

Given: n≥1, an entropy pair (η,q), a constant vector C, a bounded measurable u ⁣:ΠT→R, and a test function φ∈Cc∞(ΠT).

[F1]

The distributional divergence is defined by duality, ⟨div⁡xw,φ⟩=−∫w⋅∇φ, and a distribution is determined by its pairings with test functions; the test-function space is D(ΠT)=Cc∞(ΠT) (Distribution, Distributional derivative, Test function space d of an open set).

[F2]

Entropy pairs and entropy inequalities, including the dependence on the normalisation of the entropy flux, are as in Convex entropy--entropy flux pairs and Kruzhkov entropy solutions; since q is locally Lipschitz and u is bounded, q(u) and qC(u) are locally integrable and their divergences are defined by [F1].

Proof

technique · direct
1.1F1given

For a test function φ, [F1] and the definition of qC give ⟨div⁡xqC(u),φ⟩=−∫qC(u)⋅∇φ=−∫q(u)⋅∇φ−∑i=1nCi∫∂xiφ.

1.2given

Each ∫ΠT∂xiφ dx dt=0: the inner spatial integral vanishes because φ is compactly supported in x, so the function xi↦φ(t,xi,x′) is smooth compactly supported and the fundamental theorem of calculus applies, and the remaining integral over the other variables is finite as φ has compact support.

2.1step 1.1step 1.2F1

By steps 1.1 and 1.2, ⟨div⁡xqC(u),φ⟩=−∫q(u)⋅∇φ=⟨div⁡xq(u),φ⟩ for every test function, hence div⁡xqC(u)=div⁡xq(u) in D′(ΠT).

3.1step 2.1F2∎

Adding the common distribution ∂tη(u) to both sides of step 2.1, the two inequalities ∂tη(u)+div⁡xq(u)≤0 and ∂tη(u)+div⁡xqC(u)≤0 are literally the same distributional inequality, so they are equivalent; in particular the entropy condition does not depend on the additive normalisation of the entropy flux.

Remarks

Consequently the entropy inequality depends only on the pair (η,q) up to the normalisation of q, and statements such as The convex entropy condition for a single shock is the chord condition are independent of the chosen constant.

CorollaryStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The L∞ maximum bound for entropy solutions

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)). Let n≥1, T>0, let f be locally Lipschitz and C1, and let u be a bounded Kruzhkov entropy solution with initial datum u0∈L∞. Then, with essential extrema taken with respect to Lebesgue measure, ess inf⁡Rnu0 ≤ u(t,x) ≤ ess sup⁡Rnu0for a.e. (t,x)∈ΠT; in particular ∥u(t,⋅)∥∞≤∥u0∥∞ for almost every t. For a representative continuous in local L1, the same bound holds at every t: local L1 convergence from times in the full-measure set preserves the range bound (The essential supremum of a measurable function with respect to a measure, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable Choice, n≥1, T>0, a locally Lipschitz C1 flux f, a bounded Kruzhkov entropy solution u on ΠT with datum u0∈L∞(Rn), and the essential bounds m0=ess inf⁡u0, M0=ess sup⁡u0, both finite.

[F1]

Constant functions on ΠT are Kruzhkov entropy solutions with their own constant value as initial datum, for every flux: for u≡c the weak equation is the equality ∂tc+div⁡xf(c)=0, and for every k the functions ηk(c)=∣c−k∣ and qk(c)=sgn⁡(c−k)(f(c)−f(k)) are constant in (t,x), so ∂tηk(c)+div⁡xqk(c)=0≤0 in distributions; the strong local L1 trace of the constant c is the constant c, with ∫K∣c−c∣ dx=0 for every compact K (Kruzhkov entropy solutions).

[F2]

Order preservation: if two bounded Kruzhkov entropy solutions v,w on ΠT have ∣v∣,∣w∣≤M and v0≤w0 almost everywhere, then v≤w almost everywhere on ΠT (Uniqueness, comparison and order preservation of entropy solutions).

[F3]

The cited essential-supremum definition defines ∥u0∥∞ using bounds on ∣u0∣ (The essential supremum of a measurable function with respect to a measure). Here define the signed extrema explicitly by M0=inf⁡{b∈R:u0≤b a.e.} and m0=sup⁡{a∈R:a≤u0 a.e.}. Since u0 is essentially bounded on the nonnull space Rn, these are finite. For each integer j≥1, the infimum property gives an essential upper bound below M0+1/j, so u0≤M0+1/j a.e.; the supremum property similarly gives m0−1/j≤u0 a.e. Discarding the countable union of exceptional null sets and letting j→∞ yields m0≤u0≤M0 a.e. Thus ∥u0∥∞≤max⁡{∣m0∣,∣M0∣}. Conversely every essential absolute bound B gives m0≥−B and M0≤B, so max⁡{∣m0∣,∣M0∣}≤B; taking its infimum proves equality. Inequalities between L∞ classes are a.e. (The space Lp(μ) as the quotient by null functions). Fubini transfers null sets to spatial slices for a.e. time (Fubini's theorem for L^1 functions on a sigma-finite product).

Proof

technique · direct
1.1F1F2F3

Comparison with the constant ceilings and floors. By [F1] the constants c=M0 and c′=m0 are bounded Kruzhkov entropy solutions. Since u0≤M0 almost everywhere and m0≤u0 almost everywhere by [F3], choose a finite common bound for u, m0, and M0. Then [F2] applied to the pairs (u,M0) and (m0,u) gives u≤M0 almost everywhere and m0≤u almost everywhere on ΠT, that is, m0≤u(t,x)≤M0 for almost every (t,x).

2.1F3step 1.1

The almost-everywhere L∞ bound. Integrating the pointwise almost-everywhere bound of step 1.1 over spatial slices and using Fubini, for almost every t∈(0,T) one has m0≤u(t,x)≤M0 for almost every x, hence ∥u(t,⋅)∥∞≤max⁡{∣m0∣,∣M0∣}=∥u0∥∞ for almost every t.

3.1F3step 2.1∎

Every time for a continuous representative. Suppose u has a representative on [0,T] continuous into Lloc1(Rn): for tj→t and every compact K, u(tj)→u(t) in L1(K). Fix t∈[0,T] and choose tj→t with tj in the full-measure set of step 2.1. For each ball BR, the bound ∥(u(t)−M0)+∥L1(BR)+∥(m0−u(t))+∥L1(BR)≤2∥u(t)−u(tj)∥L1(BR)→0 preserves the range directly; exhausting Rn by countably many balls, the bound holds for almost every x∈Rn at this time t.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The entropy solution semigroup on L1∩L∞

Statement

Assume Countable Choice and Dependent Choice (The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) for the heat-kernel, L1 completeness and vanishing-viscosity extraction interfaces used below. Let n≥1 and let f ⁣:R→Rn be locally Lipschitz and C1. For u0∈L1(Rn)∩L∞(Rn) and t≥0 let Stu0 be the value at time t of the unique Kruzhkov entropy solution with datum u0 (Existence of bounded Kruzhkov entropy solutions, Uniqueness, comparison and order preservation of entropy solutions), with S0u0=u0. Then: (i) St+s=St∘Ss for all s,t≥0 (semigroup law); (ii) each St is order-preserving and an L1 contraction, ∥Stu0−Stv0∥1≤∥u0−v0∥1; (iii) ∥Stu0∥∞≤∥u0∥∞. If in addition f(0)=0 and f is globally Lipschitz on R, then each St extends uniquely to a map on L1(Rn) that is order-preserving, an L1 contraction and satisfies the same semigroup law; the extension agrees with the classical flow on L1∩L∞ (Kruzhkov entropy solutions, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable and Dependent Choice, n≥1, a locally Lipschitz C1 flux f, and data u0,v0∈L1∩L∞ with the associated unique bounded Kruzhkov entropy solutions Stu0, Stv0 on any finite time horizon.

[F1]

Existence and uniqueness: for every datum in L1∩L∞ there is a bounded Kruzhkov entropy solution, unique in the bounded Kruzhkov class, with a representative continuous in Lloc1 on [0,T] and attaining the datum in the strong local L1 sense (Existence of bounded Kruzhkov entropy solutions, Uniqueness, comparison and order preservation of entropy solutions, Kruzhkov entropy solutions).

[F2]

Comparison and contraction: if u0≤v0 almost everywhere then Stu0≤Stv0 almost everywhere, and ∥Stu0−Stv0∥1≤∥u0−v0∥1 for every t (Uniqueness, comparison and order preservation of entropy solutions, Global L1 contraction from the local estimate).

[F3]

L∞ bound: ∥Stu0∥∞≤∥u0∥∞ for every t; in particular the range of each solution is contained in a bounded interval on which f is Lipschitz (The L∞ maximum bound for entropy solutions).

[F4]

Truncation and dominated convergence: for u0∈L1, the truncations u0m=(−m)∨(u0∧m) lie in L1∩L∞ and converge to u0 in L1; limits of sequences of equivalence classes are taken in L1 and are independent of the pointwise representatives (Dominated convergence, Monotone convergence for the integral, The space Lp(μ) as the quotient by null functions). The Cauchy limits exist by Riesz-Fischer completeness of Lp for 1≤p≤∞.

[F5]

If an initial datum is supported in B(0,R), finite propagation gives support of its entropy solution in B(0,R+Lt) for almost every t, where L is a Lipschitz constant of f on the common range (Finite propagation for scalar conservation laws, Open ball, closed ball and sphere in a metric space). The representative is continuous in Lloc1 by [F1].

Proof

technique · direct
1.1F1F2F3

Semigroup law. Fix s,t≥0 and a horizon T>s+t. The solution Sτu0 is in L1 for every τ∈[0,T]: compare it with the zero solution in [F2] to get ∥Sτu0∥1≤∥u0∥1; its L∞ bound follows from [F3]. Thus Ssu0∈L1∩L∞ and [F1] supplies the entropy solution z(τ)=Sτ(Ssu0). Define w(τ,x)=Ss+τu0(x) on ΠT−s: its entropy inequalities are those of the original solution with time shifted, and its strong local L1 trace at τ=0 is Ssu0 by the representative's continuity in Lloc1. Both w and z are bounded entropy solutions with this same datum, so uniqueness [F1] gives w=z almost everywhere on ΠT−s. Their time-continuous representatives then agree at every time in Lloc1, so evaluating at τ=t<T−s gives Ss+tu0=St(Ssu0); as T is arbitrary, this holds for all s,t≥0.

1.2F2F3

Order, contraction and the maximum bound. Let u0≤v0 almost everywhere in L1∩L∞; by [F2] and [F3], Stu0≤Stv0 almost everywhere, ∥Stu0−Stv0∥1≤∥u0−v0∥1, and ∥Stu0∥∞≤∥u0∥∞ for every t≥0. This proves (ii) and (iii).

2.1F2F4step 1.2

Extension to L1: construction. Assume f(0)=0 and f globally Lipschitz, and let u0∈L1. Put u0m=(−m)∨(u0∧m) as in [F4]. For m,ℓ≥1 and every t≥0, step 1.2 gives ∥Stu0m−Stu0ℓ∥1≤∥u0m−u0ℓ∥1, so (Stu0m)m is Cauchy in L1, uniformly in t; define Stu0=lim⁡mStu0m in L1. The definition is independent of the approximating sequence: if wm∈L1∩L∞ with wm→u0 in L1, then ∥Stwm−Stu0m∥1≤∥wm−u0m∥1→0, so both sequences have the same limit.

3.1F1F2F4F5step 2.1

Time continuity of the L1 flow. First fix w0∈L1∩L∞ and T>0, and choose R>0 with w0R:=w01B(0,R). Let M=∥w0∥∞ and let L be a Lipschitz constant of f on [−M,M]. By [F5], for almost every t∈(0,T) the orbit Stw0R is supported in B(0,R+Lt). Fix t0∈[0,T] and a compact set K⊆{x:∣x∣>R+Lt0}. Its positive distance from B‾(0,R+Lt0) lets us choose times tj from that full-measure set tending to t0 with K⊆{x:∣x∣>R+Ltj}. Then Stjw0R=0 in L1(K), and the Lloc1 continuity [F1] gives St0w0R=0 in L1(K). A countable exhaustion of the strict exterior by compact sets shows that every slice is supported in B‾(0,R+Lt0); hence all slices on [0,T] are supported in the fixed ball B‾(0,R+LT). Local L1 continuity is therefore global L1 continuity for this truncated orbit. By [F2], sup⁡t∈[0,T]∥Stw0−Stw0R∥1≤∥w0−w0R∥1, which tends to 0 as R→∞. Thus the L1-continuous truncated orbits converge uniformly on [0,T] to t↦Stw0, proving continuity for every datum in L1∩L∞. For u0∈L1 in step 2.1, the extension orbit is the uniform limit of the continuous orbits t↦Stu0m, since sup⁡t≥0∥Stu0−Stu0m∥1≤∥u0−u0m∥1→0; hence the extension is continuous as well.

3.2F2F4step 1.1step 2.1

Extension: properties. The extended maps preserve order: if u0≤v0 in L1, then the truncated sequences satisfy u0m≤v0m and hence Stu0m≤Stv0m almost everywhere; passing to the L1 limit gives Stu0≤Stv0 almost everywhere. They are contractions: ∥Stu0−Stv0∥1≤lim inf⁡m∥Stu0m−Stv0m∥1≤lim inf⁡m∥u0m−v0m∥1=∥u0−v0∥1, using that truncation is a contraction in L1. The semigroup law passes to the limit: St+su0=lim⁡mSt+su0m=lim⁡mSt(Ssu0m)=St(Ssu0), the last step by the contraction property just proved applied to Ssu0m→Ssu0. Finally, the extension agrees with the original flow on L1∩L∞, because for such u0 the estimate of step 2.1 with wm=u0 gives Stu0=lim⁡mStu0m in the original sense as well. An order-preserving L1 contraction agreeing on the dense subset L1∩L∞ is unique, so the extension is unique.

4.1step 1.1step 1.2step 2.1step 3.1step 3.2∎

Conclusion. Steps 1.1–1.2 prove (i)–(iii) for data in L1∩L∞, and steps 2.1 and 3.2 construct and characterise the unique order-preserving L1 contraction extension to L1 when f(0)=0 and f is globally Lipschitz, agreeing with the classical flow on L1∩L∞ and satisfying the semigroup law. Step 3.1 proves strong L1 continuity of these orbits. This completes the proof.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Entropy solution orbits are strongly continuous in L1

Statement

Assume Countable Choice and Dependent Choice (The Axiom of Countable Choice (ACω), The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) for the heat-kernel, L1 completeness and vanishing-viscosity extraction interfaces used below. Let n≥1, let f ⁣:R→Rn be C1, and let u0∈L1(Rn)∩L∞(Rn). Then the orbit t↦Stu0 is continuous from [0,∞) to L1(Rn): for every t≥0, ∥Stu0−Ssu0∥1→0 as s→t, and in particular lim⁡t↓0∥Stu0−u0∥L1(Rn)=0. If additionally f(0)=0 and f is globally Lipschitz, its extension is a strongly continuous semigroup of L1 contractions on all of L1 (The entropy solution semigroup on L1∩L∞, Kruzhkov entropy solutions, The space Lp(μ) as the quotient by null functions).

Facts & Assumptions

Given: Countable and Dependent Choice, n≥1, a C1 flux, with the extra global Lipschitz and f(0)=0 hypotheses only for the extension, a datum u0∈L1∩L∞ with M=∥u0∥∞, and the entropy solution semigroup St of The entropy solution semigroup on L1∩L∞.

[F1]

Semigroup and contraction: St+s=StSs, each St is order-preserving, and ∥Stu0−Stv0∥1≤∥u0−v0∥1 for all t≥0 and u0,v0∈L1∩L∞; when f(0)=0 and f is globally Lipschitz the maps extend to order-preserving L1 contractions on all of L1, agreeing with the flow on L1∩L∞ (The entropy solution semigroup on L1∩L∞).

[F2]

Every any-time contraction, for data in L1∩L∞, holds for every t because the solutions have Lloc1-continuous representatives: ∥Stu0−Stv0∥1≤∥u0−v0∥1 (Global L1 contraction from the local estimate, Existence of bounded Kruzhkov entropy solutions).

[F3]

Finite propagation: if u0 vanishes almost everywhere outside B(0,R), the entropy solution vanishes almost everywhere outside B(0,R+Lt) for almost every t, where L is a Lipschitz constant of the flux on the common range; with the Lloc1-continuous representative this support statement upgrades to every t (Finite propagation for scalar conservation laws, Open ball, closed ball and sphere in a metric space).

[F4]

The strong local L1 trace: for every compact K, ess sup⁡0<t<δ∫K∣u(t,x)−u0(x)∣ dx→0 as δ↓0; and monotone convergence controls the tails of an L1 function over increasing balls (Kruzhkov entropy solutions, Monotone convergence for the integral, The space Lp(μ) as the quotient by null functions).

Proof

technique · direct
1.1F2F3F4

Continuity at time zero. Fix R>0, put u0R=u01B(0,R) and vR(t)=Stu0R; the datum u0R is compactly supported, so by [F3] the solution vR is supported in B(0,R+Lt) for every t, where L is a Lipschitz constant of f on the common range [−∥u0∥∞,∥u0∥∞]. By [F2], ∥Stu0−vR(t)∥1≤∥u0−u0R∥1=∥u0∥L1(B(0,R)c) for every t. Hence for 0<t<1/L (any t>0 when L=0) the exterior B(0,R+1)c is contained in B(0,R+Lt)c and ∥Stu0−u0∥1≤∥Stu0−u0∥L1(B(0,R+1))+2∥u0∥L1(B(0,R)c), because on B(0,R+1)c both Stu0 and u0 have L1 norms bounded by ∥u0∥L1(B(0,R)c) (for Stu0 combine the contraction bound with vR(t)=0 there). The first term tends to 0 as t↓0 by the local L1 continuity of the chosen representative [F2] and its trace [F4], and the tail term tends to 0 as R→∞ by monotone convergence. Therefore lim⁡t↓0∥Stu0−u0∥1=0.

2.1F1F2step 1.1

Continuity at every time. Let s,t≥0 with t>s. By the semigroup law and the contraction estimate of [F1], ∥Stu0−Ssu0∥1=∥Ss(St−su0)−Ssu0∥1≤∥St−su0−u0∥1, and the right side tends to 0 as t↓s by step 1.1 applied to the fixed datum u0. The case s<t is symmetric, so the orbit is continuous at every t≥0.

3.1F1step 2.1∎

Strong continuity on the closure. If f(0)=0 and f is globally Lipschitz, the extension of [F1] is defined on all of L1, and L1∩L∞ is dense in L1 (the closure appearing in the statement). For u0∈L1 and u0m=(−m)∨(u0∧m)∈L1∩L∞ with u0m→u0 in L1, the contraction property gives ∥Stu0−Stu0m∥1≤∥u0−u0m∥1 for every t, so ∥Stu0−Ssu0∥1≤2∥u0−u0m∥1+∥Stu0m−Ssu0m∥1→0 by first making the two approximation errors small with a fixed large m and then taking s→t, by step 2.1 applied to each u0m. Hence the extended semigroup is strongly continuous on all of L1, which is the closure of L1∩L∞.

5 · Examples, counterexamples and false statements

None yet.

Sources