Alphabeta Math
Physics theoremStatement: AI-adaptedProof: AI-adaptedPipeline-generatedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)
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.

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.

Depends on

Dependency tree · 2 levels

Sources