Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-generatedPipeline-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.

First-order finite-volume convergence for compatible smooth renewal data

Statement

Fix T>0 and a cutoff A>T. Let nA be the finite-age solution from Characteristic solution and mass bound for the finite-age renewal model with initial density n0A∈C2([0,A]) and assume μ,β∈C2([0,A]) are nonnegative and bounded. Assume nA extends to C2([0,T]×[0,A]). Its smoothness at the renewal corner includes n0A(0)=BA(0) and −(n0A)′(0)−μ(0)n0A(0)=BA′(0); no boundary value is imposed at the absorbing outflow. Take integers J,N≥1, set Δa=A/J and Δt=T/N, and use the corresponding uniform age and time grids. Initialize with exact cell masses and use the one-panel midpoint rule in each cell from Positive composite or capped-adaptive quadrature for renewal births. Assume Δt/Δa∈[c−,c+] for fixed 0<c−≤c+<∞ and the transport-plus-loss CFL. Define

ϵQ:=max⁡0≤k<N∣Bh(M‾(tk))−∫0Aβ(a)nA(tk,a) da∣,

where Bh(M‾(tk)) is the birth-rate quadrature applied to the exact cell-average masses M‾j(tk)=∫CjnA(tk,a) da. Then the cell-mass error Ek=∑j∣Mjk−∫CjnA(tk,a) da∣ obeys

max⁡0≤k≤NEk≤CT(Δa+Δt+ϵQ).

The piecewise-constant reconstruction uhk=Mjk/Δa on cell Cj and zero above A satisfies the same L1 error order. For the midpoint rule, ϵQ≤CTΔa under the stated regularity, so the method is first order in age mesh and time step. If the capped adaptive rule is used instead, the estimate requires an independently established true ϵQ; its reported indicator is not substituted for that hypothesis.

Facts & Assumptions

Given: The smooth finite-cutoff solution and finite-volume scheme from Finite-age renewal model with absorbing outflow, Characteristic solution and mass bound for the finite-age renewal model and Finite-volume upwind scheme for age renewal transport, under the hypotheses in the Statement.

[F1]

Finite-age renewal model with absorbing outflow: nt+na=−μn on 0<a<A, the renewal inflow is BA(t)=∫0Aβ(a)nA(t,a) da, and the upper boundary is absorbing outflow.

[F2]

Characteristic solution and mass bound for the finite-age renewal model: nA follows characteristics, satisfies the renewal relation, and its exact initial-cohort outflow is absorbing.

[F3]

Finite-volume upwind scheme for age renewal transport: the cell-mass step is Mjk+1=Mjk−Δt(Fjk−Fj−1k)−ΔtμjMjk with upwind face fluxes and absorbing top outflow.

[F4]

Positivity, L1 growth and data stability of the renewal finite-volume step: for the same coefficients and positive quadrature, a one-step mass difference is amplified by at most 1+Δtβ∞.

[F5]

Composite midpoint error bound for a twice continuously differentiable function: for f∈C2[a,b], the composite midpoint error on width h is at most (b−a)h2max⁡∣f′′∣/24.

[F6]

Finite geometric error bound from a one-step recurrence: if e0≤C0h and ek+1≤(1+Lh)ek+Ch2 on a fixed grid Nh=T, then max⁡kek≤CTh.

Proof

1.1F1F2constructalgebra

Let M‾j(t)=∫CjnA(t,a) da and u‾j(t)=M‾j(t)/Δa. The smooth finite-age characteristic solution [F2] satisfies the transport equation and renewal boundary in [F1]. Integrating [F1] over one cell and one time step gives M‾j(tk+1)=M‾j(tk)−∫tktk+1[nA(s,aj)−nA(s,aj−1)]ds−∫tktk+1∫Cjμ(a)nA(s,a) da ds, with nA(s,a0)=BA(s) and the outgoing trace at aJ=A.

2.1givenF5step 1.1algebra

Taylor's formula in age gives u‾j(t)−nA(t,aj)=−Δa2(nA)a(t,aj)+O(Δa2) uniformly. For interior cells, subtracting adjacent face-flux errors cancels the leading term and leaves O((Δa)2) per cell because (nA)aa is bounded. In the first and last cells there is only one uncancelled endpoint face error, of size O(Δa); these two cells contribute O(ΔtΔa) in total. Taylor's formula in time gives a face-integral remainder O(Δt2Δa) per cell after adjacent differences, using bounded (nA)ta, hence O(Δt2) after summing. At the inlet, the smooth trace gives (nA)t(t,a0)=BA′(t), so the time derivative of nA(t,a1)−BA(t) is also bounded by Δa∥(nA)ta∥∞; thus the boundary birth flux has no additional O(\Delta t^2) term beyond this face remainder and ΔtϵQ. For the reaction term, first note that u‾j(t)=nA(t,ajc)+O((Δa)2) by the symmetric cell average. Thus replacing μ(a)nA(t,a) by μ(ajc)u‾j(t) has per-cell error O((Δa)3): the midpoint formula [F5] bounds the error for μnA, and the cell-average replacement changes the midpoint value by O((Δa)2) times the cell width. Summing and multiplying by Δt gives O(Δt(Δa)2). Freezing the reaction term at tk instead of integrating in time contributes O(Δt2) after summing cells, since (μnA)t is bounded. Therefore ∥ρk∥ℓ1≤CΔt(Δa+Δt+ϵQ) for sufficiently fine meshes with Δa≤1.

3.1F3F4step 2.1algebra

Set ek=∑j∣Mjk−M‾j(tk)∣. Apply the one-step stability [F4] to the numerical and exact cell-average states and add the local defect from step 2.1; this yields ek+1≤(1+β∞Δt)ek+CΔt(Δa+Δt+ϵQ), with e0=0.

4.1F5F6step 3.1algebra∎

Put D=Δa+Δt+ϵQ. Iterating the recurrence in step 3.1 gives ek≤CΔtD∑r=0k−1(1+β∞Δt)r≤CTD, where the bound also holds when β∞=0, in which case the sum is k≤T/Δt. For the midpoint birth rule, split each cell defect into u‾j[QM(β;Cj,1)−∫Cjβ] and ∫Cjβ(a)(u‾j−nA(tk,a)) da. The first term is O((Δa)3) per cell by [F5] and bounded nA; in the second, subtract β at the cell midpoint because ∫Cj(u‾j−nA)=0, then use ∣β(a)−β(ajc)∣≤Lip⁡(β)Δa and ∫Cj∣u‾j−nA∣=O((Δa)2) from bounded (nA)a. Summing gives ϵQ=O((Δa)2), hence in particular O(Δa). With the two-sided mesh ratio, D=O(Δt); the recurrence then has additive term C′Δt2, and [F6] gives the stated first-order estimate. Finally, on each cell ∫Cj∣nA(tk,a)−u‾j(tk)∣da≤(Δa)2∥(nA)a∥∞; summing over A/Δa cells adds only O(Δa) to the cell-mass error, proving the reconstruction claim.

Depends on

Used by

Dependency tree · two levels

27 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources