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.

✓ 10 results · all verified · 7 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 3 not AI-judged were verified by owner audit (typically over a confirmed judge false positive), not failures.

Hamilton Jacobi Equations and Viscosity Solutions — Examples

1 · Prerequisites

2 · Summary

These companions compute the viscosity theory of the main page on explicit equations and mark its scope boundaries. The eikonal equation is tested at its ridge: the positive norm has no upper contact at the tip but fails the supersolution test, while the negative norm has upper contacts with slopes in [−1,1] and no lower contact; distance to the boundary of the unit ball is verified as the nonsmooth solution of the unit eikonal Dirichlet problem, and the counterexample of three distinct eikonal solutions on an interval shows that omitting endpoint data destroys uniqueness. The quadratic Hopf--Lax formula is identified with an infimal convolution and the quadratic Moreau envelope is computed separately for the unbounded datum. Two scope counterexamples show that nonconvexity breaks the Hopf--Lax equation and that nonsuperlinearity restricts the Lagrangian to a bounded velocity set, and that minima of subsolutions of a proper equation need not be subsolutions. A smooth datum whose characteristics cross at t=1/2 produces a Hopf--Lax solution with a forming corner, and the Cole--Hopf family for sin⁡ is shown to converge to the Hopf--Lax solution uniformly on finite time strips. None of the examples uses a choice principle.

3 · Logical flowchart

4 · Definitions, theorems and proofs

None yet.

5 · Examples, counterexamples and false statements

ExampleConstruction: Literature-sourcedVerification: AI-adaptedprecheck passOpen item page →

The eikonal equation as a viscosity equation at a tip

Example

Let n≥1 and consider the eikonal equation ∣Du∣=1 on Rn, i.e. the first-order equation ut+H(Du)=0 with H(p)=∣p∣−1 read on functions independent of t. The Euclidean norm u(x)=∣x∣ is a classical solution on Rn∖{0}, and at its tip a=0 it exhibits the difference between the subsolution inequality ∣Dϕ∣≤1 and the reverse supersolution inequality ∣Dϕ∣≥1. There is no C1 test function ϕ with u−ϕ having a local maximum at 0: such a contact would force ϕ(x)≥ϕ(0)+∣x∣ near 0, hence ⟨Dϕ(0),h⟩≥∣h∣ for every direction h, which is impossible for a linear functional. Writing x0 for the coordinate indexed by 0<n, the C1 functions ϕ(x)=εx0 with 0<ε<1 satisfy ϕ≤∣x∣ near 0 with equality at 0, so u−ϕ has a local minimum at 0 and the supersolution test would require ∣Dϕ(0)∣=ε≥1, which fails. Hence u is a subsolution of ∣Du∣=1 but not a supersolution at the tip: it satisfies the viscosity subsolution inequality ∣Du∣≤1 at the tip and is a viscosity solution of ∣Du∣=1 on Rn∖{0}.

Verification

Given: The tip a=0∈Rn, the function u(x)=∣x∣, the test-function definition of viscosity sub- and supersolutions for ut+H(Du)=0 with H(p)=∣p∣−1 (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem), and the Euclidean inner product ⟨⋅,⋅⟩ (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn) with gradient as in The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case.

[F1] u−ϕ has a local maximum (respectively minimum) at 0 for ϕ∈C1 precisely when u(x)−ϕ(x)≤u(0)−ϕ(0) (respectively ≥) for all x near 0, and the viscosity inequalities are tested at such contacts (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem). For the stationary extension U(x,t)=u(x) on Rn×(0,1), a space--time contact with Φ has Φt=0: restrict to the time line, where the differentiable function t↦Φ(x,t) has a local extremum. Restricting to the spatial slice therefore gives the inequalities ∣Dϕ∣≤1 and ∣Dϕ∣≥1 used here; conversely, a spatial test extends to a time-independent space--time test. The norm is continuous by [F2], so the required semicontinuity holds.

[F2] The Euclidean norm satisfies ∣th∣=∣t∣ ∣h∣ for t∈R, and for x≠0 it is differentiable at x with gradient x/∣x∣, of norm 1: The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn gives ∣x+th∣2=∣x∣2+2t⟨x,h⟩+t2∣h∣2, the norm axioms and the triangle inequality in clause 2 of Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation give ∣x+th∣≥∣x∣−∣t∣ ∣h∣ and ∣∣x+th∣−∣x∣∣≤∣t∣ ∣h∣, Cauchy--Schwarz (clause 1 there) gives ∣⟨x,h⟩∣≤∣x∣ ∣h∣, and since ∣x+th∣+∣x∣≥∣x∣>0 the identity ∣x+th∣−∣x∣=(2t⟨x,h⟩+t2∣h∣2)/(∣x+th∣+∣x∣) differs from t⟨x,h⟩/∣x∣ by at most 2t2∣h∣2/∣x∣, so the gradient is x/∣x∣ (The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case).

Proof technique: direct test-function computation at the tip and the classical-consistency proposition away from it.

1.1F1F2algebra

No upper test exists at the tip. Suppose ϕ∈C1(Rn) with u−ϕ having a local maximum at 0; by [F1], ∣x∣−ϕ(x)≤−ϕ(0) near 0, that is ϕ(x)≥ϕ(0)+∣x∣ there. Substituting x=th with ∥h∥=1 and t↓0 and dividing by t gives ⟨Dϕ(0),h⟩≥1 for every unit vector h; testing h and −h gives both ⟨Dϕ(0),h⟩≥1 and −⟨Dϕ(0),h⟩≥1, an impossibility. Hence there is no upper test at the tip and the subsolution inequality holds vacuously there.

1.2F1F2algebra

A lower test with a failing supersolution inequality. Every lower contact at the tip has slope of norm at most one: if ϕ is C1 with u−ϕ having a local minimum at 0, then after translating ϕ we have ∣x∣≥⟨Dϕ(0),x⟩+o(∣x∣); substituting x=th with ∥h∥=1 gives ⟨Dϕ(0),h⟩≤1 for t>0 and ⟨Dϕ(0),h⟩≥−1 for t<0, hence ∣Dϕ(0)∣≤1. Now take ϕ(x)=εx0, where x0 is the coordinate indexed by 0<n, with 0<ε<1: then ϕ(x)≤∣x∣ near 0 with equality at 0, so u−ϕ has a local minimum at 0 by [F1], while the supersolution condition requires ∣Dϕ(0)∣=ε≥1 and fails. This single lower contact shows that u is not a viscosity supersolution of ∣Du∣=1 at the tip.

1.3F2algebra

Away from the tip the equation holds classically. On the open set Rn∖{0} the function u is C1 with ∣Du∣=1 by [F2]. Its stationary extension on (Rn∖{0})×(0,1) extends continuously to the closed cylinder with initial datum ∣x∣, so Classical solutions are viscosity solutions and differentiable viscosity solutions solve the equation pointwise applies and it is a viscosity solution of ∣Du∣=1 there; in particular it is both a subsolution and a supersolution at every x≠0.

2.1step 1.1step 1.2step 1.3∎

Conclusion. By step 1.1 the subsolution test at the tip is vacuous, by step 1.2 the supersolution test fails there, and by step 1.3 both tests hold away from the tip. Hence u solves the subsolution inequality ∣Du∣≤1 everywhere and solves ∣Du∣=1 exactly on Rn∖{0}, while it is not a viscosity solution of ∣Du∣=1 on any neighbourhood of the tip.

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

The quadratic Hopf--Lax formula as an infimal convolution

Example

Let n≥1 and H(p)=12∣p∣2. Its Legendre transform is L(v)=12∣v∣2. For every bounded uniformly continuous u0:Rn→R, the Hopf--Lax operator of The Hopf--Lax operator and the Hopf--Lax formula is, for t>0, Qtu0(x)=inf⁡y∈Rn{u0(y)+∣x−y∣22t}, the infimal convolution with the quadratic kernel qt(z)=∣z∣2/(2t). The infimum is attained, and for every minimiser y the Euler relation Du0(y)=(x−y)/t holds whenever u0 is differentiable at y (A Hopf--Lax minimiser satisfies the characteristic Euler relation at differentiability points). Separately, the quadratic datum w0(y)=12∣y∣2 is unbounded and so is outside the datum class in The Hopf--Lax operator and the Hopf--Lax formula. Its algebraic infimal convolution It(x):=inf⁡y∈Rn{12∣y∣2+∣x−y∣22t} has the unique minimiser y=x/(1+t) and value It(x)=∣x∣2/(2(1+t)); this separate calculation is the Moreau envelope of the quadratic function and does not apply the bounded-data Hopf--Lax theorem to w0.

Verification

Given: The Hamiltonian H(p)=12∣p∣2 on Rn, its Legendre transform L, a bounded uniformly continuous datum u0, the operators Qt of The Hopf--Lax operator and the Hopf--Lax formula, and the unbounded quadratic datum w0(y)=12∣y∣2.

[F1] L(v)=sup⁡p∈Rn(p⋅v−H(p)), the supremum taken in R‾ (The Legendre transform of a finite-valued convex Hamiltonian).

[F2] For t>0, Qtu0(x)=inf⁡y{u0(y)+tL((x−y)/t)}, and under convexity and superlinearity of H the infimum is finite and attained for bounded uniformly continuous u0 (The Hopf--Lax operator and the Hopf--Lax formula, Finiteness, superlinearity of the Lagrangian and localisation of Hopf--Lax near-minimisers).

[F3] At a minimiser y of u0(y)+tL((x−y)/t), if u0 is differentiable at y and L at (x−y)/t, then Du0(y)=DL((x−y)/t) (A Hopf--Lax minimiser satisfies the characteristic Euler relation at differentiability points).

Proof technique: compute the quadratic conjugate, use the in-class Hopf--Lax suppliers only for bounded uniformly continuous data, and evaluate the separate quadratic infimum by completing the square.

1.1F1algebra

The conjugate of the quadratic Hamiltonian. Fix v∈Rn and complete the square: p⋅v−12∣p∣2=12∣v∣2−12∣p−v∣2, whose supremum over p is attained at p=v with value 12∣v∣2. Hence L(v)=12∣v∣2 by [F1]; in particular L is finite, convex and superlinear.

2.1step 1.1F2F3algebra

The formula, attainment and the Euler relation. Substituting L(v)=∣v∣2/2 into the definition of Qt gives the displayed infimal convolution. The infimum is attained by [F2], and for every minimiser y the conditional Euler relation Du0(y)=DL((x−y)/t)=(x−y)/t is [F3]; both suppliers use only the bounded uniformly continuous data class.

3.1step 1.1algebra∎

The separate quadratic infimum. For t>0 and all x,y, completing the square gives 12∣y∣2+∣x−y∣22t=1+t2t∣y−x1+t∣2+∣x∣22(1+t). Since the coefficient (1+t)/(2t) is positive, the infimum over y is attained uniquely at y=x/(1+t) with value ∣x∣2/(2(1+t)). This is a direct computation for the unbounded datum w0 and makes no assertion that w0 lies in the domain of the Hopf--Lax operator.

Remarks

  • What is and is not applied. The bounded-data statements are applied only to bounded uniformly continuous u0; the quadratic datum is treated by the displayed algebraic computation, which is the Moreau envelope of 12∣⋅∣2 and does not claim a Hopf--Lax solution for it.
ExampleConstruction: Literature-sourcedVerification: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

A Hopf--Lax solution with a forming corner from smooth data

Example

Let H(p)=p2/2, L(v)=v2/2, and u0(x)=e−x2−1. Then u0 is smooth, bounded and uniformly continuous, with p0(y):=u0′(y)=−2ye−y2 and u0′′(y)=(4y2−2)e−y2. The characteristic projection is Xt(y)=y+tp0(y)=y−2tye−y2, and its lifted value is Zt(y)=u0(y)+t2p0(y)2. For 0≤t<12, Xt is a diffeomorphism of R, and ucl(Xt(y),t)=Zt(y) is the classical characteristic solution. At t=12, Xt′(0)=0; for t>12, putting st=log⁡(2t) gives Xt(−st)=Xt(0)=Xt(st)=0, so the characteristic projection is no longer injective and its single-valued classical graph breaks down. The Hopf--Lax formula u(x,t)=inf⁡y∈R{e−y2−1+∣x−y∣22t},t>0, is finite, satisfies −1≤u(x,t)≤0, is uniformly continuous in x, and is a viscosity solution of ut+12ux2=0 with initial datum u0 (The Hopf--Lax formula solves the Hamilton--Jacobi Cauchy problem, The Hopf--Lax operator preserves a modulus of continuity). For each fixed t>12, the minimisers at x=0 are exactly y=±st, and u(0,t)=1+log⁡(2t)2t−1. For x>0 sufficiently close to 0, the unique minimiser tends to st as x↓0; for x<0 sufficiently close to 0, it tends to −st as x↑0. Thus the one-sided spatial derivatives tend to −st/t from the right and st/t from the left, so u(⋅,t) is continuous but has a corner at x=0.

Verification

Given: The Hamiltonian H(p)=p2/2 with Lagrangian L(v)=v2/2, the datum u0(y)=e−y2−1, its derivatives p0=u0′, u0′′, the characteristic data Xt(y)=y+tp0(y), Zt(y)=u0(y)+t2p0(y)2, and the Hopf--Lax function u=Qtu0 (The Hopf--Lax operator and the Hopf--Lax formula, Characteristic crossing and caustic for a first-order PDE).

[F1] Qtu0 satisfies the bounds −1≤Qtu0≤0, is spatially uniformly continuous, and is a viscosity solution with datum u0 (The Hopf--Lax operator preserves a modulus of continuity, The Hopf--Lax formula solves the Hamilton--Jacobi Cauchy problem); the defining infimum is attained (Finiteness, superlinearity of the Lagrangian and localisation of Hopf--Lax near-minimisers).

[F2] Differentiation rules for the elementary functions give p0(y)=−2ye−y2, p0′(y)=(4y2−2)e−y2, and Xt′(y)=1+tp0′(y)=1−2te−y2+4ty2e−y2 (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).

Proof technique: explicit characteristic and minimiser computations.

1.1F2F3algebra

Classical solution before the first singular time. For 0≤t<1/2, [F2] gives Xt′(y)≥1−2t>0. The mean value theorem [F3] makes Xt strictly increasing, and Xt(y)−y=−2tye−y2→0 at both infinities, so the intermediate value theorem gives a unique inverse Y(x,t). This inverse is continuous jointly: near any fixed (x,t), Xt′ has a positive lower bound, and the mean value theorem bounds changes in Y by changes in x and in Xt(Y). Differentiating the identity x=Y+tp0(Y) by difference quotients then gives Yx=(1+tp0′(Y))−1 and Yt=−p0(Y)/(1+tp0′(Y)), continuously. Thus ucl(x,t)=u0(Y)+(t/2)p0(Y)2 is C1, with ucl,x=p0(Y) and ucl,t=−p0(Y)2/2, by substitution of these derivatives. It solves the equation and has the initial datum u0.

1.2F2algebra

Breakdown of the projection. At t=12 we have Xt′(0)=1−2t=0. For t>12 put st=log⁡(2t)>0; then e−st2=1/(2t), so Xt(±st)=±st−2t(±st)/(2t)=0=Xt(0), and the projection is not injective; it is locally decreasing near y=0 because Xt′(0)=1−2t<0. The loss of rank at t=1/2, y=0, is the caustic of Characteristic crossing and caustic for a first-order PDE; the later equal projections show global folding, without asserting local noninjectivity at y=0 for t>1/2.

2.1step 1.1step 1.2F1F2F3algebra∎

Bounds, minimisers at x=0, and the corner. By [F1] the Hopf--Lax function is finite, −1≤u≤0 and uniformly continuous in x. For gt(y):=e−y2−1+y2/(2t) we have gt′(y)=y(1/t−2e−y2), so for t>12 the critical points are 0,±st, with gt′′(0)=1/t−2<0 and gt′′(±st)=2st2/t>0; since gt(y)→∞ as ∣y∣→∞ and gt(±st)=−1+(1+log⁡(2t))/(2t)<0=gt(0) (since gt′(y)<0 for 0<y<st and [F3] makes gt strictly decreasing there), the global minimisers of gt are exactly ±st, and the value is (1+log⁡(2t))/(2t)−1. For Gt(x,y):=u0(y)+(x−y)2/(2t) any minimiser obeys ∣x−y∣≤2t, so all minimisers for ∣x∣≤1 lie in a fixed compact interval. For every neighbourhood of {−st,st}, the complement in this interval has a positive gap above min⁡gt by continuity and compact attainment [F1]; uniform convergence Gt(x,⋅)→gt on the interval forces every minimiser into that neighbourhood for all sufficiently small ∣x∣. since Gt(x,y)−Gt(x,−y)=−2xy/t, a minimiser cannot be negative when x>0 nor positive when x<0, and y=0 is not a minimiser for x≠0 because ∂yGt(x,0)=−x/t≠0. Hence all minimisers have the sign of x and, as x↓0, they converge to st (and to −st as x↑0). The stationarity equation is x=Ft(y):=y−2tye−y2 with Ft′(±st)=2st2>0, so on small intervals around ±st, Ft′ is bounded below by a positive constant. The mean value and intermediate value theorems [F3] give a unique local inverse there, and its difference quotient has derivative 1/Ft′(y), which is continuous. Since every minimiser is on the corresponding interval for small ∣x∣, this inverse is the unique C1 minimising branch on each punctured side. Along a branch the envelope derivative is ux=(x−y(x))/t, whose one-sided limits are −st/t (from the right) and st/t (from the left); these unequal finite limits show that u(⋅,t) has a corner at x=0 while remaining continuous.

Remarks

  • What is claimed. The computation identifies the minimisers and the one-sided derivatives at the corner; it does not assert local noninjectivity of Xt near y=0, where Xt′(0)<0 and the map is locally decreasing, and it does not claim that the classical solution extends past t=12.
ExampleConstruction: Literature-sourcedVerification: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

The negative absolute value solves the eikonal equation in the viscosity sense

Example

Let u(x)=−∣x∣ on R. Then ∣u′∣=1 for x≠0, and u is a viscosity solution of the stationary eikonal equation ∣u′∣−1=0. Equivalently, for any T>0 its evolutionary extension U(x,t):=−∣x∣−t is a viscosity solution of Ut+∣Ux∣=0 on R×(0,T) (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem). At every point x≠0 the equation holds classically. At the cusp (0,t0), every C1 upper test has ϕt(0,t0)=−1 and satisfies ∣ϕx(0,t0)∣≤1, so the subsolution inequality holds, and there is no C1 test function for which U−ϕ has a local minimum at (0,t0); hence the supersolution test is vacuous. This is the complementary cusp to The eikonal equation as a viscosity equation at a tip, where the positive absolute value fails the supersolution test because it has lower tests with slopes of modulus less than one.

Verification

Given: The function U(x,t)=−∣x∣−t on U=R×(0,T), the equation Ut+∣Ux∣=0, and the test-function definition of viscosity sub- and supersolutions.

[F1] The subsolution inequality is tested at local maxima of U−ϕ and the supersolution inequality at local minima, for C1 test functions ϕ (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem).

[F2] At a local extremum of a differentiable function of two variables both partial derivatives vanish (Fermat's theorem: an interior differentiable local extremum has zero gradient); away from x=0 the function U is C1 with Ut=−1 and Ux=∓1, so it solves the equation classically there and is a viscosity solution on the open set R∖{0}×(0,T) by Classical solutions are viscosity solutions and differentiable viscosity solutions solve the equation pointwise.

Proof technique: direct contact computation at the cusp.

1.1F2

Away from the cusp. On the open set {x≠0} the function U is C1 with Ut=−1 and ∣Ux∣=1, so it is a viscosity solution of the equation there by [F2].

1.2F1F2algebra

Upper contacts at the cusp. Fix t0 and let ϕ∈C1 with U−ϕ having a local maximum at (0,t0); normalize ϕ(0,t0)=U(0,t0)=−t0. Restricting to the line x=0 and using [F2] gives ϕt(0,t0)=−1. Writing p:=ϕx(0,t0) and testing x=h, x=−h with h↓0 in the inequality −∣x∣≤ϕ(x,t0)−ϕ(0,t0)=px+o(∣x∣) gives p≥−1 and p≤1, that is ∣p∣≤1. Hence ϕt+∣ϕx∣=−1+∣p∣≤0, which is the subsolution inequality.

1.3F1F2algebra

No lower contact at the cusp. If U−ϕ had a local minimum at (0,t0), the same computation with the inequality reversed would give ϕt(0,t0)=−1 and p≤−1 (from h>0) together with p≥1 (from h<0), an impossibility; hence the set of lower contacts at the cusp is empty and the supersolution inequality holds vacuously.

2.1step 1.1step 1.2step 1.3∎

Conclusion. Steps 1.2 and 1.3 show that U is a subsolution everywhere on R×(0,T) and a supersolution everywhere, hence a viscosity solution; step 1.1 identifies the classical region. The cusp supports the subsolution inequality but admits no lower test, which is the complementary behaviour to the positive absolute value at its ridge point.

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

Vanishing viscosity selects the Hopf--Lax solution for bounded data

Example

Let u0(x)=sin⁡x on R, a bounded uniformly continuous datum, and let H(p)=p2/2. For ε>0 and t>0 define uε(x,t)=−2εlog⁡ ⁣[(4πεt)−1/2∫Rexp⁡ ⁣(−sin⁡y+(x−y)2/(2t)2ε) dy], and set uε(x,0)=sin⁡x. Then uε is a bounded classical solution for t>0 of utε+12∣uxε∣2=εuxxε, and it attains u0 uniformly as t↓0. The Hopf--Lax function Qtu0(x):=inf⁡y∈R{sin⁡y+(x−y)22t},Q0u0=u0, is a bounded uniformly continuous viscosity solution of ut+12∣ux∣2=0 with initial trace u0, and for every T<∞ one has sup⁡x∈R, 0≤t≤T∣uε(x,t)−Qtu0(x)∣≤max⁡{εlog⁡(1+T), 2εlog⁡(8πε+2)}. Both terms on the right tend to zero as ε↓0. Thus the viscous solutions converge uniformly on every finite time strip. The estimate permits an O(ε∣log⁡ε∣) error and does not assert a uniform O(ε) rate.

Verification

Given: The datum u0(x)=sin⁡x, the Hamiltonian H(p)=p2/2 with Legendre transform L(v)=v2/2, all displayed integrals interpreted as absolutely convergent improper integrals of continuous functions, the heat kernel Γ(z,s)=(4πs)−1/2e−z2/(4s) (The heat kernel on Rn and its causal extension), the viscous equations and the Hopf--Lax function Qtu0 (The Hopf--Lax operator and the Hopf--Lax formula).

[F1] sin⁡′=cos⁡, ∣cos⁡∣≤1, ∣sin⁡x−sin⁡y∣≤∣x−y∣, and sin⁡ is bounded (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine are 1-Lipschitz on R, Parity and the Pythagorean identity for sine and cosine).

[F2] The Gaussian integral and change of variable for improper integrals are The Gaussian integral ∫−∞∞e−x2 dx=π and Change of variable in an improper integral. The logarithm derivative is 1/x for x>0 (The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t), and the kernel is the explicit function of The heat kernel on Rn and its causal extension; the exponential derivative, chain and algebra rules are The exponential function is smooth and (exp⁡)′=exp⁡, 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. On a compact (x,t)-neighbourhood with t>0, every required kernel derivative is bounded by a fixed polynomial-times-Gaussian function of y, integrable by comparison with a Gaussian. On each compact y-interval its difference quotients converge uniformly; the mean value theorem bounds their tails by that same integrable majorant. Splitting into this interval and its tail justifies differentiation under the improper integral without a Lebesgue change-of-variable theorem or a choice assumption (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)).

[F3] The quadratic conjugate is L(v)=v2/2 by completing the square in The Legendre transform of a finite-valued convex Hamiltonian. For the bounded uniformly continuous datum sin⁡, the infimum is attained (Finiteness, superlinearity of the Lagrangian and localisation of Hopf--Lax near-minimisers) and The Hopf--Lax formula solves the Hamilton--Jacobi Cauchy problem supplies the viscosity solution, finite-strip bounded uniform continuity, initial trace and uniqueness. At a differentiable minimizer the derivative is zero by Fermat's theorem: an interior differentiable local extremum has zero gradient. The semigroup is The Hopf--Lax operators form a semigroup (dynamic programming).

Proof technique: Cole--Hopf representation, an approximate-identity bound, and the variational reduction of the limit.

1.1F1F2algebra

The Cole--Hopf family solves the viscous equation. Write Wε(x,t):=∫RΓ(x−y,εt)e−sin⁡y/(2ε) dy. The substitution z=2εtr and ∫e−r2dr=π show ∫Γ(z,εt)dz=1, and the same substitution with the Gaussian tail shows that Γ(⋅,εt) is an approximate identity as t↓0. Differentiating the kernel gives ∂tΓ(⋅,εt)=ε∂xxΓ(⋅,εt), so by [F2] Wtε=εWxxε, Wε>0, and uε=−2εlog⁡Wε satisfies utε=εuxxε−∣uxε∣2/2, that is utε+∣uxε∣2/2=εuxxε, on t>0. The approximate identity applied to the bounded uniformly continuous function e−sin⁡(⋅)/(2ε) gives Wε→e−sin⁡x/(2ε) uniformly as t↓0, hence uε(x,t)→sin⁡x uniformly; and −1≤uε≤1 because e−1/(2ε)≤e−sin⁡y/(2ε)≤e1/(2ε) and the kernel has unit mass.

1.2F1F2F3algebra

Uniform comparison with Hopf--Lax. Fix x, t>0, put g(y)=sin⁡y+(x−y)2/(2t) and m=min⁡g=Qtu0(x), and take one minimiser y∗. Then −1≤m≤1, g′(y∗)=0, and g′′(y)≤A:=1+1/t. Applying the mean value theorem [F2] to g′ shows that the derivative of g(y∗+r)−m−Ar2/2 is nonpositive for r>0 and nonnegative for r<0; a second application gives g(y∗+r)≤m+Ar2/2. The full Gaussian integral therefore gives Wε(x,t)≥e−m/(2ε)/1+t, hence uε−m≤εlog⁡(1+t). Conversely, g−m≥0 everywhere, and g−m≥(x−y)2/(2t)−2. If ∣y−x∣≥8t, the latter is at least (x−y)2/(4t). Splitting the integral at this radius and using [F2] bounds its normalized ratio by em/(2ε)Wε(x,t)≤(4πεt)−1/2(28t+∫Re−(x−y)2/(8εt) dy)=8πε+2. Thus uε−m≥−2εlog⁡(8/(πε)+2). These two bounds hold for every x and t>0; at t=0 the functions agree. Taking their maximum for 0≤t≤T gives the stated absolute error bound. Its right-hand side tends to zero: writing s=1/ε≥1, the integral formula for the logarithm in [F2] gives log⁡s=∫1sdr/r≤∫1sdr/r=2(s−1), hence (log⁡s)/s→0.

2.1step 1.1step 1.2F1F3∎

The limit is the viscosity solution. By [F3] the function Qtu0 is a viscosity solution of ut+∣ux∣2/2=0; independently, since u0 is 1-Lipschitz with ∣u0∣≤1, the bounds u0(x)−t/2≤Qtu0(x)≤u0(x) hold (competitor y=x and the Lipschitz bound for sin⁡), so Qtu0 has the initial trace u0 and is bounded and uniformly continuous. The uniform estimate of step 1.2, whose right-hand side tends to zero, then gives locally uniform convergence of the viscous family to this viscosity solution without invoking a general vanishing-viscosity theorem and without a momentum-Lipschitz hypothesis.

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

Nonconvexity can break the equation; nonsuperlinearity can limit the Lagrangian domain

Statement refuted

False claim: the Hopf--Lax formula of The Hopf--Lax operator and the Hopf--Lax formula solves the equation for every finite superlinear Hamiltonian, and the Legendre transform of a finite Hamiltonian is real-valued on all of Rn.

Two separate scope boundaries occur. (a) Convexity is essential if the Hopf--Lax formula is claimed to solve the equation for the original Hamiltonian: for H(p)=p2−2∣p∣ on R, which is finite and superlinear but not convex, the formula applied to u0≡0 produces u(x,t)=t, and at every interior point the smooth function itself is an upper test whose residual ut+H(ux)=1+H(0)=1>0 violates the viscosity subsolution inequality (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem). (b) Superlinearity is needed for the full-domain real-valued Lagrangian conclusion: for the finite convex Hamiltonian H(p)=1+p2, which is not superlinear, the Legendre transform equals −1−v2 for ∣v∣≤1 and +∞ for ∣v∣>1, so finite action is available only when ∣x−y∣≤t; this is exactly the full-domain finiteness conclusion of Finiteness, superlinearity of the Lagrangian and localisation of Hopf--Lax near-minimisers that fails without superlinearity. The second clause does not show failure of the value or the solution property: for u0≡0 the formula is Qtu0≡−t, whose residual is −1+H(0)=0.

Facts & Assumptions

Given: The Hamiltonians H(p)=p2−2∣p∣ and H(p)=1+p2 on R, their Legendre transforms L, the bounded uniformly continuous datum u0≡0, and the extended infimum formula Qt defined in [F2].

[F1]

L(v)=sup⁡p∈R(pv−H(p)) for a finite Hamiltonian on R, the supremum being taken in R‾, and L(v)=+∞ is possible (The Legendre transform of a finite-valued convex Hamiltonian).

[F2]

For the two Hamiltonians considered here, define the extended infimum formula Qtu0(x):=inf⁡y∈R{u0(y)+tL((x−y)/t)} for t>0 and Q0u0:=u0. This extends the same expression in The Hopf--Lax operator and the Hopf--Lax formula beyond that definition's convex-superlinear hypotheses. Here L(v)≥−H(0) and u0=0, so the infimum is well defined in R∪{+∞}, with t(+∞)=+∞. The finiteness and attainment theorem Finiteness, superlinearity of the Lagrangian and localisation of Hopf--Lax near-minimisers applies under convexity and superlinearity only; it is not invoked for these two Hamiltonians.

[F3]

A function u is a viscosity subsolution of ut+H(ux)=0 only if ϕt+H(ϕx)≤0 at every local maximum of u−ϕ with ϕ∈C1 (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem).

Counterexample

technique · split each conjugate supremum at the sign of $p$ and evaluate the resulting formula
1.1F1algebra

Conjugate of the nonconvex Hamiltonian. For H(p)=p2−2∣p∣ and p≥0 one has pv−H(p)=p(v+2)−p2, whose supremum over p≥0 is (v+2)+2/4; for p≤0, writing p=−q with q≥0 gives pv−H(p)=−qv−q2+2q=q(2−v)−q2, whose supremum is (2−v)+2/4. Hence L(v)=max⁡{(v+2)+2/4,(2−v)+2/4}, where r+=max⁡(r,0); in particular L(0)=1, L(v)=(2+∣v∣)2/4≥1 for every v and L(v)≥4 for ∣v∣≥2, so inf⁡vL(v)=1, attained at v=0.

1.2F1F2algebra

Conjugate of the non-superlinear Hamiltonian. For H(p)=1+p2 and ∣v∣<1, the function p↦pv−1+p2 has derivative v−p/1+p2, which vanishes exactly at p=v/1−v2, where the value is −1−v2; at v=±1 the supremum is 0, approached along the infinite tail p→±∞; and for ∣v∣>1 the expression tends to +∞ along p=sgn⁡(v)r as r→∞. Hence L(v)=−1−v2 for ∣v∣≤1 and L(v)=+∞ for ∣v∣>1: the Legendre transform is not real-valued on all of R, finite action being available only when ∣x−y∣/t≤1, that is ∣x−y∣≤t. For u0≡0 the formula is Qtu0=tinf⁡vL(v)=−t; the residual of u(x,t)=−t is ut+H(ux)=−1+H(0)=−1+1=0, so as a formal expression it does satisfy the equation, and no failure of the solution property is claimed in this clause.

2.1step 1.1F2F3algebra

The formula is not a solution for the nonconvex Hamiltonian. With u0≡0, the Hopf--Lax formula of [F2] is Qtu0(x)=tinf⁡v∈RL(v)=t for every x and t>0 (the substitution v=(x−y)/t turns the infimum over y into the infimum over v). The function u(x,t)=t is smooth, and for a smooth function every point is both an upper and a lower contact with the test function itself; as an upper test, [F3] requires ut+H(ux)≤0, but ut=1 and ux=0 give 1+H(0)=1>0. Hence the Hopf--Lax output is not a viscosity subsolution, therefore not a viscosity solution, of ut+H(ux)=0.

3.1step 1.1step 1.2step 2.1∎

Conclusion. Part (a) exhibits a finite superlinear nonconvex Hamiltonian whose Hopf--Lax output fails the subsolution inequality pointwise, and part (b) exhibits a finite convex non-superlinear Hamiltonian whose Legendre transform is finite only on a bounded velocity set, so superlinearity cannot be dropped from the full-domain finiteness conclusion. The two failures are of different kinds: convexity is needed for the equation to hold, superlinearity for the Lagrangian to be finite everywhere.

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

A strict subsolution can fail the supersolution lower-test condition

Statement refuted

False claim: for a first-order equation ut+H(x,t,Du)=0, passing the subsolution test at every upper contact forces the supersolution test at every lower contact, so that a viscosity subsolution which is differentiable is automatically a viscosity solution.

The claim fails for the equation ut+H(ux)=0 on U=R×(0,∞) with H(p)=p (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem): the function u(x,t)=−t is a viscosity subsolution but not a viscosity supersolution. The numerical residual ut+H(ux)=−1 is the same at upper and lower contacts; what differs is the direction of the inequality that the definition requires.

Facts & Assumptions

Given: The open set U=R×(0,∞), the Hamiltonian H(p)=p, the equation ut+H(ux)=0, and u(x,t)=−t on U.

[F1]

A viscosity subsolution is tested at local maxima of u−ϕ and must satisfy ϕt+H(z0,Dϕ)≤0 there; a viscosity supersolution is tested at local minima of v−ϕ and must satisfy ϕt+H(z0,Dϕ)≥0 there; the test function ϕ is required to be C1, whereas u and v need only have their respective semicontinuity (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem).

[F2]

At a point where a differentiable function of several variables has a local maximum or a local minimum, its total derivative vanishes (Fermat's theorem: an interior differentiable local extremum has zero gradient).

Counterexample

technique · direct evaluation of the two contact inequalities
1.1F1F2algebra

u is a viscosity subsolution. Let ϕ∈C1(U) and let u−ϕ have a local maximum at z0=(x0,t0)∈U. The function u−ϕ is differentiable with total derivative Du−Dϕ, so [F2] gives Dϕ(z0)=Du(z0)=(0,−1), that is ϕx(z0)=0 and ϕt(z0)=−1. Hence ϕt(z0)+H(ϕx(z0))=−1+0=−1≤0, which is the subsolution inequality of [F1] at the upper contact.

1.2F1algebra

u is not a viscosity supersolution. Take the test function ϕ:=u, which belongs to C1(U). Then u−ϕ≡0 has a local minimum at every point of U, so the supersolution test of [F1] applies at, say, z0=(0,1); but ϕt(z0)+H(ϕx(z0))=−1+0=−1<0, and the required inequality is ≥0. Hence the supersolution condition fails, and u is not a viscosity solution of the equation on U.

2.1step 1.1step 1.2∎

Conclusion. Step 1.1 verifies every upper contact of u and step 1.2 exhibits a lower contact at which the opposite inequality fails, so the subsolution property does not imply the supersolution property; the residual is −1 in both computations, and only the required direction of the inequality changes between them.

Remarks

  • What the example isolates. The sign asymmetry of the test-function definition is not a matter of the value of the residual but of the direction of the inequality at the two kinds of contact. This is the reason the page defines the two one-sided notions separately and why the reverse implication is false for a monotone-in-time function.
CounterexampleConstruction: Literature-sourcedVerification: AI-adaptedprecheck passjudge pass (gpt-6.1-sol)Open item page →

Minima of viscosity subsolutions need not be subsolutions

Statement refuted

False claim: a finite minimum of finitely many viscosity subsolutions of a first-order equation is again a viscosity subsolution.

Take the stationary proper equation F(x,u,u′):=u+1−(u′)2=0on Ω=(−1,1), whose operator F(r,p)=r+1−p2 is continuous and strictly increasing in r (proper), and the two affine functions u1(x)=1−2x and u2(x)=1+2x. Both are classical subsolutions: F(ui,ui′)=ui+1−4=ui−3≤0 on Ω because ui≤3 there. Their pointwise minimum is u(x)=1−2∣x∣, a peak at x0=0 with value 1. The C1 test function ϕ≡1 satisfies ϕ≥u on Ω with equality at 0, so u−ϕ has a local maximum at 0; but F(0,u(0),ϕ′(0))=1+1−0=2>0, so u fails the subsolution test at the corner. Hence the finite-minimum operation is not admissible for subsolutions of a proper equation, while the finite maximum of subsolutions and the finite minimum of supersolutions are, by the same active-index contact argument as Finite maxima of subsolutions and finite minima of supersolutions: the active function has the same value and the same test at the contact. The zero-order term contributes to this residual but is not essential to failure of the minimum operation: for F(p)=1−p2 the same two affine functions have residual −3, while their minimum has the upper-test residual F(0)=1>0.

Facts & Assumptions

Given: The interval Ω=(−1,1), the proper operator F(r,p)=r+1−p2, the functions u1=1−2x, u2=1+2x, their pointwise minimum u=min⁡(u1,u2)=1−2∣x∣, and the constant test ϕ≡1.

[F1]

For a proper operator F(x,r,p) the viscosity subsolution inequality is F(x0,u(x0),Dϕ(x0))≤0 at every point x0 where u−ϕ has a local maximum, ϕ∈C1 (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem for the test-function scheme; here this displayed inequality defines the zero-order stationary adaptation; the cited item itself treats only the evolutionary operator with no zero-order term).

[F2]

A real-valued function has a local maximum at x0 when its value near x0 is at most its value at x0; the total derivative of a C1 function is computed from its partial derivatives; at a differentiable local extremum the derivative vanishes (Fermat's theorem: an interior differentiable local extremum has zero gradient, Local and strict local extrema for scalar fields on Euclidean open sets, 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).

Counterexample

technique · explicit touching test at the corner of the minimum
1.1F1F2algebra

The two affine functions are subsolutions. For i=1,2 and x∈(−1,1) we have ui(x)=1∓2x≤3, and ui is C1 with ui′=∓2, so F(ui,ui′)=ui+1−4=ui−3≤0 pointwise on Ω. Since at a local maximum of ui−ϕ the differentiability of ui and ϕ forces Dϕ(x0)=Dui(x0) by [F2], we get F(x0,ui(x0),Dϕ(x0))=ui(x0)−3≤0; hence u1 and u2 are viscosity subsolutions.

2.1step 1.1F1F2algebra

The minimum fails the test. The pointwise minimum is u=1−2∣x∣ with u(0)=1 and u(x)<1 for x≠0; the constant ϕ≡1 is C1 with ϕ′=0 and satisfies ϕ≥u on Ω with equality exactly at 0, so u−ϕ has a local maximum at 0 by [F2]. The subsolution test of [F1] at x0=0 would require F(0,u(0),ϕ′(0))=1+1−0=2≤0, which is false. Hence the finite minimum of the two subsolutions is not a subsolution.

3.1step 1.1step 2.1∎

Conclusion. Step 1.1 exhibits two viscosity subsolutions of the proper equation u+1−(u′)2=0 and step 2.1 shows their pointwise minimum fails the subsolution test at the peak; For this stationary operator, the finite-maximum and dual finite-minimum rules follow directly by choosing an active index at the contact: its function value is unchanged there and its test is the same. This is the argument of Finite maxima of subsolutions and finite minima of supersolutions, whose stated operator has no zero-order term. The computation F(p)=1−p2 with the same slopes and constant upper test also shows failure without a zero-order term.

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

The eikonal equation on an interval has many solutions when endpoint data are omitted

Statement refuted

False claim: the stationary eikonal equation ∣u′∣=1 on Ω=(0,1) has a unique viscosity solution when no boundary data are prescribed.

The functions u1(x)=x, u2(x)=1−x and u3(x)=12−∣x−12∣=min⁡{x,1−x} are viscosity solutions of ∣u′∣=1 on (0,1) (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem, Viscosity testing by first-order jets, and closure of the jet inequality). They are pairwise distinct, and their continuous boundary traces are respectively (0,1), (1,0) and (0,0); in particular the equation without prescribed endpoint data does not select a unique solution, and the examples are distinguished by their boundary traces.

Facts & Assumptions

Given: The interval Ω=(0,1), the equation ∣u′∣−1=0 written as Ut+F(Ux)=0 on U=Ω×(0,∞) with F(p)=∣p∣−1, the profiles u1(x)=x, u2(x)=1−x, u3(x)=min⁡{x,1−x}, and their time-independent lifts Ui(x,t)=ui(x).

[F1]

A viscosity subsolution is defined by ϕt+F(ϕx)≤0 at every local maximum of U−ϕ, and a supersolution by ϕt+F(ϕx)≥0 at every local minimum, for C1 test functions ϕ (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem).

[F2]

A function U is a viscosity subsolution if and only if every p=(px,pt)∈D+U satisfies pt+F(px)≤0, and a supersolution if and only if every p∈D−U satisfies pt+F(px)≥0 (Viscosity testing by first-order jets, and closure of the jet inequality).

Counterexample

technique · direct test-function verification at the peak and explicit evaluation of the boundary traces
1.1F2algebra

Time-independent profiles have zero test time-derivative. Let U(x,t)=u(x) for a continuous u:(0,1)→R and let ϕ∈C1(U) with U−ϕ having a local extremum at (x0,t0). Then, testing along the time line x=x0 and along the space line t=t0, the jet characterisation [F2] gives pt=0 for every p∈D+U(x0,t0) or p∈D−U(x0,t0): the inequality u(x0)−ϕ(x0,t)≤u(x0)−ϕ(x0,t0) near t0 forces ∂tϕ(x0,t0)=0 in the upper case (and dually in the lower case), and the remaining spatial inequality is the one-sided differentiability statement for u at x0. Consequently the viscosity conditions for U reduce to: every upper-test spatial slope px satisfies F(px)≤0, and every lower-test spatial slope px satisfies F(px)≥0.

2.1step 1.1F1F2algebra

The three profiles are solutions. For x0≠12 each ui is C1 near x0 with ∣ui′∣=1, so at any upper or lower test the spatial slope equals ui′(x0)=±1, and [F1] gives F(±1)=0 in both directions. At x0=12, only u3 is non-differentiable: writing h=x−12 we have u3(12+h)=12−∣h∣. If ϕ is an upper test with slope p=ϕx(12,t0), then u3≤ϕ near the point gives for h>0 the inequality −h+o(h)≤ph+o(h), hence p≥−1, and for h<0 the inequality h+o(∣h∣)≤ph+o(∣h∣), hence p≤1; thus every upper-test slope lies in [−1,1] and F(p)≤0. If ϕ were a lower test, u3≥ϕ would give for h>0 that p≤−1 and for h<0 that p≥1, which is impossible; hence there is no lower test at the peak and the supersolution condition is vacuous. By step 1.1 the same reduction applies to u1 and u2 at every point. Hence all three profiles are viscosity solutions of ∣u′∣=1 on (0,1).

3.1step 2.1algebra∎

Distinctness and boundary traces. At x=14 the values are (u1,u2,u3)=(14,34,14) and at x=34 they are (34,14,14), which distinguishes every pair; also u1−u2=2x−1 changes sign on (0,1). The continuous extensions to [0,1] have boundary values u1(0)=0, u1(1)=1; u2(0)=1, u2(1)=0; and u3(0)=u3(1)=0. Thus the equation with no prescribed boundary data admits at least three distinct viscosity solutions, and the displayed candidates carry different endpoint traces.

Remarks

  • What this shows. Comparison and uniqueness on a bounded domain require the boundary condition to be imposed; without it the solution class is not a singleton even for the simplest non-smooth first-order equation. This is the reason the page states comparison and uniqueness on Rn or, on a bounded cylinder, with boundary data on the parabolic boundary.
ExampleConstruction: Literature-sourcedVerification: AI-adaptedprecheck passOpen item page →

Distance to the boundary solves the unit eikonal Dirichlet problem on the ball

Example

Let n≥1 and let Ω=B(0,1)⊆Rn be the open unit ball and let d(x):=1−∣x∣=dist⁡(x,∂Ω). Then d is Lipschitz with ∣Dd∣=1 at every x≠0, continuous on Ω‾ with d=0 on ∂Ω, and d is a viscosity solution of the stationary eikonal equation with zero boundary data: ∣Du∣=1 in Ω,u=0 on ∂Ω. At x≠0 the equation holds classically; at the centre the function has a peak, every C1 test function ϕ with d−ϕ having a local maximum at 0 satisfies ∣Dϕ(0)∣≤1 (the subsolution inequality), and no C1 test function has d−ϕ locally minimal at 0: a lower contact would require ⟨Dϕ(0),h⟩≤−∣h∣ for all h, which is impossible, so no such contact exists and the supersolution test is vacuous at the centre. Hence d is a viscosity solution. The example is the canonical nonsmooth boundary-value solution of the eikonal equation, obtained by cone comparisons at the centre rather than by an abstract existence theorem.

Verification

Given: n≥1, the open unit ball Ω=B(0,1), its boundary the unit sphere (Euclidean spheres and closed balls as subspaces of Rn), the distance function d(x)=1−∣x∣, the norm ∣⋅∣ (A norm on a real vector space, the induced metric, and the dictionary with the metric axioms) and the inner product (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[F1] The test-function definition of viscosity sub- and supersolutions of ∣Du∣−1=0, and its equivalent jet formulation (Viscosity subsolutions and supersolutions of a first-order equation and of the Cauchy problem, Viscosity testing by first-order jets, and closure of the jet inequality).

[F2] For x≠0 the map x↦∣x∣ is C1 with gradient x/∣x∣ of norm 1, as computed in The eikonal equation as a viscosity equation at a tip; hence d is C1 on Ω∖{0} with Dd(x)=−x/∣x∣ and ∣Dd(x)∣=1; also d is Lipschitz with constant 1 by the triangle inequality (A norm on a real vector space, the induced metric, and the dictionary with the metric axioms, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

Proof technique: classical verification away from the centre and cone comparisons at the centre.

1.1F2

Classical region. For x≠0 the function d is C1 near x with ∣Dd∣=1 by [F2], so it is a viscosity solution of ∣Du∣=1 on the punctured ball by the classical-consistency argument of The eikonal equation as a viscosity equation at a tip.

1.2F1F2algebra

Upper contacts at the centre. Let ϕ∈C1 with d−ϕ having a local maximum at 0; normalize ϕ(0)=d(0)=1. Then 1−∣x∣≤ϕ(x)=1+⟨Dϕ(0),x⟩+o(∣x∣) near 0, that is ⟨Dϕ(0),x⟩≥−∣x∣+o(∣x∣); substitute x=ae and x=−ae for a unit vector e, divide by a>0, and let a↓0 to obtain ∣⟨Dϕ(0),e⟩∣≤1. Taking e in the direction of Dϕ(0), when this gradient is nonzero, gives ∣Dϕ(0)∣≤1. Hence every upper test satisfies the subsolution inequality ∣Dϕ(0)∣−1≤0 at the centre.

1.3F1F2algebra

No lower contact at the centre. If d−ϕ had a local minimum at 0, the reversed inequality would give ⟨Dϕ(0),x⟩≤−∣x∣+o(∣x∣) for all x near 0; substituting x=ae and x=−ae, dividing by a>0 and taking a↓0 gives ⟨Dϕ(0),e⟩≤−1 and ⟨Dϕ(0),e⟩≥1, an impossibility. Hence no lower C1 test exists at the centre and the supersolution inequality holds vacuously.

2.1step 1.1step 1.2step 1.3F2∎

Boundary values and conclusion. For x∈Ω and y∈∂Ω, the reverse triangle inequality gives ∣x−y∣≥1−∣x∣. Equality is achieved at y=x/∣x∣ when x≠0, and at any unit vector when x=0, so d(x)=dist⁡(x,∂Ω). Since ∣x∣→1 along sequences approaching the unit sphere, the continuous extension of d to Ω‾ vanishes exactly on ∂Ω. Steps 1.2 and 1.3 give the subsolution inequality everywhere and the supersolution inequality everywhere (vacuously at the centre, classically elsewhere by step 1.1), so d is a viscosity solution of the unit eikonal equation with zero boundary data on the ball.

Sources