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.

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

Positivity, invariant domains and entropy-stable systems fluxes — Examples

1 · Prerequisites

2 · Summary

These worked cases exercise the exact assumptions of the companion page: near-vacuum Euler states on the symmetric HLL speed branch, a two-point Burgers entropy residual, a finite-node conservative positivity rescaling, and a coarse Burgers shock update that shows why semidiscrete entropy stability does not promise sharp shock resolution. The last example reports its forward-Euler state calculation without asserting a fully discrete entropy inequality.

3 · Logical flowchart

4 · Definitions, theorems and proofs

None yet.

5 · Examples, counterexamples and false statements

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Exact Tadmor entropy residual for a dissipative Burgers flux

Example

For Burgers' flux f(u)=u2/2 and entropy η(u)=u2/2, the entropy flux is q(u)=u3/3 and the entropy potential is ψ(u)=u3/6. The symmetric two-point flux fec(uL,uR)=uL2+uLuR+uR26 is entropy-conservative. Adding unit dissipation in entropy variables gives fes=fec−(uR−uL)/2; at the jump 1→0 its Tadmor entropy residual is exactly −1/2.

Facts & Assumptions

Given: Burgers' flux, the quadratic entropy pair, and the two-point flux formula. Verify compatibility, the entropy-conservative identity, and the dissipative residual at (uL,uR)=(1,0).

[F1]

Convex entropy pairs, entropy variables, and entropy potentials for systems: an entropy pair satisfies q′=η′f′, with v=η′ and ψ=vf−q.

[F2]

Tadmor entropy-conservative two-point flux: a continuous, consistent, symmetric two-point flux is entropy-conservative when [v]fec=[ψ] for every pair of states, with right-minus-left jumps.

[F3]

Entropy-stable flux from positive-semidefinite entropy-variable dissipation: fes=fec−D[v]/2 for symmetric positive-semidefinite D.

Verification

1.1givenF1algebra

Here η′′(u)=1, η′(u)=u, and q′(u)=u2=η′(u)f′(u), so this is a convex entropy pair. Thus v=u and ψ=uf−q=u3/6 by [F1].

2.1givenF1F2step 1.1algebra

The displayed fec is polynomial, hence continuous, is symmetric, and satisfies fec(u,u)=u2/2=f(u). Moreover, (uR−uL)fec(uL,uR)=uR3−uL36=ψ(uR)−ψ(uL), so [F2] applies using ψ=u3/6 from step 1.1.

3.1givenF2F3step 2.1algebra∎

Taking D=1 in [F3], at (1,0) we have fec=1/6, [v]=−1, and [ψ]=−1/6. Therefore [v]fes−[ψ]=(−1)(1/6+1/2)−(−1/6)=−1/2. In general the same substitution gives residual −12(uR−uL)2≤0, using the formula in step 2.1.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

Two Euler quadrature nodes rescaled to a positive margin

Example

With weights w1=w2=1/2, cell mean Uˉ=(1,0,1), and candidate conserved states V1=(1,4,1) and V2=(1,−4,1), apply the limiter to K7/5,1/2. It returns θ∗=1/4, giving states (1,1,1) and (1,−1,1) with positive pressure and the original cell mean.

Facts & Assumptions

Given: The weights, mean, nodes, and set K7/5,1/2. Compute the maximal scaling and verify the mean and positive margin.

[F1]

Maximal convex scaling toward an admissible cell mean: the limiter chooses the largest θ∈[0,1] keeping each Uˉ+θ(Vj−Uˉ) in K.

[F2]

Positive-margin admissible states for ideal-gas Euler: Kγ,ε requires ρ≥ε and E−m2/(2ρ)≥ε.

[F3]

Admissible states and pressure for one-dimensional ideal-gas Euler: p=(γ−1)(E−m2/(2ρ)).

[F4]

Convex scaling preserves the cell mean and positive-margin node states: the maximal scaling exists, preserves the weighted mean, and returns nodes in the positive-margin set.

Verification

1.1givenF2algebra

The weighted mean is (V1+V2)/2=(1,0,1). Its density is 1 and its internal-energy density is 1, so it belongs to K7/5,1/2 by [F2]. The scaled states are (1,4θ,1) and (1,−4θ,1), each with internal energy 1−8θ2. Hence they lie in K7/5,1/2 exactly when 1−8θ2≥1/2, or 0≤θ≤1/4.

2.1givenF1F2F3F4algebra∎

The largest feasible parameter is θ∗=1/4 by [F1]. The outputs are (1,1,1) and (1,−1,1); each has internal-energy density 1/2 and pressure (2/5)(1/2)=1/5>0 by [F3]. Their weighted mean is (1,0,1), also as guaranteed by [F4].

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

An entropy-stable Burgers flux still smears an under-resolved shock

Example

For the Burgers flux and quadratic entropy, the unit-dissipation entropy-stable flux has a negative Tadmor residual at the entropy shock 1→0. On the four-cell periodic data (1,1,0,0) and with ν=Δt/Δx=1/2, one explicit forward-Euler update replaces the two states adjacent to that jump by 11/12 and 1/3. Thus semidiscrete entropy stability coexists with visible shock smearing at this coarse resolution. This computation gives no fully-discrete entropy inequality for forward Euler.

Facts & Assumptions

Given: Burgers' flux f(u)=u2/2, the entropy-stable flux from the preceding example, the four periodic cell states, and ν=1/2. Verify entropy admissibility of the exact shock, the semidiscrete flux residual, and the displayed finite update.

[F1]

Convex entropy pairs, entropy variables, and entropy potentials for systems: for η=u2/2 and q=u3/3, v=u and ψ=u3/6.

[F2]

Tadmor entropy-conservative two-point flux: consistency and [v]fec=[ψ] define entropy conservation.

[F3]
[F4]

Semidiscrete cell entropy inequality for Tadmor-dissipative fluxes: this flux gives the semidiscrete cell entropy inequality, not a fully discrete time-step claim.

Verification

1.1givenF1algebra

The Rankine–Hugoniot speed from states uL=1,uR=0 is (f(0)−f(1))/(0−1)=1/2. Since f′(0)=0<1/2<1=f′(1), it is a compressive shock. Its entropy production is q(0)−q(1)−12(η(0)−η(1))=−1/3+1/4=−1/12<0, so it satisfies the quadratic entropy inequality.

2.1givenF1F2F3step 1.1algebra

The entropy-conservative Burgers flux is fec=(uL2+uLuR+uR2)/6: it is symmetric and diagonal-consistent, and (uR−uL)fec=(uR3−uL3)/6=[ψ], so [F2] applies. Hence fec(1,0)=1/6 and [F3] gives fes(1,0)=2/3. With [u]=−1 and [ψ]=−1/6 from step 1.1, the Tadmor residual is [u]fes−[ψ]=−2/3+1/6=−1/2<0.

3.1givenF1F2F3F4step 2.1algebra∎

Equal-state faces have flux fes(1,1)=1/2 and fes(0,0)=0 by [F2]–[F3]. Thus the left state adjacent to the jump updates to 1−12(2/3−1/2)=11/12, and the right state updates to 0−12(0−2/3)=1/3. The jump is spread over two intermediate values. This is a finite computation only; [F4] does not certify the forward-Euler step's entropy inequality. The interface flux value was obtained in step 2.1.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

A near-vacuum Euler state remains positive under symmetric HLL speeds

Example

Take γ=7/5, density ρ=10−8, pressure p=10−12, and velocities uL=−1/100, uR=1/100. The conserved states have mL,R=∓10−10 and EL,R=3×10−12. Choose the symmetric HLL speed a=11/500 at every face and Courant ratio ν=Δt/h=125/11. For the four-cell periodic data (UL,UL,UR,UR), one HLL forward-Euler update preserves positive density and pressure despite the very small initial values. The calculation illustrates a first-order positivity branch; it does not claim accuracy of the shock profile.

Facts & Assumptions

Given: The ideal-gas Euler flux and HLL update, with the parameters above. Verify the states, speed branch, CFL condition, and positivity after one update.

[F1]

Admissible states and pressure for one-dimensional ideal-gas Euler: p=(γ−1)(E−m2/(2ρ)) and c2=γp/ρ.

[F2]

Ideal-gas Euler states remain positive under Lax–Friedrichs flux splitting: if a≥∣u∣+c, then U±F(U)/a∈Gγ.

[F4]

First-order HLL updates are averages of adjacent half-cell fan averages: under 2νa≤1, the HLL update is the mean of adjacent half-cell fan averages.

[F5]

Symmetric HLL speeds preserve ideal-gas density and internal-energy positivity: the symmetric HLL branch with a≥∣u∣+c and 2νa≤1 preserves Gγ.

[F6]

First-order finite-volume update using supplied HLL speeds: Ui+=Ui−ν(F^i+1/2−F^i−1/2).

[F7]

HLL two-wave approximate flux: the symmetric middle-branch HLL flux is used at every interface, with equal-state data substituted into its displayed formula.

Verification

1.1givenF1algebra

The kinetic energy density is ρu2/2=5×10−13 and p/(γ−1)=5/2×10−12, so E=3×10−12. Hence [F1] gives the specified pressure and positive internal energy, so both states lie in Gγ.

2.1givenF1F2F5step 1.1algebra

Their sound speed satisfies c2=7/50000<(3/250)2, hence ∣u∣+c<1/100+3/250=11/500=a. Also νa=1/4, so 2νa=1/2≤1. The symmetric-speed hypotheses in [F2]–[F5] hold, with the positive admissibility verified in step 1.1.

3.1givenF1F2F3F4F6F7step 2.1algebra∎

At the discontinuity, substitution of SL=−a,SR=a into [F7] and the HLL intermediate-state formula of [F4] gives U∗=(UL+FL/a+UR−FR/a)/2. Both split states are in Gγ by [F2], so U∗∈Gγ by [F3]. At equal-state faces, the same middle-state formula gives U∗=UL or UR. Substituting the shared face fluxes in [F6], or equivalently using [F4], yields Ui+=(3/4)UL+(1/4)U∗ and Ui+1+=(3/4)UR+(1/4)U∗. Convexity keeps both states in Gγ, and [F1] gives positive density and pressure, using the CFL verification in step 2.1.

Sources