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.

Continuum Heat Model

1 · Prerequisites

2 · Summary

This page develops two conditional descriptions of heat transport on explicitly stated domains. The local continuum branch fixes a rigid rectangle and regular positive-temperature fields. Its energy balance, Fourier flux and caloric entropy density are adopted physical assumptions. They imply a local nonnegative entropy-production identity; the constant-capacity, constant-conductivity specialization has a Gaussian Neumann solution for the stated even-periodic C2 initial data, and its mean-zero variance decays at an explicit exponential rate.

The kinetic branch is a separate periodic, single-species model. Incoming molecular chaos and the weak elastic collision equation are postulates. Subject to the theorem's explicit regularity, entropy and transport integrability, spatial collision-production, pointwise weak-identity and transport-equation conditions, and its current mathematical collision-entropy supplier, the kinetic entropy is nondecreasing. This conditional result does not derive the kinetic equation from many-particle mechanics, identify kinetic entropy with the local continuum entropy, or derive Fourier's law from collisions.

The page prerequisites remain those declared by the plan: td-entropy-foundation, td-math-local-pfaff-geometry, and td-math-mixture-dynamics-and-shells. The entropy foundation supplies the macroscopic entropy context. The geometry and mixture pages are page-level track prerequisites; this page adds no item-level claim from them that is not used in its proofs. Mathematical assumptions needed locally, including the rectangular Poincare estimate and the specified heat-kernel interface, are stated and cited at their use sites.

The order above follows the in-run supplier graph: field and collision data, adopted balance and constitutive assumptions, the local estimate, then their consuming results. The extra local caloric-entropy postulate and rectangular Poincare lemma are explicit local prerequisites rather than omitted theory. The constant-coefficient existence result is deliberately limited to its stated initial-data class; the final remark records that boundary of scope.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Rectangular local-equilibrium field data

Physical scope

Smooth local-equilibrium continuum data on a fixed three-dimensional Euclidean rectangle; this definition adopts no balance law or constitutive law.

Definition

Let I⊆R be an open interval measured in seconds, let L1,L2,L3>0 be lengths in metres, and set Ω=∏i=13(0,Li)⊆R3 with closure Ω‾. A rectangular local-equilibrium field record consists of an open temperature interval T⊂(0,∞), a temperature field T:Ω‾×I→T whose time derivative through order one and spatial derivatives through order two extend continuously to Ω‾×I, a function u∈C1(T,R) of internal-energy density in J m−3, and its volumetric heat capacity c(T):=u′(T)>0 in J m−3 K−1. The coordinates xi are lengths in metres, and T is measured in kelvins, using SI quantities and unit conventions.

The record may also include a thermal-energy flux-density field jQ:Ω‾×I→R3 with components in W m−2 and a conductivity field K:Ω‾×T→R3×3 with entries in W m−1 K−1. For any unit vector n, jQ⋅n is the signed thermal-energy transfer rate per unit area in direction n, positive when energy flows along n. The flux and tensor are field data here; no relation between them is imposed by this Definition. In Euclidean coordinates, ∇T=(∂x1T,∂x2T,∂x3T). The classical divergence div⁡jQ(x,t):=∑i=13∂xi(jQ)i(x,t) is defined at interior points where these three partial derivatives exist. No spatial differentiability of an optional flux is assumed here; any balance or result using its divergence throughout a domain must impose the needed regularity. On a flat face, n denotes its outward unit normal.

This is a specified regular field class on a fixed rigid body. It does not assert that every physical medium has local-equilibrium fields, nor does it adopt a heat equation, a boundary condition, Fourier's law, or an entropy-production law.

Remarks

The label “local-equilibrium field record” names this model's field class; it is not a claim about microscopic equilibration.

DefinitionDefinition: AI-adaptedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Regular single-species elastic collision data

Physical scope

A three-dimensional single-species classical kinetic collision model with elastic binary events; no hard-sphere limit or microscopic derivation is included.

Definition

Fix a particle mass mp>0 in kilograms and velocity space V=R3 with components in m s−1. The elastic collision set is Qel={(v,w,v′,w′)∈V4:v+w=v′+w′, ∣v∣2+∣w∣2=∣v′∣2+∣w′∣2}. Thus each allowed event conserves momentum mp(v+w) and kinetic energy mp(∣v∣2+∣w∣2)/2. Equip Qel with the Borel sets generated by its relative topology.

Collision-rate data consist of a nonnegative, countably additive, σ-finite Borel measure ν on these Borel sets. It is invariant under exchanging v with w, exchanging v′ with w′, and exchanging the incoming pair (v,w) with the outgoing pair (v′,w′): for each of these maps S and Borel set E, ν(S−1E)=ν(E). This is the three-dimensional physical restriction of the collision set and exchange conventions in Elastic collision model and weak collision identity, with σ-finiteness stated explicitly for the product integrals used by the entropy theorem. For physical applications of that weak identity, at a fixed position and time let f:V→(0,∞) be the one-particle number density per spatial volume and velocity volume: f(v) d3v is number per spatial volume in the velocity element. Particle number is dimensionless, d3v has units m3 s−3, and f has units s3 m−6. A collision term Qν(f), when a representative satisfying the absolutely convergent weak identity exists, has units s2 m−6, namely f per second. Test functions φ:V→R are measurable and have one common unit Uφ for their values (dimensionless for the tests 1 and log⁡(f/fref), where fref>0 has the units of f). Assign the measure ν itself units m9 s−7: its numerical values are taken in these SI units, not as a density relative to an unstated measure on Qel. Thus ∫VφQν(f) d3v and the collision integral of f2 times the test bracket both have units Uφ m−3 s−1. The magnitude and distribution of ν remain model inputs; dimensional consistency fixes neither a numerical collision rate nor existence of Qν(f). The measure is a declared model input; no differential cross section or microscopic limit is inferred from this Definition.

PostulateStatement: AI-adaptedProof: Not applicable‡ sources checked 2026-10-15Open item page →

Fourier heat-flux constitutive model

Physical scope

An idealized local-equilibrium continuum regime in which the first-order Fourier relation is adopted. No quantitative gradient threshold or empirically validated coefficient range is specified; the conductivity may be anisotropic and state dependent.

Postulate

For the field record of Rectangular local-equilibrium field data, adopt jQ(x,t)=−K(x,T(x,t))∇T(x,t). Here K has entries in W m−1 K−1 and jQ in W m−2, so the equation is dimensionally consistent. Write Ks=(K+KT)/2 for the symmetric part and assume zTKs(x,T)z≥0 for every z∈R3. Define the thermal quadratic density by σQ(x,t):=−jQ(x,t)⋅∇T(x,t)T(x,t)2=(∇T(x,t))TKs(x,T(x,t))∇T(x,t)T(x,t)2, with units W m−3 K−1. For the antisymmetric part Ka=(K−KT)/2, the scalar zTKaz equals its own negative by transposition, so it is zero. Thus Ka may contribute to jQ but not to σQ, and σQ≥0. Interpreting this density as entropy production additionally requires a local entropy map s∈C1(T,R) in J m−3 K−1 with s′(T)=u′(T)/T, the energy balance ∂tu+div⁡jQ=0, and regularity permitting the product and chain rules; under those additional assumptions the entropy balance is ∂ts(T)+div⁡(jQ/T)=σQ. This postulate alone adopts neither that caloric entropy relation nor the energy balance.

This is an adopted constitutive model in the stated continuum regime. It is neither a consequence of equilibrium entropy alone nor a universal law for all media. In particular, it does not assert a finite propagation speed or an empirical value for K.

PostulateStatement: AI-adaptedProof: Not applicable‡ sources checked 2026-10-15Open item page →

Local energy balance for an insulated continuum body

Physical scope

A rigid, fixed-composition, nonrelativistic continuum body on the rectangular domain of def-td-rectangular-local-equilibrium-fields, with no volumetric energy source, no matter flux, and insulated outer faces.

Postulate

For a field record on Ω×I from Rectangular local-equilibrium field data, require the supplied heat-flux field jQ and each of its first spatial partial derivatives to extend continuously to Ω‾×[t0,t1] for every compact interval [t0,t1]⊂I. Thus div⁡jQ is the classical coordinate divergence in Ω×I, and jQ⋅n denotes its continuous trace on each closed face. For these records, adopt the local balance ∂tu(T(x,t))+div⁡jQ(x,t)=0 and the insulated-boundary condition jQ(x,t)⋅n(x)=0 on every face of ∂Ω. Both terms in the local equation have units W m−3. The boundary condition says that no heat energy crosses the body's outer faces in this model.

This local continuum balance is an additional adopted postulate, with the sign and energy convention of First-law energy balance for closed process records. It is not deduced from that global process-record balance. Matter transfer, volume work, internal sources, moving boundaries, and external-field energy require separately specified balance models and are outside this item.

PostulateStatement: AI-adaptedProof: Not applicable‡ sources checked 2026-10-15Open item page →

Local-equilibrium caloric entropy density

Physical scope

Fixed-composition local-equilibrium patches of the rigid continuum model; this caloric relation is not asserted for arbitrary nonequilibrium states or every thermodynamic material.

Postulate

For the local-equilibrium medium of Rectangular local-equilibrium field data, adopt a local entropy-density map s:T→R measured in J m−3 K−1 with dsdT(T)=c(T)T, where c(T)=u′(T)>0 is the volumetric heat capacity and T>0 is measured in kelvins. Equivalently, for a fixed reference Tref∈T and reference density sref∈R in J m−3 K−1, set s(T)=sref+∫TrefTc(θ)θ dθ. The volume integral of this field is the entropy assigned by this continuum model. This adopts a local thermodynamic relation extending the macroscopic entropy framework of Concave macroscopic fundamental entropy; it does not prove that every evolving medium has such a local entropy density or identify it with an arbitrary finite-system entropy function.

For constant c(T)=c0>0, the adopted relation gives s(T)=sref+c0log⁡(T/Tref). This is a restricted constitutive specialization, not a universal equation of state.

PostulateStatement: AI-adaptedProof: Not applicable‡ sources checked 2026-10-15Open item page →

Periodic Boltzmann model with incoming molecular chaos

Physical scope

A single-species periodic Boltzmann kinetic model with elastic binary collision data and the stated incoming-pair molecular-chaos closure. Dilute-gas applicability is qualitative; no density threshold, hard-sphere or Boltzmann–Grad limit, or many-particle derivation is asserted.

Postulate

Let Ωper=∏i=13(0,Li) with opposite faces identified, let I⊆R be an open time interval in seconds, and let f(x,v,t)≥0 be a one-particle number density for x∈Ωper, v∈R3, and t∈I. With x in metres and v in m s−1, f has units s3 m−6 relative to d3x d3v; its moments used below are assumed finite. For this postulate, require f to be jointly measurable and periodic in x, with (x,t)↦f(x,v,t) continuously differentiable for each fixed v and its derivatives jointly measurable. Require f, ∂tf, and v⋅∇xf to be locally integrable with respect to product Lebesgue measure. Periodic boundary conditions are imposed in x.

Adopt the incoming molecular-chaos closure: at a collision location and time, the incoming two-particle velocity density used in the collision term is f(x,v,t)f(x,w,t). This is a one-sided factorization assumption; no factorization of the outgoing pair and no derivation from reversible many-particle mechanics is asserted.

For the elastic data Qel,ν of Regular single-species elastic collision data, adopt a collision operator Qν(f), with units of f per second, as a jointly measurable locally integrable representative whose velocity-wise weak action satisfies, for almost every (x,t) and every measurable test function φ for which both sides converge absolutely, ∫R3φ(v)Qν(f)(x,v,t) d3v=14∫Qel(f(x,v′,t)f(x,w′,t)−f(x,v,t)f(x,w,t))(φ(v)+φ(w)−φ(v′)−φ(w′)) dν. The adopted kinetic equation is ∂tf+v⋅∇xf=Qν(f), as an equality of locally integrable functions almost everywhere on the product domain. These are admissibility conditions for solutions of the adopted model, not an existence assertion. The equation and weak collision identity are model postulates. This item does not assert a hard-sphere Boltzmann–Grad limit, propagation of chaos, or universal equilibration.

LemmaStatement: AI-adaptedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Mean-zero Poincare inequality on a rectangle

Statement

Let d≥1, let L1,…,Ld>0, and let Ω=∏i=1d(0,Li). If w∈C1(Ω‾,R) satisfies ∫Ωw(x) dx=0, write ∇w=(∂1w,…,∂dw) and ∣∇w∣2=∑i=1d∣∂iw∣2. Then ∫Ω∣w(x)∣2 dx≤d2∑i=1dLi2∫Ω∣∂iw(x)∣2 dx≤CΩ∫Ω∣∇w(x)∣2 dx, where CΩ=d2max⁡iLi2. All integrals are the multidimensional Riemann integrals of continuous functions on the indicated closed rectangles, with the usual extension to their interiors.

Facts & Assumptions

Given: A dimension d≥1, positive side lengths Li, the open rectangle Ω=∏i(0,Li), and w∈C1(Ω‾,R) with ∫Ωw=0.

[A1]

The zero-mean condition is ∫Ωw(x) dx=0.

[F1]

Newton–Leibniz needs only continuity on [a,b], differentiability on (a,b), and a Riemann-integrable extension of the interior derivative: if G is continuous on [a,b], differentiable on (a,b), and its interior derivative has a Riemann-integrable extension f, then ∫abf=G(b)−G(a).

[F2]

Riemann--Fubini on product rectangles, with lower and upper section integrals and content-zero exceptional sections: when every section of a Riemann-integrable function on a product of nondegenerate closed rectangles is integrable, the ordinary iterated integral equals the multiple integral. In particular, this applies to the continuous integrands used here.

Proof

1.1A1F2algebra

Expand the square in the product integral and use [A1]: ∫Ω∫Ω(w(x)−w(y))2 dy dx=2∣Ω∣∫Ωw2−2(∫Ωw)2=2∣Ω∣∫Ωw2. The integrand is continuous, so [F2] justifies the repeated integrals.

1.2givenalgebra

For x,y∈Ω, define z(0)=x and z(i)=(y1,…,yi,xi+1,…,xd) for 1≤i≤d, and set δi=w(z(i−1))−w(z(i)). Then w(x)−w(y)=∑iδi and d∑iδi2−(∑iδi)2=∑i<j(δi−δj)2≥0.

1.3F1algebra

Fix i and the other coordinates. By [F1], δi is the integral of h=∂iw over the interval between xi and yi. If the endpoints differ, let hˉ be the average of h on that interval. Nonnegativity of ∫(h−hˉ)2 gives (∫h)2≤∣xi−yi∣∫h2; if they agree, δi=0. Since ∣xi−yi∣≤Li, in both cases δi2≤Li∫0Li∣∂iw(z1,…,zi−1,s,zi+1,…,zd)∣2 ds.

2.1F2step 1.3algebra

Integrate the last bound over x,y. The xi,yi integrations contribute Li2, the displayed segment bound contributes Li, each unused variable from a pair xj,yj with j≠i contributes Lj, and the remaining variables and s give ∫Ω∣∂iw∣2. By [F2], this yields ∫Ω∫Ωδi2 dy dx≤∣Ω∣Li2∫Ω∣∂iw∣2.

3.1step 1.1step 1.2step 2.1algebra∎

By step 1.2 and the finite squared-sum inequality, the double integral in step 1.1 is at most d∑i∫Ω∫Ωδi2 dy dx. Apply step 2.1, divide by 2∣Ω∣>0, and use ∑iLi2∣∂iw∣2≤(max⁡iLi2)∣∇w∣2 pointwise. This gives both stated inequalities with the displayed CΩ.

Physics theoremStatement: AI-adaptedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Conditional entropy production in the periodic Boltzmann model

Physical scope

Conditional kinetic-entropy production for a positive periodic solution satisfying the stated regularity, entropy and transport integrability, spatial collision-production, and pointwise weak-identity hypotheses; no many-particle derivation or Maxwellian equality characterization is asserted.

Statement

Let Ωper=∏i=13(0,Li) with opposite faces identified, and let I⊆R be an open time interval. Assume the current conditional mathematical result Conditional Boltzmann entropy inequality and equality form at its exact draft interface [L1] below. Let f:Ωper×R3×I→(0,∞) satisfy the periodic Boltzmann postulate Periodic Boltzmann model with incoming molecular chaos and the regularity, integrability, and domination hypotheses in [A2]–[A3]. For the all-time conclusion below, require the kinetic equation ∂tf+v⋅∇xf=Qν(f) to hold at every (x,v,t) for the chosen collision representative, in addition to the postulate's product-almost-everywhere equation. Fix a positive constant reference density fref with the same units as f, and define the kinetic entropy Skin(t)=−kB∫Ωper∫R3f(x,v,t)log⁡ ⁣(f(x,v,t)fref) d3v d3x. Then dSkindt=kB4∫Ωper∫Qel(A−B)log⁡ ⁣(AB) dν d3x≥0, where A=f(x,v,t)f(x,w,t) and B=f(x,v′,t)f(x,w′,t); log⁡(A/B) is dimensionless because A and B have the same units. Equality holds exactly when A=B for ν-almost every collision at almost every spatial point. This is a conditional deduction inside the adopted incoming-molecular-chaos kinetic model. The cited mathematical carrier is currently draft and has not received independent proof acceptance; no microscopic derivation, experimental claim, or Maxwellian characterization is included.

Facts & Assumptions

Given: The positive regular periodic solution and collision data described in [A1]–[A3], Boltzmann's constant kB>0, and the exact current mathematical supplier interface [L1].

[A1]

Periodic Boltzmann model with incoming molecular chaos: on the periodic three-torus, the adopted one-sided incoming-pair closure gives an admissible locally integrable collision term Qν(f), the kinetic equation ∂tf+v⋅∇xf=Qν(f) almost everywhere, and the weak collision identity with factor 1/4 for absolutely convergent test functions; it is a model postulate, not a many-particle theorem.

[A2]

The solution satisfies f>0, is C1 in (x,t) and C2 in v, and is periodic in each spatial coordinate. The reference fref>0 is constant in x,v,t, so ℓ=log⁡(f/fref) and log⁡(A/B) are dimensionless.

[A3]

For each t∈I, f, f∣v∣2, and f∣ℓ∣ are integrable over Ωper×R3. For every compact interval I0⊂I, assume a common GI0∈L1(Ωper×R3) bounds ∣(1+ℓ)∂tf∣ almost everywhere for every t∈I0. Also assume (1+∣ℓ∣)∣v⋅∇xf∣∈L1(Ωper×R3) for each t∈I. For every t∈I, the map x⟼∫Qel∣(A−B)log⁡(A/B)∣ dν is measurable and belongs to L1(Ωper). In addition to the almost-everywhere weak identity in [A1], assume that the identities for tests 1 and ℓ(x,⋅,t) are absolutely convergent and hold at every (x,t), as required by the exact current interface [L1]. The collision measure is σ-finite. Also require ∂tf+v⋅∇xf=Qν(f) at every (x,v,t) for this representative, so substitution at each fixed time is licensed.

[F1]

Regular single-species elastic collision data: in dimension three, Qel consists of collisions conserving pair momentum and squared speed; ν is a nonnegative σ-finite Borel measure invariant under the two particle exchanges and incoming/outgoing exchange. The σ-finiteness here is an explicit physical restriction that makes the spatial/collision product integrals well-defined.

[F2]

SI quantities and unit conventions: kB is the positive Boltzmann constant with SI units J K−1.

[L1]

Conditional Boltzmann entropy inequality and equality form: the current draft supplier assumes a positive C1-in-space-and-time, C2-in-velocity density on a periodic rectangle; finite moments and entropy; a compact-time L1 majorant for (1+ℓ)∂tf; transport integrability at each time; an L1 spatial collision-production map; and absolutely convergent weak identities for tests 1 and ℓ at every point. It states ddtH[f](t)=−14∫Ω∫Qel(A−B)log⁡(A/B) dν dr≤0, with equality exactly when A=B for ν-almost every collision at almost every spatial point. The carrier remains draft without independent proof acceptance; this item uses only that inequality and equality condition, not the supplier's additional Maxwellian branch.

Proof

Given: [A1]–[A3], [F1]–[F2], and [L1] with its current stated proof status.

1.1A1A2A3F1L1

Set H(t)=∫Ωper∫R3fℓ d3v d3x. The hypotheses match [L1]: the spatial dimension is d=3, the opposite-face identification is periodic, the density is positive with the required regularity and moments, [A3] supplies the compact-time derivative majorant, transport integrability, pointwise weak identities, and spatial collision-production integrability, and [A3] supplies the kinetic equation at every time for the chosen representative. The collision measure is symmetric by [F1] and σ-finite by [A3]. Thus [L1] gives H′(t)=−14∫Ωper∫Qel(A−B)log⁡ ⁣(AB) dν d3x≤0, and its equality condition applies at each time where the derivative and production integral are defined.

2.1A1F2L1step 1.1algebra∎

By definition Skin=−kBH, and kB>0 by [F2]. Multiplying the identity in step 1.1 by −kB reverses its sign and yields the displayed entropy-production formula with nonnegative right-hand side. Since the multiplier is strictly positive, Skin′(t)=0 holds exactly when the production integral vanishes; the equality characterization in [L1] gives A=B for ν-almost every collision at almost every spatial point. This remains a conditional model deduction and uses no outgoing-pair molecular-chaos assumption.

Physics theoremStatement: AI-adaptedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Entropy production in the local Fourier heat model

Physical scope

Conditional local entropy and energy balance for a sufficiently regular, fixed-composition local-equilibrium continuum field on a rectangular insulated body with a positive-semidefinite symmetric conductivity tensor.

Statement

Let I⊆R be an open time interval, Ω=∏i=13(0,Li) a fixed rectangle, and (T,u,jQ,K) a field record from Rectangular local-equilibrium field data satisfying the adopted local-energy, Fourier, and local-entropy postulates Local energy balance for an insulated continuum body, Fourier heat-flux constitutive model, and Local-equilibrium caloric entropy density. Assume additionally that jQ and its first spatial derivatives are continuous on every compact subset of Ω‾×I. Then, throughout Ω×I, ∂ts(T)+div⁡ ⁣(jQT)=(∇T)TKs(x,T)∇TT2≥0, where Ks=(K+KT)/2. For every compact time interval [t0,t1]⊂I, the model's total entropy S(t)=∫Ωs(T(x,t)) dx is differentiable and nondecreasing, while ∫Ωu(T(x,t)) dx is constant. These are conditional deductions of the declared continuum model, not observations or a claim that every medium obeys Fourier's law.

Facts & Assumptions

Given: A rectangular field record on Ω‾×I with T>0, u=u(T), c(T)=u′(T)>0, jQ=−K(x,T)∇T, Ks⪰0, ∂tu+div⁡jQ=0, and jQ⋅n=0 on every face.

[F1]

Rectangular local-equilibrium field data: T is C1 in time and C2 in space up to the rectangular faces; u and its heat capacity c=u′ are the stated local functions, and ∇ and div⁡ are Euclidean coordinate derivatives.

[A1]

Local energy balance for an insulated continuum body: ∂tu+div⁡jQ=0 and jQ⋅n=0 on each insulated face.

[A2]

Fourier heat-flux constitutive model: jQ=−K∇T and zTKsz≥0 for every z∈R3.

[A3]

Local-equilibrium caloric entropy density: s′(T)=c(T)/T for the declared local-equilibrium entropy density.

[A4]

The flux jQ and its first spatial derivatives are continuous on each compact time slab through the boundary, so the coordinate divergence and face traces in the proof are defined and integrable.

[F2]

Newton–Leibniz needs only continuity on [a,b], differentiability on (a,b), and a Riemann-integrable extension of the interior derivative: a continuously differentiable coordinate function on a closed interval has its derivative integral equal to its endpoint increment.

[F3]

Riemann--Fubini on product rectangles, with lower and upper section integrals and content-zero exceptional sections: iterated Riemann integrals of continuous functions on product rectangles equal the multidimensional integral.

Proof

1.1F1A3algebra

By [A3] and the chain rule, ∂ts(T)=s′(T)∂tT=(c(T)/T)∂tT=(∂tu)/T, since u′(T)=c(T).

1.2A2F1algebra

Substituting jQ=−K∇T from [A2] makes the right side (∇T)TK∇T/T2. For any real vector z, zTKz=zTKsz because the scalar zT(K−KT)z equals its own negative; hence [A2] and T>0 make this quantity nonnegative.

1.3F2F3A1A4

For each coordinate i, apply [F2] to the ith component of jQ/T on [0,Li]: its integral is the endpoint difference [(jQ)i/T]xi=0xi=Li. By [F3] these one-dimensional integrals may be integrated over the other two coordinates and summed. The result is ∫Ωdiv⁡(jQ/T) dx=0 because jQ⋅n=0 on every face and T>0. The same calculation gives ∫Ωdiv⁡jQ dx=0.

2.1A1step 1.1algebra

The coordinate product rule and ∇(1/T)=−∇T/T2 give div⁡(jQ/T)=(div⁡jQ)/T−jQ⋅∇T/T2. Substituting [A1] and step 1.1 yields ∂ts(T)+div⁡(jQ/T)=−jQ⋅∇T/T2.

3.1F1F2A1step 1.2step 1.3step 2.1algebra∎

On Ω‾×[t0,t1], ∂ts(T) and ∂tu(T) are continuous. The one-dimensional fundamental theorem [F2] expresses each time difference quotient as the average of its continuous time derivative; uniform continuity on this compact set makes those quotients converge uniformly in x. Therefore differentiation passes through the bounded-box integrals. Integrating the local balance in step 2.1 and using the zero face integral in step 1.3 gives S′(t)=∫Ω(∇T)TKs∇T/T2 dx. Its integrand is nonnegative by step 1.2. Integrating [A1] and using step 1.3 gives ddt∫Ωu(T) dx=0.

Physics theoremStatement: AI-adaptedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Constant-conductivity heat evolution on an insulated rectangle

Physical scope

Existence and uniqueness for the three-dimensional rigid Fourier heat model with constant positive volumetric heat capacity and scalar conductivity, homogeneous Neumann walls, and even-periodic uniformly regular positive initial data; no general rough-data or arbitrary-domain PDE claim is made.

Statement

Assume Countable Choice. Let Li>0 be lengths measured in metres, set Ω=∏i=13(0,Li), and take the open temperature interval T=(0,∞) with temperatures measured in kelvins. Fix a constant volumetric heat capacity c>0 in J m−3 K−1, a constant scalar conductivity k>0 in W m−1 K−1, and an arbitrary reference energy density u⋆ in J m−3. Define u:T→R by u(θ)=u⋆+cθ, so u′(θ)=c, and define K:Ω‾×T→R3×3 by K(x,θ)=kI3 as in Fourier heat-flux constitutive model. Put a=k/c in m2 s−1.

Let g:R3→T be C2 and satisfy, for every i∈{1,2,3} and x∈R3, g(x1,…,−xi,…,x3)=g(x1,…,xi,…,x3),g(x1,…,xi+2Li,…,x3)=g(x1,…,xi,…,x3). Assume every coordinate derivative Dαg with ∣α∣≤2 is bounded and uniformly continuous, and that 0<Tmin⁡≤g(y)≤Tmax⁡<∞ for all y, with these bounds measured in kelvins.

For x∈Ω‾ and t>0 define T(x,t)=∫R3Γ(x−y,at)g(y) dy,Γ(z,θ)=(4πθ)−3/2e−∣z∣2/(4θ), and set T(x,0)=g(x). Here y is measured in metres, θ=at in the kernel has units m2, and the normalized kernel Γ(⋅,at) has units m−3, so the integral has units of kelvins. Also ∇T=(∂x1T,∂x2T,∂x3T), ΔT:=∑i=13∂xi2T, and on a face with outward unit normal n, ∂nT:=n⋅∇T. For every finite τ>0, this is the unique function T continuous on Ω‾×[0,τ] whose first time derivative for t>0 and spatial derivatives through order two extend continuously to Ω‾×[0,τ], solving c∂tT=kΔT on Ω×(0,τ], with ∂nT=0 on each face and the given initial value. The time derivative at t=0 is interpreted from the right. The field satisfies Tmin⁡≤T(x,t)≤Tmax⁡ for all t≥0, its spatial derivatives through order two are continuous to t=0, and its right time derivative there is aΔg.

On I=(0,∞), set jQ(x,t)=−k∇T(x,t). Then (T,u,jQ,K) is a field record of Rectangular local-equilibrium field data satisfying Local energy balance for an insulated continuum body and Fourier heat-flux constitutive model, so the solution is global in time for this restricted physical continuum model.

Facts & Assumptions

Given: Countable Choice and the geometric, caloric, conductivity, and initial-profile data specified in the Statement.

[A1]

The Axiom of Countable Choice (ACω): for every sequence of nonempty sets there is a function selecting one element from each set. Countable Choice is the only choice principle assumed here.

[A2]

The stated energy-density map is u(θ)=u⋆+cθ on T=(0,∞), so its heat capacity is the constant u′(θ)=c>0 on the full temperature range.

[A3]

Local energy balance for an insulated continuum body: ∂tu+div⁡jQ=0 and jQ⋅n=0 on each insulated face.

[A4]

Fourier heat-flux constitutive model: jQ=−K∇T; the constant scalar case is K=kI3.

[F1]

The heat Cauchy problem for bounded uniformly continuous data: assuming Countable Choice, convolution of a bounded uniformly continuous g:R3→R with the Gaussian heat kernel gives a bounded function, smooth for positive time, solving Uθ=ΔU, and converging locally uniformly to g as θ↓0.

[F2]

Normalisation, parabolic scaling, heat equation and derivative bounds for the heat kernel: Γ(z,θ)>0, has unit integral, and its derivatives have Gaussian tail bounds; Γ(⋅,θ) is an L1 approximate identity as θ↓0.

[F3]

Spatial and time derivatives pass through heat convolution for positive time: for bounded data the heat convolution is smooth for positive time, derivatives pass through its integral, and Uθ=ΔU.

[F4]

L1 approximate identities converge uniformly on compacta for bounded continuous functions: convolution by an L1 approximate identity converges uniformly on every compact set for bounded continuous data.

[F5]

Rectangular local-equilibrium field data: T,u,c,jQ,K have the domains, regularity classes, Euclidean derivatives, and SI units stated in its Definition.

[F8]

Newton–Leibniz needs only continuity on [a,b], differentiability on (a,b), and a Riemann-integrable extension of the interior derivative: a continuous function on a closed interval with an integrable interior derivative has endpoint increment equal to the derivative integral.

[F9]

The heat kernel on Rn and its causal extension: for θ>0, Γ(z,θ)=(4πθ)−3/2exp⁡(−∣z∣2/(4θ)) on R3.

[F10]

Fubini's theorem for L^1 functions on a sigma-finite product: absolutely integrable functions on sigma-finite product measure spaces have equal iterated and product integrals.

[F11]

Ck maps and multi-index derivative notation in Euclidean space: Ck means that every coordinate derivative through order k exists and is continuous on the open domain.

Proof

1.1A1A2A3A4F1givenalgebra

The assumptions on g make it bounded and uniformly continuous, so [F1] applies. Write U(x,θ)=Γθ∗g(x) and set T(x,t)=U(x,at). The chain rule gives ∂tT=aUθ(x,at)=aΔU(x,at)=aΔT for t>0; [F1] also gives the formula, smoothness, and local uniform initial convergence. Since k=ac, this is c∂tT=kΔT. With u′(T)=c, ut=cTt, and jQ=−k∇T from [A3, A4], the same equality gives ∂tu+div⁡jQ=c∂tT−kΔT=0.

1.2F1F2F5F6F7F9A4givenalgebra

Positivity and unit mass in [F2] imply Tmin⁡∫Γ≤T(x,t)≤Tmax⁡∫Γ, hence Tmin⁡≤T≤Tmax⁡. Translation by 2Li in the convolution variable preserves g and Lebesgue measure by [F6], so T is 2Li-periodic in xi. Reflection yi↦−yi preserves Lebesgue measure by [F7] because its matrix has determinant of absolute value one; evenness of the kernel [F9] and of g then makes T even in each coordinate. For t>0, evenness gives ∂xiT=0 on the face xi=0; periodicity gives ∂xiT=∂xiT at opposite faces xi=−Li,Li, while evenness makes these two values negatives, so both vanish. The same boundary derivative holds at t=0 for g, and [A4] gives jQ⋅n=0 on every face.

1.3F2F3F4F8F10F11givenalgebra

By [F11], the stated C2 hypothesis gives continuous derivatives through spatial order two. For ∣α∣≤2 and t>0, differentiating under the integral and then integrating by parts in each coordinate gives DαT(x,t)=∫Γ(x−y,at)Dαg(y) dy. The products of kernel derivatives with bounded derivatives of g are absolutely integrable by the Gaussian bound [F2], so [F10] licenses the iterated integrals; on each finite coordinate interval Newton–Leibniz [F8] gives the integration-by-parts identity, and the boundary terms vanish as the interval expands because of the Gaussian tails. Each Dαg is bounded and continuous, so [F4] makes this expression converge uniformly on Ω‾ to Dαg as t↓0. Thus the spatial derivatives through order two extend continuously to the initial time.

2.1F1F3F4F8step 1.1step 1.3

For t>0, the heat equation and the same integration-by-parts argument give ∂tT(x,t)=a∫Γ(x−y,at)Δg(y) dy. By [F4] this converges uniformly on Ω‾ to aΔg(x). For 0<ε<t, [F8] gives T(x,t)−T(x,ε)=∫εt∂sT(x,s) ds. Letting ε↓0 is valid uniformly in x by step 1.1 and the uniform limit just established; consequently T(x,t)−g(x)=∫0t∂sT(x,s) ds, and dividing by t proves the right derivative ∂tT(x,0)=aΔg(x).

3.1F5F8F10step 1.2step 1.3step 2.1algebra∎

For uniqueness, let T1,T2 be two solutions in the stated regularity class and set w=T1−T2. Then wt=aΔw, ∂nw=0, and w(⋅,0)=0. At each fixed time, w and its spatial derivatives through order two are continuous on the compact rectangle, so the integrands w ∂i2w and ∣∂iw∣2 are bounded and hence L1 for the finite Lebesgue product measure. Apply [F10] successively in the three coordinate factors to identify the rectangle integrals with the iterated one-dimensional integrals. The product rule and [F8] then give ∫0Liw ∂i2w=[w∂iw]0Li−∫0Li∣∂iw∣2. The boundary terms vanish by the homogeneous Neumann condition, so E′(t)=2a∫ΩwΔw=−2a∫Ω∣∇w∣2≤0. Thus E(0)=0 and E(t)≥0 imply E(t)=0; continuity of w gives w=0 everywhere. The formula therefore gives the unique solution in the stated class.

Physics theoremStatement: AI-adaptedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Exponential relaxation of the constant-conductivity Neumann heat model

Physical scope

Conditional L2 relaxation for the compatible positive-temperature solution of the constant-capacity, constant-conductivity model on a three-dimensional insulated rectangle.

Statement

Assume Countable Choice and the hypotheses and solution T of Constant-conductivity heat evolution on an insulated rectangle, let V=∣Ω∣=L1L2L3, define the spatial mean T‾(t)=V−1∫ΩT(x,t) dx, and put E(t)=∫Ω∣T(x,t)−T‾(t)∣2 dx. With a=k/c and CΩ=32max⁡iLi2, the mean is constant and E(t)≤e−2at/CΩE(0)(t≥0). Thus this model solution converges exponentially to its initial spatial mean in L2(Ω); no pointwise-uniform convergence rate is asserted.

Facts & Assumptions

Given: Countable Choice, the smooth positive solution and homogeneous Neumann boundary data from Constant-conductivity heat evolution on an insulated rectangle, the rectangular mean-zero inequality on Ω, and a=k/c>0.

[A1]

The Axiom of Countable Choice (ACω): every sequence of nonempty sets has a choice function.

[F1]

Constant-conductivity heat evolution on an insulated rectangle: Tt=aΔT on Ω, ∂nT=0 on every face, and T is C1 in time and C2 in space up to the boundary on finite intervals.

[F2]

Mean-zero Poincare inequality on a rectangle: if w has zero spatial mean, then ∫Ωw2≤CΩ∫Ω∣∇w∣2 with ∇w=(∂x1w,∂x2w,∂x3w) and CΩ=32max⁡iLi2.

[F3]

Newton–Leibniz needs only continuity on [a,b], differentiability on (a,b), and a Riemann-integrable extension of the interior derivative: a continuous function on a closed interval with an integrable interior derivative has endpoint increment equal to the derivative integral.

[F4]

Riemann--Fubini on product rectangles, with lower and upper section integrals and content-zero exceptional sections: iterated Riemann integrals of continuous functions on product rectangles equal the multidimensional integral.

Proof

1.1A1F1F3F4given

On each finite time interval the time derivative is continuous on the compact box, so its difference quotients converge uniformly under the spatial integral. By [F1] and [F4], coordinatewise application of [F3] gives T‾′(t)=aV−1∫ΩΔT dx=0: the integral of each second coordinate derivative is the difference of the corresponding first derivative on its two faces, and both normal derivatives vanish. Hence T‾(t)=T‾(0).

2.1F1F3F4step 1.1algebra

Set w(x,t)=T(x,t)−T‾(0). Then [F4] and step 1.1 give ∫Ωw=0, while wt=aΔw and ∂nw=0. Differentiating E(t)=∫Ωw2 and using [F4], the product rule, and [F3] in each coordinate gives E′(t)=2a∫ΩwΔw=−2a∫Ω∣∇w∣2; each endpoint term is w∂nw=0.

3.1F2step 2.1

Apply [F2] to w(⋅,t) in step 2.1. It gives E(t)≤CΩ∫Ω∣∇w∣2, so E′(t)≤−(2a/CΩ)E(t).

4.1F3step 3.1algebra∎

The product rule gives (e2at/CΩE(t))′=e2at/CΩ(E′(t)+(2a/CΩ)E(t))≤0 by step 3.1. Integrating this derivative with [F3] from 0 to t yields e2at/CΩE(t)≤E(0), which is the asserted estimate.

RemarkRemark: AI-adaptedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Scope of the continuum heat solution results

Physical scope

Source-backed scope and limitation note for the local rectangular Fourier heat model; it reports no experiment and asserts no universal continuum closure.

Remarks

The existence and relaxation results on this page apply to a fixed rectangle, constant positive c and k, homogeneous Neumann boundary data, and the specified even-periodic C2 initial profiles. They do not prove existence for rough initial data, arbitrary domains or boundary conditions, state-dependent conductivity, or a microscopic Boltzmann-to-Fourier limit. For the constant-coefficient heat equation of Constant-conductivity heat evolution on an insulated rectangle, the Gaussian-kernel solution gives an immediate nonzero response throughout the rectangle for the compactly supported perturbation in The rectangular Fourier model has instantaneous spatial support ↗. That is a mathematical property of this adopted diffusion model; it is not a measurement or a claim about finite-speed physical signals.

5 · Examples, experiments and counterexamples

None yet.

Sources