Documentation

LeanPool.EllipticPDE.Regularity.Interior.EnergyBound

Master interior difference-quotient energy estimate #

The interior second-derivative estimate (Evans, Partial Differential Equations (2nd ed.), §6.3.1; Gilbarg-Trudinger, Elliptic PDE of Second Order, Theorem 8.8) proceeds by testing the weak formulation of L u = f with the difference-quotient test element v_h = -Dₖ^{-h}(ξ² Dₖ^h u), using discrete integration by parts to move the outer difference quotient onto the coefficient factor, uniform ellipticity from below to control the leading term, and Cauchy-Schwarz together with the Peter-Paul (Young) inequality to absorb the commutator, cross, zeroth-order, and right-hand terms.

Main declarations #

Admissible Evans test element #

noncomputable def EllipticPdes.Regularity.evansTest {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} {ξ θ : EuclideanSpace ℝ (Fin d) → ℝ} (hΩm : MeasurableSet Ω) (hξ : Sobolev.IsTestFn Ω ξ) (hθ : Sobolev.IsTestFn Ω θ) (k : Fin d) (h : ℝ) (hsm_in : ∀ x ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y, x + hshift k h ∈ Ω) (hsm_out : ∀ x ∈ tsupport θ, x + hshift k (-h) ∈ Ω) (u : ↥(Sobolev.H01 Ω)) :
↥(Sobolev.H01 Ω)

Admissible Evans test element v_h = -Dₖ^{-h}(ξ² Dₖ^h u) ∈ H₀¹(Ω). The inner cutoff ξ² and the outer cutoff θ (which is ≡ 1 on tsupport ξ) localise the two difference quotients so that the composite stays inside H₀¹(Ω); membership is two applications of the crux admissibility lemma cutoffMul_diffQuotG_mem_H01, together with closure of the submodule under negation. This is the single admissible test vector that the weak formulation consumes in the difference-quotient energy method (Evans, Partial Differential Equations (2nd ed.), §6.3.1).

Equations
  • One or more equations did not get rendered due to their size.
Instances For
    theorem EllipticPdes.Regularity.evansTest_coe {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} {ξ θ : EuclideanSpace ℝ (Fin d) → ℝ} (hΩm : MeasurableSet Ω) (hξ : Sobolev.IsTestFn Ω ξ) (hθ : Sobolev.IsTestFn Ω θ) (k : Fin d) (h : ℝ) (hsm_in : ∀ x ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y, x + hshift k h ∈ Ω) (hsm_out : ∀ x ∈ tsupport θ, x + hshift k (-h) ∈ Ω) (u : ↥(Sobolev.H01 Ω)) :
    ↑(evansTest hΩm hξ hθ k h hsm_in hsm_out u) = -(cutoffMul hθ) ((diffQuotG k (-h) hΩm) ((cutoffMul ⋯) ((diffQuotG k h hΩm) ↑u)))

    The ambient-graph value of evansTest is the negated cutoff of the outer difference quotient of ξ² Dₖ^h u.

    Support control and θ-invisibility #

    Discrete integration by parts #

    Support of the inner cutoff data and the Evans coordinate reduction #

    Extension-by-zero weak derivative and the first-order global energy #

    theorem EllipticPdes.Regularity.hasWeakDeriv_extendL2_of_mem_H01 {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} (hΩm : MeasurableSet Ω) (k : Fin d) {U : Sobolev.H1amb Ω} (hU : U ∈ Sobolev.H01 Ω) :
    HasWeakDeriv k ((extendL2 hΩm) (U.ofLp 0)) ((extendL2 hΩm) (U.ofLp k.succ))

    Extension by zero of an H₀¹ element preserves the weak gradient. For u ∈ H₀¹(Ω), the whole-space extension by zero of the function value u₀ has whole-space L² weak k-derivative equal to the extension by zero of the gradient component u_{k+1}. Because u vanishes at the boundary (it lies in the closure of the compactly supported test functions), no boundary term appears when integrating against an arbitrary whole-space test function φ: the identity ∫ (extendL2 u₀) ∂ₖφ = -∫ (extendL2 u_{k+1}) φ is closed under L² limits and holds on every test-function graph by classical integration by parts, hence on all of H₀¹(Ω) (Evans, Partial Differential Equations (2nd ed.), §5.8.2).

    theorem EllipticPdes.Regularity.firstOrder_energy_le {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} (Op : Sobolev.FullEllipticOp d) (u : ↥(Sobolev.H01 Ω)) (f : Sobolev.L2D Ω) (hu : ∀ (v : ↥(Sobolev.H01 Ω)), ((Op.fullBilin Ω) u) v = ∫ (x : EuclideanSpace ℝ (Fin d)) in Ω, ↑↑f x * ↑↑((↑v).ofLp 0) x) :
    Op.lam / 2 * ∑ i : Fin d, ‖(↑u).ofLp i.succ‖ ^ 2 ≤ ‖f‖ * ‖(↑u).ofLp 0‖ + Op.gardingγ * ‖(↑u).ofLp 0‖ ^ 2

    First-order global energy estimate. For a weak solution u ∈ H₀¹(Ω) of L u = f whose transport field b vanishes and whose zeroth-order coefficient c is nonnegative (a.e. on Ω), the full gradient energy is bounded by the data: λ ∑ᵢ ‖u_{i+1}‖² ≤ ‖f‖ · ‖u₀‖. Testing the weak formulation with u itself, ellipticity bounds the principal part from below, the transport term drops (b = 0) and the zeroth-order term has a sign (c ≥ 0), so only the right-hand pairing ⟪f, u₀⟫ survives (Evans, Partial Differential Equations (2nd ed.), §6.2.2).

    theorem EllipticPdes.Regularity.norm_diffQuotD_u0_le_ambient {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} (hΩm : MeasurableSet Ω) (k : Fin d) (h : ℝ) (U : Sobolev.H1amb Ω) (hW : HasWeakDeriv k ((extendL2 hΩm) (U.ofLp 0)) ((extendL2 hΩm) (U.ofLp k.succ))) :
    ‖(diffQuotD k h hΩm) (U.ofLp 0)‖ ≤ ‖U.ofLp k.succ‖

    Ambient-space generalisation of norm_diffQuotD_u0_le (proof of concept). H01 is not needed for this inequality itself: it is needed only to manufacture the one whole-space weak-derivative fact hasWeakDeriv_extendL2_of_mem_H01 supplies. Factoring that fact out as a hypothesis hW shows the rest of norm_diffQuotD_u0_le's proof (the identity diffQuotD = restrict ∘ diffQuot ∘ extendL2, the non-expansive restriction, and the difference-quotient/weak-derivative bound norm_diffQuot_le_of_hasWeakDeriv) survives unchanged for an arbitrary ambient graph U : H1amb Ω, with no reference to H01 at all. Evans, Partial Differential Equations (2nd ed.), §6.3.1, Remark (i), states that the zero-trace hypothesis this lemma's H01-typed sibling states is not required for the interior estimate; this declaration isolates that the only place it was doing work in this particular step was in supplying hW, not in the inequality.

    Kept as that record. The chain runs on the H01-typed sibling, so nothing consumes this one.

    Master assembly toolkit #

    Master interior difference-quotient energy estimate #

    theorem EllipticPdes.Regularity.firstOrder_gradNorm_le {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} (Op : Sobolev.FullEllipticOp d) (u : ↥(Sobolev.H01 Ω)) (f : Sobolev.L2D Ω) (hu : ∀ (w : ↥(Sobolev.H01 Ω)), ((Op.fullBilin Ω) u) w = ∫ (x : EuclideanSpace ℝ (Fin d)) in Ω, ↑↑f x * ↑↑((↑w).ofLp 0) x) (i : Fin d) :
    ‖(↑u).ofLp i.succ‖ ≤ √((1 + 4 * Op.gardingγ) / (2 * Op.lam)) * (‖f‖ + ‖(↑u).ofLp 0‖)

    First-order gradient bound. Each gradient component of a weak solution is bounded in L² by the data: ‖∂ᵢu‖ ≤ √((1 + 4γ) / (2λ)) (‖f‖ + ‖u₀‖), where γ is the Gårding shift constant, through which the transport and zeroth-order coefficients enter. This is the first-order energy estimate firstOrder_energy_le combined with the arithmetic-geometric mean inequality.

    theorem EllipticPdes.Regularity.interior_diffQuot_energy_bound {d : ℕ} {Ω : Set (EuclideanSpace ℝ (Fin d))} {ξ θ : EuclideanSpace ℝ (Fin d) → ℝ} (Op : Sobolev.FullEllipticOp d) (hΩm : MeasurableSet Ω) (hA : IsLipCoeff Op.toEllipticCoeff) (hξ : Sobolev.IsTestFn Ω ξ) (hθ : Sobolev.IsTestFn Ω θ) (k : Fin d) :
    ∃ (C : ℝ), 0 ≤ C ∧ ∀ (u : ↥(Sobolev.H01 Ω)) (f : Sobolev.L2D Ω), (∀ (w : ↥(Sobolev.H01 Ω)), ((Op.fullBilin Ω) u) w = ∫ (x : EuclideanSpace ℝ (Fin d)) in Ω, ↑↑f x * ↑↑((↑w).ofLp 0) x) → ∀ (h : ℝ), h ≠ 0 → (∀ x ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y, x + hshift k h ∈ Ω) → (∀ x ∈ tsupport θ, x + hshift k (-h) ∈ Ω) → (∀ x ∈ Ω, ((x ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y) ∨ x + hshift k (-h) ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y) → θ x = 1) → (∀ (j : Fin d), ∀ x ∈ Ω, ((x ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y) ∨ x + hshift k (-h) ∈ tsupport fun (y : EuclideanSpace ℝ (Fin d)) => ξ y * ξ y) → Sobolev.partialD j θ x = 0) → Op.lam / 2 * ∑ i : Fin d, ‖(extendL2 hΩm) ((mulTest hξ) ((diffQuotD k h hΩm) ((↑u).ofLp i.succ)))‖ ^ 2 ≤ C * (‖f‖ ^ 2 + ‖(↑u).ofLp 0‖ ^ 2)

    Master interior difference-quotient energy estimate. For a W^{1,∞}-coefficient weak solution u ∈ H₀¹(Ω) of L u = f, an inner cutoff ξ and an outer cutoff θ ≡ 1 on the shift-reachable part of tsupport ξ², the cutoff-weighted energy of the interior difference quotient of the gradient is bounded by the data, uniformly in the step h: (λ/2) ∑ᵢ ‖ξ · Dₖ^h ∂ᵢu‖² ≤ C (‖f‖² + ‖u₀‖²). The constant is quantified before the solution and the datum, so it depends only on λ, Λ, A₁, d, γ, ‖b‖∞, ‖c‖∞, ‖ξ‖∞, ‖∂ξ‖∞, and on none of u, f, h. Testing the weak formulation with the admissible Evans element v_h = -Dₖ^{-h}(ξ² Dₖ^h u), discrete integration by parts (evansTest_bilin_L2D) moves the outer difference quotient onto the coefficient action; the discrete Leibniz split (norm_diffQuotD_actL_sub_le) exposes the translated-coefficient leading term, controlled from below by ellipticity (energy_ge for the translate, same λ); Cauchy-Schwarz and the Peter-Paul inequality absorb the commutator, cross, transport, zeroth-order and right-hand terms, five families each spending an eighth of the ellipticity lower bound, with the first-order energy bound firstOrder_energy_le supplying all gradient data (Evans, Partial Differential Equations (2nd ed.), §6.3.1; Gilbarg-Trudinger, Elliptic PDE of Second Order, Theorem 8.8).