Documentation

TauCeti.Analysis.PDE.Regularity.HarmonicReplacement

Harmonic replacement #

Let u ∈ H¹(Ω). Its harmonic replacement on Ω is the weak solution w of the Dirichlet problem Δw = 0 in Ω, w = u on ∂Ω: the function w ∈ H¹(Ω) with u - w ∈ H¹₀(Ω) and ∫_Ω ⟨∇w, ∇v⟩ = 0 for every v ∈ H¹₀(Ω). Replacing a function by the solution of the frozen problem with the same boundary values on each ball is the comparison step of Campanato's approach to Schauder estimates. Two facts make it work:

The energy estimate holds for every measurable, uniformly elliptic a with constants 0 < λ ≤ Λ, in the form λ ‖∇(u - w)‖ ≤ ‖G‖ for any w with u - w ∈ H¹₀(Ω) and a(w, v) = 0 for every v ∈ H¹₀(Ω), as is needed when the frozen problem has a constant coefficient matrix other than the identity. Its proof is one line: test the equation for u - w ∈ H¹₀(Ω) against itself and use ellipticity.

Harmonicity comes from Weyl's lemma. A function w ∈ H¹(Ω) with ∫_Ω ⟨∇w, ∇φ⟩ = 0 for every test function φ on Ω is weakly harmonic, ∫ Δφ • w = -∫ ⟪∇φ, ∇w⟫ = 0, so it agrees almost everywhere on Ω with a harmonic function. Existence and uniqueness of the replacement on a bounded domain are the Dirichlet problem with boundary values in H¹(Ω), UniformlyEllipticOn.existsUnique_sub_mem_w1p0Submodule_forall_energyFormH1_eq_of_subset_ball.

Main declarations #

References #

theorem TauCeti.PDE.UniformlyEllipticOn.mul_norm_gradient_sub_le {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} [DecidableEq ι] {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {lam Lam : ℝ} (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) {u w : ↥(W1p mu Omega 2)} (huw : u - w ∈ w1p0Submodule mu Omega 2) (hw : ∀ (v : ↥(W1p0 mu Omega 2)), energyFormH1 a 0 0 w ↑v = 0) {G : ↥(MeasureTheory.Lp (EuclideanSpace ℝ ι) 2 (mu.restrict ↑Omega))} (hG : ∀ (v : ↥(W1p0 mu Omega 2)), energyFormH1 a 0 0 u ↑v = inner ℝ G (W1p.gradient ↑v)) :

The energy estimate for the replacement. Let a be measurable and uniformly elliptic on Ω with constants 0 < λ ≤ Λ, and let u, w ∈ H¹(Ω) have the same boundary values, u - w ∈ H¹₀(Ω), with w a weak solution of -∂ᵢ(aⁱʲ ∂ⱼw) = 0, that is a(w, v) = 0 for every v ∈ H¹₀(Ω). If the functional v ↦ a(u, v) on H¹₀(Ω) is represented by a field G ∈ L²(Ω; ℝⁿ) as a(u, v) = ⟪G, ∇v⟫_{L²}, then λ ‖∇(u - w)‖_{L²} ≤ ‖G‖_{L²}.

theorem TauCeti.PDE.exists_harmonicOnNhd_ae_eq_of_forall_energyFormH1_eq_zero {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} [DecidableEq ι] {w : ↥(W1p mu Omega 2)} (hw : ∀ (φ : TestFunction Omega ℝ ⊤), energyFormH1 (fun (x : EuclideanSpace ℝ ι) => 1) 0 0 w ((W1p.ofTestFunctionₗ mu Omega 2) φ) = 0) :
∃ (h : EuclideanSpace ℝ ι → ℝ), InnerProductSpace.HarmonicOnNhd h ↑Omega ∧ ↑↑(W1p.value w) =ᵐ[mu.restrict ↑Omega] h

Weakly harmonic Sobolev functions are harmonic. If w ∈ H¹(Ω) solves Δw = 0 weakly, in the sense that ∫_Ω ⟨∇w, ∇φ⟩ = 0 for every test function φ on Ω, then w agrees almost everywhere on Ω with a function harmonic on Ω. No boundary condition on w and no regularity of Ω is assumed.

theorem TauCeti.PDE.W1p.exists_harmonicOnNhd_ae_eq_gradient {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {p : ENNReal} [Fact (1 ≤ p)] {w : ↥(W1p mu Omega p)} (hw : WeaklyHarmonicOn mu Omega ↑↑(W1p.value w)) :
∃ (H : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι), InnerProductSpace.HarmonicOnNhd H ↑Omega ∧ ↑↑(W1p.gradient w) =ᵐ[mu.restrict ↑Omega] H

The weak gradient of a weakly harmonic Sobolev function is harmonic. If the value of w ∈ W^{1,p}(Ω) is weakly harmonic on Ω, then the weak gradient of w agrees almost everywhere on Ω with a function Ω → ℝⁿ harmonic on Ω.

theorem TauCeti.PDE.exists_sub_mem_w1p0Submodule_harmonicOnNhd_ae_eq_of_subset_ball {n : ℕ} {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ (Fin (n + 1)))} {z : EuclideanSpace ℝ (Fin (n + 1))} {R : ℝ} (hball : ↑Omega ⊆ Metric.ball z R) (u : ↥(W1p MeasureTheory.volume Omega 2)) :
∃ (w : ↥(W1p MeasureTheory.volume Omega 2)), u - w ∈ w1p0Submodule MeasureTheory.volume Omega 2 ∧ (∀ (v : ↥(W1p0 MeasureTheory.volume Omega 2)), energyFormH1 (fun (x : EuclideanSpace ℝ (Fin (n + 1))) => 1) 0 0 w ↑v = 0) ∧ (∃ (h : EuclideanSpace ℝ (Fin (n + 1)) → ℝ), InnerProductSpace.HarmonicOnNhd h ↑Omega ∧ ↑↑(W1p.value w) =ᵐ[MeasureTheory.volume.restrict ↑Omega] h) ∧ ∀ (G : ↥(MeasureTheory.Lp (EuclideanSpace ℝ (Fin (n + 1))) 2 (MeasureTheory.volume.restrict ↑Omega))), (∀ (v : ↥(W1p0 MeasureTheory.volume Omega 2)), energyFormH1 (fun (x : EuclideanSpace ℝ (Fin (n + 1))) => 1) 0 0 u ↑v = inner ℝ G (W1p.gradient ↑v)) → ‖W1p.gradient (u - w)‖ ≤ ‖G‖

Harmonic replacement. Let Ω ⊆ ℝ^{n+1} be an open set contained in a ball, and let u ∈ H¹(Ω). There is w ∈ H¹(Ω) with the same boundary values, u - w ∈ H¹₀(Ω), solving Δw = 0 weakly, ∫_Ω ⟨∇w, ∇v⟩ = 0 for every v ∈ H¹₀(Ω), such that:

  • w agrees almost everywhere on Ω with a function harmonic on Ω;
  • if ∫_Ω ⟨∇u, ∇v⟩ = ⟪G, ∇v⟫_{L²} for every v ∈ H¹₀(Ω), then ‖∇(u - w)‖_{L²(Ω)} ≤ ‖G‖_{L²(Ω)}.

The replacement is unique, by UniformlyEllipticOn.existsUnique_sub_mem_w1p0Submodule_forall_energyFormH1_eq_of_subset_ball.