Alphabeta Math
TheoremStatement: AI-adaptedProof: AI-adaptedPipeline-generatedprecheck passjudge pass (gpt-6-luna)
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.

TPFA flux-energy convergence for smooth face-aligned interfaces

Statement

Let Ω⊂Rd be a bounded connected polyhedral set and let (Th) be a refining, uniformly regular, boundary-anchored family of TPFA-admissible meshes. Assume λK∈[λ∗,λ∗] with 0<λ∗≤λ∗<∞, constant on each cell, and that coefficient jumps occur only across mesh faces. Suppose u is continuous, vanishes on ∂Ω, and is C2 on each cell closure with the uniform family bound M=sup⁡hmax⁡K∈Kh∥D2u∥L∞(K)<∞. For every interior face F=K∣L fix an orientation and define its exact outward flux from K by JK,F=−∫FλK∇u∣K⋅nK,F. On a boundary face define JK,F by the same formula. Assume the exact fluxes satisfy the cell balances ∑F⊂∂KJK,F=∫Kf and JK,F=−JL,F on every interior face. Let Uh solve the TPFA cell balances with source f and zero boundary values, and let Ihu be the cell-center sample u(xK) with zero boundary face values. Then, for a constant C depending only on Ω, λ∗,λ∗, and the uniform mesh regularity bounds, ∥Uh−Ihu∥Eh≤ChM,(∑F∈Fint∣QK,F(Uh)−JK,F∣2TF+∑F∈Fbdy∣QK,F(Uh)−JK,F∣2TK,F)1/2≤2ChM. The second estimate is convergence of the reconstructed numerical face flux in the transmissibility-dual energy norm. No pointwise operator consistency on arbitrary nonorthogonal meshes or error rate for nonsmooth weak solutions is asserted.

Facts & Assumptions

Given: The exact cell and face geometry, a refining uniformly regular orthogonal mesh family, λ∗≤λK≤λ∗ with mesh-aligned jumps, the piecewise C2 function u, the exact face fluxes and their cell balances, and the TPFA scheme.

[F1]

Finite-volume cell-face meshes and geometric regularity: in a uniformly regular family, cell volumes and face measures scale as hd and hd−1, center-to-face distances scale as h, and the number of faces per cell is uniformly bounded.

[F2]

Orthogonal TPFA diffusion with harmonic face transmissibility: TF=∣F∣/(dK,F/λK+dL,F/λL) on interior faces, TK,F=∣F∣λK/dK,F on boundary faces, and ∥v∥Eh2=∑FTF(ΔFv)2.

[F3]

TPFA discrete Green identity, coercivity, and maximum principle: the TPFA stiffness matrix is symmetric positive definite, and its bilinear form is the face energy in [F2].

[A1]

Each cell restriction of u extends C2 to its closure with Hessian norm at most the family-wide bound M. The exact integrated fluxes are opposite on each interior face, and the exact cell flux balances hold as stated.

Proof

1.1F1F2A1givenalgebra

For each face fix an incident cell K and set δF=QK,F(Ihu)−JK,F. On an interior face F=K∣L, write n=nK,F, let yF be the center-normal crossing point, and put sF=−JK,F/∣F∣. Integrated conservation says that sF is the face mean of both λK∂nuK and λL∂nuL. On each cell closure the Hessian bound gives ∣λK∂nuK(yF)−sF∣≤Cλ∗Mh, and the same bound holds for L. Continuity of u and one-sided normal Taylor expansion therefore give u(xK)=u(yF)−dK,FsF/λK+RK and u(xL)=u(yF)+dL,FsF/λL+RL, with ∣RK∣+∣RL∣≤CMh2. Substitution into the harmonic formula yields QK,F(Ihu)=−∣F∣sF+TF(RK−RL), while JK,F=−∣F∣sF. Since TF≤C∣F∣/h, ∣δF∣≤CM∣F∣h. On a boundary face take sF=−JK,F/∣F∣ and yF=xF. The same face-mean estimate and u(yF)=0 give u(xK)=−dK,FsF/λK+RK, ∣RK∣≤CMh2, and the identical defect bound.

2.1F1F2step 1.1algebra

By uniform regularity, every transmissibility satisfies TF≥c∣F∣/h. There are at most N0 faces per cell, ∣F∣≤Chd−1, and ∣K∣≥chd, so #Kh≤∣Ω∣/(chd) and ∑F∣F∣≤Ch−1. Step 1.1 therefore gives ∑FδF2TF≤CM2h3∑F∣F∣≤CM2h2. Thus the vector of local flux defects has transmissibility-dual norm at most ChM.

3.1A1F3step 2.1algebra

Let e=Uh−Ihu. The exact cell balances in [A1] and the finite-volume balances for Uh imply, for every cell vector v, ah(e,v)=−∑F∈FintδF(vK−vL)−∑F∈FbdyδFvK, because the two exact interior fluxes are oppositely oriented. Take v=e. Weighted Cauchy–Schwarz gives ∥e∥Eh2≤(∑FδF2/TF)1/2∥e∥Eh. If e=0 the first estimate is immediate; otherwise divide by ∥e∥Eh and use step 2.1 to obtain ∥e∥Eh≤ChM. The bilinear identity used here is [F3].

4.1F2step 2.1step 3.1algebra∎

On every face, QK,F(Uh)−QK,F(Ihu) equals TF(eK−eL) for an interior face and TK,FeK for a boundary face. Its transmissibility-dual norm is therefore exactly ∥e∥Eh. By the triangle inequality, the numerical-flux error norm is at most ∥e∥Eh+(∑FδF2/TF)1/2≤2ChM, using steps 2.1 and 3.1. This proves both estimates.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · one level

3 results within one dependency step 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