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

First-order time error under smooth fixed-mesh solutions

Statement

Fix one finite-dimensional realization in No-flux finite-volume and conforming realizations for phase fields and T>0. Let Gc,h(u)=κAhu+∇hPc,h(u) and Ge,h(u)=∇hPe,h(u). The exact semidiscrete Allen–Cahn trajectory is the solution of y′=−M(Gc,h(y)−Ge,h(y)); the exact semidiscrete Cahn–Hilliard trajectory is the fixed-mean solution of y′=−MAh(Gc,h(y)−Ge,h(y)). Assume the applicable exact trajectory is C2 on [0,T] and remains in a compact set. Let the numerical solution start from u0=y(0) and use the convex-splitting step. Then there are constants C>0 and τ0>0, depending on the fixed spatial realization, ε,κ,M, the exact trajectory, and T, such that, for 0<τ≤τ0:

  • for Allen–Cahn, max⁡nτ≤T∥un−y(tn)∥h≤Cτ;
  • for Cahn–Hilliard, with y and un having the same mean and ∥z∥−1,h2=⟨Ah†z,z⟩h on mean-zero vectors, max⁡nτ≤T∥un−y(tn)∥−1,h≤Cτ.

This is a fixed-mesh time estimate; it is not a uniform spatial error bound.

Facts & Assumptions

Given: A fixed finite-dimensional realization, the convex energy split, the exact semidiscrete solution, its stated smoothness and compact range, the exact initial value, and a step size 0<τ≤τ0.

[F1]

Backward-Euler convex splitting for Allen–Cahn and Cahn–Hilliard: the AC update is δu/τ=−M(Gc,h(un+1)−Ge,h(un)), and the CH update is δu/τ=−MAh(Gc,h(un+1)−Ge,h(un)).

[F2]

Existence, uniqueness, mass conservation, and energy decay of a convex-splitting step: each step exists uniquely, the physical discrete energy is nonincreasing, and the convex implicit term has the monotonicity used in its proof.

[F3]

Kernel, range, and mean-zero inverse for the no-flux Laplacian: Ah is self-adjoint, nonnegative, and invertible on mean-zero vectors, so ∥⋅∥−1,h is a norm there and Ah† is available.

[F4]

Allen–Cahn and Cahn–Hilliard double-well gradient flows and No-flux finite-volume and conforming realizations for phase fields: the scalar split is Fc(s)=s4/(4ε2)+s2/(2ε2), Fe(s)=s2/ε2−1/(4ε2), so ∇hPe,h(u)=2u/ε2 is globally Lipschitz and Pc,h is convex.

[A1]

The relevant exact semidiscrete trajectory is C2 and its image lies in a compact set on which Gc,h has bounded derivative.

Proof

1.1A1F1F3givenalgebra

Substitute the exact trajectory into the backward-Euler split. For Allen–Cahn define RnAC:=y(tn+1)−y(tn)+Mτ(Gc,h(y(tn+1))−Ge,h(y(tn))). For Cahn–Hilliard define RnCH:=y(tn+1)−y(tn)+MτAh(Gc,h(y(tn+1))−Ge,h(y(tn))). The exact integral equations express these residuals as M times the quadrature difference τGc,h(y(tn+1))−∫tntn+1Gc,h(y(s))ds+∫tntn+1Ge,h(y(s))ds−τGe,h(y(tn)), with an additional Ah in the Cahn–Hilliard branch. The bounded derivatives on the compact exact trajectory and ∥y(t)−y(s)∥≤C∣t−s∣ imply ∥RnAC∥h≤Cτ2 and ∥RnCH∥−1,h≤Cτ2; in the latter case RnCH lies in the mean-zero range of Ah by [F3].

2.1F2F4step 1.1algebra

Convexity of Pc,h and nonnegativity of Ah make Gc,h monotone. By [F4], Ge,h is globally Lipschitz with constant Le=2/ε2. In the quartic realizations, the lower bound Fε(s)≥s4/(8ε2)−1/(4ε2) and the energy inequality in [F2] also bound every numerical state on the fixed finite-dimensional space.

3.1F1F4step 1.1step 2.1algebra

In the Allen–Cahn branch put en=un−y(tn). Subtract the exact residual equation from the numerical recurrence and pair with en+1. The identity 2⟨en+1−en,en+1⟩h=∥en+1∥h2−∥en∥h2+∥en+1−en∥h2 and monotonicity from step 2.1 leave only the explicit term and residual on the right. Cauchy–Schwarz and Young give MτLe∥en∥h∥en+1∥h≤MτLe2(∥en∥h2+∥en+1∥h2) and ∥RnAC∥h∥en+1∥h≤τ8∥en+1∥h2+2τ−1∥RnAC∥h2. For τ small enough that MτLe≤1/8 and τ≤1, step 1.1 then gives ∥en+1∥h2≤(1+Cτ)∥en∥h2+Cτ3.

3.2F1F3F4step 1.1step 2.1algebra

In the Cahn–Hilliard branch both errors have mean zero. Pair the original difference equation with Ah†en+1. Self-adjointness and AhAh†en+1=en+1 give 2⟨en+1−en,Ah†en+1⟩h=∥en+1∥−1,h2−∥en∥−1,h2+∥en+1−en∥−1,h2. The implicit term becomes Mτ⟨Gc,h(un+1)−Gc,h(yn+1),en+1⟩h≥Mτκ∥en+1∥1,h2, because Gc,h is the sum of κAh and the gradient of a convex functional. The explicit term obeys ∣⟨Ge,h(un)−Ge,h(yn),en+1⟩h∣≤Le∥en∥−1,h∥en+1∥1,h, where ∥v∥1,h2=⟨Ahv,v⟩h. Indeed, on the mean-zero space (x,y)Ah:=⟨Ahx,y⟩h is an inner product and ⟨en,en+1⟩h=(Ah†en,en+1)Ah, with ∥Ah†en∥Ah=∥en∥−1,h and ∥en+1∥Ah=∥en+1∥1,h. Cauchy–Schwarz in this inner product gives the bound. Young's inequality absorbs half of the positive Mτκ∥en+1∥1,h2 term and leaves Cτ∥en∥−1,h2. The residual from step 1.1 satisfies ∣⟨RnCH,Ah†en+1⟩h∣=∣⟨RnCH,en+1⟩−1,h∣≤Cτ2∥en+1∥−1,h≤τ8∥en+1∥−1,h2+C′τ3. Hence, for small τ, ∥en+1∥−1,h2≤(1+Cτ)∥en∥−1,h2+Cτ3.

4.1step 3.1step 3.2algebra∎

In either recurrence, e0=0. Iterating zn+1≤(1+Cτ)zn+Cτ3 for nτ≤T gives zn≤Cτ3∑j=0n−1(1+Cτ)j≤CTτ2, since (1+Cτ)n≤eCT. Taking square roots gives the two asserted O(τ) bounds. The constants include the fixed-mesh operator norms, the fixed parameters M,κ,ε, and the trajectory bounds.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

13 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