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

Smooth-branch velocity error for backward Euler MINI

Statement

Assume Countable Choice. Assume the MINI mesh and time-partition hypotheses of The mixed mini-element backward Euler Navier–Stokes scheme. Suppose (u,p) is a strong solution of the no-slip Navier–Stokes equations with forcing f∈L2(0,T;L2(D)d) and u∈W1,∞(0,T;H2(D)d∩Vσ) and p∈W1,∞(0,T;H1(D)∩L02(D)). Set M:=1+∥u∥W1,∞H2+∥p∥W1,∞H1,Mu:=∥u∥L∞V, M∞:=(∥u0∥22+CP(D)2ν−1∥f∥L2L22)1/2,ME:=ν−1/2∥u0∥2+CP(D)ν−1∥f∥L2L2, and Λ:=1+CDν−3Mu4, where CD depends only on the domain, dimension and the trilinear-form constant. For sufficiently fine time partitions satisfying τmax⁡Λ≤12, initialize by uh0=PZhu0 and choose any backward Euler MINI solution. Then max⁡0≤n≤N∥u(tn)−uhn∥2+ν1/2(∑n=1Nτn∥u(tn)−uhn∥V2)1/2≤C∗(hM+τmax⁡M2), where one may take C∗:=CD,shape(1+ν)2(1+ν−4)(1+T)2(1+M∞+ME+M)4exp⁡(2TΛ). Here f is the forcing of the strong solution and CD,shape is independent of h, all time steps, and ν. The polynomial factor displays the positive and inverse viscosity powers produced by Young's inequalities and static Stokes projection. The explicit factor ν−3 in Λ and exponential dependence on TΛ are retained; no Reynolds-uniform estimate is asserted. The separate instantaneous pressure assertion is addressed in Remarks.

Facts & Assumptions

Given: Countable Choice, a bounded connected Lipschitz polyhedral domain, a shape-regular MINI mesh family, a partition of [0,T], viscosity ν>0, and the strong solution and initialization in the Statement.

[F1]

The exact solution satisfies the time-integrated mixed weak equation; for tests zh∈Zh its pressure contribution is b(zh,pˉn). (The mixed mini-element backward Euler Navier–Stokes scheme, Energy weak solutions of incompressible Navier–Stokes)

[F2]

The linear solenoidal energy projection Sh satisfies ∥u−Shu∥V≤C0(1+ν−1)h(∥u∥H2+∥p∥H1) and the same estimate applies to (ut,pt). (Approximation by the MINI solenoidal energy projection)

[F3]

Mean-adjusted face-average interpolation maps H1(D) to Qh with ∥p−rh∥2≤Ch∥p∥H1. (Boundary-preserving face-average P1 interpolation, The conforming P1-bubble/P1 mini element)

[F4]

The skew form is trilinear and continuous, obeys ∣c(a;b,w)∣≤CD∥a∥21/2∥a∥V1/2∥b∥V∥w∥V, and c(a;w,w)=0. (Continuity, energy cancellation, and Sobolev bounds for skew convection, The skew convective form for incompressible flow)

[F5]

The discrete solution obeys max⁡n∥uhn∥2≤M∞ and ∑nτn∥uhn∥V2≤ME2. (Backward Euler MINI step existence and discrete energy bound)

[F6]

The L2 projector PZh onto the weak MINI kernel defined in Weak and pointwise divergence-free discrete velocities exists and is a best approximation in L2. (Discrete Helmholtz projection and pressure Schur identities)

[F7]

On the Stokes velocity space V=H01(D)d, Poincare, Korn, Holder and Young give ∥z∥2≤CP∥z∥V and aν(z,z)≥ν∥z∥V2. (Stokes velocity and pressure spaces and weak forms, The Poincare inequality for zero-boundary Sobolev closures on domains bounded in one direction, The homogeneous Dirichlet strain identity, Holder's inequality for integrals, including the endpoint cases, Young's inequality for conjugate real exponents)

[A1]

Countable Choice is assumed through the finite-dimensional projection and scheme-existence inputs. (The Axiom of Countable Choice (ACω), Approximation by the MINI solenoidal energy projection, Backward Euler MINI step existence and discrete energy bound)

[F8]

The static mixed projection contributes at most 1+ν−1, pressure and trilinear Young estimates contribute inverse powers no worse than ν−3, and the viscous consistency term contributes one factor ν before it is balanced by coercivity. Thus all residual constants are bounded by CD,shape(1+ν+ν−3); all such constants depend otherwise only on the domain, mesh-shape bound and dimension. (Mixed Galerkin error bound retaining beta_h, Approximation by the MINI solenoidal energy projection, Young's inequality for conjugate real exponents)

Proof

1.1A1F2F6F7given

Let Sh be the projection of [F2] and set η(t)=u(t)−Shu(t) and Mν=(1+ν−1)M. Linearity and time independence of Sh give ηt=ut−Shut. Applying [F2] to (u(t),p(t)) and (ut(t),pt(t)) yields ∥η∥L∞V+∥ηt∥L∞V+∥η∥L∞L2≤ChMν, where the last term follows from Poincare [F7]. Define θn=Shu(tn)−uhn∈Zh. Since PZh is a best L2 approximation and Shu0∈Zh, [F2], [F6] and Poincare give ∥θ0∥2≤2∥u0−Shu0∥2≤ChMν.

2.1F1F2step 1.1algebra

Write un=u(tn) and uˉn=τn−1∫tn−1tnu(t) dt, with analogous notation pˉn and cˉn=τn−1∫c(u(t);u(t),⋅) dt. Average the exact equation in time and subtract the discrete equation. For every zh∈Zh this gives (θn−θn−1τn,zh)+aν(θn,zh)=−(η(tn)−η(tn−1)τn,zh)−aν(η(tn)+uˉn−un,zh)−b(zh,pˉn)+c(uhn;uhn,zh)−cˉn(zh). The identity uses aν(Shun−un,zh)=0 from the defining projection; all forcing averages cancel exactly.

3.1F2F3F7step 2.1

Test step 2.1 with zh=θn and multiply by τn. For the pressure choose the mean-adjusted interpolant rh(t) from [F3] and set rˉhn=τn−1∫tn−1tnrh(t) dt∈Qh; then b(θn,rˉhn)=0. Linearity and the triangle inequality for the Bochner average give ∥pˉn−rˉhn∥2≤ChM, and ∣div⁡θn∣≤d∥θn∥V, so ∣b(θn,pˉn)∣≤CDhM∥θn∥V. Also ∥uˉn−un∥V≤τnM. The time derivative estimate [F2] gives ∥η(tn)−η(tn−1)∥2≤CτnhMν. Cauchy–Schwarz and Young, with viscous fractions ντn/16 for each of the viscous and pressure terms, yield ∣(η(tn)−η(tn−1),θn)∣+τn∣aν(η(tn)+uˉn−un,θn)∣+τn∣b(θn,pˉn−rˉhn)∣≤τn8∥θn∥22+ντn8∥θn∥V2+CD,ντn(h2Mν2+τn2M2). Here CD,ν is bounded by CD,shape(1+ν+ν−1) for these linear residuals.

3.2F4F5F7step 1.1step 2.1

The identity c(uhn;uhn,θn)−c(un;un,θn)=−c(uhn;η(tn),θn)−c(η(tn);un,θn)−c(θn;un,θn) follows by bilinearity in the first two slots and skew cancellation c(uhn;θn,θn)=0 from [F4]. Apply [F4], Poincare and Young after multiplication by τn. Each of these three terms is bounded by ντn∥θn∥V2/64 plus respectively CD,ντnM∞∥uhn∥Vh2Mν2,CD,ντn∥un∥V2h2Mν2,CDν−3τnMu4∥θn∥22. For the time-averaging remainder, [F4] and ∥u(t)−un∥V≤τnM give τn∣c(un;un,θn)−cˉn(θn)∣≤CDτn2M2∥θn∥V≤ντn∥θn∥V2/64+CD,ντn3M4.

4.1F4F5F8step 2.1step 3.1step 3.2algebra

The Hilbert identity 2(x−y,x)=∥x∥2−∥y∥2+∥x−y∥2 applied to the time term in step 2.1 converts the tested equation to an energy recurrence. The viscous Young fractions in steps 3.1–3.2 sum to 3/16 of ντn∥θn∥V2 before doubling, so at least 5/8 remains and in particular the displayed half remains. Doubling also changes the time residual in step 3.1 into a coefficient at most τn/4 of En; the convection coefficient is at most CDν−3τnMu4. Thus, enlarging CD if necessary, En−En−1+ν2τn∥θn∥V2≤τnΛEn+CD,ντn{h2Mν2(1+M∞∥uhn∥V+Mu2)+τn2M4},En:=∥θn∥22, with Λ=1+CDν−3Mu4. The residual coefficient satisfies CD,ν≤CD,shape(1+ν+ν−3) by [F8].

5.1F5F7F8step 1.1step 4.1algebra

Since τnΛ≤1/2, rearrange step 4.1 and use (1−x)−1≤e2x for 0≤x≤1/2. Iterating gives max⁡nEn≤e2TΛ(E0+CD,ν∑nτn{h2Mν2(1+M∞∥uhn∥V+Mu2)+τn2M4}). By Cauchy–Schwarz and [F5], ∑nτn∥uhn∥V≤TME and ∑nτn3≤Tτmax⁡2. Summing step 4.1 and dropping the nonnegative final EN also gives ν2∑nτn∥θn∥V2≤E0+TΛmax⁡nEn+CD,ν∑nτn{h2Mν2(1+M∞∥uhn∥V+Mu2)+τn2M4}. Since Mu≤M, these bounds and E0≤Ch2Mν2 imply max⁡n∥θn∥2+ν1/2(∑nτn∥θn∥V2)1/2≤CD,shape(1+ν)2(1+ν−4)(1+T)2(1+M∞+ME+M)4e2TΛ(hM+τmax⁡M2). Here (1+ν)2(1+ν−4) dominates 1+ν−1 from Mν and the residual inverse powers through ν−3 from [F8]; the extra Λ in the energy norm combines with those factors to give no power worse than ν−4.

6.1F2step 5.1algebra∎

Finally un−uhn=η(tn)+θn. The projection estimate [F2] bounds max⁡n∥η(tn)∥2 by ChMν and ∑nτn∥η(tn)∥V2 by CTh2Mν2. Combine with step 5.1; the explicit viscosity factor in C∗ absorbs Mν/M=1+ν−1. This proves the velocity estimate. Its ν dependence is explicit and contains an exponential in ν−3Mu4T; it makes no Reynolds-uniform claim.

Remarks

The positive velocity estimate above retains the original mesh, data, mass initialization, arbitrary-step and any-root scope. The originally paired unrestricted instantaneous Lt4/3Lx2 pressure-error rate is false: the contextual Graded MINI mass initialization can destroy the smooth instantaneous pressure rate ↗ constructs an exact-scope graded-mesh counterexample with smooth data and genuine backward Euler roots. That counterexample is not a logical premise of the velocity proof. Pressure-primitive estimates remain compatible with it.

Depends on

Used by

Dependency tree · two levels

72 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